混沌信号的自相关

%chao signal autocorrelation
%
x(1)=0.47;
for i = 1:9999
   x(i+1)=4*x(i)^3-3*x(i);
end

y=x(2:10000);
y(10000)=4*x(10000)^3-3*x(10000);

plot(x,y,':')
z=xcorr(x,'unbiased');
%计算x的自相关性。'unbiased'选项表示使用无偏自相关性估计
figure 
plot(z)

代码来自香港理工大学——刘重明教授Simulation Programs (polyu.edu.hk)