计算声学 非线性方程的数值解法.pptx

《计算声学 非线性方程的数值解法.pptx》由会员分享，可在线阅读，更多相关《计算声学 非线性方程的数值解法.pptx（115页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第三章 非线性方程的数值解法 对于代数方程，根的数目与方程的次数相同，不高于4次的代数方程已有求根公式，高于4次的代数方程没有精确的求根公式；超越方程则复杂得多，如果有解，可能是一个或多个，或无穷多个，没有精确的求解公式。因此，需要研究如何采用一定的数值算法求得满足一定精度要求的近近似根似根。内容提要内容提要：隔根区间的确定、二分法、迭代法、牛顿迭代法、弦截法重点内容重点内容：迭代法、牛顿迭代法第1页/共115页第三章 非线性方程的数值解法数值方法求方程 根的近似值，要解决三个问题：1、根的存在性存在性：方程有没有根，如果有根，有几个；2、根的隔离隔离：找出有根区间，把有根区间分成较小的子区

2、间，每个子区间只有一个根（隔根区间）；3、根的精确化精确化：确定了隔根区间后，可以用各种方法将某 一近似根 逐步精确化。按照一定的方法产生一个序列 ，此序列在一定的条件下收敛于方程 的根 。产生序列 的不同方法就构成了不 同的方程求根方法。第2页/共115页第三章 非线性方程的数值解法应用举例：本征声线求解本征声线求解本征声线本征声线：从声源出发经过一定的传播路径到达接收点的声线，接收点处的声场是所有本征声线能量叠加的结果。折射定律折射定律(Snell)：第3页/共115页第三章 非线性方程的数值解法对于给定的海洋环境，每条声线的轨迹由声线的起始掠射角唯一决定。设声源位于 处，接收点位于 处，

3、声线轨迹方程：需要导出声线以掠射角 从声源出发，经过一定水平距离 到达深度 的计算公式，有改变 ，给定精度要求，进行求解。第4页/共115页第三章 非线性方程的数值解法第5页/共115页1 根的搜索与二分法 根的搜索根的搜索定义（隔根区间）：如果在区间 内只有方程的一个根，则称区间 为隔根区间。隔根区间的确定描图法逐步搜索法第6页/共115页介值定理介值定理：定理定理1 1：设函数 在区间 上连续，且 ，则方程 在 上至少有一个根。定理定理2 2：设函数 在区间 上是单调连续函数，且 ，则方程 在 有且仅有一个根。1 根的搜索与二分法 第7页/共115页1 根的搜索与二分法描图法描图法：画出

4、的简图，由曲线与 轴的交点位置确定出隔根区间。或者将方程等价变形为 ，画出函数 和 的简图，从两条曲线交点的横坐标位置确定隔根区间。例：求方程 的隔根区间。解：由图可知，方程仅有一个实根，隔根区间为 。第8页/共115页1 根的搜索与二分法逐步搜索法：逐步搜索法：首先确定方程 的实根所在区间 再按照选定的步长 （n为正整数），逐点计算 处的函数值 ，当 与异号时，则 即为方程的一个隔根区间。对于 次代数方程其根的绝对值的上下界有如下结论：（1）如果 ，则方程的根的绝对值小于 ；（2）如果 ，则方程根的绝对值大 于 。第9页/共115页1 根的搜索与二分法例：求方程 的隔根区间。解：利用逐步搜索

5、法设方程的根为 ，则 即根的所在区间为 和 。取 计算 ，隔根区间为-4.2-0.7-0.20.20.71.21.72.24.2-137.72-2.440.281.060.910.20-0.310.1426.42第10页/共115页1 根的搜索与二分法二分法二分法基本思想基本思想：通过计算隔根区间的中点，逐步缩小隔根区间，从而得到方程的近似根序列 。过程过程：设 为连续函数，方程 的隔根区间为，则函数值 异号。（1）将区间 二分得到中点 ，计算函数值 ；（2）根据计算得到的函数值进行判断，如果 ，则方程的根为 ；否则判断 的符号，如果其与 异号，则隔根区间变为 ；如果其与 异号，则隔根区间变为

