原来是小飞 发表于 2013-11-18 16:01

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

混沌时间序列的平均周期如何求解呢?我用了一下的程序,但是求出的平均周期值等于数据的长度,不知道错误出在哪里了,希望好心人帮忙指教。data=load('b.txt');
Y = fft(data);   %快速FFT变换
N = length(Y);    %FFT变换后数据长度
Y(1) = [ ];         %去掉Y的第一个数据,它是data所有数据的和
power = abs(Y(1:N/2)).^2;%求功率谱
nyquist = 1/2;
freq = (1:N/2)/(N/2)*nyquist; %求频率
period = 1./freq;                %计算周期
plot(period,power); grid on%绘制周期-功率谱曲线
xlabel('周期')
ylabel('功率')
title('周期—功率谱图')
= max(power);       %求最高谱线所对应的下标
T_mean=period(index)      %由下标求出平均周期数据:

kingsir 发表于 2014-7-7 15:29

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

窃以为洒脱 发表于 2014-8-30 11:02

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

窃以为洒脱 发表于 2014-8-30 11:04

其次其他的计算平均周期的方法还有很多,这里用的是最大值法,有存在误差的风险,其他方法(比如加权法)你可以在百度文库里搜索“时间序列 平均周期”找到相关论文。
页: [1]
查看完整版本: 混沌时间序列的平均周期的求解