高斯(Gauss)求积公式汇总.

Gauss型积分公式摘要求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。

这时我们可以通过数值方法求出函数积分的近似值。

当然再用近似值代替真实值时，误差精度是我们需要考虑因素，但是除了误差精度以外，还可以用代数精度来判断其精度的高低。

已知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）多项式如下定义的多项式称作勒让德多项式。

由于是次多项式，所以是n次多项式，其最高次幂的系数与多项式的系数相同。

也就是说n次勒让德多项式具有正交性即勒让德多项式是在上带的n次正交多项式，而且这时Gauss型积分公式的节点就取为上述多项式的零点，相应的Gauss型积分公式为此积分公式即成为高斯-勒让德积分公式。

其中Gauss-Legendre求积公式的系数1其中k的取值范围为Gauss点和系数不容易计算，但是在实际计算中精度要求不是很高，所以给出如下表所示的部分Gauss点和系数，在实际应用中只需查表即可。

摘要求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。

这时我们可以通过数值方法求出函数积分的近似值。

当然再用近似值代替真实值时，误差精度是我们需要考虑因素，但是除了误差精度以外，还可以用代数精度来判断其精度的高低。

已知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）多项式如下定义的多项式L n(x)=12n n!d ndx n(x2−1)n,x∈[−1,1],n=0,1,2⋯称作勒让德多项式。

由于(x2−1)n是2n次多项式，所以L n(x)是n次多项式，其最高次幂的系数A n与多项式1 2n n!d ndx n(x(2n))=12n n!2n(2n−1)(2n−2)⋯(n+1)x n的系数相同。

