当前位置：首页 > news >正文

【MATLAB源码-第446期】基于MATLAB的水声时变多径信道OFDM系统仿真对比：LS、LMMSE、LMS与RLS

news 2026/5/24 1:18:19

操作环境：

MATLAB 2024a

1、算法描述

摘要

水声通信信道具有带宽受限、传播时延大、多径扩展明显和时变性强等特点，这些因素会直接破坏接收信号的稳定性，使传统单载波通信方式在高速数据传输场景中面临较大的均衡压力。OFDM通过将宽带高速数据流划分到多个低速正交子载波上，可以在一定程度上降低多径频率选择性衰落带来的影响，因此常被用于水声通信系统的仿真与算法验证。已有研究也指出，水声OFDM系统的主要难点集中在多普勒扩展、时变多径、导频资源受限以及信道状态信息难以精确获得等方面。

本文基于MATLAB构建了一个水声时变多径信道下的QPSK-OFDM通信系统仿真平台。系统包括参数配置、QPSK映射、OFDM调制、循环前缀添加、水声多径信道建模、多普勒扰动、加性复高斯噪声、导频辅助信道估计、频域均衡、误码率统计以及多类结果图输出。接收端重点比较LS、LMMSE、LMS和RLS相关处理方法在信道估计精度、均衡收敛行为和误码率性能上的差异。仿真结果表明，在当前模型参数下，LMMSE结合MMSE均衡的综合表现明显优于单纯LS结合ZF均衡；LMS和RLS自适应均衡能够体现对时变信道的跟踪能力，但在当前参数设置下，其BER改善并未明显超过LMMSE方案。这一结果说明，水声OFDM系统不能只追求算法名称上的复杂度，而必须结合导频设计、信道统计信息、步长参数、遗忘因子和噪声方差估计进行整体优化。

关键词：水声通信；OFDM系统；QPSK调制；LS信道估计；LMMSE信道估计；LMS均衡；RLS均衡；误码率

1. 引言

水声通信是海洋观测、无人潜航器、海底传感网络和水下信息传输中的关键技术之一。与无线电通信不同，水下环境中电磁波衰减严重，声波成为更可行的远距离信息传输载体。但声波传播速度远低于电磁波，水下介质又受到海面波动、海底反射、温盐深剖面变化、收发平台运动等因素影响，因此水声信道往往呈现长时延扩展、频率选择性衰落和时间选择性衰落并存的特点。对于高速通信系统而言，这些特征会造成严重码间干扰和子载波间干扰，使接收端必须依赖有效的信道估计和均衡算法。

OFDM的优势在于能够把频率选择性宽带信道拆分为多个近似平坦的窄带子信道，从而降低均衡复杂度。对于多径明显的水声场景，循环前缀可以在一定范围内吸收多径时延扩展，频域单抽头均衡则可以简化接收机结构。但水声OFDM并不是简单把无线通信中的OFDM直接搬到水下。水声信道的多普勒影响更强，路径间多普勒可能不一致，信道频率响应会随OFDM符号时间变化而波动。如果多普勒估计和信道跟踪不足，OFDM子载波正交性会受到破坏，接收星座会扩散，误码率会出现误差平台。相关研究明确指出，多普勒估计与补偿是水声OFDM移动通信设计中的关键环节。

本文代码围绕“水声时变多径信道下的OFDM信道估计与均衡”展开。系统采用128点FFT，设置循环前缀、导频间隔、QPSK调制、多径时延、路径功率衰减和多普勒频移，通过不同信噪比下的仿真比较LS、LMMSE、LMS和RLS方法。该设计的价值不在于堆叠复杂模块，而在于把一个水声OFDM接收链路的核心问题清楚呈现出来：第一，信道估计是否准确；第二，均衡后星座是否收敛；第三，BER是否随SNR提升而有效下降；第四，自适应算法是否真正带来稳定收益。

2. 系统模型与仿真流程

本系统以MATLAB为仿真平台，主程序为main.m，辅助函数放置在function文件夹中。主程序首先完成环境清理、路径添加、随机种子固定和图像目录初始化，然后调用参数配置函数生成系统参数。参数包括FFT点数、循环前缀长度、子载波间隔、采样率、OFDM符号数、训练符号数、导频间隔、信噪比扫描范围、快照信噪比、自适应均衡参数以及调制方式。这样的结构比较清楚，便于后续修改单个模块，不会让主程序堆成一大段难以维护的代码。

发射端首先产生随机比特流，并将其映射为QPSK符号。QPSK调制采用单位能量归一化处理，使不同信噪比下的噪声添加更容易控制。随后，系统将数据符号填入数据子载波，将导频符号填入固定间隔的导频子载波。导频采用恒模相位序列，可以避免导频幅度起伏过大对信道估计造成额外影响。完成频域资源映射后，系统通过IFFT生成时域OFDM符号，并在每个符号前添加循环前缀。循环前缀的作用是缓解多径延迟造成的符号间干扰，同时维持频域均衡的简化条件。

