matlab在常微分方程的应用,MATLAB在求解常微分方程中的应用

MATLAB课程论文

dy dtty 1,两边积分可得通解为y ce 1.其中c为任意常数.有些常微分方程可用一些技

巧,如分离变量法,积分因子法,常数变异法,降阶法等可化为可积分的方程而求得解析解.

线性常微分方程的解满足叠加原理,从而他们的求解可归结为求一个特解和相应齐次

微分方程的通解.一阶变系数线性微分方程总可用这一思路求得显式解。高阶线性常系数微

分方程可用特征根法求得相应齐次微分方程的基本解，再用常数变异法求特解。

一阶常微分方程与高阶微分方程可以互化，已给一个n阶方程

(n)(n 1)y f(t,y',y", ,y) (2)

(n 1)y y,y y', ,y y12n设，可将上式化为一阶方程组

y1' y2 y' y23 y' yn n 1 yn' f(t,y1,y2, ,yn) (3)

反过来，在许多情况下，一阶微分方程组也可化为高阶方程。所以一阶微分方程组与

高阶常微分方程的理论与方法在许多方面是相通的，一阶常系数线性微分方程组也可用特

征根法求解。

3、微分方程的数值解法

除常系数线性微分方程可用特征根法求解，少数特殊方程可用初等积分法求解外，大

部分微分方程无限世界，应用中主要依靠数值解法。考虑一阶常微分方程初值问题

y'(t) f(t,y(t))t,0 t tf

y(t0) y0 (4)

其中y (y1,y2, ,ym)',f (f1,f2, ,fm)',y0 (y10,y20, ,ym0)'.所谓数值解法，就是寻求

y(t)在一系列离散节点t0 t1 tn tf上的近似值yk,k 0,1, ,n称hk tk 1 tk为步

长，通常取为常量h。最简单的数值解法是Euler法。

Euler法的思路极其简单：在节点出用差商近似代替导数

y'(tk) y(tk 1) y(tk)

h (5)

这样导出计算公式(称为Euler格式)

yk 1 yk hf(tk,yk),k 0,1,2, (6)

他能求解各种形式的微分方程。Euler法也称折线法。

4、解微分方程的MATLAB命令

MATLAB中主要用dsolve求符号解析解,0de45,0de23,ode15s求数值解 。

在MATLAB中，由函数dsolve()解决常微分方程(组)的求解问题，其具体格式如下：

r = dsolve('eq1,eq2,...', 'cond1,cond2,...', 'v')