6、 ；第11页/共115页1 根的搜索与二分法（3）把新的隔根区间记为 ，继续按照上述规则对隔根区间 进行二分。满足精度要求则停止。通过不断的二分，就得到一系列隔根区间有 。由于后一区间的长度都是前一区间长度的一半，所以 的长度为当 时，区间 的长度趋于零，即区间最终会收敛于方程 的解 。实际计算中，只要二分的次数足够多，可取最后区间中点 作为方程根的近似值，即第12页/共115页1 根的搜索与二分法误差误差：如果事先给定精度要求为 ，则满足时，停止计算。第13页/共115页1 根的搜索与二分法例例：用二分法求方程 在 内的根的近似值，要求绝对误差不超过 。解解：即 严格单调增加，又 ，所以方程

7、在 上有唯一实根。令 ，得到 ，取 ，即至少二分7次。计算过程如下：第14页/共115页1 根的搜索与二分法function f=f3(x)f=x3+4*x*x-10;a0=1;b0=2;k=0;while 1 k=k+1;x=(a0+b0)/2;if f3(x)*f3(a0)0 b0=x;else a0=x;end if(b0-a0)/20.005 break;endendx=(b0+a0)/2;第15页/共115页1 根的搜索与二分法隔根区间1 11.51.52.3752.3751,1.51,1.52 21.251.25-1.7969-1.79691.25,1.51.25,1.53 31.

8、3751.3750.16210.16211.25,1.3751.25,1.3754 41.31251.3125-0.8484-0.84841.3125,1.3751.3125,1.3755 51.343751.34375-0.3510-0.35101.34375,1.3751.34375,1.3756 61.3593751.359375-0.0964-0.09641.359375,1.3751.359375,1.3757 71.36718751.36718750.03240.03241.359375,1.36718751.359375,1.3671875 所以最后得 第16页/共115页1 根

9、的搜索与二分法例例：用二分法求方程 在 内的根，要求精确到小数点后第三位小数，需要二分多少次？解解：设 ，由于 ，所以在区间 内方程 有唯一实根。令 ，求得所需对分次数至少是10次。特点特点：运算简单，方法可靠，对函数只要求在区间上连续；但收敛速度慢，不能用来求复数根及偶数重根。常用于为其它求根方法提供较好的近似初始值。第17页/共115页2 迭代法及其迭代收敛的加速方法 迭代法（逐次逼近）迭代法（逐次逼近）基本思想基本思想：利用某种递推算式，使某个预知的近似根（初值）逐次精确化，直到得到满足精度要求的近似根。算法算法：给定方程 ，其中 在有根区间 上连续。设 是方程的一个近似根，将方程 改写

10、为，代入初值 构造迭代序列第18页/共115页2 迭代法及其迭代收敛的加速方法这种求方程近似根的方法称为迭代法迭代法（或简单迭代法），称为迭代函数迭代函数。如果由迭代法产生的迭代序列 极限存在，即 ，则称 收敛，否则称 发散。称为迭代函数 的不动点不动点，或是方程 的解。由 转化为 时，迭代函数 不是唯一的，不同，会产生不同的序列 ，从而收敛情况也不一样。第19页/共115页2 迭代法及其迭代收敛的加速方法几何意义几何意义：求方程 的根 ，在几何上就是求直线 与曲线 交点 的横坐标，如图所示。从图中可以看出，当迭代函数迭代函数 的导数导数 在根 处满足不同条件时，迭代过程的收敛情况也有所不同。

11、所以迭代过程的收敛依赖于迭代函数 的构造，为使迭代法有效，必须保证其收敛性。第20页/共115页2 迭代法及其迭代收敛的加速方法第21页/共115页2 迭代法及其迭代收敛的加速方法第22页/共115页2 迭代法及其迭代收敛的加速方法例：已知方程 在 内有一个根，用两种不同的迭代公式（1）；（2）进行迭代，观察所得序列的收敛性。解：计算结果如表所示k （1）（2）00.30.31-0.00470.36172-1.01080.37323-1.90250.37534-1.98750.3757第23页/共115页2 迭代法及其迭代收敛的加速方法x0=0.3;k=0;while 1 k=k+1;xk=1

