V-BLAST基本信号检测.rar_v-blast_信号检测

clear all; close all; clc; tic Nt=4; %发射天线数 Nr=4; %接收天线数 L=10000; %仿真次数 共L组 EbNo=0:2:20; %信噪比dB 区间 M=4; %映射星座点个数 x=round(rand(Nt,L)*(M-1)); %均匀分布(并行数据) 4*10000 一次取一列4*1发送 s=qammod(x,M);% 4QAM调制 s=s/sqrt(2); %信号功率归一化 %%%%% for index=1:length(EbNo); x1=[]; x2=[]; x3=[]; x4=[]; x5=[]; for index1=1:L h=randn(Nr,Nt)+1i*randn(Nr,Nt); % 瑞利衰落信道4*4 h=h./sqrt(2); %信道系数归一化 sigma=sqrt(1/(10.^(EbNo(index)/10))); %每根接收天线的噪声标准差 n=sigma/sqrt(2)*(randn(Nr,1)+1i*randn(Nr,1)); %Nr个接收天线的高斯白噪声4*1 %每根接收天线的噪声功率都是sigma^2 %%%%%ZF和MMSE滤波矩阵%%%%% w_zf=inv(h'*h)*h'; %%%%%%%%%%%%%%%%ZF%%%%%%%%%%%%%%%%%% w_mmse=h'*inv(h*h'+sigma.^2*diag(ones(1,Nt)));%%%MMSE%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% y=h*s(:,index1)+n; %信号经过信道 一列4*1 yy=y; %%%信号检测%%% y1=w_zf*yy; %4*1 y2=w_mmse*yy; %4*1 y3=osic_zf(yy,h); %4*1 OSIC y4=osic_mmse(yy,h,sigma); %4*1 MMSE y5=ML_algorithm(yy,h,M); %%%%%%%%%%%%% temp1=qamdemod(y1,M); temp2=qamdemod(y2,M); temp3=y3; %4*1 temp4=y4; %4*1 temp5=qamdemod(y5,M); %4*1 x1=[x1;temp1']; %10000*4 xx1=x1'; x2=[x2;temp2']; xx2=x2'; x3=[x3;temp3']; xx3=x3'; x4=[x4;temp4']; xx4=x4'; x5=[x5;temp5']; xx5=x5'; end %求误码率 [~,ber1(index)]=biterr(x,xx1,log2(M)); [~,ber2(index)]=biterr(x,xx2,log2(M)); [~,ber3(index)]=biterr(x,xx3,log2(M)); [~,ber4(index)]=biterr(x,xx4,log2(M)); [~,ber5(index)]=biterr(x,xx5,log2(M)); end semilogy(EbNo,ber1,'*-k',EbNo,ber2,'o-r',EbNo,ber3,'d-b',EbNo,ber4,'+-m',EbNo,ber5,'x-g'); legend('ZF','MMSE','OSICZF','OSICMMSE','ML'); xlabel('SNR(dB)'); ylabel('BER'); grid on; toc