【信道估计】基于LS和MMSE算法导频信道估计（误差率对比）附Matlab代码

✅作者简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，matlab项目合作可私信。

🍎个人主页：Matlab科研工作室

🍊个人信条：格物致知。

⛄ 内容介绍

基于最小二乘(Least Squares, LS)和最小均方误差(Minimum Mean Square Error, MMSE)的导频信道估计是常用的信道估计方法之一，用于估计无线通信系统中的信道响应。下面将介绍LS和MMSE算法的基本原理和步骤。

1. 最小二乘(LS)算法： a. 发送导频序列：在发送端，周期性地插入已知导频序列到待传输的数据序列中。 b. 接收导频序列：接收端接收到包含导频序列的信号，并进行采样和解调。 c. 估计信道响应：使用接收到的导频序列和已知导频序列进行相关运算，得到信道响应的估计值。 d. 插值和外推：对于非导频位置，可以通过插值或外推方法获得信道响应的估计值。
2. 最小均方误差(MMSE)算法： a. 发送导频序列：同样在发送端周期性地插入已知导频序列到待传输的数据序列中。 b. 接收导频序列：接收端接收到包含导频序列的信号，并进行采样和解调。 c. 估计协方差矩阵：使用接收到的导频序列计算信道响应的协方差矩阵。 d. 计算MMSE估计：通过计算信道响应的协方差矩阵的逆矩阵与接收导频序列的乘积，得到信道响应的MMSE估计值。

LS算法是一种简单直接的信道估计方法，但在噪声较大或导频序列有限的情况下，其估计精度较低。而MMSE算法考虑了信道噪声的影响，可以提供更准确的信道估计结果。

需要注意的是，LS和MMSE算法都需要已知导频序列，因此在实际应用中，需要事先设计好导频序列并在发送端进行插入。此外，为了提高信道估计的准确性，可以采用多径信道模型、时频域插值等技术进行进一步优化。

⛄ 代码

%====================================================
%  Cobayashi Laboratory , Mie University, Japan
%
%                   Last Modified by  tao,2001/05/17 
%====================================================
function out = f_GetPoints(M)

BPSKTable = f_EnGray(2);

BPSKTable1 = [1+0i -1+0i];

QPSKTable = f_EnGray(4);
QPSKTable1 = [1 -1i 1i -1];

QAM8Table7 = f_EnGray(8);
QAM8Table = [    ...
      -1 + 1j*-1    ...
      -3 + 1j*0     ...
      -1 + 1j*1     ...
       0 + 1j*3     ...
       1 + 1j*-1    ...
       0 + 1j*-3    ...
       1 + 1j*1     ...
       3 + 1j*0  ]* sqrt(4/11);

QAM8Table0 = [    ...
      3  + 1j   ...
      1  + 1j    ...
      -3 + 1j   ...
      -1 + 1j   ...
      3  + 1j*-1 ...
      1  + 1j*-1 ...
      -3 + 1j*-1   ...
      -1 + 1j*-1];

QAM8Table6 = [    ...
      2          ...
      2+2j       ...
      -2+2j  ...
      2j         ...
      2+1j*-2 ...
      1j*-2      ...
      -2         ...
      -2 + 1j*-2];

QAM8Table5 = [    ...
      3+1j*1      ...
      1+1j*3      ...
      -3+1j*1      ...
      -1+1j*3      ...
      3+1j*-1      ...
      1+1j*-3      ...
      -3+1j*-1      ...
      -1+1j*-3 ];

QAM8Table4 = [    ...
      1+1j*1      ...
      1+1j*-1      ...
      -1+1j*1      ...
      3+1j*-1      ...
      -3+1j*1      ...
      -1+1j*-3      ...
      1+1j*3      ...
      -1+1j*-1 ];


QAM8Table3 = [    ...
      3          ...
      1+1j       ...
      -1 + 1j  ...
      3j         ...
      1 + 1j*-1 ...
      1j*-3      ...
      -3         ...
      -1 + 1j*-1];


