电磁波对流层传播抛物方程分步Pad(e)方法.pdfVIP

第22卷增刊 电波科学学报 Vbl．22，Sup CHINESE oFRADloSCIENCE 2007年9月 JOURNAL September’2007 电磁波对流层传播抛物方程的分步Pad6方法 刘艳丽1，2李廉林1李芳1 (1．中国科学院电子学研究所；2．中国科学院研究生院，北京lO0080) 摘要介绍抛物方程中对指数伪微分算子作高阶Pad6近似的方法，推导了基于 Pad6级数展开的宽角抛物方程的分步Pa彪解．该方法误差小，且易于数值实现， 保证解的稳定性。推导了采用分步Pa必方法第二辅助函数的wAPE的有限差分 形式，并用该辅助函数对阻抗边界条件建模．用基于5项Pad6级数展开的wAPE 求解对流层的电波传播问题，与几何光学双射线模型以及窄角抛物方程(NAPE) 的计算结果比较，证明了分步Pad6方法在求解对流层的电磁波犬角度传播问题的 有效性． 关键词抛物方程(PE)，分步Pad6方法(ssP)，L七ontovich边界条件 1 引 言 2抛物波动方程的分步Pa蛇求解法 20世纪80年代中期起，抛物方程(PE)方法 Z 广泛地用于研究电波传播问题。当电磁波的主要 能量在近轴方向的锥形区域内传播时，可以忽略 波的后向散射，此时椭圆型波动方程近似为PE。 19世纪40年代Leontovich和Fock在地球表面绕 射场计算中引入抛物近似【l】。Kuttler和Dockery 理论分析了PE的分步Fourier方法【21，并引入混 合Fourier变换方法来处理阻抗边界条件13】。PE ．)( 方法可以充分考虑传播区域的地形起伏、地表电 磁特性变化以及传播空间媒质折射率变化的影 圈1电磁波在对流层中的近轴传播，其主要能量在 响。 近轴方向的锥形区域内传播 目前，PE的数值解法主要有分步傅里叶变换 方法‘2】【4】、隐式有限差分法佟石】和有限元方法阴等。 考虑二维电磁问题，直角坐标系下，场分 常用的一种窄角抛物方程(NAPE)对伪微分算子量y(工，z)满足标量波动方程 作一阶Taylor近似，适用于计算大场景、远距离、 r’(1) 萼+警搿咒2(砒)y=o 孤2。az2…””。，r。 低仰角的对流层电波传播问题，当存在地形起伏 其中， 七是真空中的波数， 令 时，NAPE误差较大。 将海洋声学中的分步Pa酯方法(SsP)隅】应用 “(x，z)=已一妇少(工，z)，表示近轴方向为懈，则“ 满足的标量波动方程为 于对流层的电波传播问题，对PE指数伪微分算 子作Pad6级数展开，用牛顿迭代法求解Pad6系{a，+玻(卜√而)】×{a，+玻(1+√而)】“=o(2) 数，推导了WAJPE在阻抗边界条件下采用第二辅 助函数的SSP数值解法，最后以垂直极化的电磁 其中，z=七以a：+，12(x，z)一1，a，=a／孤， a．=a／az。 波在地面上的传播为例，比较了SSP与几何光学 当电磁波的能量主要在以帆为近轴方向的 双射线模型以及窄角抛物方程(NAPE)的计算结 锥形区域内传播时，忽略后向散