logpn_digitize.m

来自「3种pn序列的产生代码」· M 代码 · 共 26 行

n=1024;
x0=0.2;
x=zeros(1,n);
x(1,1)=x0;
for i=1:n-1
    x(1,i+1)=4*x(1,i)*(1-x(1,i));
end;
N=length(x);
A=zeros(1);
for i=1:N
    A=A+x(i);
end;
A=A/N;
for i=1:N,
    if x(i)>=A
        y(i)=1;
    else
        y(i)=-1;
    end;
end; 
subplot(2,1,1)
[c,lags]=xcorr(y(1,:),500,'unbiased');
plot(lags,(c));


figure,plot(y);