混沌时间序列的平均周期的求解

原来是小飞

混沌时间序列的平均周期如何求解呢？我用了一下的程序，但是求出的平均周期值等于数据的长度，不知道错误出在哪里了，希望好心人帮忙指教。

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

复制代码

数据： b.txt (58.6 KB, 下载次数: 0)

2013-11-18 16:01 上传

点击文件名下载附件
数据

kingsir

同问！我也遇到这个问题，但只有两组数据出现这种情况，请问楼主现在解决了吗？

窃以为洒脱

看不到图，简单分析下，这个程序出现过很多次，在大部分的论坛里大家用的都是这个方法。首先楼主要明白的是代码中没有关于采样频率的转换，因此默认的采样频率应该是1

窃以为洒脱

其次其他的计算平均周期的方法还有很多，这里用的是最大值法，有存在误差的风险，其他方法（比如加权法）你可以在百度文库里搜索“时间序列 平均周期”找到相关论文。