QAM8Table2 = [    ...
      1          ...
      0.5+1j*0.5       ...
      -0.5 + 1j*0.5  ...
      1j         ...
      0.5 + 1j*-0.5 ...
      1j*-1      ...
      -1         ...
      -0.5 + 1j*-0.5];

QAM8Table1 = [    ...
      1          ...
      1+1j       ...
      -1 + 1j  ...
      1j         ...
      1 + 1j*-1 ...
      1j*-1      ...
      -1         ...
      -1 + 1j*-1];

QAM16Table = f_EnGray(16);
QAM16Table1 = [ ...
      1  + 1j*1 ...
      1  + 1j*3 ...
      3  + 1j*1 ...
      3  + 1j*3 ...
      -1 + 1j*1 ...
      -1 + 1j*3 ...
      -3 + 1j*1 ...
      -3 + 1j*3 ...
      1 + 1j*-1 ...
      1 + 1j*-3 ...
      3  + 1j*-1 ...
      3  + 1j*-3 ...
      -1 + 1j*-1 ...
      -1 + 1j*-3 ...
      -3 + 1j*-1 ...
      -3 + 1j*-3];

QAM32Table2 = f_EnGray(32);
QAM32Table = [   ...
       1 + 1j*1  ...
       3 + 1j*1  ...
       3 + 1j*5  ...
       5 + 1j*1  ...
       1 + 1j*3  ...
       3 + 1j*3  ...
       1 + 1j*5  ...
       5 + 1j*3  ...
       1 + 1j*-1  ...
       1 + 1j*-3  ...
       5 + 1j*-3  ...
       1 + 1j*-5  ...
       3 + 1j*-1  ...
       3 + 1j*-3  ...
       5 + 1j*-1  ...
       3 + 1j*-5  ...
       -1 + 1j*1  ...
       -1 + 1j*3  ...
       -5 + 1j*3  ...
       -1 + 1j*5  ...
       -3 + 1j*1  ...
       -3 + 1j*3  ...
       -5 + 1j*1  ...
       -3 + 1j*5  ...
       -1 + 1j*-1  ...
       -3 + 1j*-1  ...
       -3 + 1j*-5  ...
       -5 + 1j*-1  ...
       -1 + 1j*-3  ...
       -3 + 1j*-3  ...
       -1 + 1j*-5  ...
       -5 + 1j*-3  ...
   ] * sqrt(1/10);

QAM32Table1 = [   ...
      -3 + 1j*5  ...
      -5 + 1j*1  ...
      1 + 1j*5  ...
      5 + 1j*3  ...
      3 + 1j*-5  ...
      5 + 1j*-1  ...
      -1 + 1j*-5  ...
      -5 + 1j*-3  ...
      1 + 1j*1  ...
      -1 + 1j*5  ...
      5 + 1j*1  ...
      1 + 1j*-1  ...
      -1 + 1j*-1  ...
      1 + 1j*-5  ...
      -5 + 1j*-1  ...
      -1 + 1j*1  ...
      -3 + 1j*-3  ...
      -1 + 1j*-3  ...
      -3 + 1j*1  ...
      -3 + 1j*3  ...
      3 + 1j*3  ...
      1 + 1j*3  ...
      3 + 1j*-1  ...
      3 + 1j*-3  ...
      5 + 1j*-3  ...
      3 + 1j*1  ...
      1 + 1j*-3  ...
      -3 + 1j*-5  ...
      -5 + 1j*3  ...
      -3 + 1j*-1  ...
      -1 + 1j*3  ...
      3 + 1j*5  ...
   ];


