复化梯形积分报告.docVIP

复化梯形积分公式 实验的目的及意义： 前边考虑到数值计算的稳定性，用增大n的方法来提高数值积分代数精度的方法是不可取的。类似于分段差值，为了减少数值积分的误差，可以把积分区间分成若干小区间，在每个小区间上采取低阶数值积分公式，然后把这些小区间上的数值积分结果加起来作为函数在整个区间上的近似，这就是复化数值积分。 将区间[a,b]n等分，取等距节点 由定积分的区间可加性，有 在每一个小区间上利用梯形积分公式有 一般记称做n+1点复化梯形积分公式。 数学公式： 算法描述： Step1:输入a,b和正整数n； Step2:置h=(b-a)/n; Step3:F=f(a)+f(b);l=0; Step4:对j=1,2,…,n循环执行5； Step5:置x=a+jh; l+=f(x); Step6:置T=h(F+2l)/2 Step7:输出T; 程序原代码如下： #include stdio.h double f(double x) { double l; l=4/(1+x*x); return l; } void main() { int n,j; double a,b,h,F,T,l,x; printf(请输入端点值a和b:); scanf(%lf%lf,a,b); printf(请输入n:); scanf(%d,n); h=(b-a)/n; F=f(a)+f(b); l=0; for(j=1;jn;j++) { x=a+j*h; l+=f(x); } T=h*(F+2*l)/2; printf(T=%lf\n,T); } 数值计算: 计算π； 请输入端点值a和b:0 1 请输入n：100 T=3.141576 Press any key to continue 对计算结果进行评价分析： 其误差可以由求出。此式说明复化梯形公式是收敛的。 参考文献： 秦新强，数值逼近，西安：西安理工大学，2010