简介：本讲义以大气波导环境为背景，系统梳理电波传播的关键理论与建模方法，适合通信工程、电磁场与微波技术方向的学生及从事超视距通信、雷达探测相关工作的技术人员学习。内容从大气折射的四种基本类型切入，逐步讲解折射率与修正折射率的计算、陷获折射与大气波导的形成条件，并重点引入抛物方程模型（PE法），结合分步傅里叶变换、粗糙海面阻抗边界条件及初始场设置，说明电波传播损耗的仿真思路。讲义还展示了MATLAB在电波传播损耗计算中的具体应用，帮助读者理解复杂大气结构下信号的传播路径与衰减规律。资源共1个PPT文件，压缩包大小472KB，整体结构紧凑，按背景介绍、大气折射、抛物方程模型、MATLAB仿真、大气波导影响与结论的顺序递进展开，便于对照公式与图示逐步理解。已有486人学习，适合希望掌握抛物方程方法、提升电波传播仿真能力的读者参考。

1. 收到这份“PPT学习教案”，先想清楚它到底在帮你解决什么问题

一台工作在 3GHz 频段的岸基雷达，标准大气条件下探测距离只有四五十公里，可某个早晨它在 120 公里外稳定地抓到了目标回波。懂行的人第一反应不是“雷达变强了”，而是“大气波导出现了”。这种超视距现象，折射率剖面里藏着看不见的导管，把电波像光纤一样夹在特定高度层里反复折射。要量化这种传播损耗、给出可预测的覆盖范围，传统的射线追踪在焦散区和绕射区集体翻车。基于抛物方程的大气波导环境下电波传播的研究，就是奔着这个黑匣子去的：它把 Helmholtz 波动方程在传播方向上做抛物近似，用步进式算法把折射率剖面、地表边界、天线方向图一层层叠进计算，最终输出一张“在某种波导条件下，电波到底往哪拐、损耗多大”的传播损耗图。

这份教案面向的读者不是来长见识的路人，而是三种人：电磁场专业的学生想把“绕射”“折射”“多径干涉”从公式变成直观图像；雷达与通信链路设计者要评估波导环境下的实际覆盖盲区；科研人员则需要一个能解释夜间超视距、海面异常传播这些怪异现象的计算工具。它的价值不在于推导过程有多严谨华丽，而在于把一条“从气象剖面到传播损耗”的计算链路讲清楚，并让学习者能自己动手复现一次。后面这些章节，我就按这条链路一层层拆给你看。

2. 大气波导环境：折射率模型与射线追踪失效的边界在哪

2.1 波导类型与折射率剖面：表面波导、抬升波导、蒸发波导怎么区分

大气波导本质上是折射率梯度异常造成的“陷获层”：当修正折射率 M 随高度下降时，电波路径向地面弯折的程度超过地球曲率，射线就被禁锢在某个高度带内反复折射，形成导管。归一化到地球曲率后，标准大气的修正折射率梯度约为 +118 M/km（教科书常取 120 M/km），只要是负梯度，就具备了波导形成的必要条件。区别在于负梯度层出现的海拔高度和厚度，这决定了它属于哪一种波导类型。

表面波导的陷获层贴在地表，M 剖面从地面开始先急剧减小再转回标准梯度，天线稍微抬升一点就可能脱离陷获区，这也是它往往只在低空起作用的原因。抬升波导则是在某个高度之上出现强负梯度层，下方是标准大气、上方也是标准大气，像三明治一样把能量夹在中间；这种波导对高度恰好处于陷获层内的目标影响最大，也是超视距传播中最难用射线法描述的形态。蒸发波导则是最特殊也最日常的一种：它只出现在海面上方几十米内，由海气界面湿度梯度形成，几乎一年四季都存在，变化规律由海面温度、风速、湿度共同决定。雷达和通信系统的工程评估里，蒸发波导的优先级往往最高，因为它稳定出现、寿命长、覆盖范围动辄上百公里。

这三类波导的 M 剖面结构可以直接画成三条曲线来对比：表面波导是“地面低凹槽”，抬升波导是“高空负梯度夹层”，蒸发波导是“海面薄层陡降”。教案里如果能配一张三剖面对照图，读者对陷获条件的理解会远好于背公式。

2.2 射线追踪在波导场景为什么会翻车：焦散、绕射与能量泄漏

