基于MATLAB微分方程与数值解实验研究.pdfVIP

实验目的

1、学会用Matlab求简单微分方程的解析解.

2、学会用Matlab求微分方程的数值解.

实验内容

1、求简单微分方程的解析解.

2、求微分方程的数值解.

3、数学建模实例

4、实验作业.

实验软件MATLAB

微分方程的解析解

求微分方程（组）的解析解命令:

dsolve(‘方程1’,‘方程2’,…‘方程n’,‘初始条件’,‘自

:DD2D3

记号在表达微分方程时，用字母表示求微分，、等

.D

表示求高阶微分任何后所跟的字母为因变量，自变量可以指

定或由系统规则选定为确省.

2

dy

例如，微分方程0应表达为：D2y=0.

2

dx

1du2

例求1+u的通解.

dt

解输入命令：dsolve(Du=1+u^2,t)

结果：utg(t-c)ToMatlab（ff1）

2

例求微分方程的特解.

2

dydy

+4+29y0

2dx

dx

y(0)0,y(0)15



解输入命令:y=dsolve(D2y+4*Dy+29*y=0,y(0)=0,Dy(0)=15,x)

ToMatlab（ff2）

结果为:y=3e-2xsin（5x）

例3求微分方程组的通解.

dx

dt2x−3y+3z



dy

4x−5y+3z

dt

dz

4x−4y+2z

dt



解输入命令：

[x,y,z]=dsolve(Dx=2*x-3*y+3*z,Dy=4*x-5*y+3*z,Dz=4*x-4*y+2*z,t)；

x=simple(x)%将x化简

y=simple(y)ToMatlab（ff3）

z=simple(z)

结果为：x=(c-c+c+ce-3t-ce-3t)e2t

12323

y-ce-4t+ce-4t+ce-3t-ce-3t+(c-c+c)e2t

1223123

z