混沌时间序列的平均周期的求解
混沌时间序列的平均周期如何求解呢?我用了一下的程序,但是求出的平均周期值等于数据的长度,不知道错误出在哪里了,希望好心人帮忙指教。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) %由下标求出平均周期数据:
同问!我也遇到这个问题,但只有两组数据出现这种情况,请问楼主现在解决了吗? 看不到图,简单分析下,这个程序出现过很多次,在大部分的论坛里大家用的都是这个方法。首先楼主要明白的是代码中没有关于采样频率的转换,因此默认的采样频率应该是1 其次其他的计算平均周期的方法还有很多,这里用的是最大值法,有存在误差的风险,其他方法(比如加权法)你可以在百度文库里搜索“时间序列 平均周期”找到相关论文。
页:
[1]