matlab 绘图求偏微分方程组程序三维程序方程组为('Dx=-y-z','Dy...

发布网友发布时间：2024-10-05 16:19

共3个回答

热心网友时间：2024-10-05 16:58

t=250;

dt=0.01;

tlen=t/dt;

x=zeros(1,tlen);

y=zeros(1,tlen);

z=zeros(1,tlen);

for k=1:tlen-1

dx=(-y(k)-z(k))*dt;

dy=(x(k)+0.2*y(k))*dt;

dz=(0.2+z(k)*(x(k)-5))*dt;

x(k+1)=x(k)+dx;

y(k+1)=y(k)+dy;

z(k+1)=z(k)+dz;

end

f_h = figure('numbertitle','off');

a_h = axes('parent',f_h);

set(a_h,'box','on','projection','perspective','dataaspectratio',[1 1 1])

line(x,y,z)

热心网友时间：2024-10-05 16:57

使用ode45

热心网友时间：2024-10-05 17:01

>> dy=inline('[-y(2)-y(3);y(1)+0.2*y(2);0.2+y(3)*(y(1)-5)]','t','y');

>> [x y]=ode45(dy,200,[0 0 0]);

>> plot3(y(:,1),y(:,2),y(:,3))

声明：本网页内容为用户发布，旨在传播知识，不代表本网认同其观点，若有侵权等问题请及时与本网联系，我们将在第一时间删除处理。E-MAIL:11247931@qq.com

matlab 绘图 求偏微分方程组程序 三维程序 方程组为('Dx=-y-z','Dy...