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

实验二： 微分方程与差分方程模型Matlab求解 一、实验目的 [1] 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； [2] 熟悉MATLAB软件关于微分方程求解的各种命令； [3] 通过范例学习建立微分方程方面的数学模型以及求解全过程； [4] 熟悉离散 Logistic模型的求解与混沌的产生过程。 二、实验原理 1. 微分方程模型与MATLAB求解 解析解 用MATLAB命令dsolve(‘eqn1’,’eqn2’, ...) 求常微分方程（组）的解析解。其中 （1） 微分方程 例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))