论文部分内容阅读
电磁计算的数值方法如矩量法(MOM),有限元法(FEM),时域有限差分方法(FDTD)可以很好地解决电小尺寸物体的散射,但在计算电大物体的散射时,对计算机的配置要求过高。而射线跟踪、物理光学等高频近似方法则只能求解电大物体的散射,对中等尺寸的目标其计算结果误差较大。抛物线方程(PE:Parabolic Equation)方法则是联系这些方法的一座桥梁,它是波动方程的一种近似形式,假设电磁波能量在沿着抛物线轴向的锥形区域内传播。本文对抛物线方程方法进行了初步的研究和计算应用,主要工作包括以下几个方面:首先介绍了PE方法的原理,如二维前向抛物线方程,二维后向抛物线方程,柱坐标系下的抛物线方程,三维抛物线方程以及矢量抛物线方程。然后对上述的方法进行数值求解,因为方程是微分形式的,所以我们采用差分方法,在二维方面,介绍了自由空间的Crank-Nicolson差分格式和pade-(1,0)差分格式。接着介绍了完全匹配层(PML)的原理,并且推导了在PML中的Crank-Nicolson差分格式和pade-(1,0)差分格式。对于抛物线方程只在近轴方向小角度范围内准确的特点,我们引入了旋转抛物线方程方法,利用旋转抛物线轴向,并对入射方向和物体边界条件进行重新设置,可以计算出各个方向的RCS。在三维物体的散射问题上,也推导了其差分格式的公式,主要是pade-(1,0)差分格式,还给出了PML中的差分公式。三维问题的一个显著特点就是未知量很大,在求解时有一定困难,对此,我们介绍了该大型稀疏矩阵的解法。最后,介绍了Double Pass方法。在本文的最后,我们还给出了一种将描述目标表面形状的三角面元数据转换为适用于抛物线方程的立方网格生成技术,利用这种技术,可以很方便地将ANSYS生成的三角面元转化为立方网格,不用通过设计图、制造图手工或半手工地去构造模型,减少了工作量;建模方法和程序均具有通用性,一次编程可完成相应情况的任何建模。抛物线方程方法是散射目标分析和雷达设计的重要理论依据,对于目标雷达探测和目标识别都有重要意义。