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

实验七 积分方程的数值计算方法 1. 实验描述 MACROBUTTON MTEditEquationSection2 SEQ MTEqn \r \h SEQ MTSec \h 计算定积分的近似值，起始容差 1.用组合梯形公式M=10计算。 2.用组合辛普生公式M=5计算。 3.用龙贝格积分计算。 4.用自适应积分方法计算。 2. 实验内容 1. 用组合梯形公式M=10计算。 Y N 子区间宽度输入区间开始结束 Y N 子区间宽度 输入区间 开始 结束 图1. 组合梯形算法流程图 将积分区间划分为宽度为的M个子区间，再将各区间的面积求和即可得到近似面积。 2.用组合辛普生公式M=5计算。 输出SY N 结束子区间宽度输入区间开始 输出S Y N 结束 子区间宽度 输入区间 开始 图2. 组合辛普森算法流程图 将积分区间划分为宽度为的2M个子区间，再由辛普森公式将各区间的面积求和即可得到近似面积。 3.用龙贝格积分计算。 结束开始输出面积近似值输入区间 结束 开始 输出面积近似值 输入区间 图3.龙贝格积分算法流程图 = 1 \* GB3 ①.由递归梯形公式序列得到递归辛普森序列序列。 = 2 \* GB3 ②.由递归辛普森序列序列得到递归布尔公式序列。 = 3 \* GB3 ③.通过理查森改进提高误差项的阶数。 4.用自适应积分方法计算 结束 结束 将辛普森公式中的区间划分为两个相等的子区间 开始 图4.自适应积分算法流程图 在辛普森公式基础上，将区间再进行划分，即为自适应积分。 3. 实验结果及分析 真实值S = 0.19971466216144 1. 用组合梯形公式M=10计算。 面积近似值S1 =0.16965032127666。误差error1=0.479。 2.用组合辛普森公式M=5计算 面积近似值S2 =0.19966805529718。误差error2=。 3.用龙贝格积分计算。 表1.龙贝格积分表 J 梯形公式 辛普森公式 布尔公式 0 -0.265 0 0 0 0 0 1 -0.023 -0.796 0 0 0 0 2 0.127 0.640 0.736 0 0 0 3 0.742 0.19817040326280 0.256 0.21227063376585 0 0 4 0.054 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%真实值和近似值的误差 2. 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%真实值和近似值的误差 3. 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