数值分析3-2课件

1牛顿法牛顿法牛顿法的局部收敛性牛顿法的局部收敛性牛顿法的大范围收敛性牛顿法的大范围收敛性简化牛顿法和牛顿下山法简化牛顿法和牛顿下山法不用导数求根不用导数求根非线性方程组非线性方程组电子工程应用电子工程应用第三章第三章 非线性方程求根非线性方程求根1牛顿法第三章 非线性方程求根2牛顿法的几何意义：牛顿法的几何意义：x0y=f(x)x1x2x*解方程解方程 f(x)=02牛顿法的几何意义：x0y=f(x)x1x2x*解方程 33.3.2 牛顿法的局部收敛性牛顿法的局部收敛性(1)当当 x*是方程的单根时，可得是方程的单根时，可得二阶局部收敛二阶局部收敛33.3.2 牛顿法的局部收敛性(1)当 x*是方程的单根4(2)当当x*是方程的是方程的m 重根时，可得重根时，可得3.3.2 牛顿法的局部收敛性牛顿法的局部收敛性4(2)当x*是方程的m 重根时，可5一阶局部收敛一阶局部收敛3.3.2 牛顿法的局部收敛性牛顿法的局部收敛性5一阶局部收敛3.3.2 牛顿法的局部收敛性63.3.2 牛顿法的局部收敛性牛顿法的局部收敛性牛顿法求牛顿法求重根重根只是只是线性收敛线性收敛则则 若用此迭代格式求解，具有若用此迭代格式求解，具有 2阶阶 收敛速度收敛速度（但需要知道根的重数但需要知道根的重数m）若取迭代函数为若取迭代函数为63.3.2 牛顿法的局部收敛性牛顿法求重根只是线性收敛则 7牛顿法牛顿法牛顿法的局部收敛性牛顿法的局部收敛性牛顿法的大范围收敛性牛顿法的大范围收敛性简化牛顿法和牛顿下山法简化牛顿法和牛顿下山法不用导数求根不用导数求根非线性方程组非线性方程组电子工程应用电子工程应用第三章第三章 非线性方程求根非线性方程求根7牛顿法第三章 非线性方程求根83.3.3 牛顿法的大范围收敛性牛顿法的大范围收敛性定理定理 3.7 设设 f(x)在在区间区间a,b 上存在二阶连续导数，上存在二阶连续导数，且满足条件：且满足条件：1)f(a)f(b)0;2)x a,b 时，时，f(x)0；3)当当 x a,b 时，时，保号；保号；则区间则区间a,b 上上存在唯一根。存在唯一根。对任意初值对任意初值 x0a,b，由牛顿迭代公式产生的序由牛顿迭代公式产生的序列列 xk 二阶收敛二阶收敛.4)83.3.3 牛顿法的大范围收敛性定理 3.7 设 f(x9条件条件1）保证方程在）保证方程在a,b 区间上存在根。区间上存在根。条件条件2）保证方程在）保证方程在a,b 区间上存在唯一根。区间上存在唯一根。条件条件3）则保证曲线的凹凸性不变，即不存在拐点。）则保证曲线的凹凸性不变，即不存在拐点。3.3.3 牛顿法的大范围收敛性牛顿法的大范围收敛性增凸增凸减凸减凸增凹增凹减凹减凹条件条件4）保证了）保证了 a,b9条件1）保证方程在a,b 区间上存在根。3.3.3 牛101011证明：设第一种情形证明：设第一种情形3.3.3 牛顿法的大范围收敛性牛顿法的大范围收敛性由于由于f(x)在区间在区间a,b上单调增加，且由条件上单调增加，且由条件(4)知，知，f(x0)0)的迭代公式，并的迭代公式，并分析算法的收敛性，计算分析算法的收敛性，计算115。解：解：x2 c=0令令 f(x)=x2 c,则则牛顿迭代法求解牛顿迭代法求解 x2 c=0 的计算格式的计算格式分析该迭代格式的收敛性。分析该迭代格式的收敛性。3.3.3 牛顿法的大范围收敛性牛顿法的大范围收敛性13例 给出用牛顿法求平方根c(c 0)的迭代公式，14条件条件1）f(0+)0,满足；满足；条件条件2）x(0,+)时，时，f(x)0,满足满足；条件条件3）x(0,+)时，时，f(x)=2 0,满足满足；条件条件4）计算计算115，若取初值，若取初值 x0=12，迭代，迭代 4次，可得结次，可得结果。果。若取若取x0(c,+)，则满足，则满足3.3.3 牛顿法的大范围收敛性牛顿法的大范围收敛性实际上，只要实际上，只要x00，迭代法都收敛，迭代法都收敛。14条件1）f(0+)0,15fi=inline(x+115/x)/2);x0=12;er=1;k=0;while er0.5e-7;x=fi(x0)er=abs(x-x0);x0=x;k=k+1;endk=4,x=10.7238053.3.3 牛顿法的大范围收敛性牛顿法的大范围收敛性15fi=inline(x+115/x)/2);k=416牛顿迭代公式牛顿迭代公式牛顿迭代格式的局部收敛性牛顿迭代格式的局部收敛性牛顿迭代格式的大范围收敛性牛顿迭代格式的大范围收敛性简化牛顿法和牛顿下山法简化牛顿法和牛顿下山法不用导数求根不用导数求根非线性方程组非线性方程组电子工程应用电子工程应用第三章第三章 非线性方程求根非线性方程求根16牛顿迭代公式第三章 非线性方程求根173.3.4 简化简化牛顿法牛顿法和牛顿下山法和牛顿下山法牛顿法每迭代一次，要计算导数值，影响计算效率牛顿法每迭代一次，要计算导数值，影响计算效率为避免计算重复计算导数值，可取其为一为避免计算重复计算导数值，可取其为一定值定值，如，如称为称为简化牛顿法简化牛顿法简化牛顿法简化牛顿法173.3.4 简化牛顿法和牛顿下山法牛顿法每迭代一次，要计183.3.4 简化简化牛顿法牛顿法和牛顿下山法和牛顿下山法x0y=f(x)x1x2x*迭代函数：迭代函数：推广的简化牛顿法：推广的简化牛顿法：把把f(xk)取为任意的常数取为任意的常数 1/C(C0)，迭代函数，迭代函数183.3.4 简化牛顿法和牛顿下山法x0y=f(x)x193.3.4 简化简化牛顿法牛顿法和牛顿下山法和牛顿下山法牛顿法的收敛性依赖于初值的选取。如果初值偏离所牛顿法的收敛性依赖于初值的选取。如果初值偏离所求的根求的根 x*较远，则牛顿法可能发散。较远，则牛顿法可能发散。为了防止迭代发散，对迭代过程附加一项要求，即：为了防止迭代发散，对迭代过程附加一项要求，即：满足此要求的算法称为满足此要求的算法称为下山法下山法将将下山法下山法与与牛顿法牛顿法结合起来使用（下山法保证函数值结合起来使用（下山法保证函数值稳定下降，牛顿法加快收敛速度）稳定下降，牛顿法加快收敛速度）牛顿下山法牛顿下山法 牛顿下山法牛顿下山法193.3.4 简化牛顿法和牛顿下山法牛顿法的收敛性依赖于初202）与前一步的近似值）与前一步的近似值 xk 适当加权平均适当加权平均其中，其中，（0 1）称为）称为下山因子下山因子。3）则有）则有试选常取试选常取 ，直到使条件，直到使条件满足，此时得迭代值满足，此时得迭代值xk+1。然后令。然后令k=k+1，将，将恢复恢复为为1，返回步骤，返回步骤3），进行下一步循环。），进行下一步循环。1）牛顿法计算结果作为中间结果：）牛顿法计算结果作为中间结果：牛顿下山法牛顿下山法的计算步骤：的计算步骤：202）与前一步的近似值 xk 适当加权平均其中，（00.00005 x=x1-f(x1)*(x1-x0)/(f(x1)-f(x0)er=abs(x-x1);x0=x1;x1=x;k=k+1;endk=4,x=0.46557割线法（弦截法）割线法（弦截法）k012345xk0.4 0.60.45745 0.46460 0.465580.4655727f=inline(x*(x+1)2-1);k=试位法位法28在给定含有根的区间在给定含有根的区间 a,b（假设假设 f(a)f(b)0），定），定义下一个点为义下一个点为因因为点点 和和 在在x轴的两的两侧，所以，所以 c点点一定在一定在a,b中中根据根据 或者或者 选择选择 a,c 或者或者 c,b作为新的有根区间作为新的有根区间试位法28在给定含有根的区间 a,b（假设 f(a)试位法位法 给定区间给定区间 a,b，使，使 f(a)*f(b)0 for i=1,2,3,.if f(c)=0 stop end if f(a)*f(c)eps&iid 时场指数衰减指数衰减 623.6 电子工程应用（TM表面波截止波数）定义两个区域的633.6 电子工程应用（电子工程应用（TM表面波截止波数表面波截止波数）通解通解上式的解，无上式的解，无论 kd和和ka是是实数数还是虚数都是是虚数都是对的，但的，但由于前由于前页截止波数的截止波数的选择确定，因此得出两者都是确定，因此得出两者都是实数。数。633.6 电子工程应用（TM表面波截止波数）通解上式的解，643.6 电子工程应用（电子工程应用（TM表面波截止波数表面波截止波数）边界条件：边界条件：连续连续 x=d643.6 电子工程应用（TM表面波截止波数）边界条件：连续653.6 电子工程应用（电子工程应用（TM表面波截止波数表面波截止波数）对于非零解，上页最后两方程的行列式必须为零，得对于非零解，上页最后两方程的行列式必须为零，得介介质空气空气构成了一构成了一组联立超越方程，立超越方程，对于于给定的定的k0和和r可利用非可利用非线性方程求根的方法解出性方程求根的方法解出传播常数播常数kd和和ka。653.6 电子工程应用（TM表面波截止波数）对于非零解，上66小小结非线性方程求根的牛顿法、局部收敛性、大范非线性方程求根的牛顿法、局部收敛性、大范围收敛性围收敛性割线法、简化牛顿法、牛顿下山法割线法、简化牛顿法、牛顿下山法牛顿法的收敛阶牛顿法的收敛阶解非线性方程组的牛顿法解非线性方程组的牛顿法电子工程应用电子工程应用练习：习题三：练习：习题三：8（1）、）、12（自愿做，自愿交，不交练习者不扣作业分）（自愿做，自愿交，不交练习者不扣作业分）66小结非线性方程求根的牛顿法、局部收敛性、大范围收敛性第三章“计算TM表面波截止波数”zd电介质0r0接地板参考文献：D.M.Pozar,微波工程（第三版），3.6节，电子工业出版社，2006。课程设计2第三章“计算TM表面波截止波数”zd电介质0r0接参考电子科技大学学报投稿格式 1.题目和摘要（中英文）2.正文：a 引言;b 原理；c 数值算例；d 结论。3.附件：程序 第八周周三上课时交打印版和电子版（发邮件：。每章约一个设计题目，每人总共做2个设计题目。格式要求格式要求参考电子科技大学学报投稿格式格式要求