第七章前半部分改完优秀PPT.ppt

《第七章前半部分改完优秀PPT.ppt》由会员分享，可在线阅读，更多相关《第七章前半部分改完优秀PPT.ppt（37页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第七章前半部分改完第七章前半部分改完现在学习的是第1页，共37页一一.研究数值解法的必要性研究数值解法的必要性 1 1 1 1 一般方程的根无法用解析表达式给出；一般方程的根无法用解析表达式给出；一般方程的根无法用解析表达式给出；一般方程的根无法用解析表达式给出；2 2 2 2 三次、四次方程的求根公式较繁。三次、四次方程的求根公式较繁。三次、四次方程的求根公式较繁。三次、四次方程的求根公式较繁。需要给出求根的近似值的方法。需要给出求根的近似值的方法。需要给出求根的近似值的方法。需要给出求根的近似值的方法。二二.方程的根方程的根 方程方程方程方程f f(x)=0 的解的解 称为方程称为方程f(

2、x x)=0)=0的根或称为的根或称为f(x x)的零点。的零点。的零点。的零点。若若f(x x)g g g g(x)(x)，其中，其中mm为正整数，为正整数，g(x x)满满足足 ，显然，显然 为为f(x x)的零点。这时，称的零点。这时，称 为为f(x)的的的的mm重零点，或称重零点，或称重零点，或称重零点，或称 为为为为f(x x)=0的的的的m重根。重根。现在学习的是第2页，共37页定理定理 若若f(x)具有具有m阶连续导数阶连续导数,则则 是是f(x)的的m重零点之充要条件为：重零点之充要条件为：证明证明 必要性必要性 设设 是是f(x)的的m重零点，则由重零点，则由 现在学习的是第

3、3页，共37页当当 时时当当 k=m时时现在学习的是第4页，共37页充分性充分性 设设 使得使得由由Taylor公式得公式得其中其中01，令，令 则有则有 且且根据定义，根据定义，为为f(x)的的m重零点重零点现在学习的是第5页，共37页三.根的搜索 求方程根的近似值之前，一般需要首先确定隔（有）根区间a,b（在a,b上方程仅有一个根）。方法：通过函数f(x)的增减性、凹凸性、变号特征等，并结合做草图来确定隔（有）根区间a,b。现在学习的是第6页，共37页 1 二分法（对分法）二分法（对分法）基本思想：基本思想：基本思想：基本思想：通过区间逐次分半，将有根区间逐通过区间逐次分半，将有根区间逐通

4、过区间逐次分半，将有根区间逐通过区间逐次分半，将有根区间逐步缩小。步缩小。步缩小。步缩小。设设设设f(x x)在在在在 a,b上连续，上连续，f(a a)f f(b b)0)0，且在，且在 a,b 内内f(x)=0 0至少有一个实根至少有一个实根至少有一个实根至少有一个实根 。记。记。记。记a,ba,b为为 ,计算计算 a1,b1 中点中点 的函数值的函数值 若若 0 0，则，则 若若若若 ，则，则，则，则 令令令令 若若 ，则，则 令令 新的有根区间新的有根区间 ,的长度的长度的长度的长度现在学习的是第7页，共37页再计算再计算 ,中点中点 的函数值的函数值若若 0 0，则，则若若 ，则，则

5、 令令若若 ，则，则 令令新的有根区间新的有根区间 ,的长度的长度如此对分下去，则得到一系列有根区间如此对分下去，则得到一系列有根区间且且现在学习的是第8页，共37页由由得得 k=1,2,3 当对分过程无限继续下去，则有根区间当对分过程无限继续下去，则有根区间必收缩为一点，即必收缩为一点，即具体做法具体做法(1)给定给定每步检查每步检查 是否成立，是否成立，若成立，取若成立，取 ，否则继续对分。，否则继续对分。(2)令令 ，先确定对分次数，先确定对分次数k，再计算，再计算(3)误差估计为误差估计为现在学习的是第9页，共37页优点 对函数性质要求不高（只要函数连续）；计算简单，且可达到任意精度。

6、缺点 计算量大；不能求复根与偶重根。现在学习的是第10页，共37页2 迭代法的算法和理论迭代法的算法和理论一不动点迭代法一不动点迭代法 对给定的方程对给定的方程对给定的方程对给定的方程f f(x x)=0，将其变为等价的方程，将其变为等价的方程，将其变为等价的方程，将其变为等价的方程 构造构造构造构造 k=1,2k=1,2k=1,2k=1,2 称为迭代序列，称为迭代序列，称为迭代序列，称为迭代序列，(x x)称为迭代函数。称为迭代函数。称为迭代函数。称为迭代函数。称为迭代格称为迭代格(公公)式或迭代过程。式或迭代过程。式或迭代过程。式或迭代过程。当当 (x)连续时，若连续时，若 则有则有 即即