很多从几何光学入门的人，第一反应是用射线追踪去算波导传播：把折射率剖面离散成一层层，让射线按 Snell 定律逐步弯折，最后统计到达某点的射线数量。这个思路在描述“波导形成”时很直观，但在定量计算传播损耗时问题非常大。首先是焦散区：弯折的射线会在空间某些位置重叠聚焦，那里场强按射线密度算会趋向无穷大，而实际电磁场是有限的；这是几何光学的奇点，也是射线法的天然缺陷。其次，波导不是无损耗的密闭管道，能量会通过绕射不断泄漏到波导外的区域，射线法对“泄漏”几乎没有建模能力，泄漏量恰恰是超视距传播评估里最核心的指标——它直接决定远区的信号还能不能解调出门限。

第三个缺陷更隐蔽：多径干涉。波导内同时存在直达、地面反射、波导壁多次反射等多条路径，它们在接收点叠加，有的频率加强、有的频率抵消，形成频率选择性衰落。射线追踪虽然能画出每条路径，但对“相位”和“幅度”的信息获取精度不够，几十条路径叠加出来往往与实测差得离谱。这些场景恰恰是抛物方程法的主场：它一步一个网格地递推场分布，天然包含绕射效应、自动描述干涉结构，不需要人为判断“这条射线要不要算”。这就是为什么业内评估超视距传播时，PE 几乎是绕不开的工具。

2.3 PE 在波导场景的优势边界：什么时候它也不是万能药

讲清楚 PE 的适用范围，比吹捧它能解决所有问题更重要。PE 对传播方向的近似让它无法处理后向散射，所以它擅长的是“前向传播”问题——雷达探测、点对点通信链路、无线覆盖评估都属于这类。但如果场景里存在强反射体、复杂的非均匀介质回波，或者需要计算双站散射，PE 就不合适了，那时需要回到全波方法如 FDTD 或矩量法。另外，PE 把大气假设为水平分层或缓慢变化，虽然可以通过“分段水平均匀 + 逐段更新”来处理非均匀路径，但每段的长度和离散策略需要人工权衡，处理不好会引入人为反射。

在实际工程里，我的判断标准很简单：目标在视距内、大气近似水平分层、只关心前向传播损耗，就优先考虑 PE；要做复杂地物绕射、多径时延精细统计、或后向散射时，再考虑混合方法或全波仿真。教案里如果能把这段“适用范围”讲明白，学习者就不会拿着 PE 到处套，也不会因为套错了而回头骂这个算法是垃圾。

3. 抛物方程法的核心脉络：从亥姆霍兹方程到分步傅里叶求解

3.1 抛物近似的物理含义：丢掉后向波，留下前向传播

抛物方程法不是从天而降的独立理论，它是从亥姆霍兹方程出发，在“传播方向占优”的前提下做近似得到的。假设电波主要沿 x 轴正向传播，把场写成缓变包络 u(x,z) 乘以快变相位 exp(ikx) 的形式，代入亥姆霍兹方程后，对传播方向的二阶导数项做近似处理，就得到窄角抛物方程：

[ \frac{\partial^2 u}{\partial z^2} + 2ik\frac{\partial u}{\partial x} + k^2(n^2-1)u = 0 ]

这个方程里，x 是传播距离，z 是高度，k 是波数，n 是折射率。丢掉的那一项是 (\partial^2 u / \partial x^2)，它对应的是后向散射和剧烈的前后向耦合。丢掉它的物理前提是：传播方向上的场变化远慢于横向变化，且不存在强烈的反向反射面。对大气波导这种缓慢变化的介质来说，这个前提基本成立；对急剧变化的城市环境或金属障碍物则完全不成立。

宽角抛物方程则是对那个被丢掉的算子做高阶展开，用 Padé 近似或 Feit-Fleck 展开保留更多项，使其支持更大的传播角度（可达 40° 甚至更大）。波导传播中，射线路径虽然会在垂直方向发生大幅度弯折，但弯折对应的角度通常不超过十几度到三十度，窄角 PE 在蒸发波导和浅表面波导场景下往往够用；抬升波导或强梯度场景则建议直接用宽角 PE 更稳妥。教案里建议把窄角和宽角的适用角度边界标清楚，这是初学者最容易忽略、也最容易导致结果翻车的地方。

3.2 分步傅里叶法（SSFT）的离散步骤：折射率相位屏和衍射传播交替进行

分步傅里叶法是求解抛物方程最经典的数值方案，核心思想是“每走一步，先让波前自由衍射传播一个步长，再叠加上折射率引起的相位修正”。之所以能这样拆，是因为在傅里叶谱域中，自由空间传播算子是对角化的，计算效率极高；折射率不均匀性则在空间域里用一个相位屏乘上去，两者的作用在局域范围内近似解耦。

