求高手指教啊

青枫林霰

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

shaox

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

2 H8 n' C/ g$ I. _$ w/ H3 |就是运用组合的4/5阶龙格库塔芬尔格算法

青枫林霰

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

shaox

ode45表示五分之四阶的龙格库塔芬尔格算法
4 K5 _9 m* Q+ B; K# Q" a* o例2-45  求解描述振荡器的经典的Ver der Pol微分方程
) l- R; d8 u! Ay(0)=1，y’(0)=0
令x1=y，x2=dy/dt，则
& }, K$ n8 S& ~# Edx1/dt = x2
7 j' n7 l0 _' B- T, E3 `4 }dx2/dt = μ(1-(x1)^2)*x2-x1
编写函数文件verderpol.m：
3 ~1 ^& K4 P- ]7 i4 Z- ifunction xprime = verderpol(t,x)
global MU
1 j- J& J- Q6 j) I  T0 Y0 qxprime = [x(2);MU*(1-x(1)^2)*x(2)-x(1)];
再在命令窗口中执行：
* X# N& K% q- g* |$ C, ?& U>>global MU
' s0 `& X3 M& L4 v& v& i$ u>>MU = 7;
>>Y0=[1;0]
>>[t,x] = ode45(‘verderpol’,0,40,Y0);
5 Q6 a# j) }$ D5 Q>>x1=x(:,1);x2=x(:,2);
2 e1 J7 b6 G1 }2 [$ m3 x) K>>plot(t,x1,t,x2)
图形结果为图2-20。

图2-20  Ver der Pol微分方程图