Matlab中龙格-库塔(Runge-Kutta)方法原理及实现_龙格库塔法 - CSDN文库

4星 · 超过85%的资源需积分: 36 123 浏览量更新于2023-03-16 评论 3 收藏 70KB DOC 举报

身份认证购VIP最低享 7 折!

领优惠券(最高得80元）

资源详情

资源评论

资源推荐

函数功能编辑本段回目录

ode 是专门用于解微分方程的功能函数，他有 ode23,ode45,ode23s 等等，采用的是

Runge-Kutta 算法。ode45 表示采用四阶，五阶 runge-kutta 单步算法,截断误差为

(Δx)³。解决的是 Nonsti(非刚性)的常微分方程.是解决数值解问题的首选方法,若长时间

没结果,应该就是刚性的,换用 ode23 来解.

使用方法编辑本段回目录

[T,Y]=ode45(odefun,tspan,y0)

odefun是函数句柄,可以是函数文件名,匿名函数句柄或内联函数名

tspan 是区间 [t0tf]或者一系列散点[t0,t1,...,tf]

y0   是初始值向量

T   返回列向量的时间点 www.iLoveMatlab.cn

Y   返回对应 T 的求解列向量

[T,Y]=ode45(odefun,tspan,y0,options) 《Simulink 与信号处理》

options 是求解参数设置,可以用 odeset 在计算前设定误差,输出参数,事件等

[T,Y,TE,YE,IE]=ode45(odefun,tspan,y0,options)

在设置了事件参数后的对应输出 www.iLoveMatlab.cn

TE   事件发生时间 book.iLoveMatlab.cn

YE   事件解决时间

IE   事件消失时间

sol=ode45(odefun,[t0tf],y0...)

sol   结构体输出结果

应用举例编辑本段回目录

1 求解一阶常微分方程

程序:

本内容试读结束，登录后可阅读更多

下载后可阅读完整内容，剩余4页未读，立即下载

评论1

ygzw_n

2013-06-20

初学常微分方程的同学可以适用，在机房用的此代码，挺简洁的

buddyyang251147

粉丝: 0
资源: 2

会员权益专享

图片转文字

全年可省5，000元立即开通

最新资源

资源上传下载、课程学习等过程中有任何疑问或建议，欢迎提出宝贵意见哦~我们会及时处理！点击此处反馈