声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 8064|回复: 23

[分形与混沌] (求助)时间序列的平均周期如何计算?

[复制链接]
发表于 2007-9-17 21:10 | 显示全部楼层 |阅读模式

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

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

x
时间序列的平均周期如何计算?我是菜鸟,还请高手赐教!!!谢谢了
回复
分享到:

使用道具 举报

发表于 2007-9-18 10:17 | 显示全部楼层
其实就是把功率谱图的横坐标f  做个1./f变换 然后画图出来对应最高峰位置的横坐标就是 周期P 了就可以了
load sunspot.dat
year = sunspot(:,1);
wolfer = sunspot(:,2);
plot(year,wolfer)
title('Sunspot Data')
Y = fft(wolfer);
N = length(Y);
Y(1) = [];
power = abs(Y(1:N/2)).^2;
nyquist = 1/2;
freq = (1:N/2)/(N/2)*nyquist;
plot(freq,power), grid on
xlabel('cycles/year')
title('Periodogram')
period = 1./freq;
plot(period,power), axis([0 40 0 2e7]), grid on
ylabel('Power')
xlabel('Period(Years/Cycle)')

评分

1

查看全部评分

 楼主| 发表于 2007-9-18 10:45 | 显示全部楼层
谢谢,让我好好看看这个程序!!!不过好像用最小数据量法求最大lyapunov时,周期P的影响不大,如果在不确定周期P的情况下,是不是可以用P=1来进行计算?(陆振波博士论文中的观点)
发表于 2007-9-18 11:23 | 显示全部楼层

回复 #3 柏莱 的帖子

这个和你要求解的数据有关,你可以先试算一下,譬如取P=1,2,3,看看结果相差大不大!

如果差别很大,那肯定不能直接设置P=1了!
发表于 2007-9-18 15:46 | 显示全部楼层

回复 #3 柏莱 的帖子

不要听信别人  自己试过才知道的
 楼主| 发表于 2007-9-19 14:35 | 显示全部楼层
也是的,陆博士在数据选取上选的数据比较多(大概几千个数),但我的数据大概也就几十到几百个,应该会有不同!谢谢大家的提醒!
 楼主| 发表于 2007-9-19 15:37 | 显示全部楼层

你给的这个程序,我有点地方不太明白意思,是不是我照葫芦画瓢就可以呢?还是有些参数需要根据实际情况而改变呢?
 楼主| 发表于 2007-9-19 16:12 | 显示全部楼层
还想问一个问题,如果求出的P是不是一个整数,那是应该向0取整呢,还是应该向正无穷取整?
发表于 2007-9-19 16:16 | 显示全部楼层
什么叫向正无穷取整啊?
P四舍五入,取整法  最多相差1   对结果没有影响的
 楼主| 发表于 2007-9-20 15:39 | 显示全部楼层
你给的这个程序,我有点地方不太明白意思,是不是我照葫芦画瓢就可以呢?还是有些参数需要根据实际情况而改变呢?
发表于 2007-9-20 16:29 | 显示全部楼层
这个程序你可以参考一下matlab中fft帮助部分的内容,应该可以找到你需要的解释的~
发表于 2007-10-19 14:28 | 显示全部楼层
这种方法求出的平均频率只是真实频率的分数倍数吧,把真实频率看作1,求出的平均频率可能是零点几,这是什么道理?
发表于 2007-10-19 20:43 | 显示全部楼层

回复 #12 JulianChin 的帖子

你这么一说还真是有问题啊,过几天有时间了再把它搞清楚!你要是有结果也请一起分享哦!
发表于 2007-10-19 22:08 | 显示全部楼层

回复 #13 octopussheng 的帖子

如果求的是真实周期的话,想想看,一个中心频率100k的信号,求出的p几乎等于0,那不如干脆设成1算了

我看原意还是求分频吧
发表于 2007-10-20 21:33 | 显示全部楼层

回复 #14 JulianChin 的帖子

昨晚又翻书了,发现用分频更解释的通:@L :@L

我再学习学习!
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-12-23 14:30 , Processed in 0.066261 second(s), 20 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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