snr_in_dB=0;
snr=10^(snr_in_dB/10);%snr_in_dB
%signal=randint(1,number);
E=1;
number=2400;
sgma=sqrt(E/snr)/2;
signal=randint(1,2400);

T=1;
M=2;Es=T/2;fc=2/T;
N=10;delta_T=T/(N-1);
t=0:delta_T:T;
carry=sqrt(2*Es/T)*cos(2*pi*fc*t)

for  i=1:number
    if signal(i)==1
        sign(i)=1;
    else
        sign(i)=-1;
    end
end
real_sign=zeros(2400,10);
for i=1:number
real_sign(i,:)=sign(i)*ones(1,size(carry,2)).*carry;%調(diào)制real_sign(i,:)
end

for  i=1:number             %加噪聲
  % N0=gngauss(sgma); 
   sign_send(i,:)=real_sign(i,:)+1*randn(1,size(carry,2));%sgma
end
  
for  i=1:number 
    sign_receive(i)=sign_send(i,:)*carry';%解調(diào)
end
    
for i=1:number
    if sign_receive(i)>0
        re_signal(i)=1;
    else
         re_signal(i)=0;
     end
 end
 
 error=0;
 for i=1:number
     if signal(i)~=re_signal(i);
         error=error+1;
     end
 end
 
 pb=error/number
 ps=pb;