非线性常微分方程组 matlab,matlab常微分方程和常微分方程组求解.doc

常微分方程和常微分方程组的求解

?

一、实验目的：

熟悉Matlab软件中关于求解常微分方程和常微分方程组的各种命令，掌握利用Matlab软件进行常微分方程和常微分方程组的求解。

?

二、相关知识

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

X=dsolve(‘eqn1’,’eqn2’,…)

函数dsolve用来解符号常微分方程、方程组，如果没有初始条件，则求出通解，如果有初始条件，则求出特解。

例1：求解常微分方程的MATLAB程序为：dsolve('Dy=1/(x+y)','x')，注意，系统缺省的自变量为t，因此这里要把自变量写明。

结果为：-lambertw(-C1*exp(-x-1))-x-1

其中：Y=lambertw(X)表示函数关系Y*exp(Y)=X。

例2：求解常微分方程的MATLAB程序为：Y2=dsolve('y*D2y-Dy^2=0’,’x’)

结果为：

Y2 =[ exp((x+C2)/C1)]

[ C2]

我们看到有两个解，其中一个是常数。

例3：求常微分方程组通解的MATLAB程序为：

[X,Y]=dsolve('Dx+5*x+y=exp(t),Dy-x-3*y=exp(2*t)','t')

例4&#x