高斯型积分公式
高斯型积分公式-CAL-FENGHAI.-(YICAI)-Company One1
2 Guass-Legendre 积分程序
1. 目的意义:
可以提高数值积分的代数精度
2. 数学公式:
)
()()(1k n k k b a x f A dx x f x ∑?=≈ρ
3. 程序:
#include<>
#include<>
#define N 10
#define f(x) (cos(x))
int main()
{
int n=0;
int k=0;
int i=0;
double x[N]={};
double A[N]={};
double s=;
n=2; switch(n)
{
case 1:
{
x[1]=0;
A[1]=2;
break;
}
case 2:
{
x[1]=;
x[2]=;
A[1]=1;
A[2]=1;
break;
}
case 3:
{
x[1]=;
x[2]=0;
x[3]=;
A[1]=;
A[2]=;
3
A[3]=;
break;
}
case 4:
{
x[1]=;
x[2]=;
x[3]=;
x[4]=;
A[1]=;
A[2]=;
A[3]=;
A[4]=;
break;
}
default:
{
printf("error! 请添加数据!\n");
return 0;
}
}
4
for(i=1;i<=n;i++)
{
s=s+A[i]*f(x[i]);
}
printf("由高斯-勒让德积分公式计算得I=%lf\n",s);
return 0;
}
4.运行结果:
5.参考文献:
[1] 谭浩强. C语言程序设计[M]. 北京:清华大学出版社,2005.
[2] 秦新强. 数值逼近, 西安,2010.
5
数值分析 高斯—勒让德积分公式
高斯—勒让德积分公式 摘要: 高斯—勒让德积分公式可以用较少节点数得到高精度的计算结果,是现在现实生活中经常运用到的数值积分法。然而,当积分区间较大时,积分精度并不理想。 T he adva ntage of Gauss-Legendre integral formula is tend to get high-precision calculational result by using fewer Gauss-points, real life is now often applied numerical integration method. But the precision is not good when the length of integral interval is longer. 关键字: … 积分计算,积分公式,高斯—勒让德积分公式,MATLAB Keyword: Integral Calculation , Integral formula ,Gauss-Legendre integral formula, Matlab 引言: 众所周知,微积分的两大部分是微分与积分。微分实际上是求一函数的导数,而积分是已知一函数的导数,求这一函数。所以,微分与积分互为逆运算。 】 实际上,积分还可以分为两部分。第一种,是单纯的积分,也就是已知导数求原函数,称为不定积分。 相对而言,另一种就是定积分了,之所以称其为定积分,是因为它积分后得出的值是确定的,是一个数,而不是一个函数。 计算定积分的方法很多,而高斯—勒让德公式就是其中之一。 高斯积分法是精度最高的插值型数值积分,具有2n+1阶精度,并且高斯积分总是稳定。而高斯求积系数,可以由Lagrange多项式插值系数进行积分得到。 高斯—勒让德求积公式是构造高精度差值积分的最好方法之一。他是通过让节点和积分系数待定让函数f(x)以此取i=0,1,2....n次多项式使其尽可能多的能够精确成立来求出积分节点和积分系数。高斯积分的代数精度是2n-1,而且是最
高斯型积分公式
高斯型积分公式-CAL-FENGHAI.-(YICAI)-Company One1
2 Guass-Legendre 积分程序 1. 目的意义: 可以提高数值积分的代数精度 2. 数学公式: ) ()()(1k n k k b a x f A dx x f x ∑?=≈ρ 3. 程序: #include<> #include<> #define N 10 #define f(x) (cos(x)) int main() { int n=0; int k=0; int i=0; double x[N]={}; double A[N]={}; double s=; n=2; switch(n)
{ case 1: { x[1]=0; A[1]=2; break; } case 2: { x[1]=; x[2]=; A[1]=1; A[2]=1; break; } case 3: { x[1]=; x[2]=0; x[3]=; A[1]=; A[2]=; 3
A[3]=; break; } case 4: { x[1]=; x[2]=; x[3]=; x[4]=; A[1]=; A[2]=; A[3]=; A[4]=; break; } default: { printf("error! 请添加数据!\n"); return 0; } } 4
for(i=1;i<=n;i++) { s=s+A[i]*f(x[i]); } printf("由高斯-勒让德积分公式计算得I=%lf\n",s); return 0; } 4.运行结果: 5.参考文献: [1] 谭浩强. C语言程序设计[M]. 北京:清华大学出版社,2005. [2] 秦新强. 数值逼近, 西安,2010. 5
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求积公式的系数
高斯-勒让德积分公式
实习论文 题目高斯勒让德积分公式 专业信息与计算科学 班级计算092 学号3090811065 学生周吉瑞 指导教师秦新强 2011 年
高斯勒让德积分公式 专 业: 信息与计算科学 学 生: 周吉瑞 指导老师: 秦新强 摘要 关于数值积分公式0 ()()b n k k k a f x dx A f x =≈∑?,除了用误差来分析其精度以外,还可以 用代数精度来判断其代数精度的高低,已知n+1点Newton-Cotes 型积分公式,当n 为奇数时,其代数精度为n ,当n 为偶数时,其代数精度达到n+1。 n+1点的Newton-Cotes 型积分公式属于插值积分型积分公式,一般地,若对随机选取的n+1个节点作插值型积分公式也仅有n 次代数精度,但是,如果求积节点选取适当,就有可能提高数值积分的代数精度,高斯型积分公式就可以实现这一目标。 关 键 词:数值积分,代数精度,高斯型积分公式
一、目的意义 构造Gaoss 型求积公式除需要求出正交多项式外,还需要求出正交多项式的零点和求积系数,当3n ≥时,这些工作均很困难,因此给出高斯-勒让德积分公式的零点和系数。 二、公式 高斯-勒让德积分公式:1 1 1 ()()n k k k f x A f x -=≈∑?; 三、算法流程 Step1:输入所用的点数n ; Step2:对i=1,2,···,n 循环执行步3; Step3:I= I+ ()i i A f x ; Step4:输出I ;结束。
四、算法程序 #include
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
高斯-拉盖尔积分公式
实习论文 题目高斯拉盖尔积分公式 专业信息与计算科学 班级计算092 学号3090811065 学生周吉瑞 指导教师秦新强 2011 年
高斯拉盖尔积分公式 专 业: 信息与计算科学 学 生: 周吉瑞 指导老师: 秦新强 摘要 关于数值积分公式0()()b n k k k a f x dx A f x =≈∑?,除了用误差来分析其精度以外,还可以 用代数精度来判断其代数精度的高低,已知n+1点Newton-Cotes 型积分公式,当n 为奇数时,其代数精度为n ,当n 为偶数时,其代数精度达到n+1。 n+1点的Newton-Cotes 型积分公式属于插值积分型积分公式,一般地,若对随机选取的n+1个节点作插值型积分公式也仅有n 次代数精度,但是,如果求积节点选取适当,就有可能提高数值积分的代数精度,高斯型积分公式就可以实现这一目标。 关 键 词:数值积分,代数精度,高斯型积分公式
一、目的意义 构造Gaoss 型求积公式除需要求出正交多项式外,还需要求出正交多项式的零点和求积系数,当3n ≥时,这些工作均很困难,因此给出高斯-拉盖尔积分公式的零点和系数。 二、公式 高斯-拉盖尔积分公式:10()()n x k k k e f x A f x ∞-=≈∑?; 三、算法流程 Step1:输入所用的点数n ; Step2:对i=1,2,···,n 循环执行步3; Step3:I= I+ ()i i A f x ; Step4:输出I ;结束。 四、算法程序 #include
高斯型积分公式
Guass-Legendre 积分程序 1. 目的意义: 可以提高数值积分的代数精度 2. 数学公式: ) ()()(1k n k k b a x f A dx x f x ∑?=≈ρ 3. 程序: #include<> #include<> #define N 10 #define f(x) (cos(x)) int main() { int n=0; int k=0; int i=0; double x[N]={}; double A[N]={}; double s=; n=2; switch(n)
{ case 1: { x[1]=0; A[1]=2; break; } case 2: { x[1]=; x[2]=; A[1]=1; A[2]=1; break; } case 3: { x[1]=; x[2]=0; x[3]=; A[1]=; A[2]=;
A[3]=; break; } case 4: { x[1]=; x[2]=; x[3]=; x[4]=; A[1]=; A[2]=; A[3]=; A[4]=; break; } default: { printf("error! 请添加数据!\n"); return 0; } }
for(i=1;i<=n;i++) { s=s+A[i]*f(x[i]); } printf("由高斯-勒让德积分公式计算得I=%lf\n",s); return 0; } 4.运行结果: 5.参考文献: [1] 谭浩强. C语言程序设计[M]. 北京:清华大学出版社,2005. [2] 秦新强. 数值逼近, 西安,2010.
Gauss-Legendre积分公式
数值逼近实习 题目二重积分的复化梯形公式 专业信息与计算科学 班级计算092 学号3090811072 学生薛藏朋 指导教师秦新强 2011 年
一、实验目的 1.利用Gauss-Legendre 公式计算积分 2.比较计算误差与实际误差 二、数学模型 ??? ????∈--=++≈?∑-=),(),(12)(][)]()(2)([2)(''211b a f h a b f R b f x f a f h dx x f n b a n k k ηη 三、算法 Step 1:输入等分数n Step2:输入积分上下限; Step3: 求出步长及对应个点; Step4: 由??? ????∈--=++≈?∑-=),(),(12)(][)]()(2)([2)(''211b a f h a b f R b f x f a f h dx x f n b a n k k ηη计算积分结果 Step5:将积分结果输出; 四、程序 #include
double f(double x) {double z; z=cos(x); return(z); } int main() { int i, n,m; double X[N],A[N],F=0; cout<<"请输入代数精度n"<
高斯积分节点表
高斯积分点表 N = 2 节点权重 -0.5773502691896250 1.0000000000000000 0.5773502691896250 1.0000000000000000 N = 3 节点权重 -0.7745966692414830 0.5555555555555550 0.0000000000000000 0.8888888888888880 0.7745966692414830 0.5555555555555550 N = 4 节点权重 -0.8611363115940520 0.3478548451374530 -0.3399810435848560 0.6521451548625460 0.3399810435848560 0.6521451548625460 0.8611363115940520 0.3478548451374530 N = 5 节点权重 -0.9061798459386640 0.2369268850561890 -0.5384693101056830 0.4786286704993660 0.0000000000000000 0.5688888888888880 0.5384693101056830 0.4786286704993660 0.9061798459386640 0.2369268850561890 N = 6 节点权重 -0.9324695142031520 0.1713244923791700 -0.6612093864662640 0.3607615730481380 -0.2386191860831960 0.4679139345726910 0.2386191860831960 0.4679139345726910 0.6612093864662640 0.3607615730481380 0.9324695142031520 0.1713244923791700 N = 7 节点权重 -0.9491079123427580 0.1294849661688690 -0.7415311855993940 0.2797053914892760 -0.4058451513773970 0.3818300505051180 0.0000000000000000 0.4179591836734690 0.4058451513773970 0.3818300505051180 0.7415311855993940 0.2797053914892760 0.9491079123427580 0.1294849661688690