12、0 x0-2;%公式1%xk=log10(x0+2)%公式2 if abs(xk-x0)=0.0005 break;end x0=xk;end第24页/共115页2 迭代法及其迭代收敛的加速方法例：对于方程 ，试采用如下迭代法求方程的根(1)；(2)；(3)选取初值 。解：计算结果见下表 (1)(2)(3)01.501.501.50110.39061.16500.151721.2585e+0061.1374-1.000131.1350-0.499941.1347-0.969751.1347-0.53846-0.95677-0.5551第25页/共115页x0=1.5;k=0;while 1 k

13、=k+1;xk=x06-1;%公式1，发散%xk=(x0+1)(1/6);%公式2，收敛%xk=1/(x05-1);%公式3，发散 if abs(xk-x0)=0.0005 break;end x0=xk;end2 迭代法及其迭代收敛的加速方法第26页/共115页2 迭代法及其迭代收敛的加速方法迭代法收敛定理迭代法收敛定理：设方程 ，若迭代函数 在有根区间 上满足：（1）当 时，；（2）在 上可导，且有 ，；则有（1）方程 在 上有唯一的根 ；（2）对任意初值 ，迭代公式 产生的数列 收敛于方程的唯一根 ，即 ；（3）误差估计第27页/共115页2 迭代法及其迭代收敛的加速方法拉格朗日中值定理

14、拉格朗日中值定理：如果函数 在闭区间 上连续，在开区间 内可导，则至少存在一点 使得第28页/共115页2 迭代法及其迭代收敛的加速方法证明证明：（1）根的存在性存在性：令 ，则 在 上连续，由条件1有由连续函数的介值定理，必存在 使得 ，即为方程的根。根的唯一性唯一性：假设另有 也是方程的根，则有两式相减并应用微分中值定理得第29页/共115页2 迭代法及其迭代收敛的加速方法式中 在 和 之间，所以 ，由条件（2）有上式不成立，所以有 。（2）收敛性收敛性：取 时，有 ，由微分中值定理可知，在 和之间存在 使得：反复应用上述关系可得由条件（2），所以即第30页/共115页2 迭代法及其迭代收

15、敛的加速方法（3）误差估计误差估计：由微分中值定理，在 和 之间存在 使得同样得到同时即因为所以反复应用 ，最后得到第31页/共115页2 迭代法及其迭代收敛的加速方法迭代次数估计迭代次数估计：给定精度要求 ，即 ，只要从而迭代次数满足 第32页/共115页2 迭代法及其迭代收敛的加速方法注意注意：由式 可以看出，当越小时，序列 收敛越快。只要相邻两次迭代的偏差足够小，就可以保证近似解 有足够的精度。所以，常采用条件 来控制迭代次数。第33页/共115页2 迭代法及其迭代收敛的加速方法例：求方程 的最大根，要求精度 。解：（1）求隔根区间方程变为 ，做 和的图形，由图可知方程的最大根在区间 内

16、。第34页/共115页2 迭代法及其迭代收敛的加速方法2）建立迭代公式，判断收敛性。方程变形为 ，迭代公式为 在区间 内可导，并且单调递增。又 ，所以，当 时，由于 ，所以 迭代法收敛。第35页/共115页2 迭代法及其迭代收敛的加速方法（3）迭代计算过程如下表所示。所以方程的最大根03.513.7720340.27203423.7882880.01625433.7892210.00093343.7892750.000054第36页/共115页x0=3.5;k=0;while 1 k=k+1;xk=(log10(x0)+7)/2;if abs(xk-x0)=0.0001 break;end x

