chap-7-微分方程问题的求解

第7章微分方程问题的计算机求解,主要内容,常系数线性微分方程的解析解方法微分方程问题的数值解法特殊微分方程的数值解边值问题的计算机求解偏微分方程求解入门微分方程的框图求解,7.1常系数线性微分方程的解析解方法,线性常系数微分方程解析解的数学描述微分方程的解析解方法Laplace变换在线性微分方程求解中的应用特殊非线性微分方程的解析解,7.1.1线性常系数微分方程解析解的数学描述,7.1.2微分方程的解析解方法,【例7-1】,【例7-2】,【例7-3】,7.1.3Laplace变换在线性微分方程求解中的应用,【例7-4】,7.1.4特殊非线性微分方程的解析解,【例7-5】,【例7-6】,7.2微分方程问题的数值解法,微分方程问题算法概述四阶定步长Runge-Kutta算法及MATLAB实现一阶微分方程组的数值解微分方程转换,7.2.1微分方程问题算法概述,微分方程求解的误差与步长问题,7.2.2四阶定步长Runge-Kutta算法及MATLAB实现,7.2.3一阶微分方程组的数值解7.2.3.1四阶五级Runge-Kutta-Felhberg算法,7.2.3.2基于MATLAB的微分方程,求解函数,【例7-7】,【例7-8】,7.2.3.3MATLAB下带有附加参数的微分方程求解,【例7-9】,7.2.4微分方程转换7.2.4.1单个高阶常微分方程处理方法,【例7-10】,7.2.4.2高阶常微分方程组的变换方法,【例7-11】,【例7-12】,【例7-13】,7.3特殊微分方程的数值解,刚性微分方程的求解隐式微分方程求解微分代数方程的求解延迟微分方程求解,7.3.1刚性微分方程的求解,【例7-14】,【例7-14】,【例7-16】,7.3.2隐式微分方程求解,【例7-17】,【例7-18】,【例7-19】,7.3.3微分代数方程的求解,【例7-20】,【例7-21】,7.3.4延迟微分方程求解,【例7-22】,【例7-23】中性延迟微分方程,7.4边值问题的计算机求解,线性方程边值问题的打靶算法非线性方程边值问题的打靶算法线性微分方程的有限差分算法,7.4.1线性方程边值问题的打靶算法,【例7-24】,7.4.2非线性方程边值问题的打靶算法,【例7-25】,7.4.3线性微分方程的有限差分算法,【例7-26】,7.5偏微分方程求解入门,偏微分方程组求解二阶偏微分方程的数学描述偏微分方程的求解界面应用举例,7.5.1偏微分方程组求解,边界条件的函数描述：,【例7-27】,7.5.2二阶偏微分方程的数学描述7.5.2.1椭圆型偏微分方程,7.5.2.2抛物线型偏微分方程,7.5.2.3双曲型偏微分方程,7.5.2.4特征值型偏微分方程,7.5.3偏微分方程的求解界面应用举例7.5.3.1偏微分方程求解程序概述,启动偏微分方程求解界面在MATLAB下键入pdetool该界面分为四个部分菜单系统工具栏集合编辑求解区域,7.5.3.2偏微分方程求解区域绘制7.5.3.3偏微分方程边界条件描述,7.5.3.4偏微分方程求解举例,【例7-28】,7.5.3.5时变解的动画显示,7.5.3.6函数参数的偏微分方程求解,【例7-29】,7.6微分方程的框图求解,Simulink简介Simulink相关模块微分方程的Simulink建模与求解,7.6.1Simulink简介,1990年前后出现最早的Simulink，当时名为SimuLAB，1992年改为SimulinkSimulink的名字有两重含义仿真(simu)与模型连接(link)odegroup命令可以打开自定义模块集,7.6.2Simulink相关模块,常用的模块：,7.6.3微分方程的Simulink建模与求解,建立起微分方程的Simulink模型可以用sim()函数对其模型直接求解得出微分方程的数值解,【例7-30】,【例7-31】,【例7-32】,本章内容简介,本章介绍了基于MATLAB符号运算工具箱dsolve()函数的线性微分方程的解析解方法，并介绍基于该函数的特殊非线性微分方程的解析解。,对一般非线性微分方程来说，解析解是不存在的，只能依赖数值解的方法对其进行研究。引入了数值解的概念，并以最简单的一阶微分方程的Euler算法为例，介绍了一般数值解法的思路并介绍了变步长求解的概念，还介绍了MATLAB下的微分方程数值求解函数ode45()，通过例子演示了该函数的使用方法。,微分方程初值函数能直接求解的方程是一阶显式微分方程组，若给出的方程不是这类函数，则需要通过本书介绍的方法选择一组状态变量，将原方程变换成一阶显式微分方程组，以便用给定的求解函数直接求解。若某微分方程模型求解速度极慢，则有可能为刚性方程，需要调用ode15s()等函数来求解，此外，其他类型的微分方程，如微分代数方程、隐式微分方程与延迟微分方程等，也可以由MATLAB语言提供的现成函数直接求解。二阶微分方程的边值问题可以由本书提供的三种算法求解。,偏微分方程可以由MATLAB提供的现成函数直接求解，而x-y平面的偏微分方程可以由MATLAB语言的偏微分方程工具箱提供的界面直接求解，而高维偏微分方程可以由该工具箱提供的现成函数直接求解。Simulink是MATLAB中的一个很重要的系统仿真平台，可以用该高阶以框图的形式建立