常微分方程(ODEs)的MATLAB数值解法.pdfVIP

造最优秀、专业和权威的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 计算比