作业帮MATLAB下的微分方程求解与画图
来源：学生作业帮 编辑：作业帮 分类：综合作业 时间：2024/05/24 07:22:30
MATLAB下的微分方程求解与画图
毕设做的身份认证,用的软件是MATLAB.有一个微分方程让做出图来,实在是才疏学浅,请各位大侠指教!
状态方程为：dx/dt=z-w;
dy/dt=2*y+z;
dz/dt=14*x-14*y;
dw/dt=100*x-100*w+100*(abs(w+1)-abs(w-1));
需要绘制出系统的各个状态和时间关系的二维曲线图,再绘出一个三维的空间曲线.
毕设做的身份认证,用的软件是MATLAB.有一个微分方程让做出图来,实在是才疏学浅,请各位大侠指教!
状态方程为：dx/dt=z-w;
dy/dt=2*y+z;
dz/dt=14*x-14*y;
dw/dt=100*x-100*w+100*(abs(w+1)-abs(w-1));
需要绘制出系统的各个状态和时间关系的二维曲线图,再绘出一个三维的空间曲线.
各个状态和时间关系可以画
但是四个变量,画三维空间曲线?
再问： 对啊！我也不知道是怎么回事，纳闷呢！你会求解吗?
再答： fun=@(t,xyzw)[ xyzw(3)-xyzw(4); 2*xyzw(2)+xyzw(3); 14*xyzw(1)-14*xyzw(2); 100*xyzw(1)-100*xyzw(4)+100*(abs(xyzw(4)+1)-abs(xyzw(4)-1)) ]; [t,data]=ode45(fun,[0,1],[1,0,0,0]); subplot(2,2,1) plot(t,data(:,1)) xlabel('t') ylabel('x') subplot(2,2,2) plot(t,data(:,2)) xlabel('t') ylabel('y') subplot(2,2,3) plot(t,data(:,3)) xlabel('t') ylabel('z') subplot(2,2,4) plot(t,data(:,4)) xlabel('t') ylabel('w')
但是四个变量,画三维空间曲线?
再问： 对啊！我也不知道是怎么回事，纳闷呢！你会求解吗?
再答： fun=@(t,xyzw)[ xyzw(3)-xyzw(4); 2*xyzw(2)+xyzw(3); 14*xyzw(1)-14*xyzw(2); 100*xyzw(1)-100*xyzw(4)+100*(abs(xyzw(4)+1)-abs(xyzw(4)-1)) ]; [t,data]=ode45(fun,[0,1],[1,0,0,0]); subplot(2,2,1) plot(t,data(:,1)) xlabel('t') ylabel('x') subplot(2,2,2) plot(t,data(:,2)) xlabel('t') ylabel('y') subplot(2,2,3) plot(t,data(:,3)) xlabel('t') ylabel('z') subplot(2,2,4) plot(t,data(:,4)) xlabel('t') ylabel('w')