Chapt-5数值微积分的数值解法幻灯片.pptVIP

§5.7 数值微分法 三、实用的通用求数值微分 前两节分别介绍了连续函数和离散函数（表格函数）的数值微分计算方法，已经能够解决现代科学或工程技术领域或科研工作中所遇见的绝大部分微分问题。但是有些数值计算过程要求的精度较高，用前面介绍的方法求出的数值导数满足不了要求，所以本章节再介绍一种精度较高的连续函数通用数值微分方法。 1、算法分析 已知连续函数f=f(x)，求任意点 处的导数 ，可用下述方法： 此时 §5.7 数值微分法 此时 此时 则有 这里m是2～10之间的正整数，一般取m=10。 2、程序框图与通用程序设计 通用的数值微分法的程序框图见图5-11，下面是根据框图编写的通用程序： §5.7 数值微分法 subroutine d=0 do x1=a,b,h comp f(x),f(x1) END DO x1 END subroutine X=x1+0.5**i D=(f(x)-f(x1))/(x-x1) 图5-11 实用数值微分子程序框图 d1=d do i=9,10 END DO i §5.7 数值微分法 10 subroutine popder(a,b,h) 20 do x1=a,b, h 30 d=0 40 do i=9,10 50 d1=d 60 x=x1+0.5**i 70 d=(f(x)-f(x1))/(x-x1) 80 end do 90 end do 100 end subroutine popder 练习题 5-1 利用9节点定步长复化梯形与Simpson求积公式计算下列积分 5-2 已知函数y=f(x)的数值如下，利用差商或插值数值微分公式计算f’(1)。 x 0.6 0.8 0.9 1.0 1.1 1.2 1.4 y 0.736 0.8365 0.9095 1.0000 1.1105 1.2446 1.6017 上机实习题 5-3 用变步长Simpson求积公式与Romberg求积公式计算5-1题中积分。 §5.4 Newton-Cotes积分法 在梯形积分法中，是用一次二项式（直线）去逼近曲线f(x);在Simpson积分法中，是用二次三项式（抛物线）去逼近曲线f(x)。而在N-C积分法中，用的是m 次多项式去逼近曲线f(x)。 设在被积区间[a, b]选取m+1个插值点，用m 次多项式pm去逼近曲线f(x),则： 利用拉格郎日插值多项式 (5-9) §5.4 Newton-Cotes积分法 其中: Ck (m)为Cotes系数，可以具体求出起其值： m=1: C0(1) =1/2; C1(1)=1/2 m=2: C0(2)=1/6; C1(2)=4/6; C2(2)=1/6 见下表 ?注：在实际应用中，一般m9，否则会产生计算不稳定。 上面所谈到的问题，是在整个区间[a, b]上用Pm去逼近f(x)。下面看看用N-C公式即(5-1)式用在小区间时情况，将[a, b]分成n个小区间。有： §5.4 Newton-Cotes积分法 K m 0 1 2 3 4 5 6 7 8 1 1/2 1/2 ? ? ? ? ? ? ? 2 1/6 4/6 1/6 ? ? ? ? ? ? 3 1/6 3/8 3/8 1/6 ? ? ? ? ? 4 7/90 32/90 12/90 32/90 7/90 ? ? ? ? 5 19/288 25/96 25/144 25/144 25/96 19/288 ? ? ? 6 41/840 9/35 9/280 34/105 9/280 9/35 41/840 ? ? 7 751/17280 3577/17280 1323/17280 2989/17280 2989/17280 1323/17280 3577/17280 751/17280 ? 8 989/28350 5838/28350 -928/28350 10496/28350 -454