Q2.zip_16-QAM MIMO_mimo 16 QAM _mimo QAM_mimo 信号检测_qam mmse

clc clear all; close all; %+++++++++++++++++++ 初始化变量 ++++++++++++++++++++++++++++++++++++++++++ NT = 4; %发送天线个数 NR = 4; %接收天线个数 No = 10000; %执行次数 Es_N0_dB = [0:2:20]; %信噪比范围 bit_errors_ZF(1,length(Es_N0_dB)) = zeros; %单独ZF误码率初始化 BER_ZF(1,length(Es_N0_dB)) = zeros; bit_errors_MMSE(1,length(Es_N0_dB)) = zeros; %单独MMSE误码率初始化 BER_MMSE(1,length(Es_N0_dB)) = zeros; bit_errors_ML(1,length(Es_N0_dB)) = zeros; %ML误码率初始化 BER_ML(1,length(Es_N0_dB)) = zeros; %==================== New Detection ===================================== disp('New Detection') for xx = 1:length(Es_N0_dB) xx for times = 1:No sigma = 10^(-Es_N0_dB(xx)/20); %信噪比转换 %++++++++++++++++++ 输入数据部分 ++++++++++++++++++++++++++++++++++++++++ x1 = randint(1,NT); x2= randint(1,NT); %++++++++++++++++++ QPSK调制部分 ++++++++++++++++++++++++++++++++++++++++ x_mapper = QPSK_mapper(x1,x2); %++++++++++++++++++ 串并转换部分，映射到NT根天线上 ++++++++++++++++++++++++ s = reshape(x_mapper,NT,1); %++++++++++++++++++ 加性高斯白噪声信道 +++++++++++++++++++++++++++++++++++ H = (randn(NR,NT) + i*randn(NR,NT))/sqrt(2); N = (randn(NR,1) + i*randn(NR,1))/sqrt(2); %++++++++++++++++++ 接收部分 ++++++++++++++++++++++++++++++++++++++++++++++ y = H*s + sigma*N; %++++++++++++++ ML检测部分 +++++++++++++++++++++ x_ML = ML_detector_QPSK(y,H); x_r=reshape(x_ML,NT,1); bit_errors_ML(xx) = bit_errors_ML(xx) + length(find(x_r-s~=0)); %++++++++++++++++++ 单独ZF检测 ++++++++++++++++++++++++++++++++++++ x_ZF = Q1(inv(H'*H)*H'*y); bit_errors_ZF(xx) = bit_errors_ZF(xx) + length(find(x_ZF - s~=0)); %++++++++++++++++++ 单独MMSE检测 ++++++++++++++++++++++++++++++++++ x_MMSE = Q1(inv(H'*H+sigma^2*eye(NT))*H'*y); bit_errors_MMSE(xx) = bit_errors_MMSE(xx) + length(find(x_MMSE-s~=0)); %++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ end BER_ML(xx) = bit_errors_ML(xx)/(NT*No); BER_ZF(xx) = bit_errors_ZF(xx)/(NT*No); BER_MMSE(xx) = bit_errors_MMSE(xx)/(NT*No); end disp('--------------------- end of file ----------------------------------'); %----------------误码率的计算---------------------------------- semilogy(Es_N0_dB,BER_ZF,'d-g',Es_N0_dB,BER_MMSE,'o-r',Es_N0_dB,BER_ML,'p-c') %粉红正方形 grid on % axis([0 20 10^-4 10^0]) legend('ZF','MMSE','ML'); xlabel('Es/No,dB'); ylabel('Bit Error Rate'); title('ZF,MMSE,ML仿真图');