信道部分构造了一个具有多条路径的水声信道。每条路径包含采样点时延、相对功率衰减、随机初始相位和多普勒频移。接收信号由不同路径的延迟分量叠加而成，并且每条路径都会随时间产生相位旋转。这个模型虽然不是完整海试信道模型，但已经抓住了水声OFDM仿真中最关键的两个因素：多径和时变。多径决定频率选择性衰落，多普勒决定信道随时间变化，两者共同影响接收端信道估计和均衡性能。

接收端首先去除循环前缀，然后通过FFT恢复频域子载波。系统根据目标信噪比添加复高斯白噪声，并保存理论频域信道响应用于估计误差评价。之后，接收端分别进行LS和LMMSE信道估计。LS方法直接利用导频位置的接收值与发送导频进行比值估计，再通过插值得到全子载波信道响应。LMMSE方法进一步利用由多径结构构造的频域相关矩阵，在噪声方差已知的条件下对全频域信道进行统计意义上的估计。已有研究也表明，水声OFDM信道估计常常需要利用信道稀疏性、时间相关性或统计先验信息来提升估计质量，而不是只依赖简单导频插值。

3. 信道估计方法分析

LS信道估计的优点是结构简单、计算量低、实现直接，对硬件资源或实时性要求较高的场景有一定吸引力。但LS方法的缺点同样明显：它不主动利用信道统计特征，也不主动抑制噪声。在导频密度有限、信噪比较低或信道频率响应剧烈变化时，LS估计容易产生较大波动。代码中的LS方案在导频点获得估计值后采用插值扩展到所有子载波，因此插值质量会受到导频间隔和信道频率选择性强弱的影响。如果某些频段存在深衰落，LS插值曲线可能无法准确跟随真实信道，从而导致后续ZF均衡放大噪声。

LMMSE信道估计相比LS更加稳健。它通过频域相关矩阵把导频位置的信息扩展到全子载波，同时考虑噪声方差的影响。直观地看，LMMSE不是简单相信每个导频点的瞬时观测，而是在“观测值”和“信道统计结构”之间做折中。因此在噪声存在时，LMMSE能够明显降低信道估计的随机波动。代码输出的信道估计图也体现了这一点：在代表性信噪比下，LMMSE曲线对真实信道幅频响应的跟踪更平滑，NMSE曲线也明显低于LS。这个结论是合理的，因为当前信道模型的多径参数已知，LMMSE能够利用这些先验结构获得估计增益。

需要指出的是，LMMSE的优势不是免费的。它依赖信道相关矩阵和噪声方差估计，如果实际信道统计特征与模型不匹配，或者噪声方差估计偏差较大，LMMSE性能会下降。在工程系统中，LMMSE的实现成本也高于LS，尤其当子载波数量增加、导频结构复杂或信道统计矩阵需要频繁更新时，矩阵运算会带来更高计算负担。深度学习和超分辨率方法近年来也被用于水声OFDM信道估计，相关研究将导频稀疏观测下的信道恢复看作类似图像超分辨率的问题，以减少导频占用并提升估计性能。但对于本文代码而言，采用LS和LMMSE作为基准方案是合理的，因为它们更适合作为传统接收机性能对比的基础。

4. 自适应均衡方法分析

本文代码还引入了LMS和RLS自适应处理。自适应均衡的核心目的，是在信道随时间变化时持续跟踪信道状态，避免固定估计在后续符号中失效。LMS算法结构简单，通过误差信号逐步调整估计值，对计算资源要求较低，适合低复杂度实现。但LMS的性能高度依赖步长参数。步长过大可能导致收敛振荡，步长过小则会导致跟踪速度不足。水声信道具有慢变和快变交替的特点，单一固定步长很难在稳定性和跟踪速度之间同时取得最优。

RLS算法通常具有更快的收敛速度和更强的跟踪能力，但计算复杂度更高，对数值稳定性也更敏感。代码中RLS通过遗忘因子和逆相关量更新导频位置的信道估计，然后再插值到全子载波进行MMSE均衡。RLS理论上可以比LMS更快利用历史观测信息，但在当前仿真参数下，BER曲线显示RLS并没有显著超过LMMSE方案。这一点不能简单理解为RLS算法“不好”，而应理解为当前模型、导频结构、遗忘因子、初始值和统计信道估计共同决定了最终效果。已有针对水声通信自适应均衡的比较研究也指出，RLS相对LMS并不总是表现出压倒性优势，实际收益与阵列规模、信道变化速度和接收结构有关。

