pmweima.asv

伪随机码

clear
t0=0.25;
ts=0.0005;
fc=200;
kf=50;
fs=1/ts;
t=[0:ts:t0];
df=0.25; %频率分辨率
%信号
x=[1 0 1 0 0 1 1];
l1=length(x);
m=8;
y=pcode1(m);
l2=length(y);
code1=zeros(1,l1*l2);
for i=1:l1
     if x(i)==1
      for j=1:l2
       code1(1,(l1-1)*l2+j)=y(1,j);
   end
    else
        for j=1:l2
        code1(1,(l1-1)*l2+j)=-y(1,j);
    end
end
end
for i=1:l1*l2
     if code1(1,i)==1
         u=cos(2*pi*fc*t);
    else
        u=-cos(2*pi*fc*t);
  end
end
figure
plot(abs(fft(u)))
figure
plot(u)