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

《计算方法》课内实验报告 学生姓名: 张 靖 2012309010111 及 学 号： 学 院: 理学院 班 级: 信计121 课程名称： 计算方法 实验题目: 插值法与函数逼近 指导教师 姓名及职称： 周 硕 教 授 朱振菊 实验师 2014年11月03日 目 录 一、实验题目 1 二、实验目的 1 三、实验内容 1 四、实验结果 2 五、实验体会或遇到问题 8 一、实验题目 1．2．3．1．2．3．1．,要求计算精度达到,给出计算结果并比较两种方法的计算节点数. 2．,使误差不超过. 3．的高斯-勒让德公式计算积分,给出计算结果. 4．) 计算二重积分 实验结果 1．问题1： 计算结果如下表  问题复合梯形求积公式-04095 复合辛普森求积公式-010001 精确解-0 复合梯形求积公式：取为迭代次数，当迭代后，，节点数为 复合辛普森求积公式取为区间数取精度为节点数为2．问题2： 计算结果如下表 2　问题210 精确解10 龙贝格数值积分：给定被积函数，积分下限，龙贝格积分表中行的最大数目362。 3．问题3： 计算结果如下表 3　问题310 精确解10 高斯勒让德积分公式 ，节点横坐标取系数取阶导数,求得高斯积分近似解为4．问题4： 计算结果如下表 4　问题40 复合辛普森公式二重积分）矩形域,复合辛普森公式的近似解为 通过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));%绝对误差 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))/2; while((wuguwucha)(km)|(k4)) k=k+1; h=h/2; s=0; for j=1:n x=a+h*(2*j-1); s=s+feval(fun,x); end RT(k+1,1)= RT(k,1)/2+h*s; n=2*n; for i=1: