数学试验及建模.PPTVIP

实验06 微分方程实验 实验目的 1 掌握用Matlab求简单微分方程的解析解. 2 掌握用Matlab求微分方程的数值解. 3 掌握利用微分方程模型解决实际问题. 实验内容 1 求简单微分方程的解析解. 2 求微分方程的数值解 3 数学建模实例 4 练习与思考 问题背景和实验目的 Matlab 解初值问题 dsolve 求解析解 dsolve 的使用 dsolve 举例 dsolve 举例 dsolve 举例 Euler 折线法 初值问题的Euler折线法 Euler 折线法举例 Euler 折线法源程序 Euler折线法举例（续） Runge-Kutta 方法 Runge-Kutta 方法 四阶 R-K 方法源程序 Runge-Kutta 方法 Euler 法与 R-K法误差比较 Matlab提供的ODE求解器 数值求解举例 数值求解举例 Matlab 求解微分方程小结 先编写函数文件 verderpol.m function xprime=verderpol(t,x) global mu; xprime=[x(2); mu*(1-x(1)^2)*x(2) - x(1)]; 再编写脚本文件 vdpl.m，在命令窗口直接运行该文件。 clear; global mu;mu=7; y0=[1;0]; [t,x]=ode45(verderpol,[0,40],y0); plot(t,x(:,1),r-); 解: 令 y1=x，y2=y1’ 1、建立m-文件vdp1000.m如下： function dy=vdp1000(t,y) dy=zeros(2,1); dy(1)=y(2); dy(2)=1000*(1-y(1)^2)*y(2)-y(1); 2、取t0=0，tf=3000，输入命令： [T,Y]=ode15s(vdp1000,[0 3000],[2 0]); plot(T,Y(:,1),-) 3、结果如图 To Matlab（ff4） 解 1、建立m-文件rigid.m如下： function dy=rigid(t,y) dy=zeros(3,1); dy(1)=y(2)*y(3); dy(2)=-y(1)*y(3); dy(3)=-0.51*y(1)*y(2); 2、取t0=0，tf=12，输入命令： [T,Y]=ode45(rigid,[0 12],[0 1 1]); plot(T,Y(:,1),-,T,Y(:,2),*,T,Y(:,3),+) 3、结果如图 To Matlab（ff5） 图中，y1的图形为实线，y2的图形为“*”线，y3的图形为“+”线. 返 回 Matlab 函数 求解析解（通解或特解），用 dsolve 求数值解（特解），用 ode45、ode23 ... Matlab 编程 Euler 折线法 Runga-Kutta 方法 1、目标跟踪问题一：导弹追踪问题 2、目标跟踪问题二：慢跑者与狗 3、地中海鲨鱼问题 返 回 数学建模实例 * * 实验目的 实验内容 2、掌握用Matlab求微分方程的数值解. 1、掌握用Matlab求简单微分方程的解析解. 1、求简单微分方程的解析解. 4、练习与思考. 2、求微分方程的数值解. 3、数学建模实例. 3、掌握利用微分方程模型解决实际问题. 自牛顿发明微积分以来，微分方程在描述事物运动规律上已发挥了重要的作用。实际应用问题通过数学建模所得到的方程，绝大多数是微分方程。 由于实际应用的需要，人们必须求解微分方程。然而能够求得解析解的微分方程十分有限，绝大多数微分方程需要利用数值方法来近似求解。 本实验主要研究如何用 Matlab 来计算微分方程（组）的数值解，并重点介绍一个求解微分方程的基本数值解法－－Euler折线法。 用 Maltab自带函数 解初值问题 求解析解：dsolve 求数值解：  ode45、ode23、 ode113、ode23t、ode15s、 ode23s、ode23tb dsolve 的使用 y=dsolve(eq1,eq2, ... ,cond1,cond2, ... ,v) 其中 y 为输出， eq1、eq2、...为微分方程，cond1、cond2、...为初值条件，v 为自变量。 例 1：求微分方程 的通解，并验证