混沌时间序列的平均周期问题
本帖最后由 原来是小飞 于 2013-11-12 16:28 编辑使用如下代码求解一时间序列的平均周期，代码为从论坛中找到的。结果如图片所示，可是不知道怎么从图片看出结果，求好心人帮忙解释一下。
代码：clear all;
data=load('b.txt');
Y = fft(data); %快速FFT变换
N = length(data); %FFT变换后数据长度
Y(1) = []; %去掉Y的第一个数据，它是data所有数据的和
power = abs(Y(1:N/2)).^2;%求功率谱
nyquist = 1/2;
freq = (1:N/2)/(N/2)*nyquist; %求频率
subplot(121)
plot(freq,power); grid on %绘制功率谱图
xlabel('频率')
ylabel('功率')
title('功率谱图')
period = 1./freq; %计算周期
subplot(122)
plot(period,power); grid on%绘制周期－功率谱曲线
ylabel('功率')
xlabel('周期')
title('周期—功率谱图')
= max(power); %求最高谱线所对应的下标
T_mean=period(index); %由下标求出平均周期运行时有个提示：Warning: Integer operands are required for colon operator when used as index图片：
页:
[1]