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放大系數(shù)
???????%公式（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.