积分方程的值计算方法.docVIP

实验七 积分方程的数值计算方法 1. 实验描述 计算定积分的近似值，起始容差 1.用组合梯形公式M=10计算。 2.用组合辛普生公式M=5计算。 3.用龙贝格积分计算。 4.用自适应积分方法计算。 2. 实验内容 1. 用组合梯形公式M=10计算。 图1. 组合梯形算法流程图 将积分区间划分为宽度为的M个子区间，再将各区间的面积求和即可得到近似面积。 2.用组合辛普生公式M=5计算。 图2. 组合辛普森算法流程图 将积分区间划分为宽度为的2M个子区间，再由辛普森公式将各区间的面积求和即可得到近似面积。 3.用龙贝格积分计算。 图3.龙贝格积分算法流程图 ①.由递归梯形公式序列得到递归辛普森序列序列。 ②.由递归辛普森序列序列得到递归布尔公式序列。 ③.通过理查森改进提高误差项的阶数。 4.用自适应积分方法计算 图4.自适应积分算法流程图 在辛普森公式基础上，将区间再进行划分，即为自适应积分。 3. 实验结果及分析 真实值S = 0.19971466216144 1. 用组合梯形公式M=10计算。 面积近似值S1 =0.16965032127666。误差error1=0.03006434088479。 2.用组合辛普森公式M=5计算 面积近似值S2 =0.19966805529718。误差error2=。 3.用龙贝格积分计算。 表1.龙贝格积分表 J 梯形公式 辛普森公式 布尔公式 0 -0.00199504693265 0 0 0 0 0 1 -0.02186444123413 -0.02848757266796 0 0 0 0 2 0.01611754400127 0.02877820574640 0.03259592430736 0 0 0 3 0.15265718844742 0.19817040326280 0.20946321643056 0.21227063376585 0 0 4 0.18800280139054 0.19978467237158 0.19989229031217 0.19974037084997 0.19969123256403 0 5 0.19679137798386 0.19972090351497 0.19971665225786 0.19971386435224 0.19971376040519 0.19971378242654 面积近似值S3 = 0.19971378242654。误差error3=。 4.用自适应积分方法计算 面积近似值S4 = 0.19971391278871。误差error4=。 4. 结论 自适应计分方法和龙贝格积分法都可以得到较高的精确度。 组合梯形和组合辛普森可以增大M值，得到较高的精确度。 附件（代码） 1. f=inline(sin(4*x)*exp(-2*x),x);%定义函数 A=0;t=0;S=0; for k=1:9 t=k*0.3; A=A+0.3*feval(f,t); end S=A+0.15*(feval(f,0)+feval(f,3))%组合梯形的近似值 RS=int(sin(4*x)*exp(-2*x),0,3)%真实值 error=RS-S%真实值与近似值的误差 f=inline(sin(4*x)*exp(-2*x),x);%定义函数 S=0;X=zeros(1,11);h=0.3; for j=1:11 X(j)=(j-1)*h; end for k=1:5 S=S+h/3*(feval(f,X(2*k-1))+4*feval(f,X(2*k))+feval(f,X(2*k+1))); end RS=int(sin(4*x)*exp(-2*x),0,3);%真实值 S error=RS-S%真实值与近似值的误差 function [R,quad,err,h]=romber(f,a,b,n,tol) %INPUT-f is the intergrand input as astring f % -a and b are upper and lower limits of interation % -n is the maxium number of rows in the table % -tol is the tolerance %OUTPT-R is the Romberg table % -quad is the quadrature value % -err is the error estimate % -h is the smallest step size uesd M=1; h=b-a;