我有函式 fi(?)=γi sin(2?sin?) for i=1,2 其中 γ1=0.01 和 γ2=0.02
φ1(0)=0.1 和 φ2(0)=0.2
?1/dt=f1(?) d?sin(?2??1)
?2/dt=f2(?) d?sin(?1??2)
其中 d=0.1
所以應該有類似這張表的東西：
t | ?1 | ?2
0.00 | 0.1 |0.2
0.01 | ... |...
0.02 | ... |...
...
100.00| ... | ...
因此，使用接收到的值需要通過坐標繪制圖形
那么問題是如何使用MATLAB在下圖上繪制函式?2(?1)？
uj5u.com熱心網友回復：
所以系統的故事可能是你從兩個不耦合且略有不同的方程開始
?1/dt=f1(?1)
?2/dt=f2(?2)
并用耦合或交換術語將它們連接起來sin(?2-?1)，
?1/dt=f1(?1) d?sin(?2??1)
?2/dt=f2(?2) d?sin(?1??2)
在 matlab 腳本中，您可以將其實作為
y0 = [ 0.1; 0.2 ];
[T,Y] = ode45(eqn,[0, 100], y0);
plot(Y(:,1),Y(:,2));
function dy_dt = eqn(t,y)
d = 0.1;
g = [ 0.01; 0.02 ];
f = g sin(2*sin(y));
exch = d*sin(y(2)-y(1));
dy_dt = f [d;-d];
end%function
這給出了幾乎一條以 結尾的對角線[pi; pi]。有了更強的耦合常數，d這會變得稍微有趣一些。
您可以將引數作為引數引數提供，然后您必須通過odeset在選項物件中宣告它們，或者使用匿名函式在求解器呼叫中系結引數。
轉載請註明出處，本文鏈接：https://www.uj5u.com/houduan/394541.html