ANC物理层网络编码

clearclc
frame=1000;
N=10000;
SNR_DB=0:1:16;
SNR=10.^（SNR_DB/10）;
VAR=1./SNR;
bpsk1=[-11];
bpsk2=[-ii];
for m=1:length（SNR）
    snr=SNR（m）;
    var=VAR（m）;
    errors（m）=0;
    s1（1）=1;
    s2（1）=1;
    for k=1:N
        h1=（randn（11）+i*randn（11））;
        h2=（randn（11）+i*randn（11））;  
        c1=randint（1（frame-1））;
        c2=randint（1（frame-1））;
        c1_mod=bpsk1（c1+1）;
        c2_mod=bpsk2（c2+1）;
        nr=sqrt（var）*（randn（1frame）+i*randn（1frame））;
   %    nr=zeros（frame）; %检错 
        yr（1）=h1*s1（1）+h2*s2（1）+nr（1）;
    for n=2:frame
        %公式（1）（2）
        s1（n）=s1（n-1）*c1_mod（n-1）;
        s2（n）=s2（n-1）*c2_mod（n-1）; 
        %公式（3）       
        yr（n）=h1*s1（n）+h2*s2（n）+nr（n）;      
    end           
       %AF放大系数
       %公式（7）
     a_demod1=sqrt（（mean（yr*yr‘））/frame）;  
       a_demod=1/a_demod1;
    
   %      a_demod=（h1*conj（h1）+h2*conj（h2）+var）^（0.5）;  %理想值  
   %      a_demod=1/a_demod;
   %     a_demod=（2+1/snr）^（-0.