17、0=xk;end2 迭代法及其迭代收敛的加速方法第37页/共115页2 迭代法及其迭代收敛的加速方法定理（局部收敛性定理）定理（局部收敛性定理）：设 是方程 的根，在 的某一邻域连续，且 ，则必存在 的一个邻域 ，对任意选取的初值 ，迭代公式 产生的数列 收敛于方程的根 ，称迭代法在 的邻域 具有局部收敛性局部收敛性。第38页/共115页注意注意：在实际应用中 事先不知道，因此条件 无法验证。但是如果已知根的初始值 在根 的附近，又根据 的连续性，可采用条件来代替 。2 迭代法及其迭代收敛的加速方法第39页/共115页2 迭代法及其迭代收敛的加速方法例例：用迭代法求方程 在隔根区间 内的根，要

18、求精确到小数点后第4位。解解：(1)构造迭代公式：，迭代公式为(2)判断收敛性：由局部收敛性定理 ，所以迭代法收敛第40页/共115页2 迭代法及其迭代收敛的加速方法(3)迭代计算过程如下表所示：01.511.48124800.018752021.47270570.008542331.46881730.003888441.46704800.001769351.46624300.000805061.46587680.000366271.46571020.000166681.46563450.000075791.46560000.0000345第41页/共115页2 迭代法及其迭代收敛的加速方法x

19、0=1.5;k=0;while 1 k=k+1;xk=(x0*x0+1)(1/3);if abs(xk-x0)=0.00005 break;end x0=xk;end第42页/共115页2 迭代法及其迭代收敛的加速方法用迭代法求方程根的近似值的计算步骤计算步骤如下：(1)准备准备：选定初值 ，确定方程 的等价方程 (2)迭代迭代：按迭代公式 计算出(3)判别判别：如果 ，则终止迭代计算，取 作为方程的近似根。否则，转到第(2)步继续迭代计算。第43页/共115页2 迭代法及其迭代收敛的加速方法迭代收敛的加速方法迭代收敛的加速方法迭代迭代-加速公式加速公式:记 ，则由微分中值定理有其中 在 和

20、之间。设 在根 附近变化不大，又设 ，由迭代收敛条件 ，有上式整理为 第44页/共115页2 迭代法及其迭代收敛的加速方法上式说明 作为根的近似值时绝对误差大致为如果把该误差作为一种补偿，可以得到更好的近似值记得到迭代-加速公式第45页/共115页2 迭代法及其迭代收敛的加速方法例例:用迭代-加速公式求 在区间 内的根。解解:，加速-迭代公式为01.511.48124801.46590550.034094521.46572331.46557420.000331331.46557261.46557130.0000029第46页/共115页2 迭代法及其迭代收敛的加速方法x0=1.5;k=0;wh

21、ile 1 k=k+1;xk1=(x0*x0+1)(1/3);xk=20*xk1/11-9*x0/11;if abs(xk-x0)=0.00005 break;end x0=xk;end第47页/共115页3 牛顿（Newton）迭代法应用范围应用范围：用于求解高次代数方程、超越方程等非线性方程；可求方程的实根、复根，也可求单根和重根。基本思想基本思想：将非线性方程 逐步转化为线性方程来进行求解。牛顿迭代法的最大优点是在方程的单根附近具有较高的收敛速度，是一种将近似根精确化的相当有效的迭代法。第48页/共115页3 牛顿（Newton）迭代法对于非线性方程 ，设 连续可微，将 在 处泰勒（泰勒

22、（TaylorTaylor）展开）展开如果 ，则可在点 附近取线性部分近似替代 ，得到 的近似方程进行变换，得第49页/共115页3 牛顿（Newton）迭代法由迭代法思想，上式写为迭代公式的一般形式称上式为牛顿（牛顿（NewtonNewton）迭代公式）迭代公式，其迭代函数为第50页/共115页3 牛顿（Newton）迭代法几何意义几何意义：第51页/共115页3 牛顿（Newton）迭代法如图所示，假设 是非线性方程 在隔根区间 内的根，在 内可导，且对于 有 。任取初值 ，过曲线 上的点 作切线以其作为曲线 的近似表达式。切线交 轴于 ，则以 作为方程根的第一次近似值。同理，过曲线 上点

