高速移动下多径信道OFDM-MIMO信道LS估计求MSE与信噪比关系图的代码

以下是基于MATLAB的OFDM-MIMO信道LS估计求MSE与信噪比关系图的代码：

clear all;
close all;

% 设置参数
M = 4; % 调制阶数
Nt = 2; % 发射天线数
Nr = 2; % 接收天线数
SNR_dB = 0:5:30; % 信噪比范围
numIter = 1000; % 模拟次数

% 生成QPSK调制信号
data = randi([0 M-1], Nt, 1);

% 构造发送信号
x = qammod(data, M);

% 构建MIMO信道
H = (randn(Nr, Nt) + 1j*randn(Nr, Nt))/sqrt(2);

% 构造接收信号
for i = 1:length(SNR_dB)
    SNR = 10^(SNR_dB(i)/10);
    noise = sqrt(1/SNR)*(randn(Nr, 1) + 1j*randn(Nr, 1))/sqrt(2);
    y(:,i) = H*x + noise;
end

% LS信道估计
H_LS = zeros(Nr, Nt, length(SNR_dB));
for i = 1:length(SNR_dB)
    H_LS(:,:,i) = y(:,i)*x'/(x*x');
end

% 计算MSE
MSE_LS = zeros(length(SNR_dB), 1);
for i = 1:length(SNR_dB)
    for j = 1:numIter
        noise = sqrt(1/SNR)*(randn(Nr, 1) + 1j*randn(Nr, 1))/sqrt(2);
        H_est = y(:,i)*x'/(x*x');
        MSE_LS(i) = MSE_LS(i) + norm(H - H_est, 'fro')^2;
    end
end
MSE_LS = MSE_LS/numIter;

% 绘制MSE与信噪比关系图
figure;
semilogy(SNR_dB, MSE_LS, 'r-o');
title('LS信道估计MSE与信噪比关系图');
xlabel('信噪比 (dB)');
ylabel('MSE');
grid on;

该代码使用LS（最小二乘）算法对OFDM-MIMO信道进行信道估计，并绘制了MSE与信噪比之间的关系图。