龙贝格积分的程序实现

1、计算方法实验报告3【课题名称】龙贝格积分的程序实现【目的和意义】函数变化有急有缓，为了照顾变化剧烈部分的误差，我们需要加密格点。对于变化缓慢的部分， 加密格点会造成计算的浪费。以此我们介绍一种算法，可以自动在变化剧烈的地方加密格点计算，而变化缓慢的地方，则取稀疏的格点。实际计算中， 由于要事先给出一个合适的步长往往很困难，所以我们往往采用变步长的计算方案， 即在步长逐步分半的过程中， 反复利用复化求积公式进行计算， 直到所求得的积分值满足精度要求为止。在步长逐步分半过程中将粗糙的积分值逐步加工为精度较高的积分值，或者说将收敛缓慢的梯形值序列加工成收敛迅速的积分值序列。这种加速方法称为龙贝格算法

2、。【计算公式】设表示复化梯形求得的积分值，其下标是等分数，由此则有递推公式其中T2 n1Tnh22n1ba , x1)hf ( x1) ，其中 h1a ( kk0k2nk22由复化梯形公式的截断误差公式可得( ba )2I ( f ) Tn( f )( b a ) h 2f ' ' () ， I ( f )hT2n( f )12f ' ' ( )。122由此可知， I4 T2n1 Tn 。这样导出的加速公式是辛普森公式：S4T21T33n3n3n同理可得 Cn16 S2n1 SnRn64 C2n _ 1 Cn 。由此便可得加速的算法：15156363龙贝格算法。

3、【龙贝格求积算法流程图】定义被积函数f ，积分上下限 a， b 和精度 c定义一个15×4的零矩阵，用于存放t 值k>=15【龙贝格求积算法Matlab 主程序】functiont=rbg(f,a,b,c)t=zeros(15,4);t(1,1)=(b-a)/2*(f(a)+f(b);%定义龙贝格积分函数，f 为待积函数， a 与为积分上下限，c 为精度控制；生成一零矩阵，用于存放t 值；由于矩阵行列值均从1 开始，所以将原本b的 t(0,0) 记为 t(1,1), 行列均加 1；for k=2:4%先算出第一列的4 个（包括t(1,1)）值，以便程后面可以直接循环计算；sum

4、=0;for i=1:2(k-2)sum=sum+f(a+(2*i-1)*(b-a)/2(k-1);endt(k,1)=*t(k-1,1)+(b-a)/2(k-1)*sum;for i=2:kt(k,i)=(4(i-1)*t(k,i-1)-t(k-1,i-1)/(4(i-1)-1);endendfor k=5:15%循环按照公式计算出t 值；sum=0;for i=1:2(k-2)sum=sum+f(a+(2*i-1)*(b-a)/2(k-1);endt(k,1)=*t(k-1,1)+(b-a)/2(k-1)*sum;for i=2:4t(k,i)=(4(i-1)*t(k,i-1)-t(k-1,i-1)/(4(i-1)-1);endif k>6%可知最小二分次数，防止假收敛；if abs(t(k,4)-t(k-1,4)<c%若此时 t 值满足精度，则输出积分值；disp('答案 ',num2str(t(k,4);break;endendendif k>=15disp('不收敛 ');%二分次数达