23、 做切线第52页/共115页3 牛顿（Newton）迭代法切线交 轴于 ，则 然后以 作为方程根的第二次近似值。反复作切线，第k+1条切线方程为它与 轴的交点横坐标为由此得到方程的近似根数列 。几何意义即依次用切线代替曲线，用线性函数 的零点作为函数零点的近似值，牛顿迭代法又称为切线法切线法。第53页/共115页3 牛顿（Newton）迭代法初始值初始值 的选取的选取：定理（局部收敛性定理局部收敛性定理）：设 是方程 的根，若（1）函数 在 的邻域内具有连续的二阶导数；（2）在 的邻域内 ，则存在 的某个邻域 ，对于任意的初值，由牛顿迭代法产生的数列收敛于根 。第54页/共115页3 牛顿（N

24、ewton）迭代法证明证明：迭代函数 的导数为由条件(1)和(2)可知，在 的邻域内可导。又 ，所以在 的某个邻域 内恒有 第55页/共115页3 牛顿（Newton）迭代法定理（非局部收敛定理非局部收敛定理）：设 是方程 在隔根区间内的根，如果（1）对于 连续且不变号；（2）选取初始值 ，使 ，则由牛顿迭代公式产生的数列收敛于根 。如果在 上 的符号不容易判断，则利用下式确定初值即如果在 处，满足上式，并且 ，则大多数情况下可以保证牛顿迭代法是收敛的。第56页/共115页3 牛顿（Newton）迭代法计算步骤计算步骤：1.选定初值 ，构造迭代公式 ；2.迭代计算，得到 ；3.迭代终止判断，由

25、迭代终止控制条件 控制迭代计算是否终止；4.迭代计算终止，得到方程根的近似值。第57页/共115页3 牛顿（Newton）迭代法例例：用牛顿迭代法求方程 在隔根区间 内的根，要求精确到小数点后第4位小数。解解:(1)令 ，牛顿迭代公式为(2)判断收敛性所以取 时满足 ，牛顿迭代法收敛。第58页/共115页3 牛顿（Newton）迭代法(3)计算过程如下表01.511.4666670.03333321.4655720.00109531.4655710.000001第59页/共115页x0=1.5;k=0;while 1 k=k+1;xk=(2*x03-x02+1)/(3*x02-2*x0);if

26、 abs(xk-x0)=0.00005 break;end x0=xk;end3 牛顿（Newton）迭代法第60页/共115页3 牛顿（Newton）迭代法例例：用牛顿迭代法求方程 在 附近的根，精确到小数点后第4位。解解：所以由此 成立，所以可取为初始值。迭代公式为：第61页/共115页3 牛顿（Newton）迭代法迭代计算过程如下表所示0-11-0.9772730.0227272-0.9604610.0168123-0.9533910.0070704-0.9525000.0008915-0.9524840.000016第62页/共115页x0=-1;k=0;while 1 k=k+1;x

27、k=(40*x041+2*x03-1)/(41*x040+3*x02);if abs(xk-x0)=0.00005 break;end x0=xk;end3 牛顿（Newton）迭代法第63页/共115页3 牛顿（Newton）迭代法例例：用牛顿迭代法建立计算 近似值的迭代公式。解解：令 ，则以上问题化为求方程的正根，牛顿迭代公式为继续证明对于任意初值 ，上面迭代公式恒收敛。证明：对于任意 有 即序列 有下界；第64页/共115页3 牛顿（Newton）迭代法又即 单调下降，所以迭代公式收敛。计算 的近似值。的值在1011之间，取初值 ，迭代公式为第65页/共115页3 牛顿（Newton）迭

28、代法x0=10;C=115;k=0;while 1 k=k+1;xk=(x0+C/x0)/2 err=abs(xk-x0)if err=0.000005 break;end x0=xk;end第66页/共115页3 牛顿（Newton）迭代法迭代计算过程如下（）：010110.750.75210.723837210.02616279310.723805290.00003192410.72380529第67页/共115页迭代过程的收敛速度定义定义：如果从任何可取得的初始值出发，由迭代法产生的迭代序列都收敛于方程的根，则称该迭代法是大范围收敛大范围收敛；如果选取的初始值充分接近于方程的根，或者说存

