混沌时间序列的平均周期问题

原来是小飞

使用如下代码求解一时间序列的平均周期，代码为从论坛中找到的。结果如图片所示，可是不知道怎么从图片看出结果，求好心人帮忙解释一下。
代码：

1. clear all;
   
2. data=load('b.txt');
   
3. Y = fft(data);       %快速FFT变换
   
4. N = length(data);    %FFT变换后数据长度
   
5. Y(1) = [];           %去掉Y的第一个数据，它是data所有数据的和
   
6. power = abs(Y(1:N/2)).^2;  %求功率谱
   
7. nyquist = 1/2;
   
8. freq = (1:N/2)/(N/2)*nyquist; %求频率
   
9. subplot(121)
   
10. plot(freq,power); grid on     %绘制功率谱图
    
11. xlabel('频率')
    
12. ylabel('功率')
    
13. title('功率谱图')
    
14. period = 1./freq;                %计算周期
    
15. subplot(122)
    
16. plot(period,power); grid on  %绘制周期－功率谱曲线
    
17. ylabel('功率')
    
18. xlabel('周期')
    
19. title('周期—功率谱图')
    
20. [mp,index] = max(power);       %求最高谱线所对应的下标
    
21. T_mean=period(index);            %由下标求出平均周期

复制代码

运行时有个提示：

1. Warning: Integer operands are required for colon operator when used as index

复制代码

图片：