梯形积分和3点Guss积分2.docVIP

实验六 一、实验名称 复合梯形积分和复化3点Gauss积分计算数值积分 二、实验目的与要求： 实验目的: 掌握复合梯形积分和复化3点Gauss积分算法。 实验要求：1.给出复合梯形积分和复化3点Gauss积分算法思路， 2.用C语言实现算法，运行环境为Microsoft Visual C++。 三、算法思路： 1．复合梯形积分 我们把整个积分区间[a,b]分成n个子区间[xi,xi+1]，i=0，1,2,…,n，其中x0=a，xn+1=b。这样求定积分问题就分解为求和问题： QUOTE abfxdx=1nxixi+1f 当这n+1个结点为等距结点时，即 SKIPIF 1 0 ，i=0,1,2,…,n，复化梯形公式的形式是 SKIPIF 1 0 算法： input n SKIPIF 1 0 for i=1 to n do SKIPIF 1 0 end do output S 2．复化3点Gauss积分 根据书上6.8节定理5，可以用递推法求正交多项式如下： SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 求 SKIPIF 1 0 的根，得到 SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 ，则积分可近似表示为如下形式 SKIPIF 1 0 ，将 SKIPIF 1 0 带入其中得到 SKIPIF 1 0 ，带入上式得到复化3点Gauss积分公式。 SKIPIF 1 0 分别在每个小区间上应用三点Gauss积分并求和，便得到复合三点Gauss积分法数值积分值。 算法： input n,a,b h=(a+b)/2 SKIPIF 1 0 for i=1 to n do SKIPIF 1 0 ， SKIPIF 1 0 ， SKIPIF 1 0 SKIPIF 1 0 ， SKIPIF 1 0 SKIPIF 1 0 end do output S 四、实验题目： 为了方便起见，我们采取复合梯形积分的方法取N充分大时（这里取N等于2的20次方）得到的积分值近似看作积分的真实值，因为题中的积分函数在积分闭区域上是一致连续的，所以由黎曼积分的定义知在N充分大的情况下，复化梯形积分的值充分接近原积分值，故我们这样做是合理的。 五、问题的解： 编写程序（程序见后面附录），输出结果如下： 为了便于看清数值积分结果与原函数积分实际结果的差异。我在运行程序时故意计算了一下原函数积分的近似真实结果。 分析并比较得到的数据可以看出，当k越来越大时，error1和error2越来越接近于零，数值积分的结果越来越靠近原函数积分实际结果，并且error2 比error1 更快地接近于零。误差比率总是大于零的，当误差比率越大，误差减小得越快。 复合梯形积分的误差项是 SKIPIF 1 0 .， SKIPIF 1 0 ，当h趋于零时，显然积分的误差项更快地趋于零，实验结果复符合这一结论。 观察复化3点Gauss积分结果，发现每个结果的误差已经早早地为零了，而三个复合梯形积分的误差还是大于零的，所以我们可以得到结论，使用复化3点Gauss积分能够更精确地得到积分值。 六、附录： 实验编程，运行环境为Microsoft Visual C++ 第一个积分的程序： #include math.h #include stdio.h #include stdlib.h double f(double x) //定义函数f(x)// { double y; y=exp(-x*x); return(y); } double R(int N,double a,double b) //建立复合梯形积分// { double s,h; int i; h=(b-a)/N; s=0.0; for(i=1;i=N;i++) { s=s+h*(f(a+(i-1)*h)+f(a+i*h))/2.0; } return(s); } double T(int N,double a,double b) //建立复化3点Gauss积分// { double s,h,A1,A2,A3,x1,x2,x3; int i; h=(b-a)/N; s=0.0; for(i=1;i=N;i++)