数值分析课 第五章.pptVIP

* * §5.5 Gauss求积公式 /*Gauss Quadrature Formula */ ?Newton—Cotes求积公式中的求积节点是等距选取的，求积系数计算方便，但代数精度要受到限制； 积分公式的一般形式： ?插值型的求积公式至少有n次代数精度，至多有多少次的代数精度？ ?如何适当选取求积节点和求积系数，使求积公式达到最高的代数精度？ 窜呢疾雾掣烫穷绪火社粱块聊丽壬瓤喳化曾宽朔囚匝不货招禁纂烘跌巷皋数值分析课 第五章数值分析课 第五章 一、 Gauss积分问题的提法 ?为了提高代数精度，需要适当选择求积节点: ①当求积节点个数确定后，不管这些求积节点如何选 取，求积公式的代数精度最高能达到多少？ ②具有最高代数精度的求积公式中求积节点如何选取？ 积分公式的一般形式： 个求积节点， 个求积系数，共 个未知量，需要 个方程，因此可以取 使公式精确成立，从而求出求积节点和系数。 膳栖呸琢嘛泊讼银浆肉锚谤咨圆纤别傍寡秩洼宰毕躁电拥冗秆对抠狮晚牵数值分析课 第五章数值分析课 第五章 只需证明：对于上述插值型求积公式，存在一个 2n+2次多项式，使得求积公式不能精确成立。 形如 的插值型求积公式的代数精度最高不超过2n+1次。 证明： 令 因为 而 故求积公式不能精确成立。 写疯毫嘱君霄残背集贾衙尚吾得梧萤住酸敲靡衡环坦敏扑呜摹肥轻爸肛碘数值分析课 第五章数值分析课 第五章 下面讨论一般积分形式： 其中 为权函数 构造积分公式（*） 具有2n+1次代数精度。 其中 求积节点 求积系数 与被积函数无关 斥驭新腮书籍才舷阜薪越驱懈傣诀歼曝盯救炒秧斟摔境孕军巷详戴拐广擎数值分析课 第五章数值分析课 第五章 如果一组节点 ，使得 上述插值型求积公式具有2n+1次代数精度，则称该组节点为Gauss点，相应的公式为Gauss型求积公式。 求积系数的特征： 五点的Gauss求积公式具有多少次代数精度？ 例1：构造下列积分的Gauss求积公式： 例题： 分析：因为n=1,所以Guass求积公式具有3次代数精度。分别取 , 得到关于 的方程组，求解非线性方程组得到求积系数和求积节点。 罗峭硅店顿陋攫浑遇膀屏忆忻允芬之共亥刘之咀谐踞篮伤死蹄帆象黍毅溶数值分析课 第五章数值分析课 第五章 2n+2个未知数， 2n+2个方程的非线性方程组 由代数精度定义，当 时， 求积公式 精确成立： 问题：如何计算Gauss点 及求积系数 ？ ? 方法一：从代数精度的定义出发，求解非线性方程组； 煽冶品烬烛狞培浆哲墟役烈隘褐供晦妆匹蟹伺斥匡傅仟庞远域息永啡株万数值分析课 第五章数值分析课 第五章 ? 方法二：两步走 问题：如何计算Gauss点 及求积系数 ？ 1. 先确定Gauss求积节点 2. 计算求积系数 ①从代数精度的定义出发，求解线性方程组； 或②用系数的表达式直接计算。 丘卓导笺杖莎嘿缄票扁顶维坑段不谰中孤狗侄勾纪管位捌央亭擞瞅重淀瘦数值分析课 第五章数值分析课 第五章 二、 Gauss求积公式的性质 ? Gauss求积公式存在的条件 插值型求积公式（*）的节点 是Gauss点的充要条件是以这些节点为零点的多项式 与任何不超过n次的多项式 带权正交: 捉锋军趴缝蔓突予萧眠鳞砌秀好慌慈妻闽粒预已蛆擎享妆懈浸改位蛰瘫表数值分析课 第五章数值分析课 第五章 证明： 必要性 设 则 因为 是Gauss点 充分性 对于 其中 胎括淆界质婿臂建运开邱痴联懂郴洁藐褂笺剥姚货芳江振凸嚣桂逐赦巡梳数值分析课 第五章数值分析课 第五章 即求积公式（*）对一切不超过2n+1次的多项式精确成立 所以节点 是Gauss点 上述定理表明： 上带权的n+1次正交多项式 的零点就是求积公式（*）的Gauss点 串练快未僚寡瞥慨微贩脓釉排峭簿哆抠悉耀彤匆衅虽昏孪侦浊单呼塔富惮数值分析课 第五章数值分析课 第五章 ? Gauss求积公式中求积系数的求法 由代数精度定义，得到n+1阶线性方程组： 设已知Gauss点 或者 缸迭懒晦亥耽雇绎错白棵他骆宵犁疹库藩墙粳富霹猾想瘪忧涸宛园泥凭寡数值分析课 第五章数值分析课 第五章 ? Gauss求积公式的余项 证明