matlab中常微分方法,Matlab中常微分方程数值解法讲解.PDF

造最优秀、专业和权威的Matlab 技术交流平台!

Matlab 中常微分方程数值解法讲解

作者：dynamic

时间：2008.12.10

版权：All Rights Reserved By

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★

Matlab Sky 联盟打造最优秀、专业和权威的Matlab 技术交流平台!

网址： /com/org/net

邮箱：matlabsky@

QQ 群 405 10634满了)满了)

论坛拥有40 多个专业版块，内容涉及资料下载、视频教学、数学建模、数学运算、程序设计、GUI 开发、simulink

仿真、统计概率、拟合优化、扩展编程、算法研究、控制系统、信号通信、图像处理、经济金融、生物化学、航

航天、人工智能、汽车设计、机械自动化、毕业设计等几十个方面!

请相信我们：1.拥有绝对优秀的技术人员，热情的版主，严谨负责的管理团队

2.免费提供技术交流和在线解答

★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★

1

造最优秀、专业和权威的Matlab 技术交流平台!

1.ODE 解算器简介3

2.微分方程转换5

3.刚性/非刚性问题8

4. 隐式微分方程(IDE) 10

5.微分代数方程(DAE) 15

6.延迟微分方程(DDE) 18

7.边值问题(BVP) 20

2

造最优秀、专业和权威的Matlab 技术交流平台!

1.ODE 解算器简介

先来认识下常微分方程(ODE)初值问题解算器(solver)

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

sxint = deval(sol,xint)

解算器(odesolver)

解算器 问题类型 精确度 说明

非 刚 性

ode45 中等 4-5 阶龙格库塔，对以所有问题的首先解算器

(nonstiff)

基于Bogacki-Shampine 2-3 阶Runge-Kutta 公式，有时对轻度的刚度方程，它

ode23 非刚性 低

可以比ode45 有更好的效果，在相 的精度是，要比ode45 更小的步长

变阶次Adams-Bashforth-Moutlon 算法，此算法使用前几次节点上的值来计算

ode113 非刚性 低到高 当前节点处的解，因此在相 jingdu 下，比ode45 和ode23 更快些，适用用

于高阶或者需要大量计算的问题，但不适合不连续的系统

刚性系统的变阶次多步解法，此算法使用最新的数值差分公式，如果使用

ode15s 刚性(stiff) 低到中

ode45 计算比