第四章数值积分与数值微分 - 2014年山东科技大学《数值分析》ppt课件.ppt

《第四章数值积分与数值微分 - 2014年山东科技大学《数值分析》ppt课件.ppt》由会员分享，可在线阅读，更多相关《第四章数值积分与数值微分 - 2014年山东科技大学《数值分析》ppt课件.ppt（86页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第四章第四章数值积分数值积分与数值微分与数值微分1 1 引引 言言一、数值积分的必要性一、数值积分的必要性本章主要讨论如下形式的一元函数积分本章主要讨论如下形式的一元函数积分在微积分里，按在微积分里，按Newton-Leibniz公式公式求定积分求定积分要求函数要求函数要求函数要求函数 的原函数的原函数的原函数的原函数 有解析表达式有解析表达式有解析表达式有解析表达式；为初等函数为初等函数为初等函数为初等函数实际问题实际问题1.1.的原函数的原函数的原函数的原函数 不能用初等函数表示不能用初等函数表示不能用初等函数表示不能用初等函数表示 建筑上用的一种铝制波纹瓦是用一种机器将一块平建筑上用的一

2、种铝制波纹瓦是用一种机器将一块平建筑上用的一种铝制波纹瓦是用一种机器将一块平建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的整的铝板压制而成的整的铝板压制而成的整的铝板压制而成的.假若要求波纹瓦长假若要求波纹瓦长假若要求波纹瓦长假若要求波纹瓦长4 4英尺英尺英尺英尺,每个波纹的高度每个波纹的高度每个波纹的高度每个波纹的高度(从中心线从中心线从中心线从中心线)为为为为1 1英寸英寸英寸英寸,且每个波纹以近似且每个波纹以近似且每个波纹以近似且每个波纹以近似 英寸为一个周英寸为一个周英寸为一个周英寸为一个周期期期期.求制做一块波纹瓦所需求制做一块波纹瓦所需求制做一块波纹瓦所需求制做一块

3、波纹瓦所需铝板的长度铝板的长度铝板的长度铝板的长度L.L.从从从从 到到到到 英寸间英寸间英寸间英寸间的的的的弧长弧长弧长弧长L L.这个问题就是要求由函数这个问题就是要求由函数给定的曲线给定的曲线给定的曲线给定的曲线,由微积分学我们知道由微积分学我们知道由微积分学我们知道由微积分学我们知道,所求的弧长可表示为所求的弧长可表示为所求的弧长可表示为所求的弧长可表示为:上述积分称为上述积分称为上述积分称为上述积分称为第二类椭圆积分第二类椭圆积分第二类椭圆积分第二类椭圆积分。Whats the Original function?!Its so complex that we can not get

4、 it.类似的，下列函数也不存在由初等函数表示的原函数类似的，下列函数也不存在由初等函数表示的原函数:2 2.有些被积函数其原函数虽然可以用初等函数表示有些被积函数其原函数虽然可以用初等函数表示有些被积函数其原函数虽然可以用初等函数表示有些被积函数其原函数虽然可以用初等函数表示,但表但表但表但表达式相当复杂达式相当复杂达式相当复杂达式相当复杂,计算极不方便计算极不方便计算极不方便计算极不方便.例如函数例如函数例如函数例如函数:并不复杂并不复杂并不复杂并不复杂,但它的原函数却但它的原函数却但它的原函数却但它的原函数却十分复杂十分复杂十分复杂十分复杂:3.3.没有解析表达式，只有数表形式没有解析表

5、达式，只有数表形式没有解析表达式，只有数表形式没有解析表达式，只有数表形式:1423454.5688.5原来通过原函数来计原来通过原函数来计算积分有它的局限性。算积分有它的局限性。那那怎么办呢？怎么办呢？呵呵呵呵这就需要积这就需要积分的数值方法来帮分的数值方法来帮忙啦。忙啦。二、数值积分的基本思想二、数值积分的基本思想1、定积分的几何意义、定积分的几何意义2、数值积分的理论依据、数值积分的理论依据依据依据积分中值定理积分中值定理,对于连续函数对于连续函数 ,在在 内存在一点内存在一点 ,使得使得称称 为为 在区间在区间 上的平均高度上的平均高度.3、求积公式的构造、求积公式的构造 若简单选取区

