求高手指教啊

青枫林霰

matlab中遇到ode45（，）是什么意思啊，可以举出例子吗？

shaox

" Q0 O4 J, f  W1 gode45表示五分之四阶的龙格库塔芬尔格算法
例2-45  求解描述振荡器的经典的Ver der Pol微分方程
y(0)=1，y’(0)=0
9 O5 e& A* T2 Q- T) _+ R  L4 P令x1=y，x2=dy/dt，则
) a' h* P0 N. E- I) Cdx1/dt = x2
! x. G" \( a. N0 _+ y& O7 Z1 }, C4 vdx2/dt = μ(1-(x1)^2)*x2-x1
编写函数文件verderpol.m：
function xprime = verderpol(t,x)
4 q/ \6 e* G  V  u0 m, f: Yglobal MU
. l, k+ b8 A/ D: Q5 wxprime = [x(2);MU*(1-x(1)^2)*x(2)-x(1)];
再在命令窗口中执行：
>>global MU
8 K" S6 E/ V: P5 W: A+ g5 {>>MU = 7;
>>Y0=[1;0]
& N& ^. j( s6 P& x( ^% m2 @3 `& m>>[t,x] = ode45(‘verderpol’,0,40,Y0);
>>x1=x(:,1);x2=x(:,2);
$ U. ^9 j0 u3 s8 \3 `4 ^2 ^) r, Q& i>>plot(t,x1,t,x2)
图形结果为图2-20。
: v1 O6 ~% _" |" }
图2-20  Ver der Pol微分方程图
! }! l  t8 N" }$ E" u$ O

青枫林霰

有参考书吗？推荐一下啊。我看得有2,3,4阶，没有看到5分之4阶啊。

shaox

青枫林霰 发表于 2012-8-6 09:55
有参考书吗？推荐一下啊。我看得有2,3,4阶，没有看到5分之4阶啊。

; x) n& C+ J6 D/ K" D. t3 G0 j# o就是运用组合的4/5阶龙格库塔芬尔格算法