数值分析3-2资料.pptVIP

* 3.6 电子工程应用（TM表面波截止波数） 通解 上式的解，无论 kd和ka是实数还是虚数都是对的，但由于前页截止波数的选择确定，因此得出两者都是实数。 * 3.6 电子工程应用（TM表面波截止波数） 边界条件： 连续 x d * 3.6 电子工程应用（TM表面波截止波数） 对于非零解，上页最后两方程的行列式必须为零，得 介质 空气 构成了一组联立超越方程，对于给定的k0和εr 可利用非线性方程求根的方法解出传播常数kd和ka。 * 小结 非线性方程求根的牛顿法、局部收敛性、大范围收敛性 割线法、简化牛顿法、牛顿下山法 牛顿法的收敛阶 解非线性方程组的牛顿法 电子工程应用 练习：习题三：8（1）、12 （自愿做，自愿交，不交练习者不扣作业分） 第三章 “计算TM表面波截止波数” z d 电介质 ε0 εrε0 接地板 参考文献：D. M. Pozar, 微波工程（第三版），3.6节，电子工业出版社，2006。 课程设计2 参考《电子科技大学学报》投稿格式 1. 题目和摘要（中英文） 2. 正文： a 引言; b 原理； c 数值算例； d 结论。 3. 附件：程序 第八周周三上课时交打印版和电子版（发邮件：xsyang@126.com。 每章约一个设计题目，每人总共做2个设计题目。 格式要求 * * Muller法 用前面估计的3个点 xi, xi+1, xi+2，画出通过它们的抛物线 y p x （抛物线的轴平行于y轴） * 抛物线一般与x轴没有交点，或者有两个交点 如果有两个交点，则距离xi+2最近的一个点被选作xi+3 如果抛物线与x轴不相交，则方程有复数解，就需要确定复根 反二次插值 Inverse Quadratic Interpolation, IQI 用前面估计的3个点xi, xi+1, xi+2，画出通过它们的抛物线，抛物线具有x p y 的形式（轴线平行于 x 轴） * 将y 0代入，整理得 这条抛物线将与x轴相交于一个点，作为xi+3 通过3个点 a, A , b, B , c, C 的二次多项式 x P y 是 其中 反二次插值 Inverse Quadratic Interpolation, IQI 给定3个初始估计，通过IQI方法，就可以计算出新的估计值 x P 0 ，即：近似根 用新的估计值去代替最早的一个估计值 * IQI的另一种实现方式，是用新的估计值代替前面3个估计值中误差最大的那个。 * 3.4.2 Brent法 Brent方法是混合方法： 把保证收敛的二分法与具有快速收敛特性的求根方法结合起来 建立了一种保留各自最有用性质的新方法 用于连续函数f，而且区间以a和b为界 。 * 算法： 首先尝试用反二次插值法 如果满足下面两条之一 （1）后向误差改进了 （2）包含根的区间至少缩短一半 则用计算结果代替xi, ai, bi中的一个 否则，尝试用割线法 如果还是失败，则用二分法进行，来保证其区间长度至少减半 3.4.2 Brent法 * MATLAB的fzero 命令实现了一种Brent方法，它包含了一个预处理步。 如果用户没有提供含有根的初始区间，则fzero就会自己找一个比较好的初始区间。 3.4.2 Brent法 如果用户提供了初始区间，则不需要预处理步。 如分别输入： fzero ‘x^3+x-1’, 1, optimset ‘Display’, ‘iter’ （在1附近求根，并显示迭代方法） 和 fzero x^3+x-1, [0, 1], optimset Display, iter （在[0,1]区间求根，并显示迭代方法） * * * 牛顿迭代公式 牛顿迭代格式的局部收敛性 牛顿迭代格式的大范围收敛性 简化牛顿法和牛顿下山法 不用导数求根 非线性方程组 电子工程应用 第三章 非线性方程求根 * 3.5 非线性方程组 当多个非线性方程构成方程组时，求解难度增加 本节主要介绍求解非线性方程组的牛顿法及其变形 3.5.1 多变量牛顿法 单变量牛顿法 来源于泰勒展开式提供的线性近似 对于非线性方程组，同样处理 * 定义向量值函数 表示对向量值函数求偏导的Jacobi矩阵为 3.5.1 多变量牛顿法 * 在点x0附近，对向量值函数进行泰勒展开，其展开式为 例如，在x0 0, 0 附近， 的泰勒展开式为 3.5.1 多变量牛顿法 * 牛顿法是忽略高阶无穷小O h2 的一个线性近似 假设x*为根，x0为假设的初始解，则 为避免进行矩阵求逆 s 故：多变量牛顿法的迭代步骤如下 Jacobi矩阵 讲义第三章p18, 3-5-2 符号有误 s是向量，可利用解线性方程组的方法，求解s * 例：取初始估计为 1, 2 ，用牛顿