计算方法第二次上机作业..docxVIP

计算方法第二次上机作业Gegebao摘要：程序基于MATLAB，包括问题陈述、算法与程序、结果与分析、讨论四个部分。一 、问题陈述 数学上已经证明了：成立，所以可以通过积分来计算的近似值。分别使用矩形、梯形和Simpson复合求积公式计算的近似值。选择不同的h，对于每种求积公式，试将误差刻画成h的函数，并比较各方面的精度。是否存在某个h值，当低于这个值之后，再继续减少h的值，计算不再有所改进？为什么？实现Romberg求积方法，并重复上面的计算。使用自适应求积方法重复上面的计算。算法与程序1．矩形求积法考虑到对于具体的某一个h，不一定能整数个地覆盖积分空间，将会严重地影响算法精度。于是我们用n代替h，即将积分区间划分为n个区间，h=1/n。函数将输出积分结果Rec和这个结果对于精确的的误差的对数值r。function [Rec r ] = rec( n )h=1/n; %求出对应的hRec=0;for i=1:n x0=h*(i-1/2); %求出每一区间的中心点的横坐标 Rec=Rec+h*4/(1+x0^2);endr=log10(abs(Rec-pi)); %取其误差的对数end2.梯形求积法 和矩形法一样，采用n来标度取点密度。同样输出积分结果Lad和这个结果对于精确的的误差的对数值r。function [Lad r ] = lad( n) h=1/n;Lad=0;