该楼层疑似违规已被系统折叠 隐藏此楼查看此楼
clear,clf
N_frame=130;%帧数
N_packet=4000;%分组数
b=2;%Set to1/2/3/4forBPSK/QPSK/16QAM/64QAM
SNRdBs=[0:2:20];sq2=sqrt(2);
for iter=1:3
if iter==1
NT=1;NR=1;gs='-kx';%SISO
elseif iter==2
NT=1;NR=2;gs='-^';%1x2
else
NT=1;NR=4;gs='-ro';%1x4
end
sq_NT=sqrt(NT);
for i_SNR=1:length(SNRdBs)
SNRdB=SNRdBs(i_SNR);sigma=sqrt(0.5/(10^(SNRdB/10)));
for i_packet=1:N_packet
msg_symbol=randint(N_frame*b,NT);
[temp,sym_tab,P]=modulator(msg_symbol.',b);
X=temp.';
Hr=(randn(N_frame,NR)+ j*randn(N_frame,NR))/sq2;
H=reshape(Hr,N_frame,NR);Habs=sum(abs(H).^2,2);Z=0;
for i=1:NR
R(:,i