代碼如下:
- 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 (90.2 KB)
下載附件
2023-6-7 17:14 上傳
1.png (38.38 KB)
下載附件
2023-6-7 17:14 上傳
【必讀】版權免責聲明
1、本主題所有言論和內容純屬會員個人意見,與本論壇立場無關。2、本站對所發內容真實性、客觀性、可用性不做任何保證也不負任何責任,網友之間僅出于學習目的進行交流。3、對提供的數字內容不擁有任何權利,其版權歸原著者擁有。請勿將該數字內容進行商業交易、轉載等行為,該內容只為學習所提供,使用后發生的一切問題與本站無關。 4、本網站不保證本站提供的下載資源的準確性、安全性和完整性;同時本網站也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的損失或傷害。 5、本網站所有軟件和資料均為網友推薦收集整理而來,僅供學習用途使用,請務必下載后兩小時內刪除,禁止商用。6、如有侵犯你版權的,請及時聯系我們(電子郵箱1370723259@qq.com)指出,本站將立即改正。
|