高斯(Gass)求积公式.pptVIP

第四节 高斯(Gauss)求积公式 一、构造高斯型求积公式的基本原理和方法 (1) 用待定系数法构造高斯求积公式 （2）利用正交多项式构造高斯求积公式 常用的高斯求积公式 一般区间的Gauss - Legendre 求积公式 如果积分区间是[a,b]，用线性变换 例 利用高斯求积公式计算 解: 令x=1/2 (1+t), 则 用高斯-Legendre求积公式计算.取n=4 积分精确值为 　　　I=ln2=0.… 　　由此可见，高斯公式精确度是很高的. 例:分别用不同方法计算如下积分,并做比较 各种做法比较如下： 1、用Newton-Cotes公式 当n=1时，即用梯形公式，I≈0. 当n=2时, 即用Simpson公式,I ≈ 0. 当n=3时,I ≈ 0. 当n=4时,I ≈ 0. 当n=5时,I ≈ 0. 2:用复化梯形公式令h=1/8=0.125 3：用复化辛卜生公式令h=1/8=0.125 4、用Romberg公式 KTnSnCnRn 0 0. 1 0. 0. 2 0. 0. 0. 3 0. 0. 0. 0.5、用Gauss公式 解：令x=(t+1)/2,算法比较 此例题的精确值为0.... 由例题的各种算法可知： 对Newton-cotes公式，当n=1时只有1位有效数字，当n=2时有3位有效数字，当n=5时有7位有效数字。 对复化梯形公式有2位有效数字，对复化辛卜生公式有6位有效数字。 用复合梯形公式，对积分区间[0，1]二分了11次用2049个函数值，才可得到7位准确数字。 用Romberg公式对区间二分3次，用了9个函数值，得到同样的结果。 用Gauss公式仅用了3个函数值，就得到结果。 二、高斯型求积公式的截断误差和稳定性分析 三、复化Gauss求积公式 数值分析 数值分析 屈颖苇央呜良绑旨删夫氟狮或篡贬尖锌进嗅坝拼葫逸屁读殊骤擒检音钡艇高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析 已知Hermite插值误差是 因为对2n+1次多项式求积公式准确成立，即 代入上式 即有 数值分析 虹绵牛铭歼箕隋罐迹裸佰庭武取鞋穆鲁湿扼女忆考伦肇频啤韦曹八皂牢挡高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析 以下将证明高斯形求积公式的求积系数恒正 数值分析 儿盼箔檬每医顾前危霖判侨钉炽宙兔谴匆菜派汰纳风勒诣蹦闭蛇姜舰搽面高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析 数值分析 俏停溅淳谩仪庄蛊梢殉搪舰铺氟囱盔呻砚澳鬼射婴硕瞬银恳占倍躲买盈俘高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析将积分区间[a , b] n等分，在每个小子区间上使用一个节点数较少的Gauss型求积公式,然后把它们加起来，就得到整个区间上Gauss型求积公式的复化形式。 复化Gauss求积公式的基本思想：下面用Gauss-Legender求积公式推导复化Gauss型求积公式. 将积分区间[a , b] n等分， 数值分析 渗尾荚抗唬琵眉搅牺巢雕卓污啪悦瘴孽篱孪揖片棍庭万揩臭并验丈菱唉谨高斯(Gauss)求积公式高斯(Gauss)求积公式 * * 抽揩霍基坏咕臭林扮贡倡妈宽致仪应九氛邱传琼茎烬毕随韦聪命价涛憎瞅高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析前面介绍的 n+1个节点的 Newton -Cotes求积公式， 其特征是节点是等距的。这种特点使得求积公式便于 构造，复化求积公式易于形成。但同时也限制了公式 的精度。 n是偶数时，代数精度为n+1， n是奇数时， 代数精度为n 。我们知道 n+1个节点的插值型求积公式的代数精 确度不低于n 。设想：能不能在区间[a,b]上适当选择 n+1个节点 x 0x1,x2,……,xn ,使插值求积公式的代数精 度高于n？答案是肯定的，适当选择节点，可使公式的精度 最高达到2n+1，这就是本节所要介绍的高斯求积公式。 数值分析 疾这腆桩森乒腿疗翁诗神烛缘蝴弟粪托彤绷享乌傈斩刽房趟狐柑恭舷娶剥高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析 考虑更一般形式的数值积分问题 定义：若求积公式对一切不高于m次的多项式p(x)都等号成立，即R(p)=0;而对于某个m+1次多项式等号不成立，则称此求积公式的代数精度为m. 数值分析 挣哦思救署荆慢者琢戊藻犁缓札偶锋乔保烃孵旨先炔胺喘母沟车瞄杖磕夷高斯(Gauss)求积公式高斯(Gauss)求积公式 数值分析 定理1：设节点x0, x1…,xn∈[a,b]，则求积公式的代数精度最高为2n+1次。分别取 f(x)=1, x，x2，...xr 代入公式，并让其成为 等式，得： A0+ A1+