QAM64Table = f_EnGray(64);
QAM64Table1 = [ ...
      1  + 1j*1 ...
      1  + 1j*3 ...
      3  + 1j*1 ...
      3  + 1j*3 ...
      7  + 1j*1 ...
      5  + 1j*1 ...
      7  + 1j*3 ...
      5  + 1j*3 ...
      1  + 1j*7 ...
      3  + 1j*7 ...
      1  + 1j*5 ...
      3  + 1j*5 ...
      7  + 1j*7 ...
      5  + 1j*7 ...
      7  + 1j*5 ...
      5  + 1j*5 ...
      1  + 1j*-1 ...
      1  + 1j*-3 ...
      3  + 1j*-1 ...
      3  + 1j*-3 ...
      7  + 1j*-1 ...
      5  + 1j*-1 ...
      7  + 1j*-3 ...
      5  + 1j*-3 ...
      1  + 1j*-7 ...
      3  + 1j*-7 ...
      1  + 1j*-5 ...
      3  + 1j*-5 ...
      7  + 1j*-7 ...
      5  + 1j*-7 ...
      7  + 1j*-5 ...
      5  + 1j*-5 ...
      -1+ 1j*1 ...
      -1+ 1j*3 ...
      -3+ 1j*1 ...
      -3+ 1j*3 ...
      -7+ 1j*1 ...
      -5+ 1j*1 ...
      -7+ 1j*3 ...
      -5+ 1j*3 ...
      -1+ 1j*7 ...
      -3+ 1j*7 ...
      -1+ 1j*5 ...
      -3+ 1j*5 ...
      -7+ 1j*7 ...
      -5+ 1j*7 ...
      -7+ 1j*5 ...
      -5+ 1j*5 ...
      -1+ 1j*-1 ...
      -1+ 1j*-3 ...
      -3+ 1j*-1 ...
      -3+ 1j*-3 ...
      -7+ 1j*-1 ...
      -5+ 1j*-1 ...
      -7+ 1j*-3 ...
      -5+ 1j*-3 ...
      -1+ 1j*-7 ...
      -3+ 1j*-7 ...
      -1+ 1j*-5 ...
      -3+ 1j*-5 ...
      -7+ 1j*-7 ...
      -5+ 1j*-7 ...
      -7+ 1j*-5 ...
      -5+ 1j*-5];

