代码如下:
- 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
复制代码
【必读】版权免责声明
1、本主题所有言论和内容纯属会员个人意见,与本论坛立场无关。2、本站对所发内容真实性、客观性、可用性不做任何保证也不负任何责任,网友之间仅出于学习目的进行交流。3、对提供的数字内容不拥有任何权利,其版权归原著者拥有。请勿将该数字内容进行商业交易、转载等行为,该内容只为学习所提供,使用后发生的一切问题与本站无关。 4、本网站不保证本站提供的下载资源的准确性、安全性和完整性;同时本网站也不承担用户因使用这些下载资源对自己和他人造成任何形式的损失或伤害。 5、本网站所有软件和资料均为网友推荐收集整理而来,仅供学习用途使用,请务必下载后两小时内删除,禁止商用。6、如有侵犯你版权的,请及时联系我们(电子邮箱1370723259@qq.com)指出,本站将立即改正。
|
|