查看: 2556|回复: 3

[代码资源] 求高手指教啊

[复制链接]

字体大小: 正常放大

青枫林霰

45 主题	7 听众	1432 积分

TA的每日心情

	奋斗 2015-6-8 12:52

签到天数: 238 天

[LV.7]常住居民III

群组: 数学建模认证项目实训

群组: 2012第二期MCM/ICM优秀

群组: 山西大学数模

群组: 英语科技论文写作实训

群组: 认证人才交流群

电梯直达

1^#

发表于 2012-8-6 09:35 |只看该作者 |倒序浏览

|招呼Ta 关注Ta |邮箱已经成功绑定

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

zan

转播0 淘帖0 分享0 收藏0 支持0 反对0 微信

使用道具举报

shaox

16 主题	42 听众	1061 积分

升级 6.1%

TA的每日心情

	开心 2013-9-15 16:21

签到天数: 266 天

[LV.8]以坛为家I

群组: 数学建摸协会

群组: 西安交大数学建模

群组: 学术交流A

群组: 数学建模认证项目实训

2^#

发表于 2012-8-6 09:47 |只看该作者 |招呼Ta 关注Ta

本帖最后由 shaox 于 2012-8-6 09:53 编辑

ode45表示五分之四阶的龙格库塔芬尔格算法
例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
编写函数文件verderpol.m：
function xprime = verderpol(t,x)
global MU
xprime = [x(2);MU*(1-x(1)^2)*x(2)-x(1)];
再在命令窗口中执行：
>>global MU
>>MU = 7;
>>Y0=[1;0]
>>[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微分方程图