从工程判断看，本文代码最值得保留的地方是同时输出均衡收敛曲线和BER曲线。只看BER容易忽略算法内部过程，只看MSE又可能无法反映最终判决质量。LMS/RLS均衡误差曲线能够显示算法在训练段和后续符号中的动态变化，BER曲线则直接反映通信可靠性。若后续要增强项目深度，可以增加多帧蒙特卡洛平均、调整LMS步长、调整RLS遗忘因子、改变导频密度，并加入多普勒补偿模块。否则，自适应均衡的性能结论会受单次随机帧影响，曲线稳定性还不够强。

5. 仿真结果分析

项目输出了七张主要图像，分别覆盖信道结构、多普勒响应、发射接收频谱、信道估计效果、均衡收敛、BER对比和星座图。信道结构图展示了多条路径的相对幅度和时延分布，可以看出系统模拟的是一个具有明显延迟扩展的多径环境。多普勒响应图进一步显示不同路径存在不同频移，说明信道不仅在频率上选择性衰落，也随时间发生变化。这类双选择性特征正是水声OFDM系统难处理的根源。

频谱图比较了发射端和接收端信号的频域分布。接收端频谱受到多径衰落和噪声影响后不再完全保持发射端形态，部分频率区域会出现起伏。这说明接收端必须进行信道估计和均衡，否则直接判决会造成严重星座扩散。信道估计图中，LS估计在部分深衰落附近出现明显偏差，而LMMSE估计整体更接近真实信道。NMSE曲线也显示，随着SNR提高，LS估计误差下降有限，而LMMSE估计误差持续降低，说明统计相关信息确实改善了估计精度。

BER曲线是本文最关键的结果之一。随着SNR从低到高增加，各方法总体误码率呈下降趋势。LS+ZF在低SNR时性能较差，并且在高SNR区域仍存在较明显误码平台，主要原因是LS估计误差和ZF噪声放大会共同影响判决。LMMSE+MMSE表现最好，误码率随SNR提升下降更明显，说明更准确的信道估计和更稳健的均衡方式能够有效改善系统性能。LMS和RLS方法虽然具备自适应跟踪能力，但在当前参数下没有明显超过LMMSE+MMSE，说明自适应算法需要精细调参，不能只靠算法名称判断优劣。

星座图能够直观反映均衡效果。均衡前星座点分布明显发散，说明多径信道和噪声已经破坏了QPSK星座结构。经过LS+ZF、LMMSE+MMSE和RLS处理后，星座点向理想判决区域收敛，其中LMMSE+MMSE通常表现更集中。这与BER曲线结论一致。需要强调的是，星座图不是装饰图，而是验证接收处理是否有效的重要证据。如果星座图没有收敛，即使BER曲线短期下降，也要警惕仿真统计不足或判决逻辑存在偏差。

6. 系统实现结构

本项目的代码结构较清楚。main.m作为主入口，只负责流程控制，不把所有算法细节堆在主程序中。cfg_ofdm.m用于统一管理系统参数，便于后续修改FFT点数、循环前缀长度、导频间隔和SNR范围。mk_ofdm.m、qpsk_mod.m和mk_pilot.m构成发射端模块，rx_ofdm.m、uwa_chan.m和true_h.m构成信道与接收预处理模块，chan_est.m和chan_cov.m构成信道估计模块，eq_adapt.m负责LMS/RLS自适应处理，get_ber.m用于误码率统计，各类plot函数负责图像输出。

这种拆分方式有利于后续扩展。例如，若要加入16QAM调制，只需要扩展调制和解调函数，并相应调整比特映射；若要加入多帧平均，只需要在run_ber.m中增加外层循环；若要加入多普勒补偿，可以在rx_ofdm.m之后、FFT之前增加重采样或相位补偿模块；若要加入深度学习信道估计，可以在chan_est.m之后增加神经网络估计分支，再与LS、LMMSE进行统一BER比较。当前结构已经具备继续扩展的基础。

代码中另一个合理点是所有图像统一保存为PNG，结果变量统一保存到result.mat。这符合仿真项目复现实验的基本要求。figure文件夹中的图像可以直接用于论文结果分析，result.mat可以用于后续复查参数和曲线数据。相比只在命令行输出结果，这种保存方式更适合反复调试和报告撰写。

7. 系统特点

本系统的第一个特点是主题聚焦明确。它不是泛泛展示OFDM调制解调，而是围绕水声时变多径信道这一具体场景展开，重点比较信道估计和均衡算法对系统性能的影响。这样的选题有实际意义，因为水声通信中的多径和多普勒问题是接收机设计的核心矛盾。

