《计算方法》课内实验报告(实验2)2014

《《计算方法》课内实验报告(实验2)2014》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、《计算方法》课内实验报告学生姓名:张靖2012309010111及学号：学院:理学院班级:信计121课程名称：计算方法实验题目:插值法与函数逼近指导教师姓名及职称：周硕教授朱振菊实验师2014年11月03日目录一、实验题目1二、实验目的1三、实验内容1四、实验结果2五、实验体会或遇到问题8一、实验题目1．熟悉matlab编写及运行数值计算程序的方法。2．进一步理解数值积分的基础理论。3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。二、实验目的1．熟悉matlab编写及运行数值计算程序的方法。2．进一步理解插值法及函数逼近方法

2、的理论基础。3．进一步掌握给定数据后应用插值法及函数逼近方法进行数据处理并给出图示结果的实际操作过程。三、实验内容1．分别用复合梯形求积公式及复合辛普森求积公式计算积分,要求计算精度达到,给出计算结果并比较两种方法的计算节点数.2．用龙贝格求积方法计算积分,使误差不超过.3．用的高斯-勒让德公式计算积分,给出计算结果.4．用辛普森公式(取)计算二重积分四、实验结果1．问题1：计算结果如下表-6-表1　问题1求解表求解值节点数复合梯形求积公式-0.444423424095复合辛普森求积公式-0.4444311310001精确解-0.44444444复合梯

3、形求积公式：取，为迭代次数，当迭代12次后，精度达到，节点数为；复合辛普森求积公式：取，为区间数，取精度为，节点数为。2．问题2：计算结果如下表表2　问题2求解表求解值龙贝格数值积分10.20759362精确解10.20759220龙贝格数值积分：给定被积函数，积分下限0，被积上限3，精度为，龙贝格积分表中行的最大数目13，计算出龙贝格数值积分近似解为10.20759362。3．问题3：计算结果如下表表3　问题3求解表求解值高斯-勒让德积分公式10.94840256精确解10.95017031-6-高斯-勒让德积分公式：取，节点横坐标取，系数取，阶导数

4、,求得高斯-勒让德积分近似解为10.94840256。4．问题4：计算结果如下表表4　问题4求解表求解值辛普森公式0.25525193复合辛普森公式（计算二重积分）：确定矩形域,求得复合辛普森公式的近似解为0.25525193。五、实验体会或遇到问题通过Matlab程序可以更为简单的运用上课所学到的，数值积分。-6-附录问题1：experiment2_1.m:clc;clear;%复合梯形求积公式[T,x1]=rctrap(@fun,eps,1,12);%T近似值，x1节点12（迭代次数）symst;fi=int(sqrt(t)*log(t),t,0,

5、1);Fs=double(fi);%精确值wT=double(abs(fi-T));%相对误差%复合辛普森求积公式Q=comsimpson(@fun,eps,1,10000);fi=int(sqrt(t)*log(t),t,0,1);Fs=double(fi);%精确解wQ=double(abs(fi-Q));%绝对误差Fs%精确解T(length(T))%复合梯形求积公式近似解Q%复合辛普森公式求积公式近似解rctrap.m:function[T,xx]=rctrap(fun,a,b,m)n=1;h=b-a;T=zeros(1,m+1);T(1)=h

6、*(feval(fun,a)+feval(fun,b))/2;xx=[];fori=1:mh=h/2;n=2*n;s=0;fork=1:n/2x=a+h*(2*k-1);s=s+feval(fun,x);xx=[xx,x];endT(i+1)=T(i)/2+h*s;endT=T(1:m);fun.m:functiony=fun(x)y=sqrt(x)*log(x);%第1题函数end问题2experiment2_2.m:clc;clear;F=inline('x*sqrt(1+x.^2)');[RT,R,err,h]=romberg(F,0,3,1.e

7、-5,13);%RT，龙贝格积分表；R，近似值；err，误差估计；h，最小步长symsx;-6-fi=int(x*sqrt(1+x.^2),x,0,3);Fs=double(fi);%精确值wR=double(abs(fi-R));%绝对误差wR1=wR-err;Fs%精确解R%龙贝格积分近似解romberg.m:function[RT,R,wugu,h]=romberg(fun,a,b,wucha,m)n=1;h=b-a;wugu=1;x=a;k=0;RT=zeros(4,4);RT(1,1)=h*(feval(fun,a)+feval(fun,b)

8、)/2;while((wugu>wucha)&(k

9、(k<4))k=k+1;h=h/2;