数学建模社区-数学中国

标题: 求高手指教啊 [打印本页]

---

作者: 青枫林霰    时间: 2012-8-6 09:35
标题: 求高手指教啊
matlab中遇到ode45（，）是什么意思啊，可以举出例子吗？

---

作者: shaox    时间: 2012-8-6 09:47


ode45表示五分之四阶的龙格库塔芬尔格算法
& \7 K/ `/ {0 T) y/ c# B) ~例2-45  求解描述振荡器的经典的Ver der Pol微分方程
y(0)=1，y’(0)=0
令x1=y，x2=dy/dt，则
dx1/dt = x2
dx2/dt = μ(1-(x1)^2)*x2-x1
( N4 D& h: W* ]编写函数文件verderpol.m：
3 w* L$ x( d# t, l/ |# @. E/ dfunction xprime = verderpol(t,x)
. J. N9 }2 u3 o9 }/ g2 @/ vglobal MU
xprime = [x(2);MU*(1-x(1)^2)*x(2)-x(1)];
7 ~" H8 a9 }8 y" L再在命令窗口中执行：
3 `. y7 X: J0 t9 i1 b7 b( W- L>>global MU
>>MU = 7;
, V- F( m& \! g) X7 L% D>>Y0=[1;0]
0 V1 l$ ^# W: k0 l) [>>[t,x] = ode45(‘verderpol’,0,40,Y0);
>>x1=x(:,1);x2=x(:,2);
>>plot(t,x1,t,x2)
图形结果为图2-20。

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

---

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

---

作者: shaox    时间: 2012-8-6 10:14

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

, X) G  ]. Y3 v" `+ y" x1 ]就是运用组合的4/5阶龙格库塔芬尔格算法

---

欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)

Powered by Discuz! X2.5