matlab微分方程画图
来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/11/18 12:01:33
matlab微分方程画图
我这儿有个matlab的代码,解四维常系数线性微分方程组的,可我不会操作,不知道该如何画图,麻烦先查下代码,帮忙画个图,
function hh
X0=[195,1500,62,200];
TSPAN=[0 6];
[t,X] = ode45(@myfun,TSPAN,X0)
plot(t,X)
legend('x','y','z','u','Location','NorthWest')
xlabel('t'),ylabel('X')
function dX=myfun(t,X)
x=X(1);y=X(2);z=X(3);u=X(4);
dxdt=-0.0027*z-0.0000067*u;
dydt=-0.0208*z-0.000033*u ;
dzdt=-0.00403*x-0.00005*y ;
dudt=-0.0284*x-0.0001*y ;
dX=[dxdt;dydt;dzdt;dudt];
我这儿有个matlab的代码,解四维常系数线性微分方程组的,可我不会操作,不知道该如何画图,麻烦先查下代码,帮忙画个图,
function hh
X0=[195,1500,62,200];
TSPAN=[0 6];
[t,X] = ode45(@myfun,TSPAN,X0)
plot(t,X)
legend('x','y','z','u','Location','NorthWest')
xlabel('t'),ylabel('X')
function dX=myfun(t,X)
x=X(1);y=X(2);z=X(3);u=X(4);
dxdt=-0.0027*z-0.0000067*u;
dydt=-0.0208*z-0.000033*u ;
dzdt=-0.00403*x-0.00005*y ;
dudt=-0.0284*x-0.0001*y ;
dX=[dxdt;dydt;dzdt;dudt];
你的方法对么,没有问题啊
我也是菜鸟,用另一个方法帮你做了
参考一下吧
key=dsolve('Dx=-0.0027*z-0.0000067*u','Dy=-0.0208*z-0.000033*u','Dz=-0.00403*x-0.00005*y','Du=-0.0284*x-0.0001*y','x(0)=195','y(0)=1500','z(0)=62','u(0)=200');
t=0:0.1:10;
x=eval(key.x);
y=eval(key.y);
z=eval(key.z);
u=eval(key.u);
plot(t,x,t,y,t,z,t,u)
祝君好运
我也是菜鸟,用另一个方法帮你做了
参考一下吧
key=dsolve('Dx=-0.0027*z-0.0000067*u','Dy=-0.0208*z-0.000033*u','Dz=-0.00403*x-0.00005*y','Du=-0.0284*x-0.0001*y','x(0)=195','y(0)=1500','z(0)=62','u(0)=200');
t=0:0.1:10;
x=eval(key.x);
y=eval(key.y);
z=eval(key.z);
u=eval(key.u);
plot(t,x,t,y,t,z,t,u)
祝君好运