自相关性MATLAB代码

clear;
clc;
x=1:1:36;
s=[0.41
0.36
0.39
0.31
0.36
0.33
0.27
0.20
0.15
0.22
0.24
0.22
0.25
0.21
0.26
0.19
0.25
0.18
0.51
0.45
0.38
0.34
0.41
0.32
0.29
0.29
0.28
0.26
0.27
0.35
0.34
0.24
0.23
0.18
0.16
0.11]';
sc=s-mean(s);
Y=sc
R0=0;
for i=1:36
R0=Y(i)^2/36+R0;
end
R0
for k=1:36
R(k)=0;
for i=k+1:36
R(k)=Y(i)*Y(i-k)/36+R(k);
end
R %自协方差函数R
end
x=R/R0 %自相关系数x
figure
plot(x);title('自相关系数分析图');