29、在根的某一邻域，使得在其中任意选取初始值，迭代序列都收敛于方程的根，则称该迭代法为局部收敛局部收敛。一种迭代法具有使用价值，不但需要是收敛收敛的，还要求收敛的比较快。所谓迭代过程的收敛速度收敛速度，是指在接近收敛时迭代误差的下降速度。第68页/共115页迭代过程的收敛速度定义定义：设迭代法 产生的数列 收敛于 的根 ，令误差 ，如果存在某个实数 以及正常数 ，使则称数列 是 阶收敛阶收敛的，相应的迭代法是 阶方阶方法法，常数 称为渐进误差常数渐进误差常数。第69页/共115页迭代过程的收敛速度 可以看出，和 为同阶无穷小量，即以 为基本无穷小量时，为 阶无穷小量，阶数 越高，收敛速度越快。收敛

30、速度是误差的收缩率，阶数越高，误差下降的越快。当 时，称数列 线性收敛线性收敛；当 时，称数列 为平方收敛平方收敛（或二阶收敛）；当 时，称数列为超线性收敛超线性收敛。显然，越大，数列收敛的越快，所以迭代法的收敛阶收敛阶是对迭代法收敛速度的一种度量。第70页/共115页迭代过程的收敛速度证明简单迭代法是线性收敛：设 是方程 的根，在 的某一邻域连续，在邻域内有 ，且 ，则简单迭代法线性收敛。证明证明：设迭代过程 收敛于 ，则式中 在 和 之间。令 ，则于是迭代过程一阶收敛（线性收敛）。第71页/共115页迭代过程的收敛速度证明牛顿迭代法平方收敛：设 是 在隔根区间 内的根，如果（1）对于 连续

31、且不变号；（2）选取初始值 ，使 ，则牛顿迭代法平方收敛。证明证明：将 在 处做泰勒展开，有式中 介于 之间。整理上式得第72页/共115页迭代过程的收敛速度即两边取极限，并利用得所以牛顿迭代法平方收敛。第73页/共115页4 弦截法如果函数 比较复杂，求导可能比较困难，有时甚至导函数不存在。尤其是当 很小时，计算中容易产生很大的误差。这时可将牛顿迭代公式中的 用差商差商来代，即就得到求解非线性方程的弦截法弦截法。第74页/共115页4 弦截法第75页/共115页4 弦截法设方程 的一个隔根区间为 ，如图所示。连结曲线 上的两点 得到弦AB，令 ，则弦AB的方程为则弦AB与 轴交点的横坐标 为

32、以 作为方程根 的一个近似值。第76页/共115页4 弦截法然后过曲线 上的两点 做弦，得到其与 轴交点的横坐标 为于是把 作为方程根 的一个新的近似值。反复做弦，得到迭代公式的一般形式为利用上述公式求方程 根的近似值的方法就叫做弦截法弦截法。第77页/共115页4 弦截法几何意义几何意义：依次用弦来代替曲线，用线性函数的零点作为函数 零点的近似值。定理（弦截法收敛定理）定理（弦截法收敛定理）：如果 在根 的某个邻域 内有二阶连续导数，且对任意 ，有 ，则当 邻域充分小时，对 邻域内任意的 ，由弦截法迭代公式得到的近似值序列 收敛到方程 的根 。并可证明弦截法是按阶 收敛的。第78页/共115

33、页4 弦截法计算步骤计算步骤：（1）准备：选定初始近似值 ，并计算 ；（2）迭代计算：依迭代公式 计算得到 ，并计算 ；（3）控制：如果 或 ，则终止迭代。否则 以 分别代替 ，转（2）继续迭代计算。弦截法和牛顿法比较：都是先把 线性化，但方式不同，牛顿法用切线方程来近似代替非线性方程，而弦截法是用弦线来近似代替非线性方程。第79页/共115页4 弦截法例例：用弦截法求方程 在区间 内根的近似值。精确到小数点后第4位。解解：令 ，取 ，弦截法迭代公式为第80页/共115页4 弦截法迭代计算过程如下表所示：00.51120.5755900.42441030.5995400.02395040.61