7、 故序列故序列 的极限为方程的极限为方程x=(x x)(或或或或f f(x)=0)=0)的根的根现在学习的是第11页，共37页若若 满足满足 ，称称 为为 (x)的不动点。的不动点。即映射关系即映射关系即映射关系即映射关系 将将 映射到映射到 本身。本身。因因求求f f(x)的零点等价求的零点等价求 的不动点。的不动点。的不动点。的不动点。也称也称 k=1,2k=1,2k=1,2k=1,2 为不动点迭代法为不动点迭代法（简单迭代法或逐次逼近法）。（简单迭代法或逐次逼近法）。迭代序列的收敛性及收敛速度依赖于迭代序列的收敛性及收敛速度依赖于迭代函数的选取。迭代函数的选取。现在学习的是第12页，共3

8、7页二 不动点迭代法的一般理论 定理（不动点定理）已知x=(x)，若 且 对 ，有 ;存在常数0L1,使对 ，有则 (x)在a,b上有唯一的不动点 ；对任意 k=0,1,2 产生的序列 必定收敛到 (x)的不动点 ；有误差估现在学习的是第13页，共37页 证明证明证明证明 作辅助函数作辅助函数作辅助函数作辅助函数(x x)=)=(x x)-x x由于由于由于由于 (x x)在在在在 a,ba,b 上连续，上连续，上连续，上连续，则则则则(x x)C C a,ba,b，且且且且 (a a)=)=(a a)-a a0 0 (b b)=)=(b b)-b b0 0故由连续函数的介值定理，故由连续函数

9、的介值定理，故由连续函数的介值定理，故由连续函数的介值定理，至少存在至少存在至少存在至少存在使使使使()=0()=0()=0()=0 。即即即即从而从而从而从而 (x x)在在在在 a,ba,b 上存在不动点上存在不动点上存在不动点上存在不动点 。又设又设又设又设 (x x)有两个不动点有两个不动点有两个不动点有两个不动点 。注意注意注意注意 且且且且 由微分中值定理由微分中值定理由微分中值定理由微分中值定理 得得得得 即即即即故故故故 即即即即 (x x)在在在在 a,ba,b 上有唯一的不动点。上有唯一的不动点。上有唯一的不动点。上有唯一的不动点。现在学习的是第14页，共37页 注意注意

10、且且由微分中值定理由微分中值定理 其中其中 介于介于 与与 之间之间得得 k=1,2k=1,2k=1,2k=1,2 因因0L0，使得，使得 （C C称为渐近误差常数）称为渐近误差常数）则称该迭代过程为则称该迭代过程为p阶收敛。阶收敛。0 0c c11称为超线性收敛；称为超线性收敛；p p2 2称为平方收敛（二次收敛）。称为平方收敛（二次收敛）。称为平方收敛（二次收敛）。称为平方收敛（二次收敛）。p越大，收敛越快。越大，收敛越快。越大，收敛越快。越大，收敛越快。现在学习的是第21页，共37页 收敛速度的判别收敛速度的判别定理定理 设迭代格式设迭代格式 中迭代函数中迭代函数的高阶导数的高阶导数 （

11、p1）在不动点）在不动点 的邻域的邻域里连续，则迭代格式为里连续，则迭代格式为p阶收敛的充要条件是阶收敛的充要条件是 且有且有证明证明 充分性充分性 当当p1，若，若 则由则由 得得 即迭代格式为线性收敛。即迭代格式为线性收敛。现在学习的是第22页，共37页 对对p1，由，由 ，知迭代格式满足，知迭代格式满足在邻域具有局部收敛的条件，故迭代式收敛。在邻域具有局部收敛的条件，故迭代式收敛。由由 介于介于 之之 间间得得故故即即 为为p阶收敛。阶收敛。现在学习的是第23页，共37页必要性必要性 设迭代格式为设迭代格式为设迭代格式为设迭代格式为p p阶收敛，则有阶收敛，则有 即即 由于由于 在在在在

12、 邻域的连续性，知邻域的连续性，知邻域的连续性，知邻域的连续性，知 下面根据下面根据 证明：证明：（反证法）设有正整数（反证法）设有正整数 ，使，使 将将 在在 展开，得展开，得 即即现在学习的是第24页，共37页显然显然当当 由由 得得故故 这与迭代格式为这与迭代格式为p阶收敛盾。阶收敛盾。当当 由由 得得故故但迭代格式为但迭代格式为p阶收敛，应有阶收敛，应有即与迭代格式为即与迭代格式为p阶收敛产生矛盾。阶收敛产生矛盾。所以所以现在学习的是第25页，共37页3 迭代加速收敛的方法迭代加速收敛的方法 对于线性收敛的迭代格式，其收敛速度较对于线性收敛的迭代格式，其收敛速度较慢（特别是慢（特别是

