Legendre-Guass求积公式的算法实现.docVIP

《 数 值 逼 近 》课 程 论 文 题目：Legendre-Guass求积公式的算法实现 学校：黑龙江科技学院 院系：计算机学院信息与计算科学系 班级：信计09-1班 姓名：杨朝旭 学号：2009025213  Legendre-Guass求积公式的算法实现 杨朝旭 （黑龙江科技学院　计算机与信息工程学院，哈尔滨　150027） 摘 要: 为更好地运用Legendre-Guass求积公式，针对给定区间[a,b]和权函数，运用计算机程序对其进行算法实现，并对Guass求积公式做了深入分析。 关键词: Guass求积公式；Legendre 中国分类号: TG174.44　　 文献标识码: A　 Graphite surfaces of research, preparation and prospect analysis Yang Zhaoxu (College of Computer and Information Engineering Heilongjiang Institute of Science and Technology, Harbin 150027 ) Abstract: For better use the Legendre-Guass for product formula, with a given interval [a, b] and weight function, using computer program of the algorithm, and the Guass product formula for further analysis. Key words: Guass for product formula, Legendre 0 一般Guass求积公式定义 定义: 对于形如： 的插值公式，当X1，X2，…Xn是[]区间上，以为权的次直交多项式的根时，公式达到最高代数精度2n-1的求积公式,称为Guass求积公式。 其中Guass求积公式的节点Xk称为Guass点,系数Ak称为Guass系数。 其中Guass余项为Rn为： 其中(( (a,b),w(x)=(x-x1)(x-x2)…..(x-xn)。 证明代数精度为2n-1 (step 1:证明代数精度达不到2n次) 证明：分别取 f(x)=1, x，x2，...xn 时代入公式，并让其成为等式得 A1 + A2 +…+ An =∫ab1dx.= b-a x1 A1 + x2 A2+ … +xn An =∫abxdx= （b2-a 2)/2 ...... x1 rA1 + x2 rA2+ …… +xn rAn =∫abxr dxr =(br+1-a r+1)/ (r+1) 上式共有 r 个 等式，2n个待定系数(变元),要想如上方程组有唯一解，应有方程组中方程的个数等于变元的个数,即 r=2n,这样求出的解答应的求积公式的代数精度至少是2n-1,下面证明代数精度只能是2n-1。 (step 2:再证明达到了2n-1次) 事实上,取 2n次多项式g(x)=(x-x1)2(x-x2)2….(x-xn)2 代入求积公式,有 左= 右= =0 左(右,故不成立等式,定理得证。 一般Gauss求积公式与古典Gauss求积公式的线性转换 对于一般有限区间[a,b]，用线性变换 ，使它变成为[-1,1]。 将积分区间从[a,b]变成[-1,1],由定积分的换元积 分法有： 这样就可以用Gauss - Legendre求积公式计算一般区间的积分。 Gauss - Legendre 求积公式的计算 Gauss - Legendre 求积公式： （1） 其中高斯点为Legendre多项式的零点，即： Ln(x)= （2） Guass点xk, Guass系数Ak都有表可以查询，如下表： 所以，古典求积 公式，可表达为： 当 ； 当