微分计算方法

1、 实 验 报 告课程名称： 计算方法 院 系： 数学科学系 专业班级： 数应1001 学 号： 学生姓名： 曹信信 指导教师： 沈 林 开课时间：2012至2013学年第一学期 一、学生撰写要求按照实验课程培养方案的要求，每门实验课程中的每一个实验项目完成后，每位参加实验的学生均须在实验教师规定的时间内独立完成一份实验报告，不得抄袭，不得缺交。学生撰写实验报告时应严格按照本实验报告规定的内容和要求填写。字迹工整，文字简练，数据齐全，图表规范，计算正确，分析充分、具体、定量。二、教师评阅与装订要求1.实验报告批改要深入细致，批改过程中要发现和纠正学生实验报告中的问题，给出评语和实验报告成绩，签名

2、并注明批改日期。实验报告批改完成后，应采用适当的形式将学生实验报告中存在的问题及时反馈给学生。2.实验报告成绩用百分制评定，并给出成绩评定的依据或评分标准（附于实验报告成绩登记表后）。对迟交实验报告的学生要酌情扣分，对缺交和抄袭实验报告的学生应及时批评教育，并对该次实验报告的分数以零分处理。对单独设课的实验课程，如学生抄袭或缺交实验报告达该课程全学期实验报告总次数三分之一以上，不得同意其参加本课程的考核。3.各实验项目的实验报告成绩登记在实验报告成绩登记表中。本学期实验项目全部完成后，给定实验报告综合成绩。4.实验报告综合成绩应按课程教学大纲规定比例（一般为10-15%）计入实验课总评成绩；实

3、验总评成绩原则上应包括考勤、实验报告、考核（操作、理论）等多方面成绩；5.实验教师每学期负责对拟存档的学生实验报告按课程、学生收齐并装订，按如下顺序装订成册：实验报告封面、实验报告成绩登记表、实验报告成绩评定依据、实验报告（按教学进度表规定的实验项目顺序排序）。装订时统一靠左侧按“两钉三等分”原则装订。实验名称MATLAB基本操作实验时间 2012年9月13日学生姓名曹信信实验地点9#405数学实验室1、实验所用软件WIN7操作系统、Matlab2、实验目的熟悉Matlab编程环境，理解和掌握Matlab执行命令的方式，会做一些简单函数的图像。通过用Matlab编程解决数值分析问题，更深一步的

4、体会计算方法这门课的重要性。3、实验内容（一）、Matlab操作界面1. 命令窗口（command window）2. 命令历史窗口（command history）3. 工作空间管理窗口（workspace）4. 当前路径窗口（current directory）（二）、具体练习1、简单矩阵的输入步骤。2、画出衰减振荡曲线及其它的包络线。的取值范围是。3、画出所表示的三维曲面。的取值范围是。4、复数矩阵的生成及运算.4、实验方法、步骤 1、简单矩阵的输入步骤。（1）在键盘上输入下列内容A = 2,6,3; 4,7,5; 3,8,9 （2）按【Enter】键，指令被执行。（3）在指令执行后，M

