数值计算方法-第7章-数值微分与数值积分-part-2.ppt

精选ppt 五. Gauss型求积公式 第7章 数值微分与数值积分 目的 求积公式: 当节点数n固定时, 选取适当的节点{xk}及系数{Ak}, 使其具有最高的代数精度. 为权函数. 对所有 精确成立. Gauss型求积公式的思想 回顾: 若具有m次代数精度, 则: 其中, 这里有m+1个方程, 未知量有2n个: xi, Ai (i=1,2,…,n) 可以证明: 当m+1≤2n, 即m≤2n-1时, 方程有解. 求积公式的最大代数精度 即, 当m=2n-1时, 可以找到一组解{xk},{Ak}, 使积分公式成立, 即代数精度可以达到2n-1. 求积公式不精确. 第7章 数值微分与数值积分 Gauss型积分公式的定义 另一方面, 当m=2n时, 取一个特殊的多项式: 求积公式: 而精确解: 结论: n个节点的积分公式最高代数精度为2n-1. 对于n个节点的积分公式: , 如果具有 2n-1次代数精度, 则称为带权函数 的Gauss型求积公式. xi 称为Gauss点. 如何求Gauss点及积分系数 直接求解上面的非线性方程组比较困难, 可采用正交多项式来求. Gauss点与正交多项式的关系 从简单问题得到的启示 第7章 数值微分与数值积分 当n=2, 要使 具有2n-1=3次代数精度, 则对任意3次多项式P3(x), 利用多项式除法, 因为积分公式具有3次精度, 对 成立. 精确成立. 从简单问题得到的启示（续） 第7章 数值微分与数值积分 分别取: 解得: Gauss点 再取两个特殊的多项式, 如: 于是, 求积公式: 具有3次代数精度. 推广到一般n节点的情况 第7章 数值微分与数值积分 将任意不超过2n-1次的多项式写为: 其中, 都是不超过n-1次的多项式. 要求积分公式具有2n-1次代数精度, 则 精确成立. 即 而 精确成立. 对 次多项式q(x)成立. 称为正交条件 分别令q(x)=1,x,x2,…,xn-1, 得到n个方程, 解出xi (i=1,2,…,n). 下面证明, xi 一定是Gauss 点. 定理7.4 Gauss积分点的充要条件 第7章 数值微分与数值积分 即, n点积分公式 中, xi (i=1,…n) 为Gauss 点 其中, 的充要条件是, 对于任意不超过 n-1 次的多项式p(x), 与 在[a,b]区间上关于权函数 正交, 即, 证明.[必要性] 若xi (i=1,…,n)为Gauss点, 则对于任意不超过n-1 次的多项式p(x), 是不超过2n-1 次的多项式, Gauss 积分精确成立, 即, 满足正交条件 必要性得证. 定理7.4 【证明】(续-1) 第7章 数值微分与数值积分 [充分性] 对于任意不超过 n-1 次的多项式q(x), 正交条件成立: 对于任意不超过2n-1次多项式 , 总可以写成: 不超过n-1次的多项式 且 充分性假设 ………. (*) 下面证明, 以xi 为积分点, 一定能找到合适的Ai, 使 精确成立. 这样, 积分公式至少具有2n-1次代数精度, xi 即为Gauss点. 为此, 取Ak 满足: 定理7.4 【证明】(续-2) 第7章 数值微分与数值积分 对于任意不超过n-1次的多项式r(x), 可写成: 这样, 积分公式至少具有2n-1次代数精度, xi 即为Gauss点. 证毕# Vandermonde矩阵, 非奇异. 有唯一解 即, 注意到, 于是, (*) 式: 精确成立. Gauss点正好是正交多项式的根 第7章 数值微分与数值积分 由第6章定理6.3: 如果正交多项式 的k个根为: xi (i=1,…,k), 则: c 是给定的常数. 对 次多项式Qk-1(x)成立. 设 是最高次系数非零的k次多项式, 则 是[a,b]上关于权函数 的正交多项式的充要条件是, 对任意次数不超过k-1次的多项式 , 都有: 由定理7.4, xi 一定是Gauss 积分点. 故Gauss点正好是正交多项式的根. 如何求Gauss积分系数Ai 第7章 数值微分与数值积分 分别取n-1