3. 水下声信道多径扩展引起的码间干扰误码分析

从实测信道冲激响应到抽头长度决策

多径扩展引起的码间干扰(ISI)本质上是同一符号经不同路径到达后,在时间轴上互相重叠。判断严重程度的关键不是“多径多不多”,而是时延扩展与符号周期的比值。本节用一组实测信道冲激响应(CIR)数据完成三件事:量化多径扩展、评估ISI程度、给出均衡器抽头长度。

一、实测CIR数据的组织方式

假设工作频段为 8–16 kHz,采样率 48 kHz,探测信号为线性调频(LFM),通过匹配滤波得到复基带冲激响应。下表为某次湖试实测的离散多径分量(以首达径为 0 μs 参考,幅度归一化到首径):

路径编号相对时延 τ (μs)归一化幅度相位 (rad)
P0(首达径)01.000.00
P13200.621.10
P27800.41-2.30
P314500.280.75
P421000.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 设计抽头会漏掉这部分拖尾能量。

多径时延扩展 vs 符号周期(Ts = 500 μs) 相对时延 (μs) 0 500 1000 1500 2000 2500 Ts 2Ts 3Ts 4Ts P0 P1 P2 P3 P4 τ_rms=427μs 线性均衡器需覆盖 0 ~ τ_max = 2100 μs(约 4.2 个符号) ISI 评估 D = 0.85 D_max = 4.2 中度偏强 ISI 建议 T/2 分数间隔 抽头 ≥ 2·4.2 ≈ 9

三、抽头长度怎么定

线性均衡器(如 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 间隔),兼顾覆盖与收敛速度。

抽头初值可用 CIR 的共轭时间反转作为匹配滤波器初始化,比全零初始化收敛快数倍;步长 μ 按 μ < 1/(N·P_in) 选取,P_in 为输入信号功率。

四、验证与判读

抽头长度是否合适,不能只看公式,要用数据验证:

  1. 残余ISI估计:用选定抽头数对实测 CIR 做迫零或MMSE均衡,计算均衡后CIR在符号采样点之外的残余能量占比。若残余能量 < -15 dB,认为可接受。
  2. 误码率对比:在相同信噪比下,分别用 N=5、8、10 跑误码率曲线。若 N=8 与 N=10 的 BER 差距小于 0.5 dB,说明尾径贡献可忽略,取 N=8 更经济。
  3. 收敛观察:记录 LMS 误差曲线达到稳态所需的符号数,抽头越多收敛越慢,需与帧长匹配。

若验证发现 N=10 仍无法把残余 ISI 压到 -15 dB 以下,说明信道存在时变或非线性,应转向判决反馈均衡(DFE)或频域均衡,而不是继续加长线性抽头。

τrms 与 τmax 必须来自同一组实测 CIR,不可混用不同时刻或不同位置的数据;水声信道时变性强,建议每帧或每若干帧重新估计 CIR 并自适应调整抽头。