多刃峰环境无线电波传播预测的双向抛物方程法.pdfVIP

第 卷 第 期. , . 电 波 科 学 学 报 年 月 . 文章编号 ? ? ? 多刃峰环境无线电波传播预测的 双向抛物方程法 王 昆 。 杨永钦 龙云亮 郭建炎。.中山大学电子与通信工程系,广东 广州 ; .广东工业大学信息工程学院, 广东 广州 ; .厦门理工学院电子与电气工程系,福建 厦门 摘 要 基于迭代双向抛物方程法求解了单刃峰和多刃峰环境下的前向和后向电波 传播。对单刃峰,将双向与单向抛物方程法的数值结果做了比较,并分析了刃峰的位 置和高度对后向反射场的影响。对多刃峰,采用迭代步进的双向抛物方程法计算相 邻刃峰间的电波传播,并将数值结果与单向抛物方程法的结果进行对比。双向抛物 方程法能计算出无线电波的前向传播场和后向传播场叠加的总场,比单向抛物方程 法的预测结果更准确。 关键词 无线电波传播;双向抛物方程;傅里叶分步步进法;刃峰 中图分类号 文献标志码 法的结果做了比较,但是,该模型的传播距 .引 言 离短,传播环境简单,也没有考虑到大气波导的影 最早由和提出的抛物方程法 响; 用双向抛物方程法计算了对流层多刃峰 是亥姆霍兹波动方程的一种前向近轴传播近似 , 环境的前向和后向传播电磁场引,考虑了大气波导 的影响,但是只讨论了窄角双向抛物方程,而且其数 其最常用的数值解法之一是傅里叶分步步进法,许多研究者采用法求解海面、平 值结果尚有令人不太满意之处,理论公式推导也不 够完善。因此,用双向抛物方程法预测刃峰环境的 地面以及不规则地形的电波传播 。 。但是,传统 的抛物方程法是一种单向算法,其缺点是只考虑了 无线电波传播的研究还有待于进一步深化。 前向传播波,而忽略了后向传播波 。通常情况 用双向抛物方程法结合傅里叶分步步进算法计 算了单个和多个刃峰环境的前向和后向传播波。算 下,前向传播波携带绝大部分能量,然而,在有障碍 物存在时,特别是在不规则大气条件下,障碍物的反 法从天线附近的初始场开始,沿距离方向前向步进, 当波遇到一个刃峰时被部分反射并沿与单向抛物方 射可能会引起较强的后向传播波。一个典型的例子 程轴向方向相反的后向前进。通过在刃峰处应用边 就是刃峰地形的电波传播问题,其中,刃峰会对入射 的电磁波产生明显的反射,如果用单向抛物方程法 界条件可以得到反射波的初始场。每次当波 前向 或后向 步进到达刃峰处时都分成两个分量,即前向 分析,则无法计算出刃峰前的后向反射场,得到的结 和后向传播波,它们分别在各自的传播方向上前进。 果与实际电磁场分布相比有较大误差,所以,有必要 在每一距离步进上,将后向传播波与前向传播波相 引入双向抛物方程法来计算刃峰前的总电磁场分 布 。 加,反复执行算法直至在整个计算域内达到预定的 收敛准则,从而能够处理刃峰 或障碍物 间的多重 双向抛物方程法最早用于计算海洋环境中声波 反射并最终得到计算域内的总场。论文第二部分给 能量的后向散射 州。后来,用双向抛物方程 出了窄角和宽角双向抛物方程法的公式推导,并总 法计算了多刃峰环境的前向和后向传播的电磁 结了刃峰地形下双向抛物方程算法的实现过程;第 场口 ,并与一致性绕射理论 和时域有限差分 收稿日期: ? ?基金项目:广州市科技计划项目 一 ;福建省教育厅科技项目联系人:王昆 :. . 第 期 王 昆等:多刃峰环境无线电波传播预测的双向抛物方程法 三部分给出了单刃峰和多刃峰环境的算例,并将双 一 向抛物方程和单向抛物方程的数 值计算结果做分析比较;最后得出结论。 一△ , 一 . .双向抛物方程法的推导与实现 一假设电磁场的时谐因子为 ,且所有的场分 式中: ? 一,引入了地球曲率的影响, 量与方位角无关, 为距离, 为高度,方位角对称, 则麦克斯韦方程组可简化为二维标量赫尔姆霍茨方 为地球半径; 和 是傅里叶变换对。边界条件 程,在直角坐标系中,方程为” 也可以用适当的傅里叶变换对表示,例如,水平极化 可以用正弦傅里叶变换对处理,而垂直极化可以用栅 一余弦傅里叶变换对处理,阻抗边界条件用混合傅里 根据微分算子理论对式进行算子的因式分 叶变换对处理 。 , 是变换域变量,通常称为垂直 解得 波数或空间频率,定义为 , 为电波到水平方 向的角度, 与所需计算的最大仰角相关。在宽√ 。 筹× 角双向抛物方程的解中,将窄角双向抛物方程 一√ 。 解中的自由空间传播算子一 用 一 一 来代替即可。 式中:第一个算子对应电波的后向传播;第二个算子 考虑如图 所示的单刃峰环境电波传播问题, 对应电波的前向传播。得到前向和后向抛物方程 单刃峰位于 ,其高度为 。从 一 处的初始场 表示前向, 表示后向 为 开始,用式将场在前向步进,直至场到达刃峰。 一 当前向场在 处遇到刃峰时被分解为前向和后向