数学建模实验报告求微分方程的解..docVIP

求微分方程的解 一、实验目的及意义 归纳和学习求解常微分方程(组)的基本原理和方法； 掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析； 熟悉MATLAB软件关于微分方程求解的各种命令。 二、实验内容 微分方程及方程组的解析求解法； 微分方程及方程组的数值求解法——欧拉、欧拉改进算法； 直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)； 利用图形对解的特征作定性分析。 三、实验步骤 开启软件平台——MATLAB，开启MATLAB编辑窗口； 根据微分方程求解步骤编写M文件 保存文件并运行； 观察运行结果(数值或图形)； 根据观察到的结果和体会写出实验报告。 四、实验要求与任务 根据实验内容和步骤，完成以下实验，要求写出实验报告 1. 求微分方程的通解． 2. 求微分方程的通解． 3. 求微分方程组 在初始条件下的特解，并画出解函数的图形． 4. 分别用 ode23、ode45 求上述第 3 题中的微分方程初值问题的数值解(近似解)，求解区间为．利用画图来比较两种求解器之间的差异． 5. 用 Euler 折线法求解微分方程初值问题 的数值解(步长h取0.1)，求解范围为区间[0,2]． 6. 用四阶 Runge-Kutta 法求解微分方程初值问题 的数值解(步长h取0.1)，求解范围为区间[0,3]． 五. 程序代码及运行结果（经调试后正确的源程序） 1. 求微分方程的通解． 程序代码： syms x y fprintf(通解为) y=dsolve((x^2-1)*Dy+2*x*y-sin(x)=0,x) 运行结果： 通解为 y = (-cos(x)+C1)/(x^2-1) 2. 求微分方程的通解． 程序代码： syms x y fprintf(通解为) y=dsolve(D2y-2*Dy+5*y=exp(x)*sin(x),x) 运行结果： 通解为 y = -1/4*exp(x)*cos(2*x)*sin(x)+1/12*exp(x)*cos(2*x)*sin(3*x)-1/12*exp(x)*sin(2*x)*cos(3*x)+1/4*exp(x)*sin(2*x)*cos(x)+C1*exp(x)*cos(2*x)+C2*exp(x)*sin(2*x) 3. 求微分方程组 在初始条件下的特解，并画出解函数的图形． 程序代码： syms x y t [x,y]=dsolve(Dx+x+y=0,Dy+x-y=0,x(0)=1,y(0)=0,t) ezplot(x,y,[0,2]); 运行结果： x = 1/2*exp(2^(1/2)*t)+1/4*2^(1/2)*exp(-2^(1/2)*t)-1/4*2^(1/2)*exp(2^(1/2)*t)+1/2*exp(-2^(1/2)*t) y = 1/4*2^(1/2)*exp(-2^(1/2)*t)-1/4*2^(1/2)*exp(2^(1/2)*t) 4. 分别用 ode23、ode45 求上述第 3 题中的微分方程初值问题的数值解(近似解)，求解区间为．利用画图来比较两种求解器之间的差异． 程序代码： M函数文件verderpol.m: function xprime=verderpol(t,x) xprime=[-x(1)-x(2); x(2)-x(1)]; 在程序中调用此函数： clear; y0=[1;0]; [t,x]=ode45(verderpol,[0,2],y0); plot(x(:,1),x(:,2),r-); hold on clear; y0=[1;0]; [t,x]=ode23(verderpol,[0,2],y0); plot(x(:,1),x(:,2),b-); 运行结果： 5. 用 Euler 折线法求解微分方程初值问题 的数值解(步长h取0.1)，求解范围为区间[0,2]． 程序代码： clear f=sym(y-(12*x^2)/y^3); a=0; b=2; h=0.1; n=(b-a)/h+1; x=0; y=1; szj=[x,y]; for i=1:n-1 y=y+h*subs(f,{x,y},{x,y}); x=x+h; szj=[szj;x,y]; end szj plot(szj(:,1),szj(:,2),or-) 运行结果： szj = 0 1.0000 0.1000 1.1000 0.2000 1.2010 0.3000 1.2934 0.4000 1.3728 0.5000