代碼如下：
clear;clc;
[T,Y] = ode45(@Lorenz,[0 300],[0.1;0.1;0.1]);
hold on
plot3(Y(:,3),Y(:,1),Y(:,2),'b','LineWidth',0.5);
view(-30,40);
xlabel('z(t)','FontName','Times New Roman','FontSize',15);
ylabel('x(t)','FontName','Times New Roman','FontSize',15);
zlabel('y(t)','FontName','Times New Roman','FontSize',15);
hold off
function dy = Lorenz(~,y)
a=10;
b=8/3;
c=28;
dy = zeros(3,1);
% a column vector
dy(1) = a*(y(2) - y(1));
dy(2) = -y(1) * y(3)+c*y(1)-y(2);
dy(3) = y(1) * y(2)-b*y(3);
end
2.png (93.84 KB, 下載次數: 59)
下載附件
2022-7-27 21:19 上傳
|
QQ好友和群
QQ空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
