空间调制代码

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））