close all
clear all
%Tc=1;
N=31;
gx='100011';
g=19;%二進制轉化為十進制
state=1;
L=310;
ts=10/3100;
td1=[0:ts:3099*ts];%3100dian時間向量
td2=[0:ts:3109*ts];%3110點時間向量
td3=[0:ts:3119*ts];%3120點時間向量
        p1=0.3*(randn(1,L*10)+j*randn(1,L*10));
        p2=0.2*(randn(1,L*10)+j*randn(1,L*10));
        p3=0.1*(randn(1,L*10)+j*randn(1,L*10));
        t1=0;
        t2=1;
        t3=2;
        %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%m序列
[pt]=mgen(g,state,L+t3); %pt是m序列 有兩個碼片的延時
pt=2*pt-1;%雙極性碼
pt1=chouyang(pt,L+t3,10);
 figure(1);
plot(td3,pt1);
axis([0 12 -1.5 1.5]);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
                      %產生信號
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

     d=sign(randn(1,10));%產生10個碼元信號
  
       for i=1:10;
          for a=1:310;
               dd(310*(i-1)+a)=d(i);
          end;
       end;

        s=conv(dd,ones(1,N*10));
        st=s(1:(L+t3)*10).*pt1(1:(L+t3)*10);%m序列
        
        
        
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%產生隨機噪聲
        z=randn(1,L*10)+j*randn(1,L*10);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
        rt=st(1:L*10).*p1+st((t2*10+1):(L+t2)*10).*p2+st((t3*10+1):(L+t3)*10).*p3+z;
        r1=rt.*conj(p1).*pt1(1:L*10);
        r2=rt.*conj(p2).*pt1((t2*10+1):(L+t2)*10);
        r3=rt.*conj(p3).*pt1((t3*10+1):(L+t3)*10);
        r1=reshape(r1,N,L*10/N);y1=sum(r1);
        r2=reshape(r2,N,L*10/N);y2=sum(r2);
        r3=reshape(r3,N,L*10/N);y3=sum(r3);
       y=y1+y2+y3;
       dc=sign(real(y));
       td2=[0:0.001:9.999];
       dc1=chouyang(dc,100,100);
       figure(3);
       subplot(211);
       
       plot(td1,dd);
       axis([0 11 -1.5 1.5])
       subplot(212);
       plot(td2,dc1);
       axis([0 11 -1.5 1.5])