close?all?
????clc??
????clear?all??
????SNR_DB=[0:1:12];?sum=1000000;??
????data=?randsrc（sum2[0?1]）;??
????[a1b1]=find（data（:1）==0&data（:2）==0）;?
????message（a1）=-1-j;??
????[a2b2]=find（data（:1）==0&data（:2）==1）;?
????message（a2）=-1+j;??
????[a3b3]=find（data（:1）==1&data（:2）==0）;
????message（a3）=1-j;??
????[a4b4]=find（data（:1）==1&data（:2）==1）;
????message（a4）=1+j;?
????scatterplot（message）?
????title（‘A點信號的星座圖‘）;
?
????A=1;?Tb=1;??Eb=A*A*Tb;
????P_signal=Eb/Tb;??
????NO=Eb./（10.^（SNR_DB/10））;
????P_noise=P_signal*NO;
????sigma=sqrt（P_noise）;?
????for?Eb_NO_id=1:length（sigma）?????
????????noise1=sigma（Eb_NO_id）*randn（1sum）;?????
????????noise2=sigma（Eb_NO_id）*randn（1sum）;???
????????receive=message+noise1+noise2*j;????
????????resum=0;?????
????????total=0;?????
????????m1=find（angle（receive）<=pi/2&angle（receive）>0）;??
????????remessage（1m1）=1+j;?
????????redata（m11）=1;???
????????redata（m12）=1;?????
????????m2=?find（?angle（receiv