matlab中函数图形洛伦茨系统中已知参数b=8、3,r=23,sigma=9.9function yp=func_nd(t,y)yp=[10*(y(2)-y(1));28*y(1)-y(2)-y(1)*y(3);-8/3*y(3)+y(1)*y(2)];要得到3个波形图 该如何做
来源:学生作业帮助网 编辑:作业帮 时间:2024/11/16 22:44:00
matlab中函数图形洛伦茨系统中已知参数b=8、3,r=23,sigma=9.9function yp=func_nd(t,y)yp=[10*(y(2)-y(1));28*y(1)-y(2)-y(1)*y(3);-8/3*y(3)+y(1)*y(2)];要得到3个波形图 该如何做
matlab中函数图形
洛伦茨系统中已知参数b=8、3,r=23,sigma=9.9
function yp=func_nd(t,y)
yp=[10*(y(2)-y(1));
28*y(1)-y(2)-y(1)*y(3);
-8/3*y(3)+y(1)*y(2)];
要得到3个波形图 该如何做
matlab中函数图形洛伦茨系统中已知参数b=8、3,r=23,sigma=9.9function yp=func_nd(t,y)yp=[10*(y(2)-y(1));28*y(1)-y(2)-y(1)*y(3);-8/3*y(3)+y(1)*y(2)];要得到3个波形图 该如何做
这里y(1),y(2),y(3)分别表示x,y,z,程序如下:
function dy=Lorenz(t,y)
dy=zeros(3,1);
dy(1)=10*(-y(1)+y(2));
dy(2)=28*y(1)-y(2)-y(1)*y(3);
dy(3)=y(1)*y(2)-8*y(3)/3;
end
%然后ode45求解
[t,y]=ode45('Lorenz',[0,30],[12,2,9]);
figure(1)
plot(t,y(:,1));
figure(2)
plot(t,y(:,2));
figure(3)
plot(t,y(:,3));
figure(4)
plot(y(:,1),y(:,2),y(:,3));
view([20,42])
%时间0-30秒,初始值(12,2,9)
得到结果图