微分方程问题的计算机求解

2020/5/30,1,高等应用数学问题的MATLAB求解东北大学信息学院,第7章微分方程问题的计算机求解,薛定宇、陈阳泉著高等应用数学问题的MATLAB求解，清华大学出版社2004CAI课件开发：刘莹莹、薛定宇,2020/5/30,2,高等应用数学问题的MATLAB求解东北大学信息学院,主要内容,常系数线性微分方程的解析解方法微分方程问题的数值解法特殊微分方程的数值解边值问题的计算机求解偏微分方程求解入门微分方程的框图求解,2020/5/30,3,高等应用数学问题的MATLAB求解东北大学信息学院,7.1常系数线性微分方程的解析解方法,线性常系数微分方程解析解的数学描述微分方程的解析解方法Laplace变换在线性微分方程求解中的应用特殊非线性微分方程的解析解,2020/5/30,4,高等应用数学问题的MATLAB求解东北大学信息学院,7.1.1线性常系数微分方程解析解的数学描述,2020/5/30,5,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,6,高等应用数学问题的MATLAB求解东北大学信息学院,7.1.2微分方程的解析解方法,2020/5/30,7,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-1】,2020/5/30,8,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,9,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,10,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,11,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,12,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-2】,2020/5/30,13,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,14,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-3】,2020/5/30,15,高等应用数学问题的MATLAB求解东北大学信息学院,7.1.3Laplace变换在线性微分方程求解中的应用,2020/5/30,16,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-4】,2020/5/30,17,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,18,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,19,高等应用数学问题的MATLAB求解东北大学信息学院,7.1.4特殊非线性微分方程的解析解,【例7-5】,2020/5/30,20,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-6】,2020/5/30,21,高等应用数学问题的MATLAB求解东北大学信息学院,7.2微分方程问题的数值解法,微分方程问题算法概述四阶定步长Runge-Kutta算法及MATLAB实现一阶微分方程组的数值解微分方程转换,2020/5/30,22,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.1微分方程问题算法概述,2020/5/30,23,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.1.1微分方程求解的误差与步长问题,2020/5/30,24,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,25,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,26,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.2四阶定步长Runge-Kutta算法及MATLAB实现,2020/5/30,27,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,28,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.3一阶微分方程组的数值解7.2.3.1四阶五级Runge-Kutta-Felhberg算法,2020/5/30,29,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.3.2基于MATLAB的微分方程,求解函数,2020/5/30,30,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,31,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,32,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-7】,2020/5/30,33,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,34,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-8】,2020/5/30,35,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,36,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.3.3MATLAB下带有附加参数的微分方程求解,【例7-9】,2020/5/30,37,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,38,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,39,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,40,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.4微分方程转换7.2.4.1单个高阶常微分方程处理方法,2020/5/30,41,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,42,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-10】,2020/5/30,43,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,44,高等应用数学问题的MATLAB求解东北大学信息学院,7.2.4.2高阶常微分方程组的变换方法,2020/5/30,45,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-11】,2020/5/30,46,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,47,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,48,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,49,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,50,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-12】,2020/5/30,51,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,52,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-13】,2020/5/30,53,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,54,高等应用数学问题的MATLAB求解东北大学信息学院,7.3特殊微分方程的数值解,刚性微分方程的求解隐式微分方程求解微分代数方程的求解延迟微分方程求解,2020/5/30,55,高等应用数学问题的MATLAB求解东北大学信息学院,7.3.1刚性微分方程的求解,【例7-14】,2020/5/30,56,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-14】,2020/5/30,57,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,58,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,59,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-16】,2020/5/30,60,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,61,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,62,高等应用数学问题的MATLAB求解东北大学信息学院,7.3.2隐式微分方程求解,【例7-17】,2020/5/30,63,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,64,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-18】,2020/5/30,65,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,66,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,67,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-19】,2020/5/30,68,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,69,高等应用数学问题的MATLAB求解东北大学信息学院,7.3.3微分代数方程的求解,2020/5/30,70,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-20】,2020/5/30,71,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,72,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,73,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,74,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,75,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,76,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-21】,2020/5/30,77,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,78,高等应用数学问题的MATLAB求解东北大学信息学院,7.3.4延迟微分方程求解,2020/5/30,79,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-22】,2020/5/30,80,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,81,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-23】中性延迟微分方程,2020/5/30,82,高等应用数学问题的MATLAB求解东北大学信息学院,7.4边值问题的计算机求解,线性方程边值问题的打靶算法非线性方程边值问题的打靶算法线性微分方程的有限差分算法,2020/5/30,83,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,84,高等应用数学问题的MATLAB求解东北大学信息学院,7.4.1线性方程边值问题的打靶算法,2020/5/30,85,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,86,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,87,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,88,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,89,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,90,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-24】,2020/5/30,91,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,92,高等应用数学问题的MATLAB求解东北大学信息学院,7.4.2非线性方程边值问题的打靶算法,2020/5/30,93,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,94,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,95,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-25】,2020/5/30,96,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,97,高等应用数学问题的MATLAB求解东北大学信息学院,7.4.3线性微分方程的有限差分算法,2020/5/30,98,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,99,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,100,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-26】,2020/5/30,101,高等应用数学问题的MATLAB求解东北大学信息学院,7.5偏微分方程求解入门,偏微分方程组求解二阶偏微分方程的数学描述偏微分方程的求解界面应用举例,2020/5/30,102,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.1偏微分方程组求解,2020/5/30,103,高等应用数学问题的MATLAB求解东北大学信息学院,边界条件的函数描述：,2020/5/30,104,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-27】,2020/5/30,105,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,106,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,107,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,108,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.2二阶偏微分方程的数学描述7.5.2.1椭圆型偏微分方程,2020/5/30,109,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,110,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.2.2抛物线型偏微分方程,2020/5/30,111,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.2.3双曲型偏微分方程,2020/5/30,112,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.2.4特征值型偏微分方程,2020/5/30,113,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.3偏微分方程的求解界面应用举例7.5.3.1偏微分方程求解程序概述,启动偏微分方程求解界面在MATLAB下键入pdetool该界面分为四个部分菜单系统工具栏集合编辑求解区域,2020/5/30,114,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.3.2偏微分方程求解区域绘制7.5.3.3偏微分方程边界条件描述,2020/5/30,115,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.3.4偏微分方程求解举例,【例7-28】,7.5.3.5时变解的动画显示,2020/5/30,116,高等应用数学问题的MATLAB求解东北大学信息学院,7.5.3.6函数参数的偏微分方程求解,【例7-29】,2020/5/30,117,高等应用数学问题的MATLAB求解东北大学信息学院,7.6微分方程的框图求解,Simulink简介Simulink相关模块微分方程的Simulink建模与求解,2020/5/30,118,高等应用数学问题的MATLAB求解东北大学信息学院,7.6.1Simulink简介,1990年前后出现最早的Simulink，当时名为SimuLAB，1992年改为SimulinkSimulink的名字有两重含义仿真(simu)与模型连接(link)odegroup命令可以打开自定义模块集,7.6.2Simulink相关模块,2020/5/30,119,高等应用数学问题的MATLAB求解东北大学信息学院,常用的模块：,2020/5/30,120,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,121,高等应用数学问题的MATLAB求解东北大学信息学院,7.6.3微分方程的Simulink建模与求解,建立起微分方程的Simulink模型可以用sim()函数对其模型直接求解得出微分方程的数值解,2020/5/30,122,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-30】,2020/5/30,123,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-31】,2020/5/30,124,高等应用数学问题的MATLAB求解东北大学信息学院,2020/5/30,125,高等应用数学问题的MATLAB求解东北大学信息学院,【例7-32】,2020/5/30,126,高等应用数学问题的MATLAB求解东北大学信息学院,本章内容简介,2020/5/30,127,高等应用数学问题的MATLAB求解东北大学信息学院,本章介绍了基于MATLAB符号运算工具箱dsolve()函数的线性微分方程的解析解方法，并介绍基于该函数的特殊非线性微分方程的解