matlab差分方程实验报告,实验二微分方程与差分方程模型matlab求解.doc

实验二微分方程与差分方程模型matlab求解

实验二： 微分方程与差分方程模型Matlab求解

一、实验目的[1] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析；[2] 熟悉MATLAB软件关于微分方程求解的各种命令；

[3] 通过范例学习建立微分方程方面的数学模型以及求解全过程；

[4]Logistic模型的求解与混沌的产生过程。

二、实验原理1 求解一阶微分方程 ???

(1) 求通解

输入：

dsolve('Dy=1+y^2')

输出：

ans =

tan(t+C1)

(2)求特解

输入：

dsolve('Dy=1+y^2','y(0)=1','x')

指定初值为1，自变量为x

输出：

ans =

tan(x+1/4*pi)

例2 求解二阶微分方程

原方程两边都除以，得

输入:

dsolve('D2y+(1/x)*Dy+(1-1/4/x^2)*y=0','y(pi/2)=2,Dy(pi/2)=-2/pi','x')

ans =

- (exp(x*i)*(pi/2)^(1/2)*i)/x^(1/2) + (exp(x*i)*exp(-x*2*i)*(pi/2)^(3/2)*2*i)/(pi*x^(1/2))

试试能不用用simplify函数化简</