clc
clear
a=3.610;
b=4.29e-2;
t=500;
x=[1:100];
r=0.410.*x./t;
p=[1 -(b+r.*t./x) a./x -a.*b./x];
roots(p)
v=roots(p);
y=v^3-(b+r.*t./x).*v^2+(a./x).*v-(a.*b./x);
plot(x,y)
xlabel P; ylabel Vm;
matlab 을 처음 사용해보는 1인입니다.
저번주부터 이리저리 검색해보고 책을 봐도 도저히 제가 원하는 그래프가 나오지 않아서 질문드립니다.
위에 보시는바와 같이 a,b,t는 고정된 값이고 x=[1:100]; <ㅡ x를 원래 p로 표현하고 싶은데 좌표상에서 x축에 p를 사용하는 방법을 몰라서
x를 그냥 넣었습니다. 이때 x는 1~100 까지의 구간을 가지고 각 구간에서 v^3-(b+r.*t./x).*v^2+(a./x).*v-(a.*b./x)=0 를 만족시키는 v값을
y축 값으로 갖도록 나타내고 싶은데 계속해서 오류만 나고 해결방법을 모르겠습니다.
p=[1 -(b+r.*t./x) a./x -a.*b./x]; 이렇게 한 이유는 p(v)=3차방정식이라 생각해서 이리 둔거고 여기서 구한 근을 roots(p)로 구하고
v=roots(p)로 두면 Inputs must be a scalar and a square matrix 이러한 에러가 나오는데 어찌 해야될지 모르겠습니다..
그리고 x가 1,2,3,4,5 ... 100 각 정수마다의 y축 값을 따로 구하고 그래프는 점이 아닌 선으로 부드럽게 연결되게 그릴수 있나요??
조언좀 꼭꼭 부탁드립니다 ㅠㅠ
글쓰기