13、），可以进行迭代加速。），可以进行迭代加速。一一.用两个迭代值进行组合用两个迭代值进行组合由由 ,得得 取取 和和则则迭代公式为迭代公式为或或当当 有有现在学习的是第26页，共37页Remark 对对 因因 收敛，在收敛，在 领域内，领域内，故改造后的迭代故改造后的迭代 有有 及及 若若 即即 这时，当这时，当L较大，加速收敛的效果明显。较大，加速收敛的效果明显。但当但当 有可能有可能 这时迭代不能加速收敛。这时迭代不能加速收敛。现在学习的是第27页，共37页对于线性迭代，有对于线性迭代，有对于线性迭代，有对于线性迭代，有通常取通常取通常取通常取 ，即，即，即，即则则则则故当故当故当故当 至少

14、为二阶收敛。至少为二阶收敛。至少为二阶收敛。至少为二阶收敛。实际中，常取实际中，常取实际中，常取实际中，常取 的近似值作为的近似值作为的近似值作为的近似值作为c c。得加速迭代公式得加速迭代公式得加速迭代公式得加速迭代公式即即即即缺点：需估计缺点：需估计缺点：需估计缺点：需估计 的近似值。的近似值。的近似值。的近似值。现在学习的是第28页，共37页二二 用三个迭代值的进行组合用三个迭代值的进行组合 设方程为设方程为 的某个预测值为的某个预测值为校正两次校正两次 由于由于 在在 与与 之间之间 在在 与与 之间之间在在 邻域中，消去导数信息，即邻域中，消去导数信息，即 即即因此可以得解下述因此可

15、以得解下述Aitken方法。方法。现在学习的是第29页，共37页AitkenAitken加速收敛迭代格式加速收敛迭代格式给出给出给出给出校正校正校正校正 再校正再校正再校正再校正改进改进或或记记 上式可写为：上式可写为：称为称为称为称为SteffensenSteffensen迭代格式迭代格式。特点：不需计算导数值；三个迭代值进行组合。特点：不需计算导数值；三个迭代值进行组合。现在学习的是第30页，共37页 对于对于对于对于 SteffensenSteffensen 迭代迭代迭代迭代，有：有：有：有：定理定理定理定理 设函数设函数 有不动点有不动点有不动点有不动点 ，且在，且在，且在，且在 邻域

16、内具有二阶连续导数，邻域内具有二阶连续导数，邻域内具有二阶连续导数，邻域内具有二阶连续导数，若若若若 且且且且则则 与与 有相同的不动点；有相同的不动点；Steffensen迭代为二阶收敛，且极限为迭代为二阶收敛，且极限为迭代为二阶收敛，且极限为迭代为二阶收敛，且极限为证明证明 设设Steffensen迭代式为迭代式为则则则则现在学习的是第31页，共37页 验证验证验证验证 与与与与 有相同的不动点有相同的不动点有相同的不动点有相同的不动点 若若若若 则则则则 即即即即 反之，若反之，若反之，若反之，若 注意注意注意注意 则则则则现在学习的是第32页，共37页即即注意注意 的连续性，知的连续性

17、，知故故 与与 的根相同。的根相同。即即 与与 有共同的不动点。有共同的不动点。现在学习的是第33页，共37页 用定义证明用定义证明Steffensen迭代为二阶收敛迭代为二阶收敛 由由 有二阶连续导数，且有二阶连续导数，且 为为 的不动点，可得关的不动点，可得关系式：系式：故故 现在学习的是第34页，共37页由由 得（假设）得（假设）即即即即 SteffensenSteffensen迭代为二阶收敛迭代为二阶收敛迭代为二阶收敛迭代为二阶收敛 现在学习的是第35页，共37页RemarksRemarks二阶收敛的证明也可以二阶收敛的证明也可以二阶收敛的证明也可以二阶收敛的证明也可以 证明。证明。上

18、述定理的证明实际上用了迭代收敛的假设。严格的说，应上述定理的证明实际上用了迭代收敛的假设。严格的说，应上述定理的证明实际上用了迭代收敛的假设。严格的说，应上述定理的证明实际上用了迭代收敛的假设。严格的说，应该首先证明收敛性（如：在证明该首先证明收敛性（如：在证明该首先证明收敛性（如：在证明该首先证明收敛性（如：在证明 的基础上，证明的基础上，证明的基础上，证明的基础上，证明 ）再证明二阶收敛。再证明二阶收敛。当当当当 线性收敛。线性收敛。线性收敛。线性收敛。当当 发散。发散。但但但但Steffensen迭代不仅能加速收敛，而且也能将发散的迭代不仅能加速收敛，而且也能将发散的迭代改进为收敛的迭代格式迭代改进为收敛的迭代格式现在学习的是第36页，共37页Steffensen迭代加速技术一般不对高阶收敛的迭代迭代加速技术一般不对高阶收敛的迭代法进行。法进行。严格的讲，严格的讲，称为称为Aitken加速迭代格式。加速迭代格式。对于具有线性收敛的任意序列对于具有线性收敛的任意序列 ，可以，可以证明证明 比比 收敛快。收敛快。只要当只要当 由不动点迭代由不动点迭代 产生，得到的公式称为产生，得到的公式称为Steffensen公式。公式。现在学习的是第37页，共37页