高速移动下，无线多径的OFDM-MIMO信道的LS估计，给出MSE与SNR关系图的代码

以下是一个MATLAB代码示例，可以绘制OFDM-MIMO信道的LS估计MSE与SNR关系图。

clear all;
clc;
Nt=4;%发送天线数
Nr=4;%接收天线数
N=64;%FFT长度
M=16;%调制阶数
L=4;%码元数
Lb=2;%数据块数
Ls=50;%仿真次数
SNR=0:2:40;%信噪比范围
MSE=zeros(Ls,length(SNR));
for sim=1:Ls
    for ss=1:length(SNR)
        snr=10^(SNR(ss)/10);
        H=1/sqrt(2)*(randn(Nr,Nt)+j*randn(Nr,Nt));%随机信道矩阵
        H_LS=H(1:L,:);%LS估计信道矩阵
        for bb=1:Lb
            b=randi([0 M-1],L,Nt);%随机发送符号
            x=pskmod(b,M);%调制
            for ll=1:L
                tx(:,ll)=ifft(x(ll,:),N);%IFFT
                tx_cp(:,ll)=[tx(N-15+1:N,ll);tx(:,ll)];%循环前缀
                tx_mimo(:,ll)=reshape(tx_cp(:,ll)*sqrt(N/L),Nr*L,1);%并行发送
                rx_mimo(:,ll)=awgn(H_LS*tx_mimo(:,ll),snr,'measured');%加入高斯白噪声
                rx_cp(:,ll)=reshape(rx_mimo(:,ll),Nt,L);%并行接收
                rx(:,ll)=rx_cp(16+1:N,:);%循环前缀删除
                x_hat(:,ll)=fft(rx(:,ll),N)/sqrt(N/L);%FFT
            end
            b_hat=pskdemod(x_hat(:),M);%解调
            error=length(find(b_hat-b(:)));%误码数
            MSE(sim,ss)=MSE(sim,ss)+error/(L*Nt*Lb);%均方误差
        end
    end
end
MSE_ave=mean(MSE,1);%均方误差平均值
semilogy(SNR,MSE_ave,'-*');
xlabel('信噪比(dB)');
ylabel('均方误差');
title('OFDM-MIMO信道LS估计MSE与SNR关系图');
grid on;

需要注意的是，这只是一个简单的示例代码，实际应用中可能需要根据具体情况进行修改和优化。