声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1576|回复: 0

[分形与混沌] 混沌时间序列的平均周期问题

[复制链接]
发表于 2013-11-12 16:25 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
本帖最后由 原来是小飞 于 2013-11-12 16:28 编辑

使用如下代码求解一时间序列的平均周期,代码为从论坛中找到的。结果如图片所示,可是不知道怎么从图片看出结果,求好心人帮忙解释一下。
代码:
  1. clear all;
  2. data=load('b.txt');
  3. Y = fft(data);       %快速FFT变换
  4. N = length(data);    %FFT变换后数据长度
  5. Y(1) = [];           %去掉Y的第一个数据,它是data所有数据的和
  6. power = abs(Y(1:N/2)).^2;  %求功率谱
  7. nyquist = 1/2;
  8. freq = (1:N/2)/(N/2)*nyquist; %求频率
  9. subplot(121)
  10. plot(freq,power); grid on     %绘制功率谱图
  11. xlabel('频率')
  12. ylabel('功率')
  13. title('功率谱图')
  14. period = 1./freq;                %计算周期
  15. subplot(122)
  16. plot(period,power); grid on  %绘制周期-功率谱曲线
  17. ylabel('功率')
  18. xlabel('周期')
  19. title('周期—功率谱图')
  20. [mp,index] = max(power);       %求最高谱线所对应的下标
  21. T_mean=period(index);            %由下标求出平均周期
复制代码
运行时有个提示:
  1. Warning: Integer operands are required for colon operator when used as index
复制代码
图片: b.bmp



回复
分享到:

使用道具 举报

您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-12-22 18:00 , Processed in 0.068937 second(s), 22 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表