具体到实现步骤可以用下面这段伪代码来描述，它表达的是算法骨架，不是某个可直接运行的程序包：

# 参数设置（教学演示用典型值，具体场景需要调整） freq = 3.0e9 # 工作频率，单位 Hz k = 2 * pi * freq / c # 波数 dx = 50.0 # 传播方向步长，单位 m，通常取 10~100 倍波长 dz = 0.05 # 垂直方向网格间距，单位 m，一般取波长的 1/4 到 1/2 zmax = 500.0 # 最大计算高度，单位 m，要高于波导陷获层和天线高度 nz = int(zmax / dz) # 垂直网格点数 # 初始化场：用高斯孔径近似天线方向图，放在指定天线高度 za 处 u = gaussian_aperture(nz, za, beam_width, freq) # 折射率剖面：从气象数据构造 M 剖面，再转换为折射率 n m_profile = build_m_profile(surface_waveguide_params) n_profile = m_to_n(m_profile) # SSFT 步进循环 for ix in range(max_range_steps): # 1. 自由空间衍射：在傅里叶谱域乘传播算子 U = fft(u) # 窄角: 相位项 exp(-i * kz^2 * dx / (2*k)) # 宽角: 相位项 exp(-i * dx * (sqrt(k^2 - kz^2) - k)) U = U * diffraction_propagator(k, kz_vector, dx, wide_angle=True) u = ifft(U) # 2. 折射率相位修正：乘折射率屏 u = u * exp(1j * k * (n_profile - 1) * dx) # 3. 吸收边界/窗函数：抑制顶部和底部的非物理反射 u = u * window_function(z, zmax, window_margin) # 4. 每隔若干步记录一条传播损耗剖面 if ix % report_step == 0: save_loss_profile(u, ix * dx, antenna_height)

这段伪代码里有几个参数值得单独说明。dx 是传播方向的步长，它决定计算效率和横向分辨率，经验上取 10 到 100 倍波长远小于场景变化尺度即可，但这只是工程经验，需要根据剖面变化的剧烈程度适当加密。dz 是垂直网格间距，理论上要满足奈奎斯特条件，通常取 0.2λ 到 0.5λ，对 3GHz 就是 0.02 到 0.05 米；这个值在很多教材里写得不直观，导致初学者把 dz 设成 1 米，然后算出来的干涉周期全是错的。zmax 要覆盖到波导上方足够远，一般取天线高度加波导厚度的 3 倍以上，否则吸收窗会截断真实场结构。

第 3 步的窗函数是 SSFT 里最容易忽视却最关键的环节。因为没有吸收边界，波场传播到计算区域顶部会被“硬边界”反射回来，和真实场叠加出假的干涉条纹。常见做法是在 zmax 上方留出 20% 到 30% 的区域作为吸收层，用 Tukey 窗或 Hanning 窗把场平滑地压到接近零。

3.3 地表边界怎么处理：从理想导电平面到阻抗边界

大地和海面对电波的影响在波导计算里不能被忽略，尤其是垂直极化波在海上传播时，海水的有限电导率会显著吸收掠射角很小的能量。最粗糙的处理是把地面当成理想导电平面，场在 z=0 处强制为零（水平极化）或导数为零（垂直极化），这在低频和干燥地表下误差不大，但在海面微波频段会高估反射系数，进而高估波导内的驻波效应。

更稳妥的做法是使用 Leontovich 阻抗边界条件，把地表阻抗与频率、极化、海水的介电常数和电导率关联起来，以复数的形式写入 PE 的边界项。这一步看起来只是数学上的一个小改动，实际对传播损耗的影响经常达到 3~5dB，在超视距链路预算里已经是决定性的误差来源。教案里如果时间有限，可以先用理想导电平面把主流程跑通，再加阻抗边界看差异；但最终交到学习者手中的版本，应该把阻抗边界作为默认配置来写。

3.4 初始场建模：天线方向图如何映射为抛物方程的起始场

PE 是初值问题，给定起始高度上的复数场分布，然后朝传播方向推进。所以初始场必须能真实表达天线在垂直面上的辐射特性，否则后续所有传播路径都会带上错误的“种子”。最常见的做法有两种：一是直接把天线垂直方向图采样后作为初始幅度分布，相位按等相位面设置；二是用高斯场分布近似主瓣，让它在传播中自然衍射展开。高斯近似的好处是解析形式简单、数值上不会出现剧烈阶跃，适合在教案中先跑通；但对多波束、赋形波束天线或低旁瓣要求高的场景，必须回到实测方向图。