6、间端点或中点的函数值作为平均高度，则若简单选取区间端点或中点的函数值作为平均高度，则可得一点求积公式如下：可得一点求积公式如下：左矩形公式：左矩形公式：中矩形公式：中矩形公式：右矩形公式：右矩形公式：左矩形公式：左矩形公式：中矩形公式：中矩形公式：右矩形公式：右矩形公式：若取若取 两点，并令两点，并令 ，则可得梯形公式（两点求积公式），则可得梯形公式（两点求积公式）则可得则可得Simpson公式公式(三点求积公式三点求积公式)若取三点，若取三点，并令并令 一般地一般地，取区间，取区间 内内 个点个点处的高度处的高度通过通过加权平均加权平均的方法近似地得出平均高度的方法近似地得出平均高度这类求积

7、方法称为这类求积方法称为机械求积机械求积:或写成或写成:数值积分公式数值积分公式求积系数求积系数 求积节点求积节点 记记称称为数值为数值求积公式求积公式称为求积公称为求积公式余项式余项(误差误差).（1）（2）构造或确定一个求积公式，要解决的问题包括构造或确定一个求积公式，要解决的问题包括:(i)确定求积系数确定求积系数 和求积节点和求积节点 (iii)求积公式的误差估计和收敛性分析求积公式的误差估计和收敛性分析.(ii)确定衡量求积公式好坏的标准；确定衡量求积公式好坏的标准；称求积公式称求积公式 具有具有m次代数精度次代数精度,如果它满足如下两个条件如果它满足如下两个条件:定义定义4.1：(

8、i)对所有次数对所有次数m次的多项式次的多项式 ,有有(ii)存在存在m+1次多项式次多项式 ,使得使得三、求积公式的代数精度三、求积公式的代数精度上述定义中的条件上述定义中的条件(i),(ii)等价于等价于:注：梯形公式与中矩形公式都只具有注：梯形公式与中矩形公式都只具有1次代数精度。次代数精度。一般的，一般的，若要使求积公式（若要使求积公式（1 1）具有）具有m m次代数精度，则只要次代数精度，则只要使求积公式对使求积公式对f(x)=1，x，x2，,xm 都准确成立，即都准确成立，即2 2 插值型求积公式插值型求积公式一、定义一、定义在积分区间在积分区间 上，上，取取 个节个节点点作作 的

9、的 次代数插值多项式次代数插值多项式（拉格朗日插值公式）（拉格朗日插值公式）:则有则有其中，其中，为插值余项。为插值余项。于是有：于是有：取取Aj由由 节点节点 决定，决定，与与 无关。无关。称为称为插值插值型求积公型求积公式式二、截断误差与代数精度二、截断误差与代数精度1、截断误差、截断误差2、代数精度、代数精度 形如形如 的求积公式至少有的求积公式至少有 n 次次代数精度代数精度 该该公式为公式为插值型插值型（即：（即：）定理定理4.1推论推论 求积系数求积系数 满足满足:3 Newton-Cotes3 Newton-Cotes公式公式一、一、Cotes系数系数取取节点为节点为等距分布等距

10、分布：由此构造的插值型求积公式称为由此构造的插值型求积公式称为Newton-Cotes公式公式,此此时时求积系数：求积系数：令令Cotes系数系数二、二、Newton-Cotes公式公式1、定义：、定义：记记则则求积公式变为求积公式变为称上式为称上式为n阶阶闭型闭型Newton-Cotes求积公式。求积公式。注意注意:由式由式确定的确定的Cotes系数只与系数只与 和和 有关有关,与与 和积分区间和积分区间无关，无关，且且满足满足:2、截断误差、截断误差Newton-Cotes公式的误差为公式的误差为:与与x有关有关3、代数精度、代数精度作为插值型求积公式，作为插值型求积公式，具有具有 次代数

11、精度，次代数精度，阶阶Newton-Cotes公式至少公式至少而实际的代数精度是否可以进一步而实际的代数精度是否可以进一步提高呢？提高呢？定理定理 4.2当阶数当阶数 为偶数时为偶数时,Newton-Cotes公式公式至少至少具有具有次代数精度。次代数精度。证明证明:只需验证当只需验证当 为偶数时为偶数时,Newton-Cotes公式公式对对的余项为零。的余项为零。由于由于 ,所以所以 即得即得引进变换引进变换 ,因为因为 为偶数为偶数,故故 为整数为整数,于是有于是有据此可断定据此可断定 ,因为上述被积函数是个奇函数因为上述被积函数是个奇函数.4、数值稳定性、数值稳定性现在讨论现在讨论舍入误

12、差舍入误差对计算结果产生的影响对计算结果产生的影响.设用公式设用公式 近似计算积分近似计算积分时时,其中计算函数值其中计算函数值 有误差有误差则在则在 的计算中的计算中,由由 引起的引起的误差为误差为没有误差没有误差,中间计算过程中的舍入误差也不考虑中间计算过程中的舍入误差也不考虑,计算计算，而，而如果如果 都是正数都是正数,并设并设则有则有故故 是有界的是有界的,即由即由 引起的误差受到控制引起的误差受到控制,的的 倍倍,不超过不超过保证了保证了数值计算的稳定性数值计算的稳定性。将出现将出现负数负数,而当而当 时时,将随将随 增大增大,因而因而不能保证数值稳定不能保证数值稳定性性.故高阶公式

13、不宜采用故高阶公式不宜采用,有实用价值的仅仅是几种有实用价值的仅仅是几种低阶的求低阶的求积公式积公式.三、几种常用的低阶求积公式三、几种常用的低阶求积公式n=1:梯形公式梯形公式/*令令 x=a+th,h=b a,用积用积分中值定理分中值定理*/代数精度代数精度=1n=2:Simpson 公式公式代数精度代数精度=3其中，其中，n=4:代数精度代数精度=5,这里这里Cotes 公式公式其中，其中，四、复化求积公式四、复化求积公式 高次插值有高次插值有Runge 现象，怎么办？现象，怎么办？可采用分段低次插值来解决可采用分段低次插值来解决高阶高阶Newton-Cotes公式会出现公式会出现数值不

14、稳定数值不稳定。而而低阶低阶Newton-Cotes公式公式有时又不能满足精度要求有时又不能满足精度要求，怎么办？，怎么办？可将积分区间可将积分区间 分成若干小分成若干小区间，在每个小区间上用区间，在每个小区间上用低阶求积公式计算，然后求和。低阶求积公式计算，然后求和。复化梯形公式：复化梯形公式：在每个在每个 上用梯形公式：上用梯形公式：=Tn/*介值定理介值定理*/复化梯形公式的几何意义复化梯形公式的几何意义 复化复化 Simpson 公式：公式：44444=Sn在每个在每个 上用上用simpson公式：公式：复化复化SimpsonSimpson公式的几何意义公式的几何意义 复化复化 Cot

15、es公式：公式：=Cn例例:利用数据表利用数据表01/83/81/25/83/47/811/422.265492.460002.876403.200003.506853.764703.938464计算积分计算积分解：解：这个问题有明显的答案这个问题有明显的答案取取n=8用复化梯形公式用复化梯形公式=3.138988494取取n=4 用辛卜生公式用辛卜生公式=3.141592502运算量基运算量基本相同本相同复化梯形公式的误差估计复化梯形公式的误差估计给定精度给定精度 ，如何取，如何取?例如：要求例如：要求 ，如何判断，如何判断 n=？1、误差先验估计式、误差先验估计式若记若记则可令则可令由由其

16、中，其中，于是，于是，?上例中若要求上例中若要求 ，则，则即：取即：取 n=409通常采取将区间通常采取将区间不断对分不断对分的方法，即取的方法，即取 n=2k上例中上例中2k 409 k=9 时，时，T512=3.141592502S4=3.141592502注意到区间再次对分时注意到区间再次对分时可用来判断迭代可用来判断迭代是否停止。是否停止。2、误差后验估计式、误差后验估计式复化复化Simpson公式的误差估计公式的误差估计1、误差先验估计式、误差先验估计式2、误差后验估计式、误差后验估计式复化复化Cotes公式的误差估计公式的误差估计1、误差先验估计式、误差先验估计式2、误差后验估计式

17、、误差后验估计式 复化求积公式的复化求积公式的收敛速度（阶）收敛速度（阶）定义定义4.2：若一个复化求积公式的误差满足若一个复化求积公式的误差满足 ，且且 ，则，则称该公式是称该公式是 p 阶收敛阶收敛的。的。【注注】根据上述定义不难验证根据上述定义不难验证 龙贝格求积公式龙贝格求积公式例例:计算计算若取若取 =10 6，则利用复化求积公式进行计算时，须将区间对分则利用复化求积公式进行计算时，须将区间对分 9 次，得到次，得到 T512=3.141592502考察考察若由若由 来计算来计算 I 效果是否好些？效果是否好些？=3.141592502=S4一般有：一般有：Romberg 求积求积公

18、式公式4 Romberg 4 Romberg 算法算法 理查德森理查德森外推加速法外推加速法利用利用低低阶公式产生阶公式产生高高精度的结果。精度的结果。由由Taylor展开得到：展开得到：i 与与 h 无关无关现将现将 对分，得：对分，得：设对于某一设对于某一 ，有公式有公式 近似计算某一未知值近似计算某一未知值 。如何将公式精度由如何将公式精度由 提高到提高到？.432112)()(23322020 =hhIhTTh 即：即：复化梯形公式的渐近展开式复化梯形公式的渐近展开式定理定理4.34.3设设则成立则成立其中，其中，系数系数 与与 无关。无关。Romberg 算法：算法：?T1=)0(0

19、T T8=)3(0T T4=)2(0T T2=)1(0T S1=)0(1T R1=)0(3T S2=)1(1T C1=)0(2T C2=)1(2T S4=)2(1T加速公式加速公式Newton-Cotes公式采用公式采用等距节点作为求积节点代等距节点作为求积节点代数精度至多可达到数精度至多可达到 。（。（为偶数）为偶数）那么，在节点个数一定的情那么，在节点个数一定的情况下，是否可以在况下，是否可以在 上自由选择节点的位置，使上自由选择节点的位置，使求积公式的精度提得更高求积公式的精度提得更高？例例：求形如求形如的两点求积公式。的两点求积公式。（1）用梯形公式（即以用梯形公式（即以x0=-1，x

20、1=1为节点的插值型为节点的插值型 求积公式）立即可得求积公式）立即可得 。只具有只具有1 次代数次代数精确度！精确度！（2）若对求积公式中的四个待定系数若对求积公式中的四个待定系数A0,A1,x0,x1适当选取，适当选取，使求积公式对使求积公式对f(x)=1，x，x2，x3都准确成立，则都准确成立，则需满足如下方程组：需满足如下方程组：具有具有2n+1次代数精度的插值型求积公式次代数精度的插值型求积公式节点称为节点称为Gauss 点点称为称为Gauss 型求积公式型求积公式。注：注：Gauss型求积公式型求积公式是代数精度最高的是代数精度最高的插值型求积公式插值型求积公式.4 4 高斯型求积

21、公式高斯型求积公式定义定义4.3事实上，对于插值型求积公式事实上，对于插值型求积公式其代数精度最高可达到其代数精度最高可达到2n+1次次(Gauss型求积公式型求积公式）。）。考虑考虑2n+2次多项式次多项式 ，其中，其中，而而故故高斯型求高斯型求积积公式的构造公式的构造将节点将节点 以及系数以及系数 都作为待定系数。都作为待定系数。并令求积公式对并令求积公式对 精确成立精确成立可得非线性方程组可得非线性方程组1、待定系数法、待定系数法求解该方程组即可得求解该方程组即可得相应的求积节点与求相应的求积节点与求积系数。积系数。例：例：求求 的的 2 点点 Gauss 公式。公式。解：解：设设 ，应

22、有，应有 3 次代数精度。次代数精度。+101100)()()(xfAxfAdxxfx令上述公式对令上述公式对f(x)=1,x,x2,x3 精确成立可得精确成立可得不是线性方程组，不是线性方程组，不易求解。不易求解。定理：定理：x0 xn 为为 Gauss 点点 与与任任意意次次数数不不大于大于n 的多项式的多项式 P(x)（带权）正交带权）正交。证明：证明：“”x0 xn 为为 Gauss 点点,则公式则公式 至少有至少有 2n+1 次代数精度。次代数精度。对任意次数对任意次数不大于不大于n 的多项式的多项式 Pm(x)，Pm(x)w(x)的次数的次数不大于不大于2n+1，则代入公式应则代入

23、公式应精确成立精确成立：=00 求求 Gauss 点点 求求w(x)2、正交多项式法、正交多项式法不大于不大于 的多项式的多项式 精确成立，即证明：精确成立，即证明：“”要证明要证明 为为 Gauss 点，点，即要证公式对任意次数即要证公式对任意次数设设0 正交多项式族正交多项式族 0,1,n,有性质：任意次数不有性质：任意次数不大于大于n 的多项式的多项式 P(x)必与必与 n+1 正交。正交。若取若取 w(x)为其中的为其中的 n+1，则则 n+1的根的根就是就是 Gauss 点。点。53=a0)(10=+dxaxx0),(10=+=+=1021102200)(53(0),(0)(0),(

24、dxcbxxxxdxcbxxx 215910=cb即：即：Step 1：构造正交多项式构造正交多项式 2设设cbxxxaxxx+=+=2210)(,)(,1)(再解上例：再解上例：+101100)()()(xfAxfAdxxfxStep 2：求求 2=0 的的 2 个根，即为个根，即为 Gauss 点点 x0，x1Step 3：代入代入 f(x)=1,x 以求解以求解 A0，A1解解线性线性方程组，方程组，简单。简单。结果与前一方法相同：结果与前一方法相同：利用此公式计算利用此公式计算 的值的值注：注：构造正交多项式也可以利用构造正交多项式也可以利用 L-S 拟合中介绍过的递推拟合中介绍过的递

25、推式进行。式进行。几种常用的几种常用的Gauss型型求积公式：求积公式：Gauss-Legendre 求积公式：求积公式：1)(xr r定义在定义在 1,1上，上，满足：满足：，递推公式：，递推公式：其中，求积节点为其中，求积节点为 Pn+1 的根（求积系数通过解线性方程组的根（求积系数通过解线性方程组得到）。得到）。Legendre多项式：多项式：Gauss-Legendre 公式：公式：区间区间a,b上的上的Gauss-Legendre 公式：公式：其中，其中，为为 n+1次次Legendre多项式多项式Pn+1 的根。的根。Gauss-Chebyshev 求积公式：求积公式：211)(x

26、x=r r定义在定义在 1,1上，上，Tn+1 的根为：的根为：k=0,n以此为节点构造公式以此为节点构造公式称为称为 Gauss-Chebyshev 公式公式。注意到积分端点注意到积分端点 1 可能是积分可能是积分的的奇点奇点，用普通，用普通Newton-Cotes公式在端点会出问题。而公式在端点会出问题。而Gauss公式可能避免此问题的发生。公式可能避免此问题的发生。Chebyshev多项式：多项式：（3）第二类第二类Gauss-Chebyshev求积公式：求积公式：以以Un+1 的零点作的零点作为求积节点构造的公式为求积节点构造的公式称为称为第二类第二类Gauss-Chebyshev 公

27、式公式。第二类第二类Chebyshev多项式：多项式：区间区间-1,1上，带权上，带权(4)Gauss-Laguerre 求积公式：求积公式：以以n+1次次Laguerre 多项式的零点多项式的零点为求积节点构造的公式为求积节点构造的公式称为称为 Gauss-Laguerre 求积求积公式公式。Laguerre多项式：多项式：(5)Gauss-Hermite 求积公式：求积公式：以以n+1次次Hermite 多项式的零点多项式的零点为求积节点构造的公式为求积节点构造的公式称为称为 Gauss-Hermite 求积求积公式公式。Hermite 多项式：多项式：Gauss 公式的余项：公式的余项：

28、插值多项式插值多项式的余项的余项/*设设P为为f 的过的过x0 xn的插值多项式的插值多项式*/*只要只要P 的阶数不大于的阶数不大于2n+1，则下一步则下一步等式成立等式成立*/Hermite 多项式！多项式！什么样的什么样的插值多项式插值多项式在在 上有上有 阶？阶？Hermite 多项式的插值条件为：多项式的插值条件为：Hermite插值余项为插值余项为其中，其中，与与有关。有关。定理定理 若若f(x)f(x)Ca,bCa,b,则则GaussGauss型求积公式所求积分型求积公式所求积分 值序列值序列 收敛于积分值收敛于积分值I(fI(f),),即即 Gauss型求积公式的收敛性型求积公

29、式的收敛性 Gauss 型求积型求积公式的公式的数值稳定性数值稳定性v定理定理 GaussGauss型求积公式的求积系数都大于型求积公式的求积系数都大于零，从而零，从而GaussGauss型求积公式是型求积公式是数值稳定的。数值稳定的。复化复化Gauss-Legendre求积公式求积公式其中，其中，为为 n+1次次Legendre多项式多项式Pn+1 的根。从而的根。从而v复化两点复化两点Gauss-Gauss-LegendreLegendre求积公式为求积公式为其中，其中，v复化三点复化三点Gauss-Gauss-LegendreLegendre求积公式为求积公式为波纹瓦长度的计算波纹瓦长度的计算