★ FEATURED ARTICLE 简介这份文档围绕短波航空移动信道的建模与仿真展开以Watterson模型为基础结合飞行器相对运动带来的多普勒效应提出了适用于远距离航空移动通信的短波信道模型并针对不同参数飞行器及典型民航场景如已知航线航迹进行了差异化与定制化仿真。内容系统介绍了Watterson模型的抽头延迟线结构、时变频响表达式、抽头增益函数及其高斯谱假设适合通信工程、航空通信领域的研究人员、高校师生以及从事短波信道仿真与系统设计的工程师参考学习。资源为单篇DOCX文档共1个文件压缩包仅526KB含完整论文式正文目前已有159人学习使用。文档重点贡献包括信道模型构建、多普勒频移与扩展的数学推导、不同飞行条件下的仿真实现思路对理解短波航空信道特性、开展相关课题研究或撰写技术报告均有直接借鉴价值。1. 短波航空移动信道的建模为什么不能直接拿 Watterson 模型用短波通信工作在 330 MHz靠电离层反射实现数千公里超视距传输在卫星链路失效、极地飞行这类场景下短波几乎是地空远程指挥控制的最后手段。现在做短波信道仿真的人大多数逃不过 Watterson 模型——它是国际电信联盟推荐的短波性能测试标准信道模型抽头延迟线结构多径、衰落、多普勒频移和扩展都覆盖到了。但问题在于Watterson 模型的成立前提是“有限时间、有限带宽下的静态信道”飞行器高速运动产生的时变多普勒效应完全不在原模型假设范围内。直接把静态模型搬到航空场景仿出来的衰落特性是失真的。这份资源就是基于 Watterson 模型把相对运动的多普勒频移拆解、建模并实现仿真的完整过程包含三种机动频率场景和固定航线定制化仿真的参数与结论适合做航空通信链路仿真、信道模拟器或短波通信系统验证的从业者直接参考。2. 先吃透 Watterson 模型抽头延迟线、复高斯增益与高斯谱叠加Watterson 模型的本质是一个抽头延迟线结构输入信号经过延迟线后产生多条路径每条路径对应电离层的一种传播模式。比如 E 层反射和 F 层反射同时到达接收机就构成了两个不同的抽头。每个抽头后面挂一个时变增益函数 $G_i(t)$用来调制这路信号在传播过程中的衰落最后所有路径的输出加上加性噪声组成信道输出。这个结构直观、实现复杂度低是它成为标准模型的核心原因。从频响角度看Watterson 模型的时变频响可以表达为$$H(f,t)\sum_{i1}^{n} e^{-j2\pi \tau_i f} G_i(t)$$其中 $i$ 是路径标号$n$ 是路径总数$\tau_i$ 是第 $i$ 条路径的时延$G_i(t)$ 是第 $i$ 条路径的抽头增益函数。公式本身不复杂真正影响仿真结果的是 $G_i(t)$ 怎么生成以及各个参数怎么定。2.1 结构输入信号先过抽头延迟线实际仿真的常见做法是先准备一组复高斯序列作为各径增益的基底然后用频移和滤波把多普勒效应加进去再把输入信号分别按 $\tau_i$ 做延时后与对应增益相乘最后把各径结果累加。抽头数不用贪多短波天波传播在大多数工程场景下 23 条路径就够用了路径太多除了拉升计算量对物理模型的解释力并没有实质帮助。需要注意的细节是Watterson 模型描述的是窄带信道带宽上限大约 12 kHz。做 330 MHz 的短波信道仿真时输入信号通常是经过数字下变频后的基带信号带宽在这个范围内才符合模型的适用条件。超过这个带宽电离层色散效应就不能忽略需要换 ITS 模型这类更复杂的统计模型。2.2 抽头增益复高斯过程才是核心$G_i(t)$ 的数学形式是两组复高斯随机过程的加权叠加$$G_i(t)G_{ia}(t)e^{j2\pi f_{ia}t}G_{ib}(t)e^{j2\pi f_{ib}t}$$$G_{ia}$ 和 $G_{ib}$ 是相互独立的零均值复高斯随机过程包络服从瑞利分布相位服从均匀分布。这里的物理含义是每条路径由大量不可分离的电磁波射线组成叠加后按中心极限定理趋于高斯统计特性所以用复高斯过程是合理的建模选择。两个高斯分量分别对应电离层磁离子分裂产生的两个传播模式。在大多数仿真场景里这两个分量的多普勒频移和扩展数值接近工程上经常简化成一个分量。原文档在仿真部分就是这么处理的少一个分量能省将近一半的计算量同时不太影响结果的可信度。2.3 自相关与功率谱里的两个高斯分量抽头增益的自相关函数为$$C_i(\Delta t)C_{ia}(0)e^{-2\pi^2\sigma_{ia}^2(\Delta t)^2j2\pi v_{ia}\Delta t}C_{ib}(0)e^{-2\pi^2\sigma_{ib}^2(\Delta t)^2j2\pi v_{ib}\Delta t}$$对应的功率谱是两个高斯函数的叠加中心位置由多普勒频移 $f_{ia}$、$f_{ib}$ 决定宽度由多普勒扩展 $\sigma_{ia}$、$\sigma_{ib}$ 决定。做仿真时一定要把这两个量分开设置多普勒频移决定频谱搬移多少多普勒扩展决定频谱展宽多少两者不是一个概念。2.4 参数从哪里来ITU 参考值与实际配置短波信道仿真参数一般参考国际电信联盟给出的典型电离层反射信道参数。原文档里给出的多普勒扩展配置为 0.5 Hz、1 Hz、1.5 Hz 三档对应三条路径由低到高的扩展程度多普勒频移则需要结合飞行器运动速度重新计算。这些数值落在短波信道仿真的常用范围内用来做链路级性能验证是够用的。仿真时还有个容易被忽略的点Watterson 模型的平稳性假设只在有限时间10 分钟内成立。单次仿真时长不能太长超过这个窗口电离层参数就可能发生明显变化抽头增益的统计特性就不再稳定。我一般把单次仿真控制在几分钟量级需要长时间统计时分成多个独立快照来做。3. 给模型装上“运动部件”多普勒频移拆解与飞行器速度模型Watterson 模型本身是静态的没有考虑收发两端的相对运动。但航空移动信道里飞行器始终在动时变多普勒效应是必须建模的核心部分。原文的思路很直接把每条路径的多普勒频移拆成两项电离层本身产生的频移加上相对运动产生的频移再结合飞行器的速度和加速度限制建立完整的时变频移模型。3.1 多普勒频移的构成两项叠加短波航空移动信道中第 $i$ 条路径的总多普勒频移为$$f_i f_{iA} f_{iB}$$$f_{iB}$ 是电离层快速运动、反射层高度变化引起的频移这部分可以直接沿用 Watterson 模型里的设定。$f_{iA}$ 是飞行器相对运动产生的频移计算公式为$$f_{iA} \frac{f_c}{c} v \cos \theta_i$$其中 $f_c$ 是载波频率$c$ 是光速$v$ 是飞行器运动速率$\theta_i$ 是接收端入射电波与飞行器运动方向的夹角。不同路径的 $\theta_i$ 通常不同工程上一般按均匀分布随机取值。这里有个工程细节$f_{iA}$ 是随时间变化的因为飞行器速度 $v$ 在变$\theta_i$ 也会随着飞行器姿态变化而变所以整套仿真里 $f_i$ 不能设成常量。3.2 飞行器速度模型匀加速、圆周运动与上限约束飞行器的瞬时速度可以统一写成$$v(t) v_0 a(t)\cos\alpha(t)\cdot t$$$v_0$ 是初速度$a(t)$ 是加速度$\alpha(t)$ 是速度方向与加速度方向的夹角。$a(t)0$ 时是匀速运动$\alpha(t)0$ 时是直线运动其余情况是圆周运动。这样一个公式就把直线、匀加速、匀速圆周都覆盖了。实际飞行器受发动机功率限制不能无限加速。发动机功率 $P Fv (kv^2 ma)v$当飞行速度接近最大值时功率保持恒定加速度会逐渐降下来最终匀速飞行。原文给出的处理办法是当 $|v/v_{max}| \geq 0.8$ 时令加速度按 $a k_0(1/x - x^2)$ 衰减其中 $x |v/v_{max}|$$k_0$ 是归一化系数。这样速度曲线就能平滑收敛到最大速度不会出现越跑越快的非物理情况。这个约束对于高机动无人载具场景尤其重要不加限幅的仿真结果会明显偏离真实飞行能力。3.3 机动频率的分级参考表飞行器轨迹不确定时多普勒频移的快变程度和机动频率直接相关。原文档给了三个典型场景的参考值我把这个表整理成更直观的形式机动频率等级典型应用场景机动频率参考值 (Hz)飞行状态持续时间 (s)低民用航空0.01≥200中私人飞机等0.11020高无人载具等1≤1低机动频率对应的飞行状态持续时间长速度曲线平稳高机动频率下飞行器速度大小和方向在短时间内快速变化多普勒频移剧烈波动。仿真时可以根据目标场景选对应的档位随机生成一段 100 秒的速度序列作为信道输入。3.4 模型汇总把时变频移代回抽头增益把 $f_i \frac{f_c}{c}[v_0 a_1(t)\cdot t]\cos\theta_i f_{iB}$ 代回 Watterson 模型的抽头增益表达式就得到了短波航空移动信道的完整模型。每条路径的多普勒频移是时变的抽头增益 $G_i(t)$ 也跟着时变信道不再平稳但这恰恰是航空移动信道需要的特性。实现时最直接的做法是用一个复指数序列 $e^{j2\pi\int f_i(t)dt}$ 去调制复高斯序列得到带时变频移的抽头增益。频移值按每个采样点重新计算而不是每段取个平均值这样才能反映出加速度带来的频谱展宽。第 4 章会把这个过程落到具体的仿真代码上。4. 动手实现三径短波航空移动信道的仿真搭建模型有了接下来就是把它变成能跑的仿真。原文档给出的实现框图里信道模型的输入端是经过数字下变频得到的 I/Q 基带信号输出端是经过多径、衰落、多普勒频移和多普勒扩展作用后的基带信号。核心模块是希尔伯特滤波器、时变多普勒频移发生器、多普勒扩展成形滤波器下面按顺序拆开讲。4.1 基带处理希尔伯特滤波器生成 I/Q 两路实际短波通信的载波在 330 MHz直接在这个频段做仿真采样率太高、计算量太大。常见做法是先把信号下变频到基带再用希尔伯特滤波器构造 I/Q 两路带通信号。原文档给出的方法是先设计一个 FIR 低通滤波器通带设为期望带通滤波器通带的一半然后通过余弦、正弦调制得到 I/Q 两路带通系数$$h_{IBP}(n) 2h_{LP}(n)\cos\left(2\pi f_0[n-\frac{N-1}{2}]T\right)$$$$h_{QBP}(n) 2h_{LP}(n)\sin\left(2\pi f_0[n-\frac{N-1}{2}]T\right)$$$h_{LP}$ 是 n 阶 FIR 低通系数$f_0$ 是通带中心频率$N$ 是滤波器阶数$T$ 是采样周期。MATLAB 里这么实现% 参数系统采样率 10 kHz基带中心频率 200 Hz fs 10000; Nf 64; % 滤波器阶数64阶足够平坦 f0 200; % 基带信号中心频率 hLP fir1(Nf, 100/(fs/2)); % 原型低通通带100 Hz nIdx 0:Nf; hIBP 2*hLP .* cos(2*pi*f0*(nIdx - Nf/2)/fs); hQBP 2*hLP .* sin(2*pi*f0*(nIdx - Nf/2)/fs);这段代码先设计通带 100 Hz 的 FIR 低通原型再通过频率搬移得到 I/Q 两路带通滤波器。滤波器阶数 64 在 10 kHz 采样率下幅频特性已经足够平滑阶数太高反而增加卷积计算量。$f_0 200$ Hz 对应原文档仿真里下变频后的单音信号中心频率如果你想仿 400 Hz 中心频率把 $f_0$ 改成 400 即可两个带通滤波器的中心频率会跟着移动。4.2 时变多普勒频移时域乘复指数一个复信号频移 $\Delta f$等价于在时域乘以 $e^{j2\pi\Delta f t}$。如果输入信号是 $e^{j2\pi f t}$频移后变成 $e^{j2\pi(f\Delta f) t}$频率就完成了叠加。当 $\Delta f$ 随时间变化时相位是频率的积分不能简单用乘法一次搞定。下面这段代码实现了带速度上限约束的时变多普勒频移生成% 输入t 为时间轴长度 N % 输出f_i 为每条路径的时变多普勒频移Git 为带频移的抽头增益 fc 15e6; % 载波频率 15 MHz c 3e8; v0 250; % 初速度 m/s a_max 80; % 最大加速度 m/s^2 v_max 600; % 最大速度 m/s theta deg2rad(30); % 入射角与运动方向夹角 f_iB 1; % 电离层本底频移 Hz % 速度序列生成接近最大速度时加速度衰减 v zeros(size(t)); v(1) v0; for k 2:N if abs(v(k-1)/v_max) 0.8 a a_max; else x abs(v(k-1)/v_max); a a_max * (1/x - x^2) / (1/0.8 - 0.8^2); end v(k) v(k-1) a/fs; if v(k) v_max, v(k) v_max; end end % 时变多普勒频移 f_iA fc/c * v .* cos(theta); f_i f_iA f_iB; % 抽头增益复高斯基底乘时变频移复指数 Gia (randn(1,N) 1j*randn(1,N))/sqrt(2); Git Gia .* exp(1j*2*pi*cumsum(f_i)/fs);关键逻辑在最后一行频率对时间做积分再乘到复高斯序列上。直接用exp(1j*2*pi*f_i.*t)的写法在频移恒定时没问题但频移随时间变化时相位会算错频谱会多出虚假的分量。cumsum(f_i)/fs本质是矩形法数值积分采样率 10 kHz 下精度足够。速度循环里对加速度做了上限限幅处理确保飞行器速度不会越过最大速度这对应原文里发动机功率受限的物理约束。4.3 多普勒扩展高斯白噪声成形与插值多普勒扩展在 Watterson 模型里假设为高斯功率谱。数字实现方法是让高斯白噪声通过一个高斯型滤波器得到功率谱呈高斯形状的噪声序列再与输入信号乘起来。时域乘法对应频域卷积等效于把信号频谱展宽。这里有个工程坑如果直接在系统采样率 10 kHz 下设计高斯滤波器多普勒扩展只有 0.5 Hz归一化截止频率只有 0.00005滤波器阶数会非常大计算量暴涨。原文档的处理办法是先低采样率生成再插值回系统采样率% 低采样率生成高斯谱噪声 fs_low 50; % 低采样率 50 Hz N_low round(T * fs_low); wn randn(1, N_low); % 高斯白噪声 sigma_i 1; % 多普勒扩展 1 Hz lpFilt designfilt(lowpassfir, FilterOrder, 20, ... CutoffFrequency, sigma_i/(fs_low/2)); noise_low filter(lpFilt, 1, wn); noise_low noise_low / std(noise_low); % 归一化到单位功率 % 线性插值回系统采样率 noise interp1(linspace(0, T, N_low), noise_low, t, linear); % 与抽头增益相乘得到扩展后的增益 G_ext Git .* noise;20 阶 FIR 滤波器在 50 Hz 采样率下就能近似出 1 Hz 带宽的高斯谱计算量几乎可以忽略。插值后的噪声序列统计特性基本保持不变多普勒扩展效果与直接在高采样率下滤波一致但计算量差了几个数量级。sigma_i对应各条路径的多普勒扩展值原文档里三条路径设为 0.5 Hz、1 Hz、1.5 Hz分别模拟不同电离层传播模式的抖动程度。归一化那行很重要滤波会改变噪声功率不归一化的话信号幅度会被整体缩放后续信噪比计算就不准了。4.4 仿真参数配置载波、时延、信噪比的匹配原文档的仿真配置是三径结构载波频率分别取了 6 MHz 和 15 MHz 两档。下面是完整的参数配置表可以直接作为复现基线仿真参数设置值说明载波频率6 MHz / 15 MHz下变频后在基带处理基带信号200 Hz 单音便于直观观察频移和扩展路径数3时延分别为 0、1 ms、2 ms多普勒扩展0.5 / 1 / 1.5 Hz对应三条路径信噪比10 dB加性高斯白噪声最大加速度80 m/s²机动上限最大速度600 m/s速度限幅确定性航迹场景下原文档给了 15 秒的运动参数分段表可以组合出匀加速直线、匀速圆周、匀速直线三种运动状态用于对特定场景做定制化仿真时间段运动方式初速度加速度05 s匀加速直线050 m/s²510 s匀速圆周 (R1000 m)250 m/s向心加速度 V²/R1015 s匀速直线250 m/s0仿真结果里有一个现象值得注意匀加速直线运动下信号频谱扩展比匀速直线更大因为加速度导致多普勒频移在积分区间内持续变化匀速圆周运动则因为运动方向连续改变会同时出现正、负频偏。这两点可以直接用来验证自己的实现是否正确。5. 常见问题与避坑记录滤波器阶数、时变频移与基带换算这套模型从理论到仿真落地有几处地方特别容易翻车。我按踩过的顺序整理成五条避坑记录每一条都是实际能复现出来的问题。5.1 高斯滤波器阶数爆炸导致计算量失控现象直接把多普勒扩展 0.5 Hz 的高斯滤波器设计在系统采样率 10 kHz 上滤波器阶数轻松上千运行时明显卡顿仿真 10 秒信号要等几十秒。原因滤波器截止频率与采样率的比值太小0.5 Hz 对应归一化截止频率 0.00005FIR 滤波器要达到足够的阻带衰减阶数必然飙升。解决采用原文的思路先在低采样率如 50 Hz下设计滤波器并完成滤波再用插值把噪声序列恢复到系统采样率。处理后的信道特性不变计算量从 O(N²) 量级降到接近 O(N)这是最值得优先落地的优化。5.2 把时变多普勒频移做成恒定值现象仿真出来的频谱只有一条干净的谱线搬移看不到理论上的频谱展宽尤其在高机动场景下结果明显偏离预期。原因代码里把 $f_i$ 当常量处理直接乘 $e^{j2\pi f_i t}$加速度带来的频移变化被完全忽略。解决用cumsum(f_i)/fs对时变频移做积分后再生成复指数让相位随频率的积分连续变化。匀加速场景下频谱会出现明显的展宽效果这正是航空移动信道区别于静态 Watterson 模型的关键特征。5.3 仿真时长超出 Watterson 模型的有效窗口现象仿真跑 30 分钟甚至更长时间信道统计特性出现漂移某些时刻的冲激响应明显异常。原因Watterson 模型的平稳性假设只在有限时间10 分钟内成立超过这个时长远超电离层的准静态变化尺度。解决单次仿真严格控制在 10 分钟以内需要跑长时间统计时分段生成信道快照每段重新生成电离层参数。原文档里 100 秒速度仿真和 5 秒频谱对比都落在这个窗口内这是有意的设计。5.4 基带仿真里的载波参数换算疏漏现象直接把 15 MHz 载波频率代进公式算出的多普勒频移是 15 MHz 下的真实物理频移但基带仿真里信号中心频率只有 200 Hz频移后信号可能超出基带带宽滤波后整个信号失真。原因多普勒频移公式 $f_{iA} f_c/c \cdot v\cos\theta$ 用的是真实载波频率这个物理频移需要对应到基带上。数字基带仿真里信号本身就代表下变频后的频谱频移量直接用物理值叠加即可不需要再乘基带与载波的频率比。解决先按真实载波算出物理多普勒频移再在基带复指数乘法里用这个物理值不要额外做一次比例换算。6 MHz 和 15 MHz 两套载波仿真结果不同反映的就是真实频移量的差异。5.5 多普勒扩展与多普勒频移混用现象输出信号的频谱整体搬移了但频谱形状没有展宽高机动频率场景下信号依然“很干净”。原因只做了频移乘法复指数漏掉了多普勒扩展高斯噪声成形后相乘。扩展和频移在 Watterson 模型里是两个独立参数分别控制频谱的位置和形状。解决检查流程里是否同时包含两个模块——频移用复指数乘法扩展用成形噪声乘法。两者缺一不可顺序上先做频移再做扩展扩展在频域上表现为卷积不会覆盖掉频移效果。6. 六个验证动作冲激响应、频谱对比与确定性运动场景仿真代码写完先别急着往项目里接按下面这套动作验证一遍确认模型行为正确再做集成。第一步看冲激响应。三径配置下冲激响应应该出现三个明显的峰值区域峰值位置对应 0、1 ms、2 ms 时延。同一机动频率下不同时刻的信号衰落幅度不一致说明信道是时变的不同载波频率6 MHz 对比 15 MHz下衰落特性不同说明信道同时具有频率选择性这两个特征是短波航空信道的基本行为。第二步看频谱搬移和展宽程度。15 MHz 载波下跑三个机动频率档位低机动场景的频谱应该只有轻微搬移高机动场景应该出现明显的频谱展宽甚至信号接近淹没在噪声里。如果高机动场景和低机动场景的频谱看不出差别基本可以断定时变多普勒频移没做进去。第三步用确定性运动场景做交叉验证。按表 2 的 15 秒运动参数跑一遍单路径频谱匀加速直线段应该比匀速直线段的频谱更宽匀速圆周段应该同时出现正、负频偏。这三段结果如果都对得上说明时变频移的积分和速度限幅逻辑没有问题。这三个动作跑完模型基本可以放心用。这套仿真代码里最值得反复检查的就是cumsum(f_i)/fs那一行。从那以后我每次搭短波信道模型都会先跑一遍匀加速对比匀速的频谱再跑一遍三径冲激响应确认时变频移和扩展都生效了才敢把结果往链路仿真里接。希望这些验证习惯能帮你少走几步弯路祝你顺利。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
