第4章节 - 数值积分.pptVIP

第四章 数值积分 经典方法(插值型数值求积方法) Gauss数值积分方法 复化数值求积公式 积分方程的数值求解 设f (xk)的近似值为 (12) 由近似值 所得数值积分为 误差E 对误差 若 称 为数值稳定的，反之为数值不稳定的。 推导见下页 所以， N-C是数值稳定的。 (13) 当 则 (b). 收敛性 当 系数Ak有正有负， N-C是数值不稳定的。 若 称 是收敛的，反之为不收敛的。 对于余项 注： 给定n+1个节点，插值型求积公式： 优点：代数精度高： 问题：代数精度最大是多少？如何寻求数值稳定的方法？ 缺点：数值不一定稳定。 由于f (x) - Ln(x)不收敛到0，因此Rn[ f ]不收敛到0，即Q[ f ]不收敛，因此，对于节点较多的情况，需要使用分段线性或者Hermite插值。 Runge现象 其中 一、 最高代数精度求积公式 问题： 结论： 本节关键 §2 Gauss型求积公式 设有n+1个节点，插值型求积公式的代数精度 m的最大值？如何确定？ 由求积系数及n+1个节点xi，i = 0,1,…, n 的分布确定。 四个未知量A0，A1，x0，x1，已知插值型求积公式的代数精度最高。可按插值型求积公式来求A0，A1。 解： 具有尽可能高的代数精度。 例 求节点 ，使插值型求积公式 (14) 分析： x0，x1待定 插值型求积的代数精度 于是，求积公式为 一般地，对于任意求积节点 任意求积系数，求积公式 对于 代数精度 结论 分析： 只需证明 使得 事实上，令 有 从而 而 前例中m=3=2×1+1=2n+1是能达到的最高代数精度。 次多项式f (x), Gauss型求积公式的构造 ——利用正交多项式的根构造代数精度最高的求积公式 分析： 引理1： (15) 令f (x)是任意次数 的代数多项式，则 其中，q(x)为商式，r(x)为余式，均为任意次数 的多项式 证明： 插值型求积公式代数精度大于n 正交多项式基本性质 ωn+1(x)是n+1次正交多项式，则 定义： 正交多项式的根一定是Gauss点，那么Gauss点是否一定是正交多项式的根？ n+1个节点(a x0… xn b)的求积公式(15)若其代数精度m=2n+1，即达到最高，称之为Gauss型求积公式, 并称其节点为Gauss点。 而 是[a,b]上关于权ρ(x)的n+1次正交多项式的根。 二、Gauss点与正交多项式的关系 定理4： 分析： Gauss点ax0…xn b 求积公式(15)是Gauss型的 充分性即是引理1的结论。以下只证必要性 只需证 关于 正交。 必要性，即Gauss点作为节点正是n+1次正交多项式的根。 证明： 取2n+1次多项式 则有 三、Gauss求积公式的余项（截断误差） 由引理1知, xi(i=0,1,…,n)是Gauss点,则m=2n+1, 由n+1个点, 确定2n+1次多项式。联想Hermite插值多项式。 定理5： ，则Gauss求积公式(15)的余项为 分析： (16) 证明： 若f(x)的Hermite插值多项式H2n+1(x) 满足插值条件 利用 由积分中值定理 Gauss型求积公式是数值稳定的 四、Gauss求积公式的数值稳定性和收敛性 1. 稳定性 首先，Gauss型求积系数Ak 0, (k = 0,1,…,n) 事实上，m = 2n+1 为L-插值基函数, 为2n次多项式 其次，取f (x)=1 故数值稳定。 证明： 上的连续函数 可以用代数多项式一致逼近, 对任意给定的 存在某个多项式 2. 收敛性 引理2： 上的任何连续函数 对于有限闭区间 (17) 从而 3. 结论： Gauss型求积公式是数值稳定的；且对有限闭区间上的连续函数，Gauss型求积公式的值随节点数目的增加而收敛到准确积分值。 收敛、稳定； 计算量小，代数精度高。 Gauss点难求（即多项式的根难求）； Gauss点是无理数， Gauss求积系数也是无理数。 定理6 优点： 缺点： 五、常用的Gauss型求积公式 Gauss型节点是多项式的根，因此与正交多项式联系起来，有以下几种常用求积公式。 1. Gauss-Legendre （勒让德）求积公式(Gauss点和积分系数可查表得到) 若区间[a, b]≠[-1,1], 可用变量替换把区间 2. Gauss- chebyshev （切比雪夫）求积公式 3. Gauss-Laguerre （拉盖尔）求积公式 4. Gauss-Hermite 求积公式 插值型求积