[問題] 一個chaotic system的問題已回收
小弟我是一個Matlab的新手
想要問一個chaotic system方程式的寫法
今天我要寫一個
x(t)''+a*x(t)'+ b*x(t) + c*x(t)^3 = k cos (wt)的方程式
我所看的原文資料有將方程式作改寫
其中
x1(t)=x(t)
x2(t)=x1(t)'
移項過後變成
x2(t)'=k*cos(wt)-a*x2(t)-b*x1(t)-c*x1(t)^3
我想要畫出它的圖形
這是我的原本程式：
close all;
a=0.4;
b=-1.1;
r=1;
w=1.8;
t0=[0,50];
t=t0*pi;
k=7
y0 = [1 0];
duffing=@(t0,x)[x(2);k*cos(w*t0)-a*x(2)-b*x(1)-r*x(1).^3];
[T,X] = ode45(duffing,t0,y0);
% plot phase plane
figure(1)
plot(X(:,1),X(:,2),'-k');
xlabel('x'); ylabel('x''');
% plot exact and numerical answers
figure(2)
plot(T,X(:,1),'-b',T,X(:,2),'--r');
xlabel('t'); legend('x','x''',0);
我的原本程式和原文中所跑出來的圖
差異有點大
如果我今天要做這樣的方程式
我應該用什麼樣的指令比較適合???
先謝謝各位的解答!!
http://www.wretch.cc/album/album.php?id=vul3coy93&book=7
不好意思，小弟不會貼圖
只好放在無名的相簿裡面
我自己所跑出來的裡面
紅色虛線的部份
很明顯和原文中的圖不太一樣
而且差異有點大
不知道是不是我的寫法不太對
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 203.68.167.87
※ 編輯: vul3coy93 來自: 203.68.167.87 (07/22 15:56)
推
07/23 11:48, , 1F
07/23 11:48, 1F
※ 編輯: vul3coy93 來自: 203.68.167.87 (07/23 13:22)
推
07/23 22:09, , 2F
07/23 22:09, 2F
→
07/23 22:09, , 3F
07/23 22:09, 3F
MATLAB 近期熱門文章
PTT數位生活區 即時熱門文章