工程上我一般会建议把天线方向图离散成与垂直网格一致的点数，然后按方向图幅度和相位直接赋给初始场。天线高度本身也决定了波导是否能把能量有效“装”进陷获层：天线在波导上方，只有部分波束进入波导；天线在波导内部，大部分能量被限制传播。错过这一步的理解，后面看到“天线升高反而信号变差”会觉得不可思议，实际是波导耦合效率在起主导作用。

4. 把教案变成可复现的仿真：参数配置、验证流程与结果判读

4.1 教案的结构组织：从现象对照到理论推导再到仿真复现的最佳授课顺序

拿到“基于抛物方程的大气波导环境下电波传播的研究”这份教案，授课顺序和章节结构直接决定学习者能不能跟上。我讲这门课的顺序几经调整，下面这个顺序是最稳妥的：先把超视距案例现场讲透，再引入修正折射率与波导分类，之后讲 PE 原理，紧接着进入数值实现，最后做验证与结果分析。

这个顺序和单纯按“数学推导在前、应用在后”的教科书路线最大的区别是：学习者先看到了“抛物方程能解什么”，带着问题去听推导就不会在符号里迷路。教案整体可以设计成六个环节，下表给出建议的内容划分和每节的教学目标：

这份教案若是做成 PPT，每一个环节控制在 15 到 20 分钟之间，加上串讲和练习，半天到一天的教学周期比较合适。流程跑不完不要紧，关键是把“从剖面到损耗输出”这条主线保住。

4.2 最小仿真配置：一张参数表让新手避免 80% 的初始翻车

给入门者一个能跑出合理结果的参数模板，比让他们自己摸索要高效得多。下面这张表是基于教学经验总结的典型配置，适用频段为微波雷达与通信系统：

这些参数不是拍脑袋的固定值，而是给初学者的“安全起点”。调参时有一条主线：先固定频率和剖面，扫描 dx 看结果是否稳定；再扫描 dz 看干涉结构是否收敛；最后调天线高度和波导厚度做物理敏感性分析。每一步只动一个变量，这是基本素质。

4.3 怎么判断 PE 算得对不对：双射线模型、标准大气、光滑地面三件套

新手最常问的问题是“我算出来的结果怎么知道对不对”。PE 程序跑出来的是一堆复数场数据，看着像模像样，但如果没有独立参照系，错误可能潜伏很久。我采用的三件套验证法，每一步看一个明确对照。

第一步，在无波导标准大气下，把 PE 计算结果与双射线模型做对比。双射线模型只考虑直射波和地面反射波，它在视距内的干涉结构和幅度分布是可以用手算验证的。如果 PE 与双射线模型在视距内前几十公里处趋势不一致，大概率是初始场或地形边界设置出了问题。第二步，把折射率剖面设为标准线性梯度（+118 M/km），计算一条传播损耗距离曲线，看它是否符合教科书上标准大气超视距衰减规律——损耗随距离单调增加且没有异常抬升。若出现长距离周期性的“假增强”，基本可以断定吸收窗或边界条件有 bug。第三步，换一个已知场景（比如均匀折射率+理想导电平面），这时可以用解析公式直接算出传播损耗的参考值，用 PE 重算一遍。

这三件套跑完，才能有底气说程序是可信的。之后再进入波导场景，才能区分“物理现象”和“数值假象”。我在实际项目中见过太多人一上来就套波导剖面，算出一堆漂亮图，结果是窗函数反射叠加了边界伪影，白白烧掉几天时间。

5. 抛物方程大气波导仿真常见问题排查：现象、原因与解决

5.1 损耗曲线在远端出现周期性抖动：吸收窗边界反射在作怪

现象：传播损耗随距离增大时，本应单调平滑或缓慢起伏的曲线在远端出现规则的波浪状抖动，周期随距离缩短而加密，幅度可达 3~6dB。原因：顶部吸收窗没有把场平滑压到零，或窗口区太短，波场在计算区域顶部被硬边界反射回来，与正向传播场干涉叠加。解决：把吸收窗宽度加长到总高度的 20% 以上，改用 Tukey 窗（过渡区余弦平方渐变），并在测试时把 zmax 增大 1.5 倍看曲线变化——如果抖动消失，说明窗口余量不足。

5.2 波导内的干涉结构与理论不符：垂直网格间距 dz 设得过大