第二个特点是评价指标完整。系统不仅输出BER曲线，还输出NMSE、均衡MSE、信道频响、频谱和星座图。BER反映最终通信可靠性，NMSE反映信道估计精度，MSE反映自适应均衡收敛过程，星座图反映判决空间中的信号聚集情况。这些图像相互支撑，使结果分析不至于单薄。

第三个特点是对比方法具有层次。LS+ZF属于低复杂度基准方案，LMMSE+MMSE属于利用统计信息的稳健方案，LMS和RLS属于自适应跟踪方案。四类结果放在同一仿真条件下比较，可以说明不同算法的适用边界，而不是孤立展示某一个算法。

第四个特点是代码结构便于二次开发。主函数、参数函数、信道函数、估计函数、均衡函数和绘图函数分工清楚，后续可以较容易地扩展调制阶数、信道模型、导频结构和算法分支。这比把所有内容写在一个脚本里更适合长期维护。

8. 结论

本文基于MATLAB实现了水声时变多径信道下的QPSK-OFDM通信系统仿真，并围绕LS、LMMSE、LMS和RLS相关方法进行了信道估计与均衡性能对比。系统通过多径时延、路径衰减和路径多普勒构造水声信道，使接收信号同时受到频率选择性衰落和时间变化影响。仿真结果表明，简单LS估计配合ZF均衡在噪声和深衰落环境下容易出现误差放大，BER性能有限；LMMSE估计配合MMSE均衡能够有效利用信道统计相关性和噪声方差信息，在NMSE和BER上均表现更优；LMS和RLS自适应均衡能够体现对时变信道的跟踪能力，但其效果依赖参数设置和导频结构，在当前配置下没有明显超过LMMSE方案。

整体来看，该项目适合作为水声OFDM接收机仿真的基础框架。它覆盖了发射、信道、接收、估计、均衡、判决和结果可视化的完整链路，也能够直观说明信道估计质量对最终误码率的决定性作用。后续若要进一步提升研究深度，应优先增加多帧平均、多普勒补偿、导频优化和参数敏感性分析，而不是盲目增加算法名称。真正有价值的改进，应体现在BER、NMSE、星座收敛和复杂度之间的综合收益上。

参考文献

Rudander, Jacob Hiden, et al. “Comparing RLS and LMS Adaptive Equalizers for Large Hydrophone Arrays in Underwater Acoustic Communication Channels.”Proceedings of OCEANS 2019 - Marseille, IEEE, 2019. doi:10.1109/OCEANSE.2019.8867313.
Wan, Lei, et al. “Fine Doppler Scale Estimations for an Underwater Acoustic CP-OFDM System.”Signal Processing, vol. 170, 2020, article 107439. doi:10.1016/j.sigpro.2019.107439.
Wang, Shuaijun, Mingliu Liu, and Deshi Li. “Bayesian Learning-Based Clustered-Sparse Channel Estimation for Time-Varying Underwater Acoustic OFDM Communication.”Sensors, vol. 21, no. 14, 2021, article 4889. doi:10.3390/s21144889.
Ouyang, Donghong, Yuzhou Li, and Zhizhan Wang. “Channel Estimation for Underwater Acoustic OFDM Communications: An Image Super-Resolution Approach.”IEEE International Conference on Communications, 2021. doi:10.1109/ICC42927.2021.9500819.
Li, Bin, et al. “Cross-Correlation Quasi-Gradient Doppler Estimation for Underwater Acoustic OFDM Mobile Communications.”Applied Acoustics, vol. 190, 2022, article 108640. doi:10.1016/j.apacoust.2022.108640.
Cho, Young-Hwan. “Fast Sparse Bayesian Learning-Based Channel Estimation for Underwater Acoustic OFDM Systems.”Applied Sciences, vol. 12, no. 19, 2022, article 10175. doi:10.3390/app121910175.
Feng, Xiao, et al. “Channel Estimation for Underwater Acoustic OFDM Communications via Temporal Sparse Bayesian Learning.”Signal Processing, vol. 207, 2023, article 108951. doi:10.1016/j.sigpro.2023.108951.
Cui, Xuerong, et al. “Channel Estimation for Underwater Acoustic OFDM Based on Super-Resolution Network.”Internet Technology Letters, 2024. doi:10.1002/itl2.496.
Guo, Jiasheng, et al. “Underwater-Acoustic-OFDM Channel Estimation Based on Deep Learning and Data Augmentation.”Electronics, vol. 13, 2024, article 689. doi:10.3390/electronics13040689.
Zhao, Hao, et al. “Deep Learning Aided Underwater Acoustic OFDM Receivers: Model-Driven or Data-Driven?”Digital Communications and Networks, 2025. doi:10.1016/j.dcan.2024.10.006.