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

《计算方法》课内实验报告 学生姓名: 张 靖 2012309010111 及 学 号： 学 院: 理学院 班 级: 信计121 课程名称： 计算方法 实验题目: 插值法与函数逼近 指导教师 姓名及职称： 周 硕 教 授 朱振菊 实验师 2014 目 录 一、实验题目 1 二、实验目的 1 三、实验内容 1 四、实验结果 2 五、实验体会或遇到问题 8 一、实验题目 1．熟悉matlab编写及运行数值计算程序的方法。 2．进一步理解数值积分的基础理论。 3．进一步掌握应用不同的数值积分方法求解给定的积分并给出数据结果及误差分析。 二、实验目的 1．熟悉matlab编写及运行数值计算程序的方法。 2．进一步理解插值法及函数逼近方法的理论基础。 3．进一步掌握给定数据后应用插值法及函数逼近方法进行数据处理并给出图示结果的实际操作过程。 三、实验内容 1．分别用复合梯形求积公式及复合辛普森求积公式计算积分,要求计算精度达到,给出计算结果并比较两种方法的计算节点数. 2．用龙贝格求积方法计算积分,使误差不超过. 3．用的高斯-勒让德公式计算积分,给出计算结果. 4．用辛普森公式 (取) 计算二重积分 实验结果 1．问题1： 计算结果如下表  表1　问题1求解表 求解值 节点数 复合梯形求积公式 -04095 复合辛普森求积公式 -010001 精确解 -0复合梯形求积公式：取，为迭代次数，当迭代12次后，精度达到，节点数为； 复合辛普森求积公式：取，为区间数，取精度为，节点数为。 2．问题2： 计算结果如下表 表2　问题2求解表 求解值 龙贝格数值积分 10精确解 10龙贝格数值积分：给定被积函数，积分下限0，被积上限3，精度为，龙贝格积分表中行的最大数目13，计算出龙贝格数值积分近似解为10 3．问题3： 计算结果如下表 表3　问题3求解表 求解值 高斯-勒让德积分公式 10精确解 10高斯-勒让德积分公式：取 ，节点横坐标取，系数取，阶导数,求得高斯-勒让德积分近似解为10 4．问题4： 计算结果如下表 表4　问题4求解表 求解值 辛普森公式 0复合辛普森公式（计算二重积分）：确定矩形域,求得复合辛普森公式的近似解为0 五、实验体会或遇到问题 通过Matlab程序可以更为简单的运用上课所学到的，数值积分。 附 录 问题1： experiment2_1.m: clc;clear; %复合梯形求积公式 [T,x1]=rctrap(@fun,eps,1,12); %T近似值，x1节点 12（迭代次数） syms t; fi=int(sqrt(t)*log(t),t,0,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*(feval(fun,a)+feval(fun,b))/2; xx=[]; for i=1:m h=h/2; n=2*n; s=0; for k=1:n/2 x=a+h*(2*k-1); s=s+feval(fun,x); xx=[xx,x]; end T(i+1)=T(i)/2+h*s; end T=T(1:m); fun.m: function y=fun(x) y=sqrt(x)*log(x);%第1题函数 end 问题2 experiment2_2.m: clc;clear; F=inline('x*sqrt(1+x.^2)'); [RT,R,err,h]=romberg(F,0,3,1.e-5,13); %RT，龙贝格积分表；R，近似值；err，误差估计；h，最小步长 syms x ; fi=int(x*sqrt(1+x.^2),x,0,3); Fs=double(fi);%精确值 wR=double(abs(fi-R));%绝对