资源描述：
《20115394王福临数学实验作业(微分方程)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、重庆大学学生实验报告实验课程名称数学实验开课实验室DS1401学院计算机年级2011专业班计科6班学生姓名王福临学号20115394开课时间至学年第学期总成绩教师签名开课学院、实验室：ds1401实验时间：2013/4/7课程名称数学实验实验项目名称微分方程实验项目类型验证演示综合设计其他指导教师何光辉成绩√一、实验目的及意义[1]归纳和学习求解常微分方程(组)的基本原理和方法；[2]掌握解析、数值解法，并学会用图形观察解的形态和进行解的定性分析；[3]熟悉MATLAB软件关于微分方程求解的各种命令；[4]通过范例学习建立微分方程方面的数学模型以及求解
2、全过程；通过该实验的学习，使学生掌握微分方程(组)求解方法（解析法、欧拉法、梯度法、改进欧拉法等），对常微分方程的数值解法有一个初步了解，同时学会使用MATLAB软件求解微分方程的基本命令，学会建立微分方程方面的数学模型。这对于学生深入理解微分、积分的数学概念，掌握数学的分析思维方法，熟悉处理大量的工程计算问题的方法是十分必要的。二、实验内容1．微分方程及方程组的解析求解法；2．微分方程及方程组的数值求解法——欧拉、欧拉改进算法；3．直接使用MATLAB命令对微分方程(组)进行求解(包括解析解、数值解)；4．利用图形对解的特征作定性分析；5．建立微分方
3、程方面的数学模型，并了解建立数学模型的全过程。三、实验步骤1．开启软件平台——MATLAB，开启MATLAB编辑窗口；2．根据微分方程求解步骤编写M文件3．保存文件并运行；4．观察运行结果(数值或图形)；5．根据观察到的结果和体会写出实验报告。四、实验要求与任务根据实验内容和步骤，完成以下实验，要求写出实验报告（实验目的→问题→数学模型→算法与编程→计算结果→分析、检验和结论）基础实验1、求微分方程的解析解,并画出它们的图形，1)y(4)=y,y(0)=y’(0)=2,y’’(0)=y’’’(0)=1代码如下：y=dsolve('D4y=y','y(0
4、)=2,Dy(0)=2,D2y(0)=1,D3y(0)=1','x')y=3/2*exp(x)+1/2*sin(x)+1/2*cos(x)ezplot('3/2*exp(x)+1/2*sin(x)+1/2*cos(x)',[-5,5]);2)x3y’’’+x2y’’-4xy’=3x2y(0.1)=y’(0.1)=2,y’’(0.1)=1,代码如下：y=dsolve('(x^3)*D3y+(x^2)*D2y-4*x*Dy=3*x^2','y(0.1)=2,Dy(0.1)=2,D2y(0.1)=1')y=1/68*(11*17^(1/2)+51)*(x+1
5、0-2*17^(1/2))/exp(1/20*(-1+17^(1/2))/x)/(-1+17^(1/2))*x*exp(1/2*(-1+17^(1/2))/x*t)+1/68*(11*17^(1/2)-51)*(2*17^(1/2)+x+10)/exp(-1/20*(1+17^(1/2))/x)/(1+17^(1/2))*x*exp(-1/2*(1+17^(1/2))/x*t)-3/4*x*t+2-7/16*x^2-17/40*xezplot('1/68*(11*17^(1/2)+51)*(x+10-2*17^(1/2))/exp(1/20*(-1+1
6、7^(1/2))/x)/(-1+17^(1/2))*x*exp(1/2*(-1+17^(1/2))/x*t)+1/68*(11*17^(1/2)-51)*(2*17^(1/2)+x+10)/exp(-1/20*(1+17^(1/2))/x)/(1+17^(1/2))*x*exp(-1/2*(1+17^(1/2))/x*t)-3/4*x*t+2-7/16*x^2-17/40*x',[-5,5])2、用向前欧拉公式和改进的欧拉公式求方程y’=y-2x/y,y(0)=1的数值解（0≤x≤1,h=0.1）要求编写程序。1、分析：解：（1）解析解法得到其精确解：
7、（2）向前欧拉法：迭代公式为 ，其中（3）改进欧拉法：迭代公式为，其中2、Matlab编码x1(1)=0;y1(1)=1;y2(1)=1;h=0.1;fork=1:10x1(k+1)=x1(k)+h;y1(k+1)=(1-h)*y1(k)+2*h*x1(k)/y1(k);y2(k+1)=(1+h)*y2(k)+(h*h)/2-h*x1(k)/y2(k)-h*(x1(k)+h)/(y2(k)+h);endx=0:0.1:1;y=(2*x+1).^(1/2);x1=x1(1:11),y=y(1:11),y1=y1(1:11),y2=y2(1:11),plo
8、t(x,y,x1,y1,'k:',x1,y2,'r--')显示图像及结果：x1=00.1000
