qpsk的matlab实现

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