现象：多径干涉的深度和间隔与双射线模型预测明显不符，干涉条纹间距看起来比理论值宽甚至完全抹平。原因：PE 输出的是网格点上的场值，而垂直方向的干涉条纹周期与波长、射线掠射角相关。dz 远大于半个波长时，场结构被空间采样抹平。解决：计算 dz 的上限公式为 λ/（2·sin(最大掠射角)），宽角传播场景下要适当加密；经验上先取 λ/4，算完看干涉结构是否与理论一致再逐步放宽。

5.3 损耗值整体偏高且伴随“假覆盖”：地表阻抗边界设成了理想导体

现象：远区损耗整体偏小，波导内覆盖范围看起来比实测乐观很多，尤其在掠射角极小的地方。原因：理想导体边界把海面反射系数模值设成了 1，实际海面在微波频段有限电导率会吸收部分能量，尤其是在低掠射角下垂直极化反射系数的幅度远小于 1。解决：换成 Leontovich 阻抗边界，把海水介电常数设为 80-60j 的数量级（具体随温度和盐度变化），并对比两种边界的差值；若差值超过 3dB，说明边界条件对结果有显著影响，必须保留阻抗边界。

5.4 抬升波导场景算出的超视距覆盖远偏弱：窄角 PE 角度限制不够

现象：用窄角 PE 计算强抬升波导或厚表面波导时，远区场强比宽角 PE 或实测低 5~10dB，且覆盖边缘呈现非物理的截断。原因：窄角抛物方程对传播方向的二阶导数做了简单抛弃，有效支持的最大射线角只有十几度。抬升波导使得射线在高度方向大幅弯折，局部传播角超过窄角限制，截断效应被放大。解决：改用宽角 PE 传播算子（Feit-Fleck 或 Padé 近似），同时把 dz 适当加密到 λ/4 以下，对比窄角与宽角在强波导场景下的结果差异，差异超过 1dB 时就以宽角结果为准。

5.5 不同机器上相同参数算出不同结果：FFT 归一化与随机相位种子没统一

现象：同一个教案里的 PE 程序在两台电脑上跑出不同的数值结果，幅度差几个百分点，干涉位置也有偏移。原因：不同版本的 FFT 库在归一化约定上不一致（有的除以 N，有的除以根号 N），导致每一步衍射传播的幅度因子不同；另外如果初始场含有随机相位（用于模拟粗糙面或随机介质），随机种子没固定也会带来差异。解决：在代码开头显式定义 FFT 归一化常数，并先用单步解析验证确认幅度不随步进衰减；所有随机过程固定 seed 并写入配置文件，保证结果可复现。

6. 进阶用法：把 PE 当作分析工具，做波导参数灵敏度扫描

跑通一个波导场景只是开始，真正有价值的工程用法是用 PE 做“波导参数灵敏度扫描”。具体做法是固定频率、天线高度和地表条件，然后系统性地改变修正折射率的陷获层强度与厚度，一次次重算传播损耗，把结果整理成一张“高度-距离-损耗”二维图或一组距离曲线。这样能回答几个实测中极其重要的问题：波导陷获层降低 10m 时，覆盖距离会缩多少？天线从 20m 升到 40m，是进入还是脱离陷获区？频段从 3GHz 换到 9GHz，波导维持超视距的能力是增还是减？

这类扫描的实用价值在于它能暴露“临界翻转”现象。很多工程师直觉以为天线越高越好，但 PE 扫描结果往往显示恰好相反：蒸发波导高度有限，天线升到波导顶之上后，大部分能量从波导上方漏走，远区信号反而断崖式下跌。又比如波导强度变大时，陷获能力增强，但能量在导管内反复折射导致损耗集中，某些距离上的覆盖空洞反而增大。这时候 PE 就不再只是“计算工具”，而是帮你在设计阶段避开盲区的决策辅助。我在实际项目里，就是把 PE 当作一种黑匣子式的“电磁仪表”在使用，先校准再扫描，用扫描结果反过来指导实验布站和天线选高。

还有一个我踩过多次的教训：灵敏度扫描时，每组参数跑完都要保存 M 剖面和损耗矩阵的原始数据，不要只存图。后期分析时你常常会想“刚才那组参数的具体损耗值是多少”，重新跑一遍几小时，有原始数据十几秒就能回答。把扫描脚本写成循环调用一个纯函数，输入是一组剖面参数，输出是损耗数组，这样既方便并行，也方便做批处理。最终你会养成一个习惯：任何 PE 结果，先画损耗-距离曲线，再看高度-距离热力图，最后对照 M 剖面确认几何关系，三步走完才算真正看懂了一组数据。希望这些经验能帮你在抛物方程这条路上少走几个来回。

本文还有配套的精品资源，点击获取