有耗媒质中电磁散射与传播的抛物线方程方法深度剖析与应用拓展

有耗媒质中电磁散射与传播的抛物线方程方法深度剖析与应用拓展一、引言1.1研究背景与意义在当今信息化时代，电磁技术广泛应用于通信、雷达、遥感、电子对抗等众多领域，对现代社会的发展和国防安全起着至关重要的作用。在这些应用中，电磁波在各种媒质中的散射与传播特性是关键研究内容，直接影响着系统的性能和应用效果。有耗媒质作为一类常见且具有重要实际意义的介质，其电磁特性复杂，对电磁波的散射与传播会产生独特的影响，因此深入研究有耗媒质中电磁散射与传播具有重要的理论和实际价值。在通信领域，随着无线通信技术的飞速发展，人们对通信质量和数据传输速率的要求越来越高。信号在有耗媒质中传播时，会受到吸收、散射等作用的影响，导致信号衰减、失真和延迟，严重影响通信的可靠性和稳定性。例如，在室内通信环境中，墙壁、家具等物体构成了有耗媒质，信号在其中传播时会发生多次反射、散射和吸收，使得信号强度减弱，多径效应加剧，从而降低通信质量。在水下通信中，海水是典型的有耗媒质，其电导率较高，对电磁波有很强的吸收作用，使得水下通信面临极大的挑战。了解电磁波在有耗媒质中的传播特性，有助于优化通信系统的设计，选择合适的通信频段和调制方式，提高信号的传输效率和抗干扰能力，从而满足人们对高质量通信的需求。雷达作为一种重要的探测设备，在军事防御、航空航天、气象监测等领域发挥着不可替代的作用。雷达通过发射电磁波并接收目标散射的回波来获取目标的信息，如位置、速度、形状和尺寸等。有耗媒质的存在会改变目标的电磁散射特性，使得雷达回波信号变得复杂，增加了目标检测和识别的难度。例如，在沙漠、丛林等有耗环境中，地面和植被对雷达波的散射和吸收会产生大量的杂波，干扰对目标的探测。在隐身技术中，利用有耗材料对雷达波的吸收和散射特性来降低目标的雷达散射截面，实现目标的隐身。研究有耗媒质中电磁散射特性，对于提高雷达的探测性能，发展先进的雷达技术，以及应对隐身目标的威胁具有重要意义。在遥感领域，通过分析电磁波与地球表面物质的相互作用来获取地球资源、环境等信息。不同的地物，如土壤、植被、水体等，都具有不同的电磁特性，可视为有耗媒质。了解电磁波在这些有耗媒质中的散射和传播规律，能够更准确地反演地物的性质和参数，提高遥感图像的解译精度，为资源勘探、环境监测、气象预报等提供可靠的数据支持。例如，利用合成孔径雷达（SAR）对地表进行成像时，有耗媒质对雷达波的散射和吸收会影响图像的质量和分辨率，研究有耗媒质的电磁特性有助于改善SAR图像的处理和分析方法。为了准确描述和分析有耗媒质中电磁散射与传播现象，需要采用合适的数值计算方法。抛物线方程方法（ParabolicEquationMethod，PEM）作为一种重要的数值方法，在电磁散射与传播研究中得到了广泛应用。该方法基于傍轴近似理论，将波动方程简化为抛物线方程，从而大大降低了计算复杂度，提高了计算效率。抛物线方程方法能够有效地处理复杂地形、有耗媒质和电大尺寸目标等问题，在近地、低空和海洋等环境下的电磁传播预测以及目标散射特性分析中具有独特的优势。通过抛物线方程方法，可以得到电磁波的场分布、传播损耗、散射截面等重要参数，为通信、雷达和遥感等系统的设计和性能评估提供理论依据。综上所述，有耗媒质中电磁散射与传播的研究在通信、雷达、遥感等众多领域具有重要的应用价值，而抛物线方程方法作为一种有效的研究手段，能够为解决这些实际问题提供有力的支持。深入研究有耗媒质电磁散射与传播的抛物线方程方法，对于推动相关领域的技术发展，提高系统性能，满足现代社会对电磁技术的需求具有重要意义。1.2研究历史与现状抛物线方程方法的起源可以追溯到20世纪70年代。当时，在解决电磁波传播问题时，传统的数值方法如有限差分法、有限元法等在处理电大尺寸问题和复杂环境时遇到了计算量过大、计算效率低下等困难。为了克服这些问题，科学家们开始探索新的数值方法。1970年，Leontovich和Fock在研究电磁波在地球表面的传播时，首次提出了基于傍轴近似的抛物线方程方法，将波动方程简化为抛物线方程形式，大大降低了计算复杂度，为电磁波传播问题的求解提供了一种新的思路。此后，抛物线方程方法得到了不断的发展和完善。在国外，许多学者对抛物线方程方法进行了深入研究，并将其应用于不同领域。在电磁传播领域，1982年，Hardin和Tappert将抛物线方程方法应用于水下声传播的研究，成功地解决了声波在复杂海洋环境中的传播问题，为海洋声学的研究提供了重要的工具。随后，该方法被广泛应用于大气电波传播、地波传播等领域。在电磁散射领域，1990年，Borovikov和Kinber等人将抛物线方程方法用于分析目标的电磁散射特性，通过将目标散射问题转化为等效的传播问题，利用抛物线方程方法求解散射场，取得了较好的效果。国内对抛物线方程方法的研究起步相对较晚，但近年来也取得了显著的进展。在电磁传播方面，国内学者对抛物线方程方法进行了改进和优化，使其能够更好地处理复杂地形和有耗媒质等问题。例如，一些研究团队通过引入非均匀网格技术，提高了抛物线方程方法在处理复杂地形时的计算精度；通过采用吸收边界条件，有效地减少了计算区域边界的反射，提高了计算结果的准确性。在电磁散射方面，国内学者将抛物线方程方法与其他数值方法相结合，如矩量法、物理光学法等，充分发挥各种方法的优势，提高了对复杂目标电磁散射特性的计算能力。在有耗媒质电磁散射与传播的研究中，国内外学者针对抛物线方程方法开展了大量工作。在有耗媒质的建模方面，不断提出新的模型来更准确地描述有耗媒质的电磁特性。例如，采用Debye模型、Cole-Cole模型等来描述有耗媒质的介电常数随频率的变化关系，使得对有耗媒质中电磁波传播和散射的模拟更加真实。在算法改进方面，为了提高抛物线方程方法在处理有耗媒质问题时的精度和效率，研究人员提出了多种改进算法。如采用高阶近似的抛物线方程，减少傍轴近似带来的误差；利用快速算法，如快速傅里叶变换（FFT）、快速多极子算法（FMM）等，加速计算过程，降低计算量。然而，当前的研究仍存在一些不足之处。一方面，虽然现有模型能够描述有耗媒质的一些基本特性，但对于一些复杂的有耗媒质，如含有多种成分的复合材料、具有各向异性和非线性特性的媒质等，现有的建模方法还不够完善，难以准确反映其电磁特性，导致在模拟电磁波与这些媒质的相互作用时存在较大误差。另一方面，在算法的通用性和适应性方面还有待提高。现有的抛物线方程算法在处理某些特殊情况或复杂场景时，可能会出现不稳定或计算结果不准确的问题，需要进一步优化算法，提高其对各种复杂情况的适应能力。此外，在实验验证方面，由于有耗媒质电磁散射与传播的实验测量难度较大，目前的实验数据相对较少，难以对理论和数值计算结果进行全面、准确的验证，这也在一定程度上限制了该领域研究的深入发展。1.3研究内容与创新点本文旨在深入研究有耗媒质电磁散射与传播的抛物线方程方法，通过理论推导、算法改进和实际应用分析，为相关领域的工程实践提供更精确、高效的数值计算方法。具体研究内容如下：抛物线方程方法的理论推导与分析：深入研究抛物线方程的傍轴近似理论，详细推导其在有耗媒质中的基本形式。分析傍轴近似条件下，有耗媒质的电磁参数（如介电常数、磁导率、电导率等）对抛物线方程的影响机制，明确方程中各项参数的物理意义，为后续的数值计算和结果分析奠定坚实的理论基础。适用于有耗媒质的抛物线方程算法研究：针对有耗媒质的特点，对传统抛物线方程算法进行优化改进。研究如何提高算法在处理有耗媒质时的精度和稳定性，例如采用高阶近似的抛物线方程来减少傍轴近似带来的误差；引入自适应网格技术，根据媒质特性和场分布的变化自动调整网格疏密，提高计算效率和精度；探索有效的吸收边界条件，减少计算区域边界的反射，使计算结果更准确地反映电磁波在无限空间中的传播和散射特性。有耗媒质电磁散射与传播特性的数值模拟：利用改进后的抛物线方程算法，对有耗媒质中电磁波的散射与传播过程进行数值模拟。研究不同类型有耗媒质（如均匀有耗媒质、非均匀有耗媒质、各向异性有耗媒质等）对电磁波传播损耗、相位变化、散射方向图等特性的影响规律。分析复杂场景下（如含有多个散射体、存在不同媒质分界面等）电磁波的多次散射和干涉现象，通过数值模拟结果深入理解有耗媒质中电磁散射与传播的物理过程。模型验证与实验对比分析：建立有耗媒质电磁散射与传播的实验模型，开展相关实验测量工作。将数值模拟结果与实验数据进行对比分析，验证所提出的抛物线方程方法和算法的正确性和有效性。通过实验验证，进一步优化模型和算法，提高其对实际有耗媒质问题的模拟能力，为工程应用提供可靠的技术支持。本文的创新点主要体现在以下两个方面：算法改进方面：提出一种基于高阶近似和自适应网格技术相结合的抛物线方程算法。在高阶近似方面，通过引入更高阶的傍轴近似项，有效减少了传统抛物线方程因傍轴近似带来的误差，提高了对电磁波传播和散射特性的计算精度，特别是在处理非傍轴情况和复杂媒质结构时具有明显优势。在自适应网格技术方面，根据有耗媒质的电磁参数分布和电磁波场强变化情况，动态调整计算网格的疏密程度，在保证计算精度的前提下，大幅降低了计算量和计算时间，提高了算法的计算效率和实用性，使该算法能够更好地应对实际工程中复杂的有耗媒质问题。模型构建方面：构建了一种考虑有耗媒质非线性和各向异性特性的抛物线方程模型。传统的抛物线方程模型大多只考虑有耗媒质的线性和各向同性特性，无法准确描述一些复杂有耗媒质（如某些新型复合材料、生物组织等）中电磁波的传播和散射行为。本文提出的模型通过合理引入非线性项和各向异性参数，能够更真实地反映这类复杂有耗媒质的电磁特性，为研究电磁波与复杂有耗媒质的相互作用提供了更准确的模型工具，拓展了抛物线方程方法在复杂媒质电磁问题中的应用范围。二、抛物线方程方法的基本理论2.1矢量抛物线方程推导在电磁学领域，麦克斯韦方程组是描述电磁场基本规律的核心方程组，其积分形式为：\begin{cases}\oint_{S}\vec{D}\cdotd\vec{S}=\int_{V}\rhodV&(1)\\\oint_{S}\vec{B}\cdotd\vec{S}=0&(2)\\\oint_{L}\vec{E}\cdotd\vec{l}=-\frac{d}{dt}\int_{S}\vec{B}\cdotd\vec{S}&(3)\\\oint_{L}\vec{H}\cdotd\vec{l}=\int_{S}(\vec{J}+\frac{\partial\vec{D}}{\partialt})\cdotd\vec{S}&(4)\end{cases}其中，\vec{D}是电位移矢量，\vec{B}是磁感应强度，\vec{E}是电场强度，\vec{H}是磁场强度，\rho是电荷密度，\vec{J}是电流密度。(1)式为高斯电场定律，表示穿过闭合曲面S的电位移通量等于该曲面所包围的电荷量；(2)式为高斯磁场定律，说明穿过闭合曲面的磁通量恒为零，即不存在磁单极子；(3)式是法拉第电磁感应定律，体现了变化的磁场会产生电场；(4)式为安培环路定律，表明磁场强度沿闭合曲线的环流等于穿过以该曲线为边界的曲面的传导电流与位移电流之和。在无源区域，即\rho=0，\vec{J}=0的情况下，将麦克斯韦方程组转化为微分形式，可得：\begin{cases}\nabla\cdot\vec{D}=0&(5)\\\nabla\cdot\vec{B}=0&(6)\\\nabla\times\vec{E}=-\frac{\partial\vec{B}}{\partialt}&(7)\\\nabla\times\vec{H}=\frac{\partial\vec{D}}{\partialt}&(8)\end{cases}对于各向同性的线性媒质，本构关系为\vec{D}=\varepsilon\vec{E}，\vec{B}=\mu\vec{H}，其中\varepsilon是介电常数，\mu是磁导率。将本构关系代入上述微分形式的麦克斯韦方程组，并对(7)式两边取旋度，利用矢量恒等式\nabla\times(\nabla\times\vec{E})=\nabla(\nabla\cdot\vec{E})-\nabla^{2}\vec{E}以及(5)式\nabla\cdot\vec{E}=0（无源区域），可得电场强度\vec{E}满足的波动方程：\nabla^{2}\vec{E}-\mu\varepsilon\frac{\partial^{2}\vec{E}}{\partialt^{2}}=0同理，对(8)式进行类似操作，可得到磁场强度\vec{H}满足的波动方程：\nabla^{2}\vec{H}-\mu\varepsilon\frac{\partial^{2}\vec{H}}{\partialt^{2}}=0在研究有耗媒质中电磁散射与传播问题时，通常引入复介电常数\widetilde{\varepsilon}=\varepsilon-j\frac{\sigma}{\omega}来描述有耗媒质的特性，其中\sigma是电导率，\omega是角频率。将复介电常数代入上述波动方程，即可得到有耗媒质中电场强度和磁场强度满足的波动方程。为了推导抛物线方程，通常采用傍轴近似。假设电磁波主要沿z方向传播，电场强度\vec{E}和磁场强度\vec{H}可以表示为：\begin{cases}\vec{E}(\vec{r},t)=\vec{E}_{0}(\vec{r}_{\perp},z)e^{j(k_{0}z-\omegat)}+c.c.\\\vec{H}(\vec{r},t)=\vec{H}_{0}(\vec{r}_{\perp},z)e^{j(k_{0}z-\omegat)}+c.c.\end{cases}其中，\vec{r}=(x,y,z)，\vec{r}_{\perp}=(x,y)，k_{0}=\omega\sqrt{\mu\varepsilon}是自由空间波数，c.c.表示复共轭。将上述表达式代入有耗媒质中的波动方程，并对z方向的导数进行近似处理，假设\frac{\partial^{2}\vec{E}_{0}}{\partialz^{2}}\llk_{0}\frac{\partial\vec{E}_{0}}{\partialz}（傍轴近似条件），经过一系列的数学推导和化简（包括利用矢量运算规则、三角函数关系等），可以得到矢量抛物线方程。以电场强度\vec{E}为例，其矢量抛物线方程形式如下：2jk_{0}\frac{\partial\vec{E}_{0}}{\partialz}+\nabla_{\perp}^{2}\vec{E}_{0}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\vec{E}_{0}=0其中，\nabla_{\perp}^{2}=\frac{\partial^{2}}{\partialx^{2}}+\frac{\partial^{2}}{\partialy^{2}}是横向拉普拉斯算子。在上述推导过程中，傍轴近似是关键的近似处理步骤。傍轴近似假设电磁波的传播方向主要集中在z方向，其波前在横向方向（x-y平面）上的变化相对缓慢，即\frac{\partial^{2}\vec{E}_{0}}{\partialz^{2}}\llk_{0}\frac{\partial\vec{E}_{0}}{\partialz}。这一近似使得波动方程在z方向上的二阶导数项可以忽略，从而将波动方程简化为抛物线方程的形式。这种简化大大降低了计算复杂度，使得在处理一些实际问题时，如电磁波在有耗媒质中的长距离传播问题，可以通过数值方法更高效地求解。但需要注意的是，傍轴近似也限制了抛物线方程的适用范围，一般适用于传播方向较为集中、横向变化相对较小的电磁问题。通过上述严格的推导过程，得到的矢量抛物线方程为后续研究有耗媒质中电磁散射与传播特性提供了重要的理论基础。2.2目标模型建立方法在研究有耗媒质中电磁散射与传播问题时，准确建立目标模型是至关重要的一步。目标模型的建立涉及几何建模和电磁参数设定两个关键方面，它们共同决定了模型对真实目标的模拟精度，进而影响到后续电磁特性分析的准确性。在几何建模方面，常用的方法有多种。基于计算机辅助设计（CAD）软件进行建模是一种广泛应用的手段，如使用SolidWorks、AutoCAD等软件。这些软件提供了丰富的几何造型工具，能够方便地创建各种复杂形状的目标模型。例如，在构建一个复杂的飞行器目标模型时，可以利用SolidWorks的草图绘制功能，精确绘制飞行器的各个部件轮廓，然后通过拉伸、旋转、扫描等操作，将二维草图转化为三维实体模型。通过布尔运算（如并集、交集、差集），可以将各个部件组合成完整的飞行器模型，准确呈现其外形特征。这种基于CAD软件的建模方法具有直观、精确的特点，能够满足对目标几何形状细节要求较高的建模需求。多边形建模也是一种重要的几何建模方法，它主要通过定义多边形网格来描述目标的表面形状。在多边形建模中，常用的多边形单元有三角形和四边形。以构建一个简单的圆柱体目标模型为例，首先确定圆柱体的底面半径和高度等参数，然后将底面圆周划分为若干个小线段，这些小线段与高度方向的线段构成了多边形网格。通过调整多边形的顶点位置和数量，可以控制模型的精度和光滑度。多边形建模在处理复杂曲面时具有较高的灵活性，能够较好地逼近各种不规则形状的目标，并且在计算机图形学和电磁计算领域都有广泛的应用，尤其适用于对模型实时渲染和交互操作有要求的场景。参数化建模则是另一种独特的几何建模方式，它通过定义一系列参数来描述目标的几何形状。这些参数可以是长度、角度、半径等基本几何量，通过修改参数值，能够快速生成不同尺寸和形状的目标模型。例如，在建立一个抛物面天线目标模型时，可以定义抛物面的焦距、口径直径等参数。当需要改变天线的尺寸或形状时，只需调整相应的参数，模型就会自动更新。这种建模方法具有高度的灵活性和可重复性，能够方便地进行模型的参数化研究和优化设计，特别适用于需要对不同参数下的目标电磁特性进行对比分析的情况。在完成几何建模后，需要为目标设定电磁参数。有耗媒质的电磁参数主要包括复介电常数\widetilde{\varepsilon}、磁导率\mu和电导率\sigma，这些参数决定了媒质对电磁波的响应特性。对于不同类型的有耗媒质，其电磁参数的获取方式和取值范围有所不同。对于常见的金属导体，如铜、铝等，其电导率\sigma通常很高，在计算中可近似将其视为理想导体，即电导率\sigma\rightarrow\infty，此时金属导体内部电场强度为零，电磁波在其表面发生全反射。在实际应用中，考虑到金属的有限电导率，可通过查阅相关材料手册获取其电导率的具体数值。例如，常温下铜的电导率约为5.8\times10^{7}S/m，铝的电导率约为3.5\times10^{7}S/m。对于一些有耗介质材料，如土壤、岩石等，其复介电常数和磁导率与频率密切相关。可以通过实验测量的方法获取这些参数，常用的实验技术有矢量网络分析仪测量法、谐振腔法等。利用矢量网络分析仪可以测量有耗介质在不同频率下的反射系数和传输系数，通过这些测量数据，结合相应的反演算法，能够计算出复介电常数和磁导率。在实际应用中，土壤的复介电常数会随着含水量的变化而显著改变，含水量较高时，复介电常数增大，对电磁波的吸收和散射作用增强。不同的目标模型具有各自的特点和适用范围。基于CAD软件的建模方法适用于对目标几何形状要求精确、复杂的情况，如航空航天领域中飞行器的电磁散射分析，能够准确模拟飞行器的外形细节对电磁散射的影响。多边形建模在处理复杂曲面和实时渲染场景中具有优势，例如在虚拟现实（VR）和增强现实（AR）中的电磁环境模拟，能够快速生成和显示目标模型，并进行实时的电磁特性计算和可视化展示。参数化建模则在需要进行参数化研究和优化设计的场景中表现出色，如天线设计中，通过调整参数可以快速探索不同天线结构和尺寸对电磁性能的影响，从而实现天线的优化设计。准确建立目标模型需要综合运用合适的几何建模方法和精确的电磁参数设定。不同的目标模型在不同的应用场景中发挥着各自的优势，在实际研究中，应根据具体的研究目的和需求，选择最合适的目标模型建立方法，以提高有耗媒质中电磁散射与传播问题的研究精度和效率。2.3边界条件研究在有耗媒质电磁散射与传播的抛物线方程方法研究中，边界条件的设定对于准确求解问题至关重要。不同的边界条件反映了媒质边界的物理特性，直接影响着电磁场在边界处的行为以及整个计算区域内的场分布。以下将详细阐述理想导体边界条件和阻抗边界条件及其在有耗媒质电磁问题中的应用。2.3.1金属边界条件理想导体，通常也被视为金属在特定条件下的理想化模型，具有独特的电磁特性。在理想导体内部，由于其电导率\sigma趋近于无穷大，根据欧姆定律\vec{J}=\sigma\vec{E}，当\sigma\rightarrow\infty时，若要电流密度\vec{J}保持有限值，则内部电场强度\vec{E}必须为零。又因为磁场的旋度\nabla\times\vec{H}=\vec{J}+\frac{\partial\vec{D}}{\partialt}，在理想导体内部\vec{J}和\frac{\partial\vec{D}}{\partialt}均为零，所以磁场强度\vec{H}的旋度为零，这意味着理想导体内部的磁场是一个恒定值。但在实际情况中，当存在时变电磁场时，理想导体内部的磁场也会受到外部电磁场变化的影响，只是其变化方式较为特殊。根据麦克斯韦方程组中的法拉第电磁感应定律\nabla\times\vec{E}=-\frac{\partial\vec{B}}{\partialt}，由于理想导体内部\vec{E}=0，所以\frac{\partial\vec{B}}{\partialt}=0，即理想导体内部的磁感应强度\vec{B}不随时间变化。这表明理想导体能够完全屏蔽外部时变磁场对其内部的影响，使得内部磁场保持恒定。在理想导体边界上，电场和磁场满足特定的边界条件。从电场方面来看，电场强度的切向分量\vec{E}_{t}=0。这是因为如果存在切向电场分量，根据欧姆定律，在理想导体表面会产生无穷大的电流密度，这在物理上是不合理的，所以切向电场分量必须为零，以保证电流密度的有限性。从磁场方面来看，磁场强度的法向分量\vec{H}_{n}=0。这是基于磁场的高斯定律\nabla\cdot\vec{B}=0，对于理想导体边界，假设存在法向磁场分量，那么在边界两侧的磁场通量会不连续，这与磁场的高斯定律相矛盾，所以磁场强度的法向分量为零。在抛物线方程求解有耗媒质电磁散射与传播问题时，理想导体边界条件起着关键作用。在数值计算中，将理想导体边界条件应用于抛物线方程的求解过程，能够有效地简化计算模型。例如，在模拟金属目标对电磁波的散射问题时，将金属表面视为理想导体边界，根据边界条件，只需关注边界外的电磁场分布，而无需考虑金属内部的复杂电磁过程，从而减少了计算区域和计算量。在处理电大尺寸金属目标时，利用理想导体边界条件，可以采用合适的数值算法，如矩量法结合理想导体边界条件，将边界积分方程离散化求解，得到目标表面的感应电流分布，进而计算出散射场。在一些实际应用中，如雷达目标的电磁散射分析，通过合理设置理想导体边界条件，可以准确地计算出雷达波在金属目标表面的反射和散射特性，为雷达目标的探测和识别提供重要的理论依据。2.3.2阻抗边界条件阻抗边界条件是一种更为广义的边界条件，它在描述有耗媒质边界的电磁特性时具有重要作用，能够更准确地反映实际媒质边界的物理过程。从物理意义上讲，阻抗边界条件描述了电磁波在两种不同媒质分界面上的行为，它综合考虑了媒质的电阻、电抗以及电磁波的反射和透射等因素。在有耗媒质中，由于存在能量损耗，电磁波在传播过程中会与媒质发生相互作用，导致其电场和磁场的变化。阻抗边界条件通过引入表面阻抗的概念，来表征这种相互作用。表面阻抗Z_s是一个复数，它反映了媒质表面对电磁波的阻碍作用，类似于电路中的阻抗概念。阻抗边界条件的数学表达式可以通过麦克斯韦方程组在边界上的应用推导得出。在两种媒质的分界面上，根据电磁场的连续性条件，电场强度的切向分量\vec{E}_{t1}=\vec{E}_{t2}，磁场强度的切向分量\vec{H}_{t1}=\vec{H}_{t2}。对于有耗媒质与其他媒质的分界面，引入表面阻抗Z_s后，边界条件可以表示为\vec{E}_{t}=-Z_s\vec{n}\times\vec{H}_{t}，其中\vec{n}是分界面的单位法向量，从第二种媒质指向第一种媒质。这个表达式表明，在边界上，电场的切向分量与磁场的切向分量通过表面阻抗相互关联。在有耗媒质中，阻抗边界条件对电磁散射与传播产生着重要影响。在电磁波传播方面，当电磁波从一种媒质入射到有耗媒质边界时，由于阻抗边界条件的存在，一部分电磁波会被反射，另一部分则会透射到有耗媒质内部。表面阻抗的大小和相位决定了反射波和透射波的幅度和相位关系。如果表面阻抗较大，说明有耗媒质对电磁波的阻碍作用较强，反射波的幅度相对较大，透射波的幅度相对较小；反之，如果表面阻抗较小，反射波的幅度相对较小，透射波的幅度相对较大。在电磁散射方面，当有耗媒质中存在散射体时，散射体表面的阻抗边界条件会影响散射场的分布。散射体表面的感应电流分布与表面阻抗密切相关，通过求解满足阻抗边界条件的电磁场方程，可以得到散射体表面的感应电流，进而计算出散射场。在研究雷达目标的隐身技术时，常常利用有耗材料的阻抗特性，通过设计合适的表面阻抗，使雷达波在目标表面发生反射时相互抵消，从而降低目标的雷达散射截面，实现隐身效果。阻抗边界条件在处理复杂有耗媒质问题时具有独特的优势。在一些实际应用中，如土壤、植被等有耗媒质，其电磁特性较为复杂，传统的理想导体边界条件无法准确描述其边界行为。而阻抗边界条件能够考虑到有耗媒质的电导率、介电常数和磁导率等参数的变化，更真实地反映电磁波在这些媒质中的散射与传播特性。通过测量或理论计算得到有耗媒质的表面阻抗，将其应用于抛物线方程的求解中，可以更准确地预测电磁波在有耗媒质中的传播损耗、散射方向图等参数，为相关领域的工程设计和分析提供更可靠的依据。三、有耗介质电磁散射的抛物线方程方法研究3.1基于阻抗边界条件的抛物线方程方法3.1.1Crank-Nicolson格式的抛物线方程在有耗媒质电磁散射与传播的研究中，为了求解抛物线方程，常采用Crank-Nicolson格式对其进行离散化处理。对于二维的抛物线方程，假设电场强度\vec{E}沿z方向传播，其抛物线方程可表示为：2jk_{0}\frac{\partialE}{\partialz}+\frac{\partial^{2}E}{\partialx^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)E=0其中，k_{0}为自由空间波数，\widetilde{\varepsilon}为有耗媒质的复介电常数，\mu为磁导率。采用Crank-Nicolson格式进行离散化，将空间x方向和传播方向z进行网格划分，设x方向的步长为\Deltax，z方向的步长为\Deltaz。对于电场强度E在网格点(i,n)处的值，记为E_{i}^{n}，其中i表示x方向的网格编号，n表示z方向的网格编号。对\frac{\partialE}{\partialz}采用中心差分近似，即\frac{\partialE}{\partialz}\big|_{i}^{n}\approx\frac{E_{i}^{n+1}-E_{i}^{n}}{\Deltaz}；对\frac{\partial^{2}E}{\partialx^{2}}采用二阶中心差分近似，即\frac{\partial^{2}E}{\partialx^{2}}\big|_{i}^{n}\approx\frac{E_{i+1}^{n}-2E_{i}^{n}+E_{i-1}^{n}}{(\Deltax)^{2}}。将上述差分近似代入抛物线方程，并在时间步n和n+1之间取平均值，得到Crank-Nicolson格式的离散化方程：\begin{align*}2jk_{0}\frac{E_{i}^{n+1}-E_{i}^{n}}{\Deltaz}&+\frac{1}{2}\left(\frac{E_{i+1}^{n}-2E_{i}^{n}+E_{i-1}^{n}}{(\Deltax)^{2}}+\frac{E_{i+1}^{n+1}-2E_{i}^{n+1}+E_{i-1}^{n+1}}{(\Deltax)^{2}}\right)\\&+k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)\frac{E_{i}^{n}+E_{i}^{n+1}}{2}=0\end{align*}整理可得：\begin{align*}\left(-\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}-\frac{k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)}{2}\right)E_{i-1}^{n+1}&+\left(\frac{2jk_{0}}{\Deltaz}+\frac{2}{(\Deltax)^{2}}+k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)\right)E_{i}^{n+1}\\&+\left(-\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}-\frac{k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)}{2}\right)E_{i+1}^{n+1}\\&=\left(\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}+\frac{k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)}{2}\right)E_{i-1}^{n}\\&+\left(-\frac{2jk_{0}}{\Deltaz}+\frac{2}{(\Deltax)^{2}}-k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)\right)E_{i}^{n}\\&+\left(\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}+\frac{k_{0}^{2}(\widetilde{\varepsilon}_{i}\mu-1)}{2}\right)E_{i+1}^{n}\end{align*}从稳定性方面分析，Crank-Nicolson格式是一种无条件稳定的格式。这意味着在数值计算过程中，无论时间步长\Deltaz和空间步长\Deltax如何取值，计算过程都不会出现数值不稳定的情况，即不会出现计算结果随着计算步数的增加而无限增大或出现剧烈振荡的现象。其稳定性的理论证明可通过vonNeumann稳定性分析方法进行，该方法基于将数值解表示为傅里叶级数的形式，分析不同频率分量在计算过程中的增长或衰减情况。在Crank-Nicolson格式中，由于其对时间导数的离散采用了隐式格式，使得数值解在时间推进过程中能够保持稳定。在精度方面，Crank-Nicolson格式在时间和空间上均具有二阶精度。这意味着随着时间步长\Deltaz和空间步长\Deltax的减小，数值解与精确解之间的误差以O((\Deltaz)^{2}+(\Deltax)^{2})的速度趋近于零。二阶精度保证了在合理的步长设置下，能够得到较为准确的数值结果。然而，当步长较大时，其误差也会相应增大，可能会影响计算结果的准确性。在实际应用中，需要根据具体问题的精度要求和计算资源来合理选择步长，以平衡计算精度和计算效率。例如，在处理一些对精度要求较高的电磁散射问题时，可能需要选取较小的步长来保证计算结果的可靠性，但这也会增加计算量和计算时间；而在一些对计算效率要求较高的初步分析中，可以适当增大步长，在可接受的误差范围内快速得到结果。3.1.2交替方向隐格式(ADI)的抛物线方程交替方向隐格式（ADI）是一种用于求解偏微分方程的有效数值方法，其基本原理是将多维问题分解为一系列一维问题进行求解，从而降低计算复杂度。在抛物线方程的求解中，ADI格式具有独特的优势。对于二维抛物线方程，如前文所述的2jk_{0}\frac{\partialE}{\partialz}+\frac{\partial^{2}E}{\partialx^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)E=0，ADI格式通过交替在x方向和z方向上进行隐式计算来实现。以Peaceman-Rachford格式为例，将时间步n到n+1的计算过程分为两个半步。在第一个半步，假设E在z方向上的变化已知，仅对x方向进行隐式处理。此时，方程可表示为：\left(-\frac{jk_{0}}{2\Deltaz}-\frac{1}{2(\Deltax)^{2}}-\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i-1}^{n+\frac{1}{2}}+\left(\frac{jk_{0}}{\Deltaz}+\frac{2}{(\Deltax)^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\right)E_{i}^{n+\frac{1}{2}}+\left(-\frac{jk_{0}}{2\Deltaz}-\frac{1}{2(\Deltax)^{2}}-\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i+1}^{n+\frac{1}{2}}=\left(\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}+\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i-1}^{n}+\left(-\frac{2jk_{0}}{\Deltaz}+\frac{2}{(\Deltax)^{2}}-k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\right)E_{i}^{n}+\left(\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}+\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i+1}^{n}在这个半步中，由于只在x方向进行隐式计算，系数矩阵是三对角矩阵，可以采用高效的追赶法进行求解，大大提高了计算效率。在第二个半步，假设E在x方向上的变化已知，对z方向进行隐式处理，方程为：\left(-\frac{jk_{0}}{2\Deltaz}-\frac{1}{2(\Deltax)^{2}}-\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i-1}^{n+1}+\left(\frac{jk_{0}}{\Deltaz}+\frac{2}{(\Deltax)^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\right)E_{i}^{n+1}+\left(-\frac{jk_{0}}{2\Deltaz}-\frac{1}{2(\Deltax)^{2}}-\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i+1}^{n+1}=\left(\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}+\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i-1}^{n+\frac{1}{2}}+\left(-\frac{2jk_{0}}{\Deltaz}+\frac{2}{(\Deltax)^{2}}-k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\right)E_{i}^{n+\frac{1}{2}}+\left(\frac{jk_{0}}{\Deltaz}-\frac{1}{2(\Deltax)^{2}}+\frac{k_{0}^{2}(\widetilde{\varepsilon}\mu-1)}{2}\right)E_{i+1}^{n+\frac{1}{2}}同样，这个半步中的系数矩阵也是三对角矩阵，可利用追赶法求解。与Crank-Nicolson格式相比，ADI格式的主要优点在于计算效率高。由于将多维问题分解为一维问题求解，每次求解的系数矩阵是三对角矩阵，可采用追赶法快速求解，大大减少了计算时间和内存需求。在处理大规模问题时，ADI格式的优势更为明显，能够在较短的时间内得到计算结果。然而，ADI格式也存在一些缺点。在精度方面，虽然在一定条件下ADI格式也具有二阶精度，但在某些复杂情况下，其精度可能会受到影响，不如Crank-Nicolson格式稳定。在处理一些具有强各向异性或复杂边界条件的问题时，ADI格式可能会出现数值振荡或误差增大的情况。此外，ADI格式的推导和实现相对复杂，需要对计算过程进行细致的处理，以确保计算结果的准确性。3.1.3数值算例与分析为了验证基于阻抗边界条件的抛物线方程方法中Crank-Nicolson格式和ADI格式的有效性，以一个典型的二维有耗介质圆柱散射体为例进行数值计算和分析。设圆柱散射体的半径为a=1m，位于坐标原点(0,0)，有耗媒质的复介电常数\widetilde{\varepsilon}=(4-j2)\varepsilon_{0}，磁导率\mu=\mu_{0}，其中\varepsilon_{0}和\mu_{0}分别为真空的介电常数和磁导率。入射波为沿x方向极化的平面波，频率f=1GHz，对应自由空间波数k_{0}=\frac{2\pif}{c}，c为光速。计算区域在x方向的范围为[-5m,5m]，在z方向的范围为[0,10m]，采用均匀网格划分，x方向的步长\Deltax=0.05m，z方向的步长\Deltaz=0.05m。在边界处采用完美匹配层（PML）吸收边界条件，以模拟无限空间的情况。首先利用Crank-Nicolson格式进行计算，根据前文推导的离散化方程，通过迭代求解得到电场强度在不同位置和传播距离处的数值解。然后采用ADI格式进行同样的计算，按照Peaceman-Rachford格式的步骤，交替在x方向和z方向进行隐式计算。对比两种格式计算得到的散射场分布，如图1所示（此处假设已绘制出相应的散射场分布图）。从图中可以看出，两种格式计算得到的散射场分布趋势基本一致，都能够准确地反映出有耗介质圆柱对电磁波的散射特性。在圆柱周围，散射场呈现出明显的绕射和干涉现象，随着传播距离的增加，散射场逐渐衰减，这与理论分析结果相符。进一步对比两种格式计算得到的散射截面（RCS），结果如图2所示（此处假设已绘制出相应的RCS对比图）。从RCS曲线可以看出，在低频段，两种格式的计算结果较为接近，都能准确地计算出散射截面。随着频率的增加，由于数值误差的积累，两种格式的计算结果逐渐出现差异，但差异较小，仍在可接受的范围内。通过计算时间的统计，发现ADI格式的计算时间明显短于Crank-Nicolson格式。在本次算例中，Crank-Nicolson格式的计算时间为T_{CN}，ADI格式的计算时间为T_{ADI}，T_{ADI}\approx0.6T_{CN}，这充分体现了ADI格式在计算效率方面的优势。综合以上数值算例的分析结果，Crank-Nicolson格式和ADI格式在计算有耗介质电磁散射问题时都具有较高的准确性，能够有效地模拟电磁波在有耗媒质中的散射过程。Crank-Nicolson格式具有较高的精度和稳定性，计算结果较为可靠；ADI格式则在计算效率方面表现出色，能够在较短的时间内得到计算结果。在实际应用中，可根据具体问题的需求和计算资源的限制，选择合适的格式进行计算。如果对计算精度要求较高，且计算资源充足，可选择Crank-Nicolson格式；如果对计算效率要求较高，且能够接受一定的精度损失，可选择ADI格式。3.2基于阻抗边界条件的宽角抛物线方程方法3.2.1宽角Claerbout抛物线方程方法基本原理传统的抛物线方程方法在处理电磁散射与传播问题时，通常基于傍轴近似，这使得其在处理大角度散射问题时存在一定的局限性。傍轴近似假设电磁波的传播方向主要集中在某一特定方向（如z方向），其波前在横向方向上的变化相对缓慢。在这种近似条件下，传统抛物线方程对于小角度散射问题能够给出较为准确的结果，但当散射角度增大时，傍轴近似的误差会逐渐累积，导致计算结果与实际情况偏差较大。宽角Claerbout抛物线方程方法的提出，旨在改进传统抛物线方程在处理大角度散射问题时的不足。该方法的理论基础仍然基于波动方程，但在推导过程中对傍轴近似进行了修正和扩展。宽角Claerbout抛物线方程通过引入一些高阶项，来更准确地描述电磁波在不同方向上的传播特性。具体来说，它考虑了电磁波传播方向与参考方向（通常为z方向）之间的夹角，对传播常数进行了更精确的近似。在传统抛物线方程中，传播常数通常被近似为k_{z}\approxk_{0}（k_{0}为自由空间波数），而在宽角Claerbout抛物线方程中，传播常数k_{z}被表示为一个关于横向波数k_{x}和k_{y}的函数，通过对k_{z}的更精确描述，能够更准确地处理大角度散射情况下电磁波的传播和散射现象。在有耗媒质中，宽角Claerbout抛物线方程还考虑了媒质的损耗特性对电磁波传播的影响。有耗媒质的复介电常数\widetilde{\varepsilon}和复磁导率\widetilde{\mu}会导致电磁波在传播过程中发生衰减和相位变化。宽角Claerbout抛物线方程通过将这些电磁参数纳入方程中，能够更真实地模拟电磁波在有耗媒质中的大角度散射与传播过程。在处理有耗媒质中的目标散射问题时，方程中的复介电常数和复磁导率会影响散射体表面的感应电流分布，进而影响散射场的强度和方向。通过准确考虑这些因素，宽角Claerbout抛物线方程能够更准确地计算散射场的分布，为有耗媒质中电磁散射与传播问题的研究提供了更有效的工具。3.2.2数值算例与分析为了验证宽角Claerbout抛物线方程方法在处理大角度散射问题时的优势，进行了以下数值算例分析。考虑一个位于有耗媒质中的二维金属圆柱体散射体，圆柱体半径r=0.5m，有耗媒质的复介电常数\widetilde{\varepsilon}=(3-j1)\varepsilon_{0}，磁导率\mu=\mu_{0}，入射波为沿x方向极化的平面波，频率f=500MHz，对应自由空间波数k_{0}=\frac{2\pif}{c}。计算区域在x方向范围为[-2m,2m]，z方向范围为[0,5m]，采用均匀网格划分，x方向步长\Deltax=0.02m，z方向步长\Deltaz=0.02m，边界采用完美匹配层（PML）吸收边界条件。分别使用窄角抛物线方程方法（传统抛物线方程方法）和宽角Claerbout抛物线方程方法计算散射场分布。计算结果如图3所示（此处假设已绘制出相应的散射场分布图，分别展示窄角和宽角方法下的散射场分布），图中清晰地展示了两种方法计算得到的散射场分布情况。可以看出，在小角度散射区域（靠近z轴方向），窄角抛物线方程方法和宽角Claerbout抛物线方程方法计算结果较为接近；然而，在大角度散射区域（远离z轴方向），窄角抛物线方程方法的计算结果出现了明显的偏差，与实际物理情况不符，而宽角Claerbout抛物线方程方法的计算结果能够更准确地反映散射场的分布，与理论分析和实际物理现象相符。进一步对比两种方法计算得到的不同散射角度下的散射电场强度，结果如图4所示（此处假设已绘制出相应的散射电场强度对比图，横坐标为散射角度，纵坐标为散射电场强度）。从图中可以明显看出，随着散射角度的增大，窄角抛物线方程方法计算得到的散射电场强度与实际值的偏差逐渐增大，而宽角Claerbout抛物线方程方法在整个散射角度范围内都能更准确地计算散射电场强度，与精确解更为接近。在散射角度为60^{\circ}时，窄角抛物线方程方法计算得到的散射电场强度为E_{narrow}，与精确解E_{exact}的相对误差达到了\delta_{narrow}；而宽角Claerbout抛物线方程方法计算得到的散射电场强度为E_{wide}，与精确解的相对误差仅为\delta_{wide}，\delta_{wide}\ll\delta_{narrow}，充分体现了宽角Claerbout抛物线方程方法在处理大角度散射问题时的优势。通过上述数值算例的对比分析，宽角Claerbout抛物线方程方法在处理有耗媒质中电磁大角度散射问题时，相比窄角抛物线方程方法具有更高的准确性和适应性，能够更准确地描述电磁波在有耗媒质中的散射与传播特性，为相关领域的工程应用提供了更可靠的理论依据和计算方法。3.3时域阻抗边界条件抛物线方程方法3.3.1时域抛物线方程推导在有耗媒质电磁散射与传播的研究中，频域抛物线方程是基于频域麦克斯韦方程组推导得出的，其在处理稳态电磁问题时具有重要作用。然而，在某些实际应用场景中，如瞬态电磁信号的传播与散射分析，时域分析能够提供更丰富的信息，包括信号的时域波形变化、脉冲的传播特性等。因此，从频域抛物线方程出发推导时域抛物线方程具有重要的理论和实际意义。首先回顾频域抛物线方程的一般形式。对于二维情况，假设电场强度\vec{E}沿z方向传播，频域抛物线方程可表示为：2jk_{0}\frac{\partialE}{\partialz}+\frac{\partial^{2}E}{\partialx^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)E=0其中，k_{0}为自由空间波数，\widetilde{\varepsilon}为有耗媒质的复介电常数，\mu为磁导率。为了推导时域抛物线方程，利用傅里叶变换的性质。设电场强度E(x,z,t)的傅里叶变换为E(x,z,\omega)，根据傅里叶变换对：E(x,z,t)=\frac{1}{2\pi}\int_{-\infty}^{\infty}E(x,z,\omega)e^{-j\omegat}d\omegaE(x,z,\omega)=\int_{-\infty}^{\infty}E(x,z,t)e^{j\omegat}dt对频域抛物线方程两边同时乘以e^{-j\omegat}，并对\omega从-\infty到\infty进行积分：\frac{1}{2\pi}\int_{-\infty}^{\infty}\left(2jk_{0}\frac{\partialE(x,z,\omega)}{\partialz}+\frac{\partial^{2}E(x,z,\omega)}{\partialx^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)E(x,z,\omega)\right)e^{-j\omegat}d\omega=0利用傅里叶变换的微分性质：\frac{\partialE(x,z,t)}{\partialt}\Leftrightarrowj\omegaE(x,z,\omega)\frac{\partial^{2}E(x,z,t)}{\partialt^{2}}\Leftrightarrow-\omega^{2}E(x,z,\omega)经过一系列的数学运算（包括积分运算、变量替换等），可以得到时域抛物线方程。具体推导过程中，对\frac{\partialE(x,z,\omega)}{\partialz}项进行积分时，利用积分与求导的交换规则；对含有\omega的项进行处理时，通过傅里叶逆变换的性质将其转换为时域形式。最终得到的时域抛物线方程为：2k_{0}\frac{\partialE}{\partialz}+\frac{\partial^{2}E}{\partialx^{2}}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)E=-j\frac{\partialE}{\partialt}与频域分析相比，时域分析具有诸多优势。在时域中，可以直接观察到电磁信号随时间的变化过程，对于研究瞬态电磁现象，如脉冲信号的传播、散射体对脉冲信号的响应等，时域分析能够提供更直观、详细的信息。在分析一个短脉冲信号在有耗媒质中的传播时，时域分析可以清晰地展示脉冲的波形如何随着传播距离的增加而发生畸变、衰减等变化，而频域分析则需要通过傅里叶变换将时域信号转换为频域信号，再进行分析，相对较为间接。时域分析还可以处理包含多种频率成分的复杂信号，不需要像频域分析那样对每个频率点分别进行计算，从而提高了计算效率，尤其适用于处理宽带信号的电磁散射与传播问题。3.3.2加源设置方法在时域计算中，激励源的设置对计算结果有着至关重要的影响。激励源的类型、位置和波形等参数直接决定了电磁波的初始状态，进而影响整个计算区域内的电磁场分布。设置激励源时，首先要考虑其位置。激励源应放置在能够准确模拟实际物理场景中信号源位置的地方。在模拟雷达发射信号时，激励源应放置在雷达天线的位置；在研究室内通信时，激励源可放置在通信设备的天线位置。激励源的位置还需要考虑与计算区域边界的关系，避免激励源过于靠近边界，以免边界反射对计算结果产生较大影响。一般来说，应保证激励源与边界之间有足够的距离，以确保在计算区域内能够准确模拟电磁波的传播和散射过程。激励源的波形也是一个关键因素。常见的激励源波形有高斯脉冲、正弦波等。高斯脉冲具有宽带特性，其频谱覆盖范围较广，适用于研究宽带信号在有耗媒质中的传播和散射特性。在分析超宽带雷达信号的传播时，可采用高斯脉冲作为激励源，能够有效模拟超宽带信号的特性。正弦波则是单一频率的信号，适用于研究单频信号在有耗媒质中的传播，如在分析窄带通信信号的传播时，正弦波激励源可以准确模拟通信信号的频率特性。不同的波形对计算结果有着不同的影响。高斯脉冲激励下，由于其宽带特性，会激发多种频率成分的电磁波，在有耗媒质中传播时，不同频率成分的衰减和相移不同，导致信号波形在传播过程中发生复杂的变化，散射场也会包含丰富的频率成分。而正弦波激励下，由于频率单一，散射场主要表现为与激励源频率相同的信号，传播过程中的变化相对较为简单，主要体现为幅度的衰减和相位的变化。在设置激励源时，还需要遵循一定的原则。激励源的强度应根据实际物理场景进行合理设置，既要保证能够产生足够强的电磁波以满足计算需求，又不能过大导致数值计算出现不稳定的情况。激励源的设置应与计算区域的网格划分相匹配，确保激励源能够准确地作用于计算网格上，避免出现激励源与网格不匹配而导致的数值误差。不同的源设置会对计算结果产生显著影响。以高斯脉冲和正弦波激励源为例，在相同的有耗媒质和计算区域条件下，高斯脉冲激励下的散射场在时域上呈现出复杂的脉冲响应，包含多个散射回波，且由于不同频率成分的衰减不同，散射回波的波形会发生明显的畸变；而正弦波激励下的散射场在时域上则呈现出稳定的正弦振荡，散射回波的频率与激励源相同，只是幅度和相位会随着传播距离和散射体的作用而发生变化。在有耗媒质中放置一个金属散射体，当采用高斯脉冲激励时，散射场中会出现多个散射峰，分别对应不同的散射机制和散射路径，通过分析这些散射峰的时间延迟和幅度变化，可以获取散射体的形状、位置等信息；当采用正弦波激励时，散射场主要表现为正弦信号的幅度调制，通过分析调制的幅度和相位变化，可以得到散射体对正弦波的散射特性。因此，在实际应用中，需要根据具体的研究目的和需求，合理选择激励源的类型、位置和波形，以获得准确的计算结果。3.3.3时域阻抗边界条件在时域抛物线方程求解有耗媒质电磁散射与传播问题中，时域阻抗边界条件起着关键作用。时域阻抗边界条件是描述电磁波在有耗媒质边界上行为的重要依据，它反映了边界处电磁场的连续性和媒质的阻抗特性。首先推导时域阻抗边界条件的表达式。从麦克斯韦方程组在时域的形式出发，对于两种媒质的分界面，根据电磁场的连续性条件，电场强度的切向分量\vec{E}_{t}和磁场强度的切向分量\vec{H}_{t}在分界面两侧是连续的。设分界面的单位法向量为\vec{n}，从第二种媒质指向第一种媒质。引入时域表面阻抗Z_s(t)，它是一个与时间相关的复数，反映了有耗媒质边界对电磁波的阻碍作用随时间的变化。根据阻抗的定义和麦克斯韦方程组的边界条件，可以得到时域阻抗边界条件的表达式为：\vec{E}_{t}(t)=-Z_s(t)\vec{n}\times\vec{H}_{t}(t)在推导过程中，利用了电场和磁场的切向分量连续性条件，以及电场和磁场之间的相互关系。通过对麦克斯韦方程组中的安培环路定律和法拉第电磁感应定律在边界上的应用，结合表面阻抗的定义，经过一系列的矢量运算和数学推导，最终得到上述时域阻抗边界条件的表达式。在时域抛物线方程求解中，时域阻抗边界条件的应用步骤如下：首先，根据实际问题确定有耗媒质边界的位置和形状，以及表面阻抗Z_s(t)的具体表达式。对于不同的有耗媒质，其表面阻抗可以通过实验测量、理论计算或经验公式等方法获得。在研究土壤等有耗媒质时，可以通过测量土壤的电磁参数，利用相关的理论模型计算其表面阻抗。然后，在数值计算过程中，将时域阻抗边界条件应用到计算区域的边界节点上。在采用有限差分法或有限元法等数值方法求解时域抛物线方程时，根据边界节点的电场和磁场分量，利用时域阻抗边界条件建立边界节点的方程，从而将边界条件融入到整个数值计算体系中。通过应用时域阻抗边界条件，可以更准确地模拟有耗媒质边界对电磁波的散射和吸收作用。在处理有耗媒质与空气的分界面时，考虑时域阻抗边界条件后，能够更真实地反映电磁波在分界面上的反射和透射情况，得到更符合实际物理现象的计算结果。它可以有效地提高计算结果的准确性，为有耗媒质电磁散射与传播问题的研究提供更可靠的理论依据。3.3.4数值算例与分析为了深入分析有耗媒质中电磁散射的时域特性，并验证时域阻抗边界条件抛物线方程方法的正确性，进行了以下数值算例研究。考虑一个位于有耗媒质中的二维金属圆柱体散射体，圆柱体半径r=0.3m，有耗媒质的复介电常数\widetilde{\varepsilon}=(2-j1)\varepsilon_{0}，磁导率\mu=\mu_{0}，其中\varepsilon_{0}和\mu_{0}分别为真空的介电常数和磁导率。入射波为高斯脉冲，其表达式为：E_{inc}(t)=E_{0}e^{-\frac{(t-t_{0})^{2}}{\tau^{2}}}其中，E_{0}为脉冲幅度，取1V/m；t_{0}为脉冲峰值时刻，取5ns；\tau为脉冲宽度，取1ns。计算区域在x方向范围为[-1m,1m]，z方向范围为[0,3m]，采用均匀网格划分，x方向步长\Deltax=0.01m，z方向步长\Deltaz=0.01m，边界采用完美匹配层（PML）吸收边界条件，以模拟无限空间的情况。利用时域阻抗边界条件抛物线方程方法进行计算，得到不同时刻的电场强度分布，结果如图5所示（此处假设已绘制出相应的电场强度分布图，展示不同时刻t_1、t_2、t_3的电场强度分布）。从图中可以清晰地观察到高斯脉冲在有耗媒质中的传播和散射过程。在初始时刻，高斯脉冲从激励源位置出发，向有耗媒质中传播。随着传播距离的增加，由于有耗媒质的吸收作用，脉冲的幅度逐渐衰减。当脉冲遇到金属圆柱体散射体时，发生散射现象，散射场向各个方向传播。在散射场中，可以看到明显的绕射和干涉现象，这与理论分析和实际物理现象相符。进一步分析散射场的时域特性，绘制出接收点处的电场强度随时间的变化曲线，如图6所示（此处假设已绘制出相应的电场强度随时间变化曲线）。从曲线中可以看出，在入射脉冲到达接收点之前，电场强度为零。当入射脉冲到达接收点时，电场强度迅速上升，达到峰值后，由于有耗媒质的吸收和散射作用，电场强度逐渐衰减。在衰减过程中，可以观察到多个散射回波，这些散射回波是由于脉冲在金属圆柱体表面的多次反射和散射形成的。通过分析散射回波的时间延迟和幅度变化，可以获取散射体的相关信息，如散射体的位置、形状等。为了验证方法的正确性，将计算结果与其他成熟的数值方法（如时域有限差分法（FDTD））进行对比。对比结果表明，两种方法计算得到的电场强度分布和散射场的时域特性基本一致，误差在可接受的范围内，这充分验证了时域阻抗边界条件抛物线方程方法的正确性和有效性。通过上述数值算例的分析，时域阻抗边界条件抛物线方程方法能够准确地模拟有耗媒质中电磁散射的时域特性，为有耗媒质电磁散射与传播问题的研究提供了一种可靠的数值计算方法。四、有耗媒质中电波传播的抛物线方程方法4.1电波传播的抛物线方程方法原理在研究有耗媒质中电波传播时，抛物线方程方法是一种基于傍轴近似理论的有效数值方法。该方法通过对波动方程进行合理的近似和简化，将其转化为抛物线方程的形式，从而降低计算复杂度，提高计算效率，能够有效地处理电波在有耗媒质中的传播问题。从波动方程的基本原理出发，在无源的有耗媒质中，电场强度\vec{E}满足的波动方程为：\nabla^{2}\vec{E}-\mu\widetilde{\varepsilon}\frac{\partial^{2}\vec{E}}{\partialt^{2}}=0其中，\mu为磁导率，\widetilde{\varepsilon}为有耗媒质的复介电常数，\widetilde{\varepsilon}=\varepsilon-j\frac{\sigma}{\omega}，\varepsilon为介电常数，\sigma为电导率，\omega为角频率。为了推导抛物线方程，通常采用傍轴近似。假设电波主要沿z方向传播，电场强度\vec{E}可表示为：\vec{E}(\vec{r},t)=\vec{E}_{0}(\vec{r}_{\perp},z)e^{j(k_{0}z-\omegat)}+c.c.其中，\vec{r}=(x,y,z)，\vec{r}_{\perp}=(x,y)，k_{0}=\omega\sqrt{\mu\varepsilon_{0}}是自由空间波数，\varepsilon_{0}为真空介电常数，c.c.表示复共轭。将上述电场强度表达式代入波动方程，并对z方向的导数进行近似处理。在傍轴近似条件下，假设\frac{\partial^{2}\vec{E}_{0}}{\partialz^{2}}\llk_{0}\frac{\partial\vec{E}_{0}}{\partialz}，经过一系列的数学推导（包括利用矢量运算规则、三角函数关系等），可以得到抛物线方程。以电场强度\vec{E}为例，其抛物线方程形式如下：2jk_{0}\frac{\partial\vec{E}_{0}}{\partialz}+\nabla_{\perp}^{2}\vec{E}_{0}+k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\vec{E}_{0}=0其中，\nabla_{\perp}^{2}=\frac{\partial^{2}}{\partialx^{2}}+\frac{\partial^{2}}{\partialy^{2}}是横向拉普拉斯算子。从物理意义上理解，抛物线方程描述了电波在有耗媒质中传播时，电场强度在横向和传播方向上的变化关系。2jk_{0}\frac{\partial\vec{E}_{0}}{\partialz}项表示电波沿z方向传播时的相位变化和幅度变化，\nabla_{\perp}^{2}\vec{E}_{0}项反映了电场强度在横向方向上的变化，k_{0}^{2}(\widetilde{\varepsilon}\mu-1)\vec{E}_{0}项则体现了有耗媒质对电波的影响，包括吸收和散射等作用。在实际应用中，抛物线方程方法对传播环境有一定的简化和假设。它假设电波的传播方向主要集中在某一特定方向（如z方向），即满足傍轴近似条件。这意味着电波的波前在横向方向上的变化相对缓慢，波的传播方向与参考方向的夹角较小。在这种假设下，抛物线方程方法能够有效地处理许多实际问题，如电波在大气、海洋等有耗媒质中的传播，以及在复杂地形环境下的传播等。然而，当传播环境不满足傍轴近似条件时，如存在大角度散射、强各向异性媒质等情况，抛物线方程方法的计算精度可能会受到影响，需要采用更精确的数值方法或对抛物线方程进行改进。4.2大气吸收的抛物线方程模型4.2.1氧气吸收复折射率模型在大气环境中，氧气对电波的吸收是影响电波传播特性的重要因素之一。氧气分子具有特定的能级结构，当电波的频率与氧气分子的能级跃迁频率相匹配时，会发生共振吸收现象，导致电波能量的损耗。从分子物理学的角度来看，氧气分子的吸收主要源于其转动和振动能级的跃迁。氧气分子是双原子分子，其转动能级的跃迁会吸收特定频率的电波。根据量子力学理论，分子的转动能级可以表示为E_{J}=hBJ(J+1)，其中h是普朗克常数，B是转动常数，J是转动量子数。当电波的能量h\nu（\nu为频率）与氧气分子的转动能级差\DeltaE_{J}相等时，即h\nu=\DeltaE_{J}，会发生转动能级跃迁吸收。为了准确描述氧气对电波的吸收特性，建立复折射率模型。复折射率n可以表示为n=n_{r}-jn_{i}，其中n_{r}是实部，代表媒质对电波的折射能力；n_{i}是虚部，代表媒质对电波的吸收能力。对于氧气吸收，其复折射率模型可以通过经验公式或理论推导得到。常用的Liebe吸收模式将氧气的吸收谱分为线状吸收和非谐振连续吸收。线状吸收由44条离散的吸收线组成，这些吸收线对应着氧气分子特定的转动能级跃迁。非谐振连续吸收则是由于氧气分子之间的相互作用等因素引起的，其吸收特性相对较为连续。在Liebe吸收模式中，氧气吸收的复折射率虚部n_{i}可以表示为：n_{i}=\sum_{k=1}^{44}\frac{S_{k}\Delta\nu_{k}}{\left(\nu-\nu_{k}\right)^{2}+\left(\frac{\Delta\nu_{k}}{2}\right)^{2}}+n_{i_{cont}}其中，S_{k}是第k条吸收线的强度，\Delta\nu_{k}是第k条吸收线的半宽度，\nu_{k}是第k条吸收线的中心频率，\nu是电波的频率，n_{i_{cont}}是非谐振连续吸收对应的复折射率虚部。从物理意义上理解，上述公式中的第一项表示线状吸收的贡献，每一项\frac{S_{k}\Delta\nu_{k}}{\left(\nu-\nu_{k}\right)^{2}+\left