相位圖是指p(t)與q(t)之間的變化。如何得到可以按下列方法來求解。
1、根據給出的微分方程組,自定義壹個函數ode_fun,其內容
function dy=ode_fun(t,x)
dy=[x(1)*(1-x(1))*(2-2*x(2));x(2)*(1-x(2))*(0.5-2*x(1))]; %p—x(1)、q—x(2)
2、利用ode45()函數,求解t、p、q 的數值解
3、利用plot()函數,繪制p(t)與q(t)相位曲線圖
4、假定本微分方程組的初始條件,p(0)=0.01;q(0)=0.05
5、按上述方法編程,可得到如下結果。