clear;?clf;
N_P=1000;?
Na=20;?Nb=4;?Ne=4;?O=0.8;
SNRdB=[-6:2:12];??
for?i_SNR=1:length（SNRdB）
????SNR=SNRdB（i_SNR）;?n=0;
???for?i_packet=1:N_P
?????????s_1=round（rand（12））?;
????????m=s_1*[2;1]+1;
????????I=eye（4）;?e=I（:m）;
?????????s_2=round（rand（14））;
?????????[Bsym_table]?=?QAM（s_2）;
?????????S=e*B;
?????????Hb=?sqrt（0.5）*（randn（NbNa）?+?j*randn（NbNa））;
?????????He=?sqrt（0.5）*（randn（NeNa）?+?j*randn（NeNa））;
?????????[ULV]?=?svd?（Hb）?;V0=V（:（Nb+1）:Na）;
?????????v=trace（inv（Hb*Hb‘））;
?????????ci=sqrt（（O*Nb）/trace（inv（Hb*Hb‘）））;%功率約束系數
?????????Or=（1-O）/（Na-Nb）;??????????????????%人工噪聲的方差?
?????????P1=ci*Hb‘*（inv（Hb*Hb‘））;P2=V0*sqrt（Or）*（randn（Na-NbNb）?+?j*randn（Na-NbNb））;
?????????Y=ci*S+sqrt（0.5/（10^（SNR/10））