Gauss型积分公式题稿.docx

摘要 求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。这时我们可以通过数值方法求出函数积分的近似值。 当然再用近似值代替真实值时，误差精度是我们需要考虑因素，但是除了误差精度以外，还可以用代数精度来判断其精度的高低。已知n+1点的Newton-Cotes型积分公式，当n为奇数时，其代数精度为n；当n为偶数时，其代数精度达到n+1。若对随机选取的n+1个节点作插值型积分公式也仅有n次代数精度。 如何选取适当的节点，能使代数精度提高？Gauss型积分公式可是实现这一点，但是Gauss型求积公式，需要被积函数满足的条件是正交，这一条件比较苛刻。因此本实验将针对三种常用的Gauss型积分公式进行讨论并编程实现。 关键词：Newton-Cotes型积分公式 正交多项式 代数精度  1、实验目的 1) 通过本次实验体会并学习Gauss型积分公式，在解决如何取节点能提高代数精度这一问题中的思想方法。 2) 通过对Gauss型积分公式的三种常见类型进行编程实现，提高自己的编程能力。 3) 用实验报告的形式展现，提高自己在写论文方面的能力。 2、算法流程 下面介绍三种常见的Gauss型积分公式 1) 高斯-勒让德（Gauss-Legendre）积分公式 勒让德（Legendre）多项式 如下定义的多项式 Lnx=12nn!dndxn(x2-1)n,x∈-1,1,n=0,1,2? 称作勒让德多项式。由于(x2-1)n是2n次多项式，所以Lnx是n次多项式，其最高次幂的系数An与多项式 12nn!dndxnx2n=12nn!2n(2n-1)(2n-2)?(n+1)xn 的系数相同。也就是说n次勒让德多项式具有正交性即勒让德多项式Lnx是在-1,1上带ρx=1的n次正交多项式，而且 Lm,Ln=-11Lm(x)Ln(x)dx=0, m≠n22n+1, m=n 这时Gauss型积分公式的节点就取为上述多项式Lnx的零点，相应的Gauss型积分公式为 -11f(x)dx≈k=1nAkf(xk) 此积分公式即成为高斯-勒让德积分公式。 其中Gauss-Legendre求积公式的系数 Ak=-11ρxωnxx-xkωnxdx=-11ρxLnxx-xkLnxdx 其中k的取值范围为k=1,2,?,n Gauss点和系数不容易计算，但是在实际计算中精度要求不是很高，所以给出如下表所示的部分Gauss点xk和系数Ak，在实际应用中只需查表即可。 n x A n x A 1 0 2 6 ±0.9324695142 ±0.6612093865 ±1.2386191816 0.171324492 0.360761573 0.467913934 2 ±0.5773502692 1 7 ±0.9491079123 ±0.7415311856 ±0.4058451514 0 0.129484966 0.279705391 0.381830050 0.417959183 3 ±0.77459666920 0 0.5555555556 0.8888888889 4 ±0.8611363116 ±0.3399810436 0.3478548451 0.6521451549 8 ±0.9602898565 ±0.7966664774 ±0.5255324099 ±0.1834346425 0.101228536 0.222381034 0.313706645 0.362683783 5 ±0.9061798459 ±0.5384693101 0 0.2369268851 0.4786286705 0.5688888889 2) 高斯-拉盖尔（Gauss-Laguerre）积分公式 拉盖尔（Laguere）多项式 Lnx=exdndxn(xne-x),0≤x+∞,n=0,1,2? 称为拉盖尔多项式。其首项系数为(-1)n,且具有性质： 正交性，在区间0,+∞上关于权函数ρx=e-x正交，而且 Lm,Ln=0∞e-xLm(x)Ln(x)dx=0, m≠n(n!)2, m=n 积分区间为0,+∞，权函数为ρx=e-x的Gauss型积分公式称为高斯-拉盖尔积分公式，其中Gauss点为拉盖尔多项式Lnx的零点，高斯-拉盖尔积分公式为 0∞e-xf(x)dx≈k=1nAkf(xk) 同样高斯-拉盖尔积分公式的Gauss点和求积系数如下表所示： n x A n x A 2 0.5857864376 3.