作业帮matlab中函数图形
来源：学生作业帮 编辑：拍题作业网作业帮 分类：综合作业 时间：2024/04/27 18:35:41
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个波形图 该如何做
洛伦茨系统中已知参数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个波形图 该如何做
这里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)
得到结果图
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)
得到结果图