Davidenko方法在搜寻分层介质表面波极点中的应用研究.pdf

下犷 Davidenko方法在搜寻分层介质 表面波极点中的应用 周桂蓉 冯宁宁方大纲 南京理工大学电光学院毫米波技术研究室，南京，210094 【摘要】:本文先简要介绍了Davidenko这种复根搜寻方法，该方法在函数具有 解析表达式时能克服 Mulle:方法等对初值依赖较人的缺点。在此基础上我们将 Davidenko方法推广到用于搜寻分层各向异性介质表面波的多个极点的情形，从 而解决了该方法只能用于搜寻单根的不足。 引言 在很多电磁场数值问题中会碰到复根搜寻的问题，诸如:求解有耗传输线本 征问题的传播常数、抽取复极点或奇异点 (例如表面波极点)从而在计算中加速 收敛等。因此，也可以据此将复根搜寻问题分为相应的两大类:所涉及的复函数 无显式表达形式和具有显式解析表示式。传统的二分法、牛顿迭代法以及Mullet 方法等都可以用来计算这两类问题[1][2]=但它们的问题是:二分法计算时间长， 尤其对于后一类问题，函数奇异性较大，如果计算间隔不够密就有可能漏根。牛 顿迭代法以及 Muller方法则对初值依赖性较大。对于后一类问题除了这些方法 以外，还可以考虑用Davidenko方法，它的最大好处就是对初值的依赖性降低。 Davidenko方法最早由文献[3]引入化学问题计算领域，后来被文献4【]用于电磁 场数值计算领域损耗波导的色散特性的求解，文献[[5]又将这一方法用于求两层 各向同性介质的表面波极点以及用于反射系数测x中，但它的所有计算结果都只 局限于单根情形 。本文将给出Davidenko方法用于多根情况的处理办法，计算出 了两层各向异性介质在不同介质厚度 卜的的多个表面波极点的位置。 2. 有关计算公式 (1)Davidenko方法的原理[4][5]: 对一个复函数方程: F(z=z+分)=Q (1) 我们可以通过求解下面的常微分方程组来解得它的复根。 adzt一尚R(eF[]ReF[,]+ImF[]，一F!,b (2) _dd,t,一尚R(eF[]ImF[,’一mF[‘jReF[,b (3) 其中，F.代表aF/az,:为辅助变量，先设 个初始值:。二((Y+JY)!从。开始 以一定的步长△t递增，月」龙格一库塔方法解此微分方程组 。当 t增加到满足 }、*:一z卜EPS，此时%十、即为所求根 其中n、代表迭代步数，EPS代表一 个很小的常数。 一 63一 (2)多根情形的处理: 假设此复函数有n个复根 ，则可以表示成为以 卜n个根乘积形式: F(z)=kn(:一Zpi) (3) 设已搜得复根个数为n，,，则每次重新搜 卜面一个根的时候，必须在原函数中除 去已有的根，同时导数也要相应变化，此时函数变为: 产 4 、 F(z)=F(z)IK(i) ‘ ， 、 沙 其中， 气 了 扮 3 . 、 K(0)=1， K(i卜K(i一1)(:一Z,)，i=1,2，二n, 、 r ，