MIMO MMSE、ZF检测

clc
clear all
close all
%%%%%%%%%%%parameter
N=10^7;
SNR_dB=15:5:50;
SNR=10.^（SNR_dB/10）;
nTx=2;
nRx=2;
%%%%%%%%%%%produce QPSK signals
for i=1:length（SNR_dB）
    S0=rand（1N）>0.5;
    S=zeros（1N/2）;
    for n=1:length（S0）/2
        p=S0（2*n-1）;
        q=S0（2*n）;
        if （p==0）&&（q==0）   %45 
            S（n）=exp（1j*pi/4）;
        end
        if （p==1）&&（q==0）   %135
            S（n）=exp（1j*3*pi/4）;
        end
        if （p==1）&&（q==1）   %225
            S（n）=exp（1j*5*pi/4）;
        end
        if （p==0）&&（q==1）   %315
            S（n）=exp（1j*7*pi/4）; 
        end
    end
   %%%%%%%%transmitter
    smod=reshape（S[nTx1N/nTx/2]）; 
    y=zeros（nRx1）;
    re=zeros（nTx1size（smod3））;
    re1=zeros（nTx1size（smod3））;
    for t=1:size（smod3）
        h=sqrt（1/2）*（randn（n