QAM128Table1 = f_EnGray(128);
QAM128Table = [ ...
      1 + 1j*1  ...
      3 + 1j*1  ...
      1 + 1j*3  ...
      3 + 1j*3  ...
      7 + 1j*1  ...
      5 + 1j*1  ...
      7 + 1j*3  ...
      5 + 1j*3  ...
      7 + 1j*9  ...
      5 + 1j*9  ...
      7 + 1j*11 ...
      5 + 1j*11 ...
      9 + 1j*1  ...
      11 + 1j*1 ...
      9 + 1j*3  ...
      11 + 1j*3 ...
      1 + 1j*7  ...
      3 + 1j*7  ...
      1 + 1j*5  ...
      3 + 1j*5  ...
      7 + 1j*7  ...
      5 + 1j*7  ...
      7 + 1j*5  ...
      5 + 1j*5  ...
      1 + 1j*9  ...
      3 + 1j*9  ...
      1 + 1j*11 ...
      3 + 1j*11 ...
      9 + 1j*7  ...
      11 + 1j*7 ...
      9 + 1j*5  ...
      11 + 1j*5 ...
      -1 + 1j*1  ...
      -1 + 1j*3  ...
      -3 + 1j*1  ...
      -3 + 1j*3  ...
      -1 + 1j*7  ...
      -1 + 1j*5  ...
      -3 + 1j*7  ...
      -3 + 1j*5  ...
      -9 + 1j*7  ...
      -9 + 1j*5  ...
      -11 + 1j*7 ...
      -11 + 1j*5 ...
      -1 + 1j*9  ...
      -1 + 1j*11 ...
      -3 + 1j*9  ...
      -3 + 1j*11 ...
      -7 + 1j*1  ...
      -7 + 1j*3  ...
      -5 + 1j*1  ...
      -5 + 1j*3  ...
      -7 + 1j*7  ...
      -7 + 1j*5  ...
      -5 + 1j*7  ...
      -5 + 1j*5  ...
      -9 + 1j*1  ...
      -9 + 1j*3  ...
      -11 + 1j*1 ...
      -11 + 1j*3 ...
      -7 + 1j*9  ...
      -7 + 1j*11 ...
      -5 + 1j*9  ...
      -5 + 1j*11 ...
      -1 + 1j*-1  ...
      -3 + 1j*-1  ...
      -1 + 1j*-3  ...
      -3 + 1j*-3  ...
      -7 + 1j*-1  ...
      -5 + 1j*-1  ...
      -7 + 1j*-3  ...
      -5 + 1j*-3  ...
      -7 + 1j*-9  ...
      -5 + 1j*-9  ...
      -7 + 1j*-11 ...
      -5 + 1j*-11 ...
      -9 + 1j*-1  ...
      -11 + 1j*-1 ...
      -9 + 1j*-3  ...
      -11 + 1j*-3 ...
      -1 + 1j*-7  ...
      -3 + 1j*-7  ...
      -1 + 1j*-5  ...
      -3 + 1j*-5  ...
      -7 + 1j*-7  ...
      -5 + 1j*-7  ...
      -7 + 1j*-5  ...
      -5 + 1j*-5  ...
      -1 + 1j*-9  ...
      -3 + 1j*-9  ...
      -1 + 1j*-11 ...
      -3 + 1j*-11 ...
      -9 + 1j*-7  ...
      -11 + 1j*-7 ...
      -9 + 1j*-5  ...
      -11 + 1j*-5 ...
      1 + 1j*-1  ...
      1 + 1j*-3  ...
      3 + 1j*-1  ...
      3 + 1j*-3  ...
      1 + 1j*-7  ...
      1 + 1j*-5  ...
      3 + 1j*-7  ...
      3 + 1j*-5  ...
      9 + 1j*-7  ...
      9 + 1j*-5  ...
      11 + 1j*-7 ...
      11 + 1j*-5 ...
      1 + 1j*-9  ...
      1 + 1j*-11 ...
      3 + 1j*-9  ...
      3 + 1j*-11 ...
      7 + 1j*-1  ...
      7 + 1j*-3  ...
      5 + 1j*-1  ...
      5 + 1j*-3  ...
      7 + 1j*-7  ...
      7 + 1j*-5  ...
      5 + 1j*-7  ...
      5 + 1j*-5  ...
      9 + 1j*-1  ...
      9 + 1j*-3  ...
      11 + 1j*-1 ...
      11 + 1j*-3 ...
      7 + 1j*-9  ...
      7 + 1j*-11 ...
      5 + 1j*-9  ...
      5 + 1j*-11] * sqrt(1/41);


QAM256Table = f_EnGray(256);

switch M,
case 0,   table = [0];
case 2,   table = BPSKTable;
case 4,   table = QPSKTable;
case 8,   table = QAM8Table;
case 16,  table = QAM16Table;
case 32,  table = QAM32Table;
case 64,  table = QAM64Table;
case 128, table = QAM128Table;
case 256, table = QAM256Table;
otherwise, error('No such MQAM!');
end

out = table;

⛄ 运行结果

⛄ 参考文献

[1] 朱竞.基于导频的MB-LSF算法的OFDM系统信道估计研究[D].华东理工大学,2014.

[2] 王宏宇.快时变信道下OFDM系统频率同步与信道估计技术研究[D].西南交通大学,2014.

⛳️ 代码获取关注我

❤️部分理论引用网络文献，若有侵权联系博主删除

❤️ 关注我领取海量matlab电子书和数学建模资料

🍅 仿真咨询

1.卷积神经网络(CNN)、LSTM、支持向量机(SVM)、最小二乘支持向量机(LSSVM)、极限学习机(ELM)、核极限学习机(KELM)、BP、RBF、宽度学习、DBN、RF、RBF、DELM实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

2.图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

3.旅行商问题(TSP)、车辆路径问题(VRP、MVRP、CVRP、VRPTW等)、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划

4.无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

5.传感器部署优化、通信协议优化、路由优化、目标定位

6.信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号

7.生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化

8.微电网优化、无功优化、配电网重构、储能配置

9.元胞自动机交通流 人群疏散 病毒扩散 晶体生长