第4章 数值积分.pptVIP

一般地,其中说明：为了既保证精度又节约时间，尽量选用n是偶数的情形。*第29页，共71页，星期日，2025年，2月5日3.Newton—Cotes公式的数值稳定性和收敛性(a).数值稳定性若某个求积公式的舍入误差，即f(xk)的误差对数值积分的结果影响较小，则称该数值求积公式是稳定的；否则，若影响较大，则称为不稳定的。由实验和观测得到，本身有误差精确值*第30页，共71页，星期日，2025年，2月5日设f(xk)的近似值为(12)由近似值所得数值积分为误差E对误差若称为数值稳定的，反之为数值不稳定的。推导见下页*第31页，共71页，星期日，2025年，2月5日所以，N-C是数值稳定的。(13)当则*第32页，共71页，星期日，2025年，2月5日(b).收敛性当系数Ak有正有负，N-C是数值不稳定的。若称是收敛的，反之为不收敛的。对于余项*第33页，共71页，星期日，2025年，2月5日注：给定n+1个节点，插值型求积公式：优点：代数精度高：问题：代数精度最大是多少？如何寻求数值稳定的方法？缺点：数值不一定稳定。由于f(x)-Ln(x)不收敛到0，因此Rn[f]不收敛到0，即Q[f]不收敛，因此，对于节点较多的情况，需要使用分段线性或者Hermite插值。Runge现象其中*第34页，共71页，星期日，2025年，2月5日一、最高代数精度求积公式问题：结论：本节关键§2Gauss型求积公式设有n+1个节点，插值型求积公式的代数精度m的最大值？如何确定？由求积系数及n+1个节点xi，i=0,1,…,n的分布确定。*第35页，共71页，星期日，2025年，2月5日四个未知量A0，A1，x0，x1，已知插值型求积公式的代数精度最高。可按插值型求积公式来求A0，A1。解：具有尽可能高的代数精度。例求节点，使插值型求积公式(14)分析：x0，x1待定插值型求积的代数精度*第36页，共71页，星期日，2025年，2月5日于是，求积公式为*第37页，共71页，星期日，2025年，2月5日一般地，对于任意求积节点任意求积系数，求积公式对于代数精度结论*第38页，共71页，星期日，2025年，2月5日分析：只需证明使得事实上，令有从而而前例中m=3=2×1+1=2n+1是能达到的最高代数精度。次多项式f(x),*第39页，共71页，星期日，2025年，2月5日Gauss型求积公式的构造——利用正交多项式的根构造代数精度最高的求积公式分析：引理1：(15)*第40页，共71页，星期日，2025年，2月5日令f(x)是任意次数的代数多项式，则其中，q(x)为商式，r(x)为余式，均为任意次数的多项式证明：插值型求积公式代数精度大于n正交多项式基本性质ωn+1(x)是n+1次正交多项式，则*第41页，共71页，星期日，2025年，2月5日定义：正交多项式的根一定是Gauss点，那么Gauss点是否一定是正交多项式的根？n+1个节点(ax0…xnb)的求积公式(15)若其代数精度m=2n+1，即达到最高，称之为Gauss型求积公式,并称其节点为Gauss点。而*第42页，共71页，星期日，2025年，2月5日是[a,b]上关于权ρ(x)的n+1次正交多项式的根。二、Gauss点与正交多项式的关系定理4：分析：Gauss点ax0…xnb求积公式(15)是Gauss型的充分性即是引理1的结论。以下只证必要性只需证关于正交。必要性，即Gauss点作为节点正是n+1次正交多项式的根。*第43页，共71页，星期日，2025年，2月5日证明：取2n+1次多项式则有*第44页，共71页，星期日，2025年，2月5日三、Gauss求积公式的余项（截断误差）由引理1知,xi(i=0,1,…,n)是Gauss点,则m=2n+1,由n+1个点,确定2n+1次多项式。联想Hermite插值多项式。定理5：，则Gauss求积公式(15)的余项为分析：(16)*第45页，共71页，星期日，2025年，2月5日证明：若f(x)的Hermite插值多项式H2n+1(x)满足插值条件利用由积分中值定理*第46页，共71页，星期日，2025年