也就是说n次勒让德多项式具有正交性即勒让德多项式L n(x)是在[−1,1]上带ρ(x)=1的n次正交多项式，而且(L m,L n)=∫L m(x)L n(x)dx1−1={0, m≠n22n+1, m=n这时Gauss型积分公式的节点就取为上述多项式L n(x)的零点，相应的Gauss型积分公式为∫f(x)dx 1−1≈∑A k f(x k) nk=1此积分公式即成为高斯-勒让德积分公式。

gauss型求积公式一、Gauss型求积公式的基本概念。

1. 定义。

- 在数值积分中，Gauss型求积公式是一种高精度的求积公式。

对于积分∫_a^bf(x)ρ(x)dx（其中ρ(x)为权函数），Gauss型求积公式的形式为∫_a^bf(x)ρ(x)dx≈∑_i = 1^nA_if(x_i)。

这里x_i称为求积节点，A_i称为求积系数，n为求积公式的节点个数。

2. 特点。

- 高精度：Gauss型求积公式具有很高的代数精度。

对于n个节点的Gauss型求积公式，其代数精度为2n - 1。

这意味着对于次数不超过2n-1的多项式f(x)，该求积公式能精确成立，即∫_a^bP_m(x)ρ(x)dx=∑_i = 1^nA_iP_m(x_i)，其中m≤slant2n - 1，P_m(x)是m次多项式。

- 节点分布：Gauss型求积公式的节点x_i不是等距分布的。

这些节点是关于权函数ρ(x)正交的多项式的零点。

例如，当ρ(x) = 1，[a,b]=[- 1,1]时，对应的正交多项式是勒让德多项式P_n(x)，Gauss型求积公式的节点就是勒让德多项式的零点。

二、求积节点与求积系数。

1. 求积节点的确定。

- 以勒让德 - Gauss求积公式为例（ρ(x)=1，[a,b]=[-1,1]），求积节点x_i是勒让德多项式P_n(x)的零点。

勒让德多项式P_n(x)可以通过递推公式(n + 1)P_n +1(x)=(2n + 1)xP_n(x)-nP_n - 1(x)，P_0(x)=1，P_1(x)=x来计算。

通过求解P_n(x)=0得到求积节点x_i。

2. 求积系数的计算。

- 求积系数A_i可以通过多种方法计算。

一种常见的方法是利用正交性条件。

对于勒让德 - Gauss求积公式，求积系数A_i可以通过公式A_i=(2)/((1 -x_i)^2)[P_{n'(x_i)]^2}计算，其中P_n'(x)是勒让德多项式P_n(x)的导数。

Gauss型积分公式摘要求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。

这时我们可以通过数值方法求出函数积分的近似值。

当然再用近似值代替真实值时，误差精度是我们需要考虑因素，但是除了误差精度以外，还可以用代数精度来判断其精度的高低。

已知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）多项式如下定义的多项式称作勒让德多项式。

由于是次多项式，所以是n次多项式，其最高次幂的系数与多项式的系数相同。

也就是说n次勒让德多项式具有正交性即勒让德多项式是在上带的n次正交多项式，而且这时Gauss型积分公式的节点就取为上述多项式的零点，相应的Gauss型积分公式为此积分公式即成为高斯-勒让德积分公式。

其中Gauss-Legendre求积公式的系数1其中k的取值范围为Gauss点和系数不容易计算，但是在实际计算中精度要求不是很高，所以给出如下表所示的部分Gauss点和系数，在实际应用中只需查表即可。

G a u s s型积分公式-CAL-FENGHAI.-(YICAI)-Company One1摘要求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。

这时我们可以通过数值方法求出函数积分的近似值。

当然再用近似值代替真实值时，误差精度是我们需要考虑因素，但是除了误差精度以外，还可以用代数精度来判断其精度的高低。

已知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）多项式如下定义的多项式称作勒让德多项式。

由于是次多项式，所以是n次多项式，其最高次幂的系数与多项式的系数相同。

也就是说n 次勒让德多项式具有正交性即勒让德多项式是在上带的n次正交多项式，而且这时Gauss 型积分公式的节点就取为上述多项式的零点，相应的Gauss型积分公式为12此积分公式即成为高斯-勒让德积分公式。

其中Gauss-Legendre 求积公式的系数其中k 的取值范围为Gauss 点和系数不容易计算，但是在实际计算中精度要求不是很高，所以给出如下表所示的部分Gauss 点，在实际应用中只需查表即可。

高斯（Gauss）求积公式的系数和确定方法如下：确定节点：首先确定求积公式所使用的节点，这些节点通常选择为高斯点。

构造高斯型求积公式：根据所选的节点，构造高斯型求积公式。

高斯型求积公式的一般形式为：∫f(x)dx≈∑(A*f(x_i))，其中A是求积系数，x_i是高斯点。

确定求积系数：通过求解线性方程组来确定求积系数。

具体地，根据高斯型求积公式的构造原理，可以建立一个线性方程组，该方程组由节点处的函数值和高斯型求积公式中的求积系数组成。

解这个线性方程组可以得到求积系数。

验证求积公式的精度：通过数值试验来验证求积公式的精度。

例如，可以选择一些已知的函数进行测试，比较使用高斯型求积公式计算的结果与真实值之间的误差。

Central South University《MATLAB 程序设计实践》《matlab 程序设计实践》课程考核1、编程实现以下科学计算算法，并举一例应用之。

（参考书籍《MATLAB 与科学计算》，王沫然著，电子工业出版社，2009年）“高斯（Gauss ）求积法”（1）算法说明：为了使求积公式得到较高的代数精度，可以使用到Guass 公式：∑⎰=-≈nkk k x f A dx x f 011)()(式中k A 为Gauss-Legendre 系数。

它可使计算达到2n-1次的代数精度。

对任意的求积区间[a,b]，通过变换x=(a-b)x/2+(a+b)/2，可以转化到区间[-1,1]上。

此时⎰⎰-++--≈badt ba t ab f a b dx x f 11)22(2)(（2）流程图（3）源程序代码： gauss1.mfunction s=gauss1(a,b,n) h=(b-a)/n; s=0.0;for m=0:(1*n/2-1)s=s+h*(gaussf(a+h*((1-1/sqrt(3))+2*m))+gaussf(a+h*((1+1/sqrt(3))+2*m)));end应用举例：求下面的积分值。

G=dx x ⎰1)cos(解：先确定n 的值。

由高斯余项可得下式：6451052701270-<=x nnh 计算得到n=6。

（1）源程序代码：（编制函数文件）gaussf.mfunction y=gaussf(x) y=cos(x)（在Matlab 命令窗口中计算得）>> gauss1(0,1,6)ans =0.8415（2）流程图：2.已知Appolo 卫星的运动轨迹（x,y ）满足下面方程：323122)()(2r x r x x dt dy dt x d λμμλ--+-+=， 3231222r y r y y dt dx dt y d μλ--+-= 其中，222221)(,)(,1,45.82/1y x r y x r ++=++=-==λμμλμ ，试在初值x(0)=1.2,x ’(0)=0, y(0)=0, y ’(0)=-1.04935371下求解，并绘制Appolo 卫星轨迹图。

gauss-legendre求积公式
高斯-勒让德求积公式是指利用勒让德多项式的零点来求解定积分的数值近似方法。

它是一种高精度、高效率的积分数值求解方法。

具体地，用勒让德多项式的根来代替单位区间上的等距节点，再用适当的权系数乘以函数值，然后把结果求和，便得到积分的数值近似值。

高斯-勒让德求积公式可以表示为：
$\int_{-1}^{1}f(x)dx≈\sum_{i=0}^{n}w_if(x_i)$。

其中，$f(x)$为被积函数，$x_i$为勒让德多项式的$n$个零点，$w_i$为相应的权系数。

当$n$越大时，数值近似的精度越高。

高斯-勒让德求积公式的主要优点在于可以在任何有限区间内求解积分。

此外，它还具有收敛速度快、精度高、有理权系数等优点。

但是，它的缺点也是显而易见的，主要在于它只适用于比较光滑的函数，对于函数间断点、壁垒等特殊情况不太适用。

Gauss型（Gaussianquadrature）求积公式和⽅法⽬录0、Gauss型积分通⽤形式1、Gauss–Legendre quadrature勒让德2、Gauss–Laguerre quadrature拉盖尔——积分区间[0，inf]3、Chebyshev–Gauss quadrature切⽐雪夫0、Gauss型积分通⽤形式The integration problem can be expressed in a slightly more general way by introducing a positive weight functionω into the integrand（被积函数）, and allowing an interval other than（除了，不同于） [−1, 1]. That is, the problem is to calculatefor some choices of a, b, and ω. For a = −1, b = 1, and ω(x) = 1, the problem is the same as that considered above（勒让德问题）. Other choices lead to other integration rules. Some of these are tabulated（列表） below.1、Gauss–Legendre quadrature勒让德——积分区间[-1，1]The most common domain of integration for such a rule is taken as [−1, 1], so the rule is stated aswhich is exact for polynomials of degree 2n − 1 or less. This exact rule is known as the Gauss-Legendre quadrature rule. The quadrature rule will only be an accurate approximation to the integral above if f(x) is well-approximated by a polynomial of degree 2n − 1 or less on [−1, 1]. The Gauss-Legendre quadrature rule is not typically used for integrable functions with endpoint singularities.（端点奇点）（1）基本概念注：P0没有根（与x轴⽆交点），P1有1个根（与x轴有⼀个交点），P2有2个根（与x轴有两个交点），。

我们要找出高斯-勒让德求积公式的复合求积公式。

首先，我们需要了解高斯-勒让德求积公式的基本形式。

高斯-勒让德求积公式是用于数值积分的一种方法，它基于多项式的零点（即，使多项式为0的x值）来选择求积点。

对于一个n次多项式P(x)，其零点为x1, x2, ..., xn，高斯-勒让德求积公式为：
∫P(x) dx = ∑ (wi × P(xi))
其中，wi和xi是高斯-勒让德求积公式的权重和零点，它们可以通过求解某些方程组来得到。

当我们需要使用高斯-勒让德求积公式进行多次积分时，我们可以使用复合求积公式。

复合求积公式的一般形式为：
∫∫ P(x, y) dxdy = ∑ (wi × P(xi, yi))
其中，P(x, y)是我们要积分的函数，wi, xi, yi是高斯-勒让德求积公式的权重和零点。

为了得到具体的复合求积公式，我们需要知道P(x, y)的形式以及它的零点。

对于给定的函数P(x, y) = x^2 + y^2，其零点为(0, 0), (1, 0), (0, 1), (-1, 0), (0, -1)。

使用高斯-勒让德求积公式，我们可以得到复合求积公式为：
∫∫ (x^2 + y^2) dxdy = 1/36 * (3pi + 4sqrt(2) + 4)
其中，pi是圆周率。

摘要求函数在给定区间上的定积分，在微积分学中已给出了许多计算方法，但是，在实际问题计算中，往往仅给出函数在一些离散点的值，它的解析表达式没有明显的给出，或者，虽然给出解析表达式，但却很难求得其原函数。

这时我们可以通过数值方法求出函数积分的近似值。

当然再用近似值代替真实值时，误差精度是我们需要考虑因素，但是除了误差精度以外，还可以用代数精度来判断其精度的高低。

已知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）多项式如下定义的多项式L n(x)=12n n!d ndx n(x2−1)n,x∈[−1,1],n=0,1,2⋯称作勒让德多项式。

由于(x2−1)n是2n次多项式，所以L n(x)是n次多项式，其最高次幂的系数A n与多项式1 2n n!d ndx n(x(2n))=12n n!2n(2n−1)(2n−2)⋯(n+1)x n的系数相同。

也就是说n次勒让德多项式具有正交性即勒让德多项式L n(x)是在[−1,1]上带ρ(x)=1的n次正交多项式，而且(L m,L n)=∫L m(x)L n(x)dx1−1={0, m≠n22n+1, m=n这时Gauss型积分公式的节点就取为上述多项式L n(x)的零点，相应的Gauss型积分公式为∫f(x)dx 1−1≈∑A k f(x k) nk=1此积分公式即成为高斯-勒让德积分公式。