34、06930.01115350.6103580.00033560.6103620.000003第81页/共115页4 弦截法x0=0.5;x1=1;k=0;while 1 k=k+1;xk=x1-(exp(2*x1)+x1-4)*(x1-x0)/(exp(2*x1)+x1-exp(2*x0)-x0);if abs(xk-x0)=0.00005 break;end x0=x1;x1=xk;end第82页/共115页课后题1、二分法求方程 在区间 内的根，要求误差不超过 。2、已知方程 在区间 内有一根，试问用二分法求根，使其具有5位有效数字至少应二分多少次？3、判断下列方程由几个根，求出隔根区间，

35、并写出收敛的迭代公式。（1）（2）5、用迭代法求的 正根，要求准确到小数点后第5位。第83页/共115页课后题4、方程 在1.5附近有根，把方程写成4种不同的等价形式，并建立相应的迭代公式：判断每种迭代公式产生的数列在1.5附近的收敛性，并用（1）的迭代公式求方程的根，要求误差不超过 。第84页/共115页课后题6、用迭代法求方程 的根，要求准确到小数点后第4位。7、用迭代-加速公式求方程 在 附近的根，要求准确到小数点后第4位。9、应用牛顿法于方程 和 ，分别导出求 的迭代公式，并求 10、用牛顿法求方程 在 附近的根，要求准确到小数点后第3位。第85页/共115页课后题11、证明计算 的牛

36、顿迭代公式为取 ，计算 ，要求 。12、用弦截法求 的根，取 ，计算到 为止。第86页/共115页1、解解：令 ，则所以 单调增；所以 在区间 内单调减。又 ，所以方程有唯一实根。有课后题隔根区间10.35-0.0350.3,0.3520.3250.0360.325,0.3530.33750.0004750.3375,0.3540.34375-0.0170.3375,0.34375第87页/共115页2、解解：根的区间为 ，近似值具有5位有效数字，则精确到小数点后4为小数，即取 。4、解解：（1），所以迭代公式收敛。课后题第88页/共115页（2），收敛（3）（4）利用公式（1）计算方程的根，

37、列表如下 课后题01.511.4812480.01875261.4658770.000366第89页/共115页5、解解：方程变形为 ，分别画出 和的图形，可知交点在 的范围；取 ，则，可知交点在 的范围，所以隔根区间可取为 。令 ，从而 单调减；单调增；，满足条件1；，满足条件2，所以迭代公式收敛。课后题第90页/共115页迭代公式为：计算过程如下：课后题0111.0371370.03713721.0434790.00634031.0445460.00106841.0447260.00017951.04475661.0447671.04476第91页/共115页6、解解：（1）求隔根区间方程

38、变形为 ，做 和 的图形。由图可知，方程的根在区间 内。（2）建立迭代公式，判断其收敛性将方程等价变形为 ，迭代函数 ，迭代公式为 。课后题第92页/共115页因为 ，当 时，所以 单调递减；又 ，所以 ；，所以 ；由定理可知，迭代公式收敛。（3）计算：取 ，列表计算课后题0010.10.120.0894830.01051750.0905260.000013第93页/共115页7、解解：取 ，迭代-加速公式为：计算过程如下：课后题00.510.566580.0665820.567130.000550230.56714第94页/共115页9、证明证明：（1）对于方程 ，牛顿迭代公式为所以 课后题

39、第95页/共115页（2）对于方程 ，牛顿迭代公式为所以课后题第96页/共115页10、解解：令 ，牛顿迭代公式为 牛顿迭代公式为满足局部收敛性定理，计算过程如下表所示课后题0211.8888890.11111121.8794520.00943731.8793850.000067第97页/共115页11、证明证明：由 ，牛顿迭代公式为取 ，计算过程如下：课后题0111.6666670.666666721.4711110.195555631.4428120.028299041.44225051.442250第98页/共115页课后题12、解解：令 ，满足收敛条件，取 ，弦截法迭代公式为计算过程如

