rake1.m

Rake接收的相关程序

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])