功率谱筛选显著周期得到这样的图，怎么这么丑，和一般的功率谱图不一样啊

zyr1995

楼主~可以贴一下你画功率谱的gs文件吗？我最近也在做这个，但我比较菜。。不会写gs文件。。

jolincai

zyr1995 发表于 2018-3-16 11:26
楼主~可以贴一下你画功率谱的gs文件吗？我最近也在做这个，但我比较菜。。不会写gs文件。。

啊 处理后得到数据了 画图就很简单了啊 类似折线图 只是要画两个分量

anndi

楼主解决问题了么 我的图也没有山峰样的

jolincai

anndi 发表于 2018-5-21 15:38
楼主解决问题了么 我的图也没有山峰样的

用距平做，时间序列最好长一些，时间在200以内的图都不会好看...

jolincai

anndi 发表于 2018-5-21 15:38
楼主解决问题了么 我的图也没有山峰样的

用距平做，时间序列最好长一些，时间在200以内的图都不会好看...

anndi

jolincai 发表于 2018-5-22 17:55
用距平做，时间序列最好长一些，时间在200以内的图都不会好看...

时间序列也只有夏季92天的，看文献也是夏季的，但是图是山峰形状的。
我用李建平老师的程序画的连续的功率谱估计方法，认真看了统计书，编的matlab程序，自己认为没有错，可是画出来的图很奇怪。
这俩天在研究matlab自带函数periodogram，psd，fft。
想请问楼主：
（1）psd和periodogram求出来的是功率谱密度，想要得到功率谱是不是需要功率谱密度乘以频率？看到文献上纵坐标是power *frequency
（2）红噪音检验时，得出的谱求平均，乘以将自相关系数r（1）带入公式一阶马尔科夫的基础上得出的，那么通过matlab算的功率谱密度乘以频率求得谱也可以这么求红噪音么？
（3）文献里有个a priori 99% confidence bounds,a posteriori  99% confidence level.都是通过99%显著性检验的线，是因为根据不同的最大滞后时间长度m求得的么？m的取值区间【n/3  n/10】
谢谢楼主

anndi

jolincai 发表于 2018-5-22 17:55
用距平做，时间序列最好长一些，时间在200以内的图都不会好看...

时间序列也只有夏季92天的，看文献也是夏季的，但是图不是山峰形状的。
我用李建平老师的程序画的连续的功率谱估计方法，认真看了统计书，编的matlab程序，自己认为没有错，可是画出来的图很奇怪。
这俩天在研究matlab自带函数periodogram，psd，fft。
想请问楼主：
（1）psd和periodogram求出来的是功率谱密度，想要得到功率谱是不是需要功率谱密度乘以频率？看到文献上纵坐标是power *frequency
（2）红噪音检验时，得出的谱求平均，乘以将自相关系数r（1）带入公式一阶马尔科夫的基础上得出的，那么通过matlab算的功率谱密度乘以频率求得谱也可以这么求红噪音么？
（3）文献里有个a priori 99% confidence bounds,a posteriori  99% confidence level.都是通过99%显著性检验的线，是因为根据不同的最大滞后时间长度m求得的么？m的取值区间【n/3  n/10】
谢谢楼主

jolincai

anndi 发表于 2018-5-22 18:32
时间序列也只有夏季92天的，看文献也是夏季的，但是图不是山峰形状的。
我用李建平老师的程序画的连续的 ...

时间这么短 为什么不用小波分析

jolincai

anndi 发表于 2018-5-22 18:32
时间序列也只有夏季92天的，看文献也是夏季的，但是图不是山峰形状的。
我用李建平老师的程序画的连续的 ...

时间这么短 为什么不用小波分析

anndi

jolincai 发表于 2018-5-23 09:34
时间这么短 为什么不用小波分析

因为文献里做的功率谱分析，想要看典型年份的情况