数值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 多变量牛顿法 * 牛顿法是忽略高阶无