【MATLAB源码-第446期】基于MATLAB的水声时变多径信道OFDM系统仿真对比:LS、LMMSE、LMS与RLS
操作环境MATLAB 2024a1、算法描述摘要水声通信信道具有带宽受限、传播时延大、多径扩展明显和时变性强等特点这些因素会直接破坏接收信号的稳定性使传统单载波通信方式在高速数据传输场景中面临较大的均衡压力。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从低到高增加各方法总体误码率呈下降趋势。LSZF在低SNR时性能较差并且在高SNR区域仍存在较明显误码平台主要原因是LS估计误差和ZF噪声放大会共同影响判决。LMMSEMMSE表现最好误码率随SNR提升下降更明显说明更准确的信道估计和更稳健的均衡方式能够有效改善系统性能。LMS和RLS方法虽然具备自适应跟踪能力但在当前参数下没有明显超过LMMSEMMSE说明自适应算法需要精细调参不能只靠算法名称判断优劣。星座图能够直观反映均衡效果。均衡前星座点分布明显发散说明多径信道和噪声已经破坏了QPSK星座结构。经过LSZF、LMMSEMMSE和RLS处理后星座点向理想判决区域收敛其中LMMSEMMSE通常表现更集中。这与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反映自适应均衡收敛过程星座图反映判决空间中的信号聚集情况。这些图像相互支撑使结果分析不至于单薄。第三个特点是对比方法具有层次。LSZF属于低复杂度基准方案LMMSEMMSE属于利用统计信息的稳健方案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.2、仿真结果演示3、关键代码展示略4、MATLAB 源码获取V点击下方名片关注公众号获取