5、ATLAB指令窗中将显示结果2、画出衰减振荡曲线及其它的包络线。的取值范围是。t=0:pi/50:4*pi;y0=exp(t/3);y=exp(t/3).*sin(2*t);plot(t,y,-r,t,y0,:b,t,-y0,:b)grid 3、画出所表示的三维曲面。的取值范围是。clear;x=-8:0.5:8;y=x;X=ones(size(y)*x;Y=y*ones(size(x);R=sqrt(X.2+Y.2)+eps;%Z=cos(R)./R;%surf(X,Y,Z);colormap(cool)xlabel(x),ylabel(y),zlabel(z)4、复数矩阵的生成及运算A=5

6、,3;2,8-7,8;3,9*iB=1+4i,2+5i;8+3*i,4+2*i C=A*B 5、实验数据记录与分析1、 2、3、4、 6、实验结论通过本节课实验，我初步熟悉Matlab的编程环境，认识了其中的命令窗口、命令历史窗口、工作空间管理窗口以及当前路径窗口，并学习了简单的运算操作。在这个过程中，我掌握了用Matlab执行命令做一些简单函数的图像。虽然是初步接触，但是我觉得自己做得还不错，下次继续努力！指导教师评语和成绩评定 指导教师签字： 年 月 日实验名称插值算法实验时间2012年9月27日学生姓名曹信信实验地点9#405数学实验室1、实验所用软件WINxp操作系统、Matlab2、

7、实验目的设计一个程序，通过键盘指令输入参数，可以实现Taylor插值方法、线性插值方法、抛物插值方法以及一般情形的Lagrange插值方法的程序。使学生掌握插值算法的基本理论，并提高学生在程序设计中发现问题、分析问题和解决问题的能力。3、实验内容1、确定方程x3-x+4=0的实根的分布情况。2、已知函数值 x=0 1 2 3 4 5 6 7 8 9 10，y=0 0.79 1.53 2.19 2.71 3.03 3.27 2.89 3.06 3.19 3.29和边界条件s(0)=0.8 s(10)=0.2， 求三次样条插值函数的图形3、用牛顿迭代法求解方程组:x12+x22=4;x12-x22

8、=1。4、已知： sin0.32=0., sin0.34=0., sin0.36=0., 用抛物插值计算 sin0.3367。1、在MATLAB工作窗口输入x=-4:0.1:4;y=x.3-x +4; plot(x,y) grid,gtext(y=x3-x+4)2、输入xx=1:0.1:10;yy=spline(x,y,xx); plot(x,y)3、输入x0=1;y0=1;n=10;tol=1e-6;x(1)=x0;y(1)=y0;i=1;u=1 1;k(1)=1;while(norm(u)tol*norm(x(i),y(i)A=2*x(i),y(i);x(i),-y(i);b=4-x(i)

9、2-y(i)2,1-x(i)2+y(i)2;u=Ab;x(i+1)=x(i)+u(1);y(i+1)=y(i)+u(2);i=i+1;k(i)=i;if(in)error(n is full);endendk,x,y4、x=0.32;0.34;0.36;y=0.;0.;0.;a=0.3367;a=polyfit(x,y,2)t=0.3367;s=t2;t;1ans5、实验数据记录与分析 1、2、3、 4、6、实验结论通过本次试验操作，我进一步了解了matlab的上机操作，通过键盘指令输入参数，实现Taylor插值方法、线性插值方法、抛物插值方法以及一般情形的Lagrange插值方法的程序。在做

10、题的过程中，不免会遇到很多问题，不过经过努力和同学们一起探讨，程序一步步都成功运行了。课后我会多加练习，争取更熟练的掌握插值算法。指导教师评语和成绩评定 指导教师签字： 年 月 日实验名称数值积分与数值微分实验时间2012年10月25日学生姓名曹信信实验地点9#405数学实验室1、实验所用软件WIN7操作系统、Matlab2、实验目的1.掌握复合梯形法求积分、复合辛普生法方法求积分公式的程序设计2.掌握Romberg积分的用法。3.掌握高斯公式的程序设计3、实验内容1) 由，所围成的平面区域,用复合梯形法求平面区域的面积。2) 估计用辛普森公式计算定积分的值I和递归次数n。3) 用两点高斯勒让

11、德积分公式计算d，并将计算结果与精确值进行比较，求出绝对误差。4、实验方法、步骤1) x=-3:0.001:3; F1= sin(x); F2=cos(x);plot(x ,F1,r:,x ,F2,b-), axis(-1,pi/4+1,-1.5,1.5),xlabel(x), ylabel(y),title(y=sinx , y=cosx 和x=-1及x=0.5所围成的平面区域的图形)syms x f1= cos(x)-sin(x); f2=-f1; S1=int(f1,x,-1,pi/4); S2=int(f2,x, pi/4,0.5); S=S1+S2,Sjg= double (S)2)

12、 syms xL= inline( exp(sin(x); QS,FCNTS =quad(L,- pi/4, pi/4,1.e-4,2)3) x1=-1/sqrt(3);x2=1/sqrt(3); y1=2/(3+x1); y2=2/(3+x2);G=y1+y2syms x f=int(2/(3+x),x,-1,1); Fjg= double(fi)Gjq= double(abs(fi-G2)5、实验数据记录与分析1) S =sin(1) - cos(1/2) - cos(1) - sin(1/2) + 2*2(1/2)Sjg =1.77262）QS,FCNTS =quad(L,- pi/4,

13、 pi/4,1.e-4,2) 9 -0. 4.e-01 0. 11 -0. 7.e-01 0. 13 0. 4.e-01 0.QS = 1.7173 FCNTS = 133) G = 1.3846 Fjg =0.4055 Gjq = 5.9703e-056、实验结论 本节实验我学会了利用复合梯形法、复合辛普生法方法求积分，在学习的时候简单地进行了程序设计，锻炼了求积分的快捷方法，也学习了更深刻的定理公式。此外，我们还学习了Romberg积分的用法，练习了高斯公式的程序设计，在以前的基础上进步许多，但是问题也暴露不少，例如一些函数的用法，特殊函数的名称等，在这些方面需要进一步提升。指导教师评语和

14、成绩评定 指导教师签字： 年 月 日实验名称常微分方程差分方法实验时间2012年11月8日学生姓名曹信信实验地点9#405数学实验室1、实验所用软件WIN7操作系统、Matlab2、实验目的1.编写程序实现Euler求解方法及改进Euler求解方法2.掌握龙格-库塔方法的用法3.掌握方程组和高阶方程的程序实现3、实验内容1) 用向后欧拉公式求解区间上的初值问题，2) 判断下列线性方程组解的情况，并求出该解。4、实验方法、步骤1) 编写程序function X,Y,n,P=Heuler1(funfcn,x0,b,y0,h,tol)n=fix(b-x0)/h); X=zeros(n+1,1); Y

15、=zeros(n+1,1);k=1; X(k)=x0; Y(k,:)=y0; Y1(k,:)=y0; clc,x0,h,y0 for i=2:n+1X(i)=x0+h; Y(i,:)=y0+h*feval(funfcn,x0,y0);Y1(i,:)=y0+h*feval(funfcn,X(i),Y(i,:); Wu=abs(Y1(i,:)-Y(i,:);while Wutol p=Y1(i,:); Y1(i,:)=y0+h*feval(funfcn,X(i),p); Y(i,:)=p; endx0=x0+h; y0=Y1(i,:); Y(i,:)=y0; plot(X,Y,ro)grid on

16、xlabel(自变量X), ylabel(因变量Y)title(用向后欧拉公式计算dy/dx=f(x,y)，y(x0)=y0在x0,b上的数值解)end X=X(1:n+1); Y=Y(1:n+1,:); n=1:n+1; P=n,X,Y输入S1=dsolve(Dy=8*x-3*y-7,y(0)=1,x) x0=0;y0=1; b=2;tol=1.e-1;subplot(2,1,1)h1=0.01; X1,Y1,n,P1=Heuler1(funfcn,x0,b,y0,h1,tol) hold onS2= 8/3*X1-29/9+38/9*exp(-3*X1), plot(X1,S2,b-)le

17、gend(h=0.01用向后欧拉公式计算dy/dx=8x-3y-7，y(0)=1在0,2上的数值解)hold off2)编写程序function RA,RB,n=jiepb(A,b)B=A b;n=length(b); RA=rank(A); RB=rank(B); cha=RB-RA;if cha0,disp(RA=RB，方程组无解)returnendif RA=RB if RA=ndisp(RA=RB=n，方程组有唯一解) else disp(RA=RBn，方程组有无穷多解)endend输入A=2 5 -1 5;3 1 2 -7;4 1 -3 6;3 -2 4 -7; b= 0; 2; 0

18、; 1; RA,RB,n=jiepb(A,b) X=Ab5、实验数据记录与分析1) x0 = 0h = 0.0000y0 = 12) RA=RB=n，方程组有唯一解RA =4RB =4n =4X =0.1647 0.1511 -0.2733 -0.27166、实验结论本次实验重点学习了求解常微分方程的差分方法，在求解过程中用到了Euler求解方法及改进的Euler求解方法，并且可以用两个方法求解比较不同之处，看到改进后的明显效果。同时，我们还学习了龙格-库塔方法的用法及解方程组和高阶方程的程序实现，通过编写、运行程序自己动手操作熟练掌握matlab解方程的应用。指导教师评语和成绩评定 指导教师

19、签字： 年 月 日实验名称方程求根的迭代法实验时间2012年11月29日学生姓名曹信信实验地点9#405数学实验室1、实验所用软件WIN7操作系统、Matlab2、实验目的1.掌握二分法与牛顿法的基本原理及应用2.掌握牛顿法、弦截法的程序设计并用具体例子加以实现3、实验内容1) 确定方程的实根的分布情况，画出函数在开区间 (-4,4)内的分布情况。2) 求x5-2x+3=0在0.5附近的根，精确到4位小数.4、实验方法、步骤1)x=-4:0.1:4;y=x.2-2*x; plot(x,y)grid,gtext(y=x2-2x);2)编写程序functionk,piancha,xdpiancha

20、,xk,yk=diedai2(x0,tol,ddmax)x(1)=x0;for i=1: ddmax x(i+1)=fun(x(i);piancha=abs(x(i+1)-x(i); xdpiancha=piancha/( abs(x(i+1)+eps);i=i+1; xk=x(i);yk=fun(x(i); (i-1) piancha xdpiancha xk yk if (pianchatol)|(xdpianchaddmax disp(迭代次数超过给定的最大值ddmax) k=i-1; xk=x(i);yk=fun(x(i); (i-1) piancha xdpiancha xk yk;

21、return; endP=(i-1),piancha,xdpiancha,xk,yk;输入k,piancha,xdpiancha,xk,yk=diedai2(0.5,1e-4,15)5、实验数据记录与分析1) 2）ans = 1.0000 0.7500 0.6000 1.2500 4.0625ans = 2.0000 2.8125 0.6923 4.0625 24.6289ans = 3.0000 20.5664 0.8351 24.6289 655.8408ans = 1.0e+005 * 0.0000 0.0063 0.0000 0.0066 4.3144ans = 1.0e+011 *

22、0.0000 0.0000 0.0000 0.0000 1.8614ans = 1.0e+022 * 0.0000 0.0000 0.0000 0.0000 3.4648ans = 1.0e+045 * 0.0000 0.0000 0.0000 0.0000 1.2005ans = 1.0e+090 * 0.0000 0.0000 0.0000 0.0000 1.4412ans = 1.0e+180 * 0.0000 0.0000 0.0000 0.0000 2.0771ans = 1.0e+180 * 0.0000 2.0771 0.0000 2.0771 Infans = 11 Inf N

23、aN Inf Infans = 12 NaN NaN Inf Infans = 13 NaN NaN Inf Infans = 14 NaN NaN Inf Infans = 15 NaN NaN Inf Inf迭代次数超过给定的最大值ddmaxk = 15piancha = NaNxdpiancha = NaNxk = Infyk = Inf6、实验结论本次实验我先学习了二分法与牛顿法的基本原理，了解其原理后，应用此原理解决微分方程，并且大大提高了其精度；学习了牛顿法、弦截法的程序设计后应用到实际例子中去，在做题中锻炼应用牛顿法的能力，加强了对该公式的熟练运用，记忆起来更为方便。由于课后练习

24、试验次数不多，做起来有点生疏，需要多加练习，熟练掌握。指导教师评语和成绩评定 指导教师签字： 年 月 日实验名称线性方程组的迭代法实验时间2012年12月13日学生姓名曹信信实验地点9#405数学实验室1、实验所用软件WIN7操作系统、Matlab2、实验目的设计程序实现雅可比迭代法，高斯-赛德尔迭代法等算法。掌握方程求根的基本理论。并用实例在计算机上计算。3、实验内容(1)用判别雅可比迭代收敛性的MATLAB主程序，判别由下列方程组的雅可比迭代产生的序列是否收敛？(2) 用高斯-塞德尔迭代定义的MATLAB主程序解下列线性方程组,取初始值，要求当时，迭代终止。4、实验方法、步骤(1)A=9

25、-1 -2;-1 10 -2;-1 -1 1.5;a=jspb(A)(2)编写程序function X=gsdddy(A,b,X0,P,wucha,max1)D=diag(diag(A);U=-triu(A,1);L=-tril(A,-1); dD=det(D);if dD=0disp(对角矩阵D奇异，此方程组无解.)elsedisp(对角矩阵D非奇异，方程组有解.)iD=inv(D-L); B2=iD*U;f2=iD*b;jX=Ab; X=X0; n m=size(A);for k=1:max1X1= B2*X+f2; djwcX=norm(X1-X,P);xdwcX=djwcX/(norm

26、(X,P)+eps);if (djwcXwucha)|(xdwcXwucha) returnelse k,X1,k=k+1;X=X1;endendif (djwcXwucha)|(xdwcX0 disp(请注意：因为RA=RB，所以此方程组无解.) endif RA=RB if RA=n disp(请注意：因为RA=RB=n，所以此方程组有唯一解.) X=zeros(n,1); C=zeros(1,n+1); for p= 1:n-1 for k=p+1:n m= B(k,p)/ B(p,p);B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1); end end b=B(1:n,n+1);A=B(1: