《数值分析第四章数值积分与数值微分ppt课件》由会员分享，可在线阅读，更多相关《数值分析第四章数值积分与数值微分ppt课件（74页珍藏版）》请在金锄头文库上搜索。
1、第四章数值积分数值积分与数值微分与数值微分:1 1 引引 言言一、数值积分的必要性一、数值积分的必要性本章主要讨论如下形式的一元函数积分本章主要讨论如下形式的一元函数积分在微积分里，按在微积分里，按Newton-Leibniz公式求定积分公式求定积分要求被积函数要求被积函数要求被积函数要求被积函数 有解析表达式；有解析表达式；有解析表达式；有解析表达式； 的原函数的原函数的原函数的原函数 为初等函数为初等函数为初等函数为初等函数:实际问题实际问题1. 1. 的原函数的原函数的原函数的原函数 不能用初等函数表示不能用初等函数表示不能用初等函数表示不能用初等函数表示例如函数例如函数:考虑一个实际问
2、题考虑一个实际问题考虑一个实际问题考虑一个实际问题: :建筑上用的一种铝制波纹瓦是用一种机器将一块平整的建筑上用的一种铝制波纹瓦是用一种机器将一块平整的建筑上用的一种铝制波纹瓦是用一种机器将一块平整的建筑上用的一种铝制波纹瓦是用一种机器将一块平整的铝板压制而成的铝板压制而成的铝板压制而成的铝板压制而成的. .:假若要求波纹瓦长假若要求波纹瓦长假若要求波纹瓦长假若要求波纹瓦长4 4英尺英尺英尺英尺, , 每个波纹的高度每个波纹的高度每个波纹的高度每个波纹的高度( (从中心线从中心线从中心线从中心线) )为为为为1 1英寸英寸英寸英寸, ,且每个波纹以近似且每个波纹以近似且每个波纹以近似且每个波纹
3、以近似 英寸为一个周期英寸为一个周期英寸为一个周期英寸为一个周期. . 求制做一块波纹瓦所需求制做一块波纹瓦所需求制做一块波纹瓦所需求制做一块波纹瓦所需铝板的长度铝板的长度铝板的长度铝板的长度L.L.从从从从 到到到到 英寸间的弧长英寸间的弧长英寸间的弧长英寸间的弧长L.L.这个问题就是要求由函数这个问题就是要求由函数给定的曲线给定的曲线给定的曲线给定的曲线, ,: 由微积分学我们知道由微积分学我们知道,所求的弧长可表示为所求的弧长可表示为:上述积分称为第二类椭圆积分。上述积分称为第二类椭圆积分。上述积分称为第二类椭圆积分。上述积分称为第二类椭圆积分。Whats the Original fu
4、nction?!Its so complex that we can not get it.:2. 2. 有些被积函数其原函数虽然可以用初等函数表示成有限有些被积函数其原函数虽然可以用初等函数表示成有限有些被积函数其原函数虽然可以用初等函数表示成有限有些被积函数其原函数虽然可以用初等函数表示成有限形式形式形式形式, ,但表达式相当复杂但表达式相当复杂但表达式相当复杂但表达式相当复杂, ,计算极不方便计算极不方便计算极不方便计算极不方便. .例如函数例如函数例如函数例如函数: :并不复杂并不复杂并不复杂并不复杂, , , ,但它的原函数却十分复杂但它的原函数却十分复杂但它的原函数却十分复杂但它的
5、原函数却十分复杂: : : :3. 3. 没有解析表达式，只有数表形式没有解析表达式，只有数表形式没有解析表达式，只有数表形式没有解析表达式，只有数表形式: :1423454.5688.5原来通过原函数来计原来通过原函数来计算积分有它的局限性。算积分有它的局限性。那那怎么办呢？怎么办呢？呵呵呵呵这就需要积这就需要积分的数值方法来帮分的数值方法来帮忙啦。忙啦。:二、数值积分的基本思想二、数值积分的基本思想1、定积分的几何意义、定积分的几何意义:2、数值积分的理论依据、数值积分的理论依据依据积分中值定理依据积分中值定理, 对于连续函数对于连续函数 ,在在 内存在一点内存在一点 ,使得使得称称 为区
6、间为区间 的平均高度的平均高度.:3、求积公式的构造、求积公式的构造 若简单选取区间端点或中点的函数值作为平均高度，则若简单选取区间端点或中点的函数值作为平均高度，则可得一点求积公式如下：可得一点求积公式如下：左矩形公式：左矩形公式：中矩形公式：中矩形公式：右矩形公式：右矩形公式：:左矩形公式：左矩形公式：:中矩形公式：中矩形公式：:右矩形公式：右矩形公式：: 若取若取 两点，并令两点，并令 ，则可得梯形公式两点求积公式），则可得梯形公式两点求积公式）:则可得则可得Simpson公式公式(三点求积公式三点求积公式) 若取三点，若取三点， 并令并令 : 一般地一般地 ，取区间，取区间 内内 个点
7、个点处的高度处的高度通过加权平均的方法近似地得出平均高度通过加权平均的方法近似地得出平均高度这类求积方法称为机械求积这类求积方法称为机械求积: 或写成或写成: :数值积分公式数值积分公式求积系数求积系数 求积节点求积节点 :记记称为数值称为数值求积公式求积公式称为求积公称为求积公式余项式余项(误差误差).:三、求积公式的代数精度三、求积公式的代数精度1、问题的提出、问题的提出构造或确定一个求积公式，要讨论解决的问题有构造或确定一个求积公式，要讨论解决的问题有:(i) 确定求积系数确定求积系数 和求积节点和求积节点 (iii)求积公式的误差估计和收敛性分析求积公式的误差估计和收敛性分析.(ii)
8、 判定求积公式精度的衡量标准；判定求积公式精度的衡量标准；: 称求积公式称求积公式 具有具有m次代数精度次代数精度,如果它满足如下两个条件如果它满足如下两个条件:2、定义、定义(i) 对所有次数对所有次数m次的多项式次的多项式 ,有有(ii)存在存在m+1次多项式次多项式 ,使得使得:上述定义中的条件上述定义中的条件(i),(ii)等等价于价于:2 2 插值型求积公式插值型求积公式一、定义一、定义在积分区间在积分区间 上，上， 取取 个节个节点点作作 的的 次代数插值多项式拉格朗日插值公式）次代数插值多项式拉格朗日插值公式）:则有则有其中，其中，为插值余项。为插值余项。:于是有：于是有：取取A
9、k由由 节点节点 决定，决定，与与 无关。无关。称为插值称为插值型求积公型求积公式式:二、截断误差与代数精度二、截断误差与代数精度1、截断误差、截断误差:2、代数精度、代数精度推论推论 求积系数求积系数 满足满足: 形如形如 的求积公式至少有的求积公式至少有 n 次次代数精度代数精度 该公式为插值型即：该公式为插值型即： ）定理定理:3 Newton-Cotes3 Newton-Cotes公式公式一、一、Cotes系数系数取节点为等距分布：取节点为等距分布：由此构造的插值型求积公式称为由此构造的插值型求积公式称为Newton-Cotes公式公式,此此时求积系数：时求积系数：令令Cotes系数系
10、数:二、二、Newton-Cotes公式公式1、定义：、定义：记记那那么么求积公式变为求积公式变为称上式为称上式为n阶闭型阶闭型Newton-Cotes求积公式。求积公式。:注意注意:由式由式确定的确定的Cotes系数只与系数只与 和和 有关有关, 与与 和积分区间和积分区间无关，无关， 且满足且满足:2、截断误差、截断误差Newton-Cotes公式的误差为公式的误差为:与与x x有关有关:3、代数精度、代数精度作为插值型求积公式，作为插值型求积公式，具有具有 次代数精度，次代数精度，阶阶Newton-Cotes公式至少公式至少而实际的代数精度是否可以进一步而实际的代数精度是否可以进一步提高
11、呢？提高呢？定理定理当阶数当阶数 为偶数时为偶数时,Newton-Cotes公式至少具有公式至少具有次代数精度。次代数精度。:证明证明:只需验证当只需验证当 为偶数时为偶数时,Newton-Cotes公式公式对对的余项为零。的余项为零。由于由于 ,所以所以 即得即得引进变换引进变换 ,因为因为 为偶数为偶数,故故 为整数为整数,于是有于是有据此可断定据此可断定 ,因为上述被积函数是个奇函数因为上述被积函数是个奇函数.:4、数值稳定性、数值稳定性现在讨论舍入误差对计算结果产生的影响现在讨论舍入误差对计算结果产生的影响.设用公式设用公式 近似计算积分近似计算积分时时, 其中计算函数值其中计算函数值
12、 有误差有误差则在则在 的计算中的计算中,由由 引起的引起的误差为误差为没有误差没有误差, 中间计算过程中的舍入误差也不考虑中间计算过程中的舍入误差也不考虑,计算计算，而，而:假如假如 都是正数都是正数,并设并设则有则有故故 是有界的是有界的, 即由即由 引起的误差受到控制引起的误差受到控制,的的 倍倍,不超过不超过保证了数值计算的稳定性。保证了数值计算的稳定性。将出现负数将出现负数,而当而当 时时,将随将随 增大增大,因而不能保证数值稳定因而不能保证数值稳定性性.故高阶公式不宜采用故高阶公式不宜采用,有实用价值的仅仅是几种低阶的求有实用价值的仅仅是几种低阶的求积公式积公式.:三、几种常用的低
13、阶求积公式三、几种常用的低阶求积公式n = 1:梯形公式梯形公式/* 令令 x = a+th, h = b a, 用中用中值定理值定理 */代数精度代数精度 = 1:n = 2:Simpson 公式公式代数精度代数精度 = 3:n = 4: Cotes 公公式式 代数精度代数精度 = 5,这里这里:四、复化求积公式四、复化求积公式 高次插值有高次插值有Runge 现象，怎么办？现象，怎么办？可采用分段低次插值来解决可采用分段低次插值来解决高阶高阶Newton-Cotes公式会出现公式会出现数值不稳定。数值不稳定。而低阶而低阶Newton-Cotes公式公式有时又不能满足精度要求，怎么办？有时又
14、不能满足精度要求，怎么办？可将积分区间可将积分区间 分成若干小分成若干小区间，在每个小区间上用区间，在每个小区间上用低阶求积公式计算，然后求和。低阶求积公式计算，然后求和。: 复化梯形公式：复化梯形公式：在每个在每个 上用梯形公式：上用梯形公式：= Tn/*中值定理中值定理*/: 复化梯形公式积分法: 复化复化 Simpson 公式：公式：44444= Sn: 复化Simpson公式积分法: 复化复化 Cotes公式：公式：= Cn: 收敛速度与误差估计：收敛速度与误差估计：定义：定义：若一个积分公式的误差满足若一个积分公式的误差满足 ，且且 ，则称该公式是，则称该公式是 p 阶收敛的。阶收敛
15、的。:例例:利用数据表利用数据表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.14159202S4 = 3.141592502注意到区间再次对分时注意到区间再次对分时可用来判断迭代可用来判断迭代是否停止。是否停止。2、误差后验估计式、误差后验估计式:复化复化Simpson公式的误差估计公式的误差估计1、误差先验估计式、误差先验估计式2、误差后验估计式、误差后验估计式:复化复化Cotes公式的误差估计公式的误差估计1、误差先验估计式、误差先验估计式2、误差后验估计式、误差
17、后验估计式:四、龙贝格积分四、龙贝格积分例例:计算计算已知对于已知对于 = 10 6 须将区间对分须将区间对分 9 次，得到次，得到 T512 = 3.14159202考察考察由由 来计算来计算 I 效果是否好些效果是否好些？= 3.141592502= S4一般有：一般有：Romberg求积求积公式公式: Romberg 算法：算法： ? ? ? T1 =)0(0T 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: 理查德森外推法理查德森外推法利用低阶公式产生高
18、精度的结果。利用低阶公式产生高精度的结果。由由Taylor展开得到：展开得到： i 与与 h 无关无关现将现将 对分，得：对分，得：设对于某一设对于某一 ， 有公式有公式 近似计算某一未知值近似计算某一未知值 。:如何将公式精度由如何将公式精度由 提高到提高到 ？.432112)()(23322020 = = hhIhTTh 即：即：:计算步骤：计算步骤：1取 ，计算2对对k = 1, 2, 计算下列各步计算下列各步:3对对n = 0, 1, 2, k = n 1, n 2, 4收敛控制收敛控制假设假设或或则输出积分值则输出积分值 ，否则转，否则转3 3。 :Newton-Cotes公式采用公
19、式采用等距节点作为求积节点代等距节点作为求积节点代数精度至多可达到数精度至多可达到 。（。（ 为偶数）为偶数）那么，在节点个数一定的情那么，在节点个数一定的情况下，是否可以在况下，是否可以在 上自由选择节点的位置，使上自由选择节点的位置，使求积公式的精度提得更高求积公式的精度提得更高 ？:例例 ：求形如求形如的两点求积公式。的两点求积公式。 （1用梯形公式即以用梯形公式即以x0 = -1，x1 = 1为节点的插值型为节点的插值型 求积公式立即可得求积公式立即可得 。只具有一次代数只具有一次代数精确度！精确度！:（2若对求积公式中的四个待定系数若对求积公式中的四个待定系数A0, A1, x0,
20、x1适当选取，适当选取，使求积公式对使求积公式对f (x) = 1f (x) = 1，x x，x2x2，x3x3都准确成立，那么都准确成立，那么需满足如下方程组：需满足如下方程组：:五、高斯型积分五、高斯型积分构造具有构造具有2n+1次代数精度的求积公式次代数精度的求积公式将节点将节点 以及系数以及系数 都作为待定系数。都作为待定系数。令令 代入可求解，代入可求解，得到的公式得到的公式具有具有 次代数精次代数精度。度。节点称为节点称为Gauss 点点此公式称为此公式称为Gauss 型求型求积公式积公式:例：求例：求 的的 2 点点 Gauss 公式。公式。解：设解：设 ，应有，应有 3 次代数
21、精度。次代数精度。 + + 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，则代入公式应精确成立：，则代入公式应精确成
22、立：= 00 求求 Gauss 点点 求求w(x):不大于不大于 的多项式的多项式 精确成立，即证明：精确成立，即证明：“”要证明要证明 为为 Gauss 点，点， 即要证公式对任意次数即要证公式对任意次数设设0 : 正交多项式族正交多项式族 0, 1, , n, 有性质：有性质：任意次数不大于任意次数不大于n 的多项式的多项式 P(x) 必与必与n+1 正交。正交。若取若取 w(x) 为其中的为其中的n+1，那么，那么n+1的根就是的根就是 Gauss 点。点。:53 = =a0)(10= =+ + dxaxx0),(10= =j jj j = =+ + + = = =+ + += =102
23、1102100)(53(0),(0)(0),(dxcbxxxxdxcbxxxj jj jj jj j215910= = = =cb即：即：Step 1：构造正交多项式：构造正交多项式2设设cbxxxaxxx+ + += =+ += = =2210)(,)(, 1)(j jj jj j再解上例：再解上例： + + 101100)()()(xfAxfAdxxfx:Step 2：求：求2 = 0 的的 2 个根，即为个根，即为 Gauss 点点 x0 ，x1Step 3：代入：代入 f (x) = 1, x 以求解以求解 A0 ，A1解线性方程组，解线性方程组，简单。简单。结果与前一方法相同：结果与
24、前一方法相同： 利用此公式计算利用此公式计算 的的值值注：构造正交多项式也可以利用注：构造正交多项式也可以利用 L-S 拟合中介绍过的递推拟合中介绍过的递推式进行。式进行。: 特殊正交多项式族：特殊正交多项式族： Legendre 多项式族：多项式族：1)( xr r定义在定义在 1, 1上，上，满足：满足：由由 有递推有递推以以 Pn+1 的根为节点的求积公式称为的根为节点的求积公式称为Gauss-Legendre 公式。公式。: Chebyshev 多项式族：多项式族：211)(xx = =r r定义在定义在 1, 1上，上，Tn+1 的根为的根为k = 0, , n以此为节点构造公式以此
25、为节点构造公式称为称为 Gauss-Chebyshev 公式。公式。注意到积分端点注意到积分端点 1 可能是积分可能是积分的奇点，用普通的奇点，用普通Newton-Cotes公公式在端点会出问题。而式在端点会出问题。而Gauss公公式可能避免此问题的发生。式可能避免此问题的发生。: Gauss 公式的余项：公式的余项：插值多项式的余项插值多项式的余项/* 设设P为为f 的过的过x0 xn的插值多项式的插值多项式 */*只要只要P 的阶数不大于的阶数不大于2n+1，则下一步，则下一步等式成立等式成立*/:Hermite 多项式！多项式！什么样的插值多项式什么样的插值多项式在在 上有上有 阶？阶？:Hermite 多项式的插值条件为：多项式的插值条件为：插值余项为插值余项为其中，其中，与与有关。有关。:Hermite求积公式的余项求积公式的余项:
