数学实验 课件 第9、10章 常微分方程实验、优化方法实验.pptx

9.1常微分方程的解析解

定义9.1含有未知函数的导数的方程称为微分方程.如果未知函数是一元函数，称为常微分方程.常微分方程的一般形式为.由几个微分方程联立而成的方程组称为微分方程组．定义9.2微分方程中出现的未知函数最高阶导数的阶数称为微分方程的阶．若方程中未知函数及其各阶导数都是一次的，称为线性常微分方程，一般表示为.若上式中的系数均与x无关，称之为常系数.二阶及二阶以上的微分方程称为高阶微分方程．

?

一阶常微分方程与高阶微分方程可以互化，已给一个n阶方程设，可将上式化为一阶方程组反过来，在许多情况下，一阶微分方程组也可化为高阶方程.所以一阶微分方程组与高阶常微分方程的理论与方法在许多方面是相通的，一阶常系数线性微分方程组也可用特征根法求解．

MATLAB中主要用dsolve求符号解析解，调用格式如下：S=dsolve(eqn,cond)求带初始条件的微分方程特解，若缺省cond，则求微分方程的通解.利用“diff”和“==”表示微分方程．例如diff(y,x)==y表示dy/dx=y．例9.1求下列微分方程的解析解（1）（2）（3）

解（1）symsy(t)abeqn=diff(y,t)==a*y+b;dsolve(eqn)ans=-(b-C1*exp(a*t))/a可知微分方程的通解为，其中，C1是任意常数．

（2）symsy(x)eqn=diff(y,x,2)==sin(2*x)-y;Dy=diff(y,x);cond=[y(0)==0,Dy(0)==1];dsolve(eqn,cond)ans=(5*sin(x))/3-sin(2*x)/3可知微分方程的特解为.

（3）symsf(t)g(t)eqns=[diff(f,t)==f+g,diff(g,t)==g-f];cond=[f(0)==1,g(0)==1];S=dsolve(eqns,cond)S=包含以下字段的struct:g:[1×1sym]f:[1×1sym]S.fans=exp(t)*cos(t)+exp(t)*sin(t)S.gans=exp(t)*cos(t)-exp(t)*sin(t)可知微分方程组的特解为.

9.2微分方程的数值解

除常系数线性微分方程可用特征根法求解，少数特殊方程可用初等积分法求解外，大部分微分方程无解析解，应用中主要依靠数值解法．考虑一阶常微分方程初值问题其中所谓数值解法，就是寻求y(x)在一系列离散节点上的近似值yk，．称为步长，通常取为常量h．最简单的数值解法是Euler法．

9.2.1欧拉法Euler法的思路极其简单：在节点处用差商近似代替导数这样导出计算公式（称为Euler格式）它能求解各种形式的微分方程.Euler法也称折线法．Euler方法只有一阶精度，改进方法有二阶Runge-Kutta法、四阶Runge-Kutta法、五阶Runge-Kutta-Felhberg法和先行多步法等，这些方法可用于解高阶常微分方程（组）初值问题.边值问题采用不同方法，如