这样求的平均周期和时间延迟对吗？

wystar

我求解方程之后提取了一组数据，然后用FFT方法求平均周期，用互信息法求时间延迟，但是感觉平均周期和时间延迟太小了点，用ＧＰ算法计算关联维也不收敛，请大家看看。

1. %function T_mean=period_mean_fft(data)
   
2. %该函数使用快速傅里叶变换FFT计算序列平均周期
   
3. %data：时间序列
   
4. %T_mean：返回快速傅里叶变换FFT计算出的序列平均周期
   
5. clc
   
6. clear
   
7. data=load('x.txt');%导入数据
   
8. 
   
9. Y = fft(data);       %快速FFT变换
   
10. N = length(data)   %FFT变换后数据长度
    
11. Y(1) = [];           %去掉Y的第一个数据，它是data所有数据的和
    
12. power = abs(Y(1:N/2)).^2 ;%求功率谱?
    
13. 
    
14. size(power)
    
15. nyquist = 1/2;
    
16. freq = (1:N/2)/(N/2)*nyquist;%求频率?
    
17. size(freq)
    
18. subplot(211)
    
19. plot(freq,power); grid on     %绘制功率谱图
    
20. xlabel('频率')
    
21. ylabel('功率')
    
22. title('功率谱图')
    
23. period = 1./freq;                %计算周期
    
24. subplot(212)
    
25. plot(period,power); grid on  %绘制周期－功率谱曲线
    
26. ylabel('功率')
    
27. xlabel('周期')
    
28. title('周期—功率谱图')
    
29. [mp,index] = max(power);       %求最高谱线所对应的下标
    
30. T_mean=period(index)            %由下标求出平均周期

复制代码

wystar

计算得到的平均周期为3.58，时间延迟等于1。

wystar

互信息法是gghhjj前辈分享的，由于现在没有找到原帖，所以再分享一份。

lzmhncs

我用FFT方法求平均周期，怎么求出来的总是data的长度