[数学]实验三 微分方程模型.doc

微分方程模型 一、实验目的及意义 掌握解析数值解法，并学会用图形观察解的形态和进行解的定性分析； 熟悉MATLAB软件关于微分方程求解的各种命令； 通过范例学习建立微分方程方面的数学模型以及求解全过程； 二、实验内容 直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)； 利用图形对解的特征作定性分析； 建立微分方程方面的数学模型，并了解建立数学模型的全过程。 三、实验步骤 开启软件平台——MATLAB，开启MATLAB编辑窗口； 根据微分方程求解步骤编写M文件 保存文件并运行； 观察运行结果(数值或图形)； 根据观察到的结果和体会写出实验报告。 四、实验要求与任务 　　根据实验内容和步骤，完成以下实验，要求写出实验报告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论） 求微分方程的解析解, 并画出它们的图形 （1）y = y + 2 x , y (0) = 1, 0 x 1; 程序： y=dsolve(Dy=y+2*x,y(0)=1,x) y = -2*x-2+3*exp(x) x=0:pi/100:1; y=-2*x-2+3*exp(x); 画出该微分方程的解析解plot(x,y,r-) 图见： （2）y + y cos( x ) = 0, y (0)=1, y (0)=0; 程序一： y=dsolve(D2y+y*cos(x)=0,y(0)=1,Dy(0)=0,x) ??? Error using == dsolve Error, (in dsolve/IC) The implicit option is not available when giving Initial Conditions. 程序二： 选择状态变量y1=y’,y2=y,则可写出题中的方程的状态方程形式为： y1’=-y2*cos(x); y2’=y1; 建立函数文件vdpol.m,程序如下： function ydot=vdpol(x,y) ydot(1)=-y(2)*cos(x); ydot(2)=y(1); ydot=ydot; 求解微分方程： [x,y]=ode45(@vdpol,[0,5],[0;1]); [x,y] ans = 0 0 1.0000 0.0001 -0.0001 1.0000 0.0001 -0.0001 1.0000 0.0002 -0.0002 1.0000 0.0002 -0.0002 1.0000 0.0005 -0.0005 1.0000 0.0007 -0.0007 1.0000 0.0010 -0.0010 1.0000 0.0012 -0.0012 1.0000 0.0025 -0.0025 1.0000 0.0037 -0.0037 1.0000 0.0050 -0.0050 1.0000 0.0062 -0.0062 1.0000 0.0125 -0.0125 0.9999 0.0188 -0.0188 0.9998 0.0251 -0.0251 0.9997 0.0313 -0.0313 0.9995 0.0627 -0.0627 0.9980 0.0941 -0.0939 0.9956 0.1255 -0.1249 0.9921 0.1569 -0.1557 0.9877 0.2819 -0.2746 0.9608 0.4069 -0.3853 0.9194 0.5319 -0.4848 0.8649 0.6569 -0.5711 0.7988 0.7819 -0.6427 0.7227 0.9069 -0.6993 0.6387 1.0319 -0.7414 0.5485 1.1569 -0.7702 0.4539 1.2819 -0.7878 0.3565 1.4069 -0.7965 0.2574 1.5319 -0.7993 0.1576 1.6569 -0.7991 0.0577 1.7819