3. 水下声信道多径扩展引起的码间干扰误码分析
从实测信道冲激响应到抽头长度决策
多径扩展引起的码间干扰(ISI)本质上是同一符号经不同路径到达后,在时间轴上互相重叠。判断严重程度的关键不是“多径多不多”,而是时延扩展与符号周期的比值。本节用一组实测信道冲激响应(CIR)数据完成三件事:量化多径扩展、评估ISI程度、给出均衡器抽头长度。
一、实测CIR数据的组织方式
假设工作频段为 8–16 kHz,采样率 48 kHz,探测信号为线性调频(LFM),通过匹配滤波得到复基带冲激响应。下表为某次湖试实测的离散多径分量(以首达径为 0 μs 参考,幅度归一化到首径):
| 路径编号 | 相对时延 τ (μs) | 归一化幅度 | 相位 (rad) |
|---|---|---|---|
| P0(首达径) | 0 | 1.00 | 0.00 |
| P1 | 320 | 0.62 | 1.10 |
| P2 | 780 | 0.41 | -2.30 |
| P3 | 1450 | 0.28 | 0.75 |
| P4 | 2100 | 0.15 | -1.60 |
由此可算出两个核心量:均方根时延扩展 τrms 与最大时延扩展 τmax。按功率加权计算(功率取幅度平方):
总功率 P = 1.00² + 0.62² + 0.41² + 0.28² + 0.15²
= 1.0000 + 0.3844 + 0.1681 + 0.0784 + 0.0225 = 1.6534
平均时延 τ̄ = Σ(Pi·τi)/P
= (1.0000·0 + 0.3844·320 + 0.1681·780 + 0.0784·1450 + 0.0225·2100)/1.6534
= (0 + 123.0 + 131.1 + 113.7 + 47.3)/1.6534 ≈ 251.2 μs
τ_rms = sqrt( ΣPi(τi-τ̄)² / P )
= sqrt( (1.0000·251.2² + 0.3844·68.8² + 0.1681·528.8²
+ 0.0784·1198.8² + 0.0225·1848.8²)/1.6534 )
≈ sqrt( (63101 + 1820 + 47000 + 112700 + 76900)/1.6534 )
≈ sqrt( 301521/1.6534 ) ≈ sqrt(182365) ≈ 427 μs
τ_max = 2100 μs
注意 τrms 明显小于 τmax,说明能量集中在 P0–P2 附近,P3/P4 是弱尾径——这直接影响抽头长度的取舍。
二、时延扩展与符号周期的关系
设符号率 Rs = 2 kBaud,则符号周期 Ts = 500 μs。定义归一化时延扩展:
D = τ_rms / T_s = 427 / 500 ≈ 0.85
D_max = τ_max / T_s = 2100 / 500 = 4.2
工程判据(经验阈值,非绝对标准):
- D < 0.1:近似平坦,ISI 可忽略;
- 0.1 ≤ D < 0.5:轻度 ISI,线性均衡即可;
- 0.5 ≤ D < 1.0:中度 ISI,需判决反馈或较长线性均衡;
- D ≥ 1.0:强 ISI,考虑 OFDM/扩频或分数间隔均衡。
本例 D ≈ 0.85,属中度偏强 ISI;而 Dmax = 4.2 意味着最弱尾径会拖尾到第 4 个符号之后,若只按 τrms 设计抽头会漏掉这部分拖尾能量。
三、抽头长度怎么定
线性均衡器(如 LMS/RLS 的 FIR 横向滤波器)要覆盖整个拖尾,抽头跨度应不小于 τmax。若采用符号间隔(T 间隔)均衡,抽头数 N 满足:
N_T ≥ ceil(τ_max / T_s) + 1 = ceil(4.2) + 1 = 6
但 T 间隔均衡对非整数倍时延的采样相位敏感,工程上更常用分数间隔(T/2)均衡,此时:
N_(T/2) ≥ ceil(2·τ_max / T_s) + 1 = ceil(8.4) + 1 = 10
考虑 P3/P4 为弱尾径,若允许少量残余 ISI,可截断到 2·τrms 附近:
N_(T/2) ≈ ceil(2·2·427/500) + 1 = ceil(3.42) + 1 = 5
两种取值的取舍:10 抽头覆盖全部多径、残余 ISI 最小但收敛慢、计算量大;5 抽头只覆盖主能量、收敛快但 P3/P4 拖尾会留下约 -12 dB 量级的残余干扰。本例建议取 N = 8~10(T/2 间隔),兼顾覆盖与收敛速度。
四、验证与判读
抽头长度是否合适,不能只看公式,要用数据验证:
- 残余ISI估计:用选定抽头数对实测 CIR 做迫零或MMSE均衡,计算均衡后CIR在符号采样点之外的残余能量占比。若残余能量 < -15 dB,认为可接受。
- 误码率对比:在相同信噪比下,分别用 N=5、8、10 跑误码率曲线。若 N=8 与 N=10 的 BER 差距小于 0.5 dB,说明尾径贡献可忽略,取 N=8 更经济。
- 收敛观察:记录 LMS 误差曲线达到稳态所需的符号数,抽头越多收敛越慢,需与帧长匹配。
若验证发现 N=10 仍无法把残余 ISI 压到 -15 dB 以下,说明信道存在时变或非线性,应转向判决反馈均衡(DFE)或频域均衡,而不是继续加长线性抽头。