40、下表所示 001120.5430440.05978930.5080930.00539540.5109860.000023第99页/共115页第三章 非线性方程的数值解法主要内容主要内容：1.确定隔根区间；2.二分法；3.简单迭代法；4.迭代收敛的加速方法；5.牛顿迭代法；6.弦截法；7.迭代法的收敛阶。第100页/共115页第三章 非线性方程的数值解法二分法二分法：第101页/共115页迭代法收敛定理迭代法收敛定理：设方程 ，若迭代函数 在有根区间 上满足：（1）当 时，；（2）在 上可导，且有 ，；则有（1）方程 在 上有唯一的根 ；（2）对任意初值 ，迭代公式 产生的数列 收敛于方程的唯一

41、根 ，即 ；（3）误差估计第三章 非线性方程的数值解法第102页/共115页第三章 非线性方程的数值解法定理（局部收敛性定理）定理（局部收敛性定理）：设 是方程 的根，在 的某一邻域连续，且 ，则必存在 的一个邻域 ，对任意选取的初值 ，迭代公式 产生的数列 收敛于方程的根 ，称迭代法在 的邻域 具有局部收敛性局部收敛性。第103页/共115页第三章 非线性方程的数值解法迭代迭代-加速公式加速公式：牛顿迭代法牛顿迭代法第104页/共115页第三章 非线性方程的数值解法初始值初始值 的选取的选取：定理（局部收敛性定理局部收敛性定理）：设 是方程 的根，若（1）函数 在 的邻域内具有连续的二阶导数

42、；（2）在 的邻域内 ，则存在 的某个邻域 ，对于任意的初值，由牛顿迭代法产生的数列收敛于根 。第105页/共115页第三章 非线性方程的数值解法定理（非局部收敛定理非局部收敛定理）：设 是方程 在隔根区间内的根，如果（1）对于 连续且不变号；（2）选取初始值 ，使 ，则由牛顿迭代公式产生的数列收敛于根 。如果在 上 的符号不容易判断，则利用下式确定初值即如果在 处，满足上式，并且 ，则大多数情况下可以保证牛顿迭代法是收敛的。第106页/共115页第三章 非线性方程的数值解法二分法二分法：由公式有取 。即采用二分法求解，得到满足精度要求的根需要二分13次。第107页/共115页第三章 非线性方

43、程的数值解法简单迭代法简单迭代法：方程 变形为而 即迭代函数单调减，又 ；所以 单调减，所以所以在区间 上迭代法收敛，取初始值 ；第108页/共115页第三章 非线性方程的数值解法迭代计算过程如下表所示：00.510.6263810.12638120.6079930.01838830.6107110.00271840.6103100.00040150.6103700.00006060.6103610.000009第109页/共115页第三章 非线性方程的数值解法迭代迭代-加速公式加速公式：取 ，公式为取初始值 ，第110页/共115页第三章 非线性方程的数值解法迭代计算过程如下表所示：00.5

44、10.6053180.10531820.6101410.00482330.6103520.00021140.6103610.000009第111页/共115页第三章 非线性方程的数值解法牛顿迭代法牛顿迭代法：当 时，不变号；又 ，所以取 ；迭代公式为第112页/共115页第三章 非线性方程的数值解法迭代计算过程如下表所示：0110.7218260.27817420.6206930.10113430.6104540.01023840.6103620.00009350.6103620.000000007第113页/共115页总结：总结：1、对非线性方程 求根，首选要对 的性态及根的近似位置有一个大致了解，使用较小的有根区间把方程的根分离出来；2、简单迭代法是一种逐次逼近的方法，需要判断收敛性，了解收敛速度，可以采用加速公式；3、牛顿迭代法是最常用的一种迭代方法，具有二阶收敛速度，但是对初值的选取比较苛刻，需要靠近方程的根，否则可能不收敛；4、弦截法是牛顿迭代法的一种修改，虽然收敛速度较牛顿法慢，但不需要计算导数，使用更加方便。同样要求初始值的选取比较靠近方程的根。第三章 非线性方程的数值解法第114页/共115页感谢您的观看。第115页/共115页