Lorzen混沌方程求解 ode45

ode45是matlab中的求解微分方程數值解的函式，使用語法為：

[t,y] = ode45(odefun,tspan,y0)

[t,y] = ode45(odefun,tspan,y0,options)

[t,y,te,ye,ie] = ode45(odefun,tspan,y0,options)

sol = ode45(___)

odefun需要在函式中編寫好，以lorzen混沌方程為例

function dydt =

lorenz_diff
(t,y)
dydt =[-
10*(y
(1)-
y(2)
);30*
y(1)
-y(1
)*y(
3)-y
(2);
-8/3
*y(3
)+y(
1)*y
(2)]
;

clear

clcy0 =[0
,2,9
];[t,y]
=ode45
(@lorenz_diff,[0
,200
],y0)
;plot(y
(:,1
),y(
:,3)
,'-');
grid on;

ode求解器的事件 Event 屬性

檢測事件 matlab微分方程如何設定變數的範圍 如dy1 y2 dy2 y1 1 其中y1的範圍為0 解 m檔案 function value,isterminal,direction events1 t,y value y 1 4 isterminal 1 direction 0 命令視窗 dy...

ode求解器的事件 Event 屬性

檢測事件 matlab微分方程如何設定變數的範圍 如dy1 y2 dy2 y1 1 其中y1的範圍為0 解 m檔案 function value,isterminal,direction events1 t,y value y 1 4 isterminal 1 direction 0 命令視窗 dy...

不定方程求解

不定方程ax by c，對於給定的係數a和b 常數項c，求出在區間p q範圍內方程的乙個整數解x y，滿足x和y不全為0且x取值最大。題目保證a b c x y p q都是 10000 10000範圍之內的整數。輸入分為兩部分，第一部分是三個整數a b c且a b不為0，第二部分包含多組輸入，每組輸...