(求助)：时间序列的功率谱密度PSD图怎么画？

ChaChing · 发表于 2011-4-10 23:25

aliu 发表于 2011-4-10 19:49

登录/注册后可看大图

ChaChing,你有没有发现样本图和你画出来的图的差别有点大呢？？呵呵，，我在好好找原因.....

并非我画的图!? 是Apollo
可能是log scale的差异吧!?

Apollo · 发表于 2011-4-11 09:15

我有一下疑问：
1、分析频率是多少？楼主提供的是一年中每小时每点的数据，难道实际采样间隔为1小时？肯定不是
2、频率能与时间在频谱中同时出现吗？像楼主要的效果图也只能如我贴的波谱图那样分别显示，在频谱图中读不出时间信息。从图形看，楼主画的并非是功率谱密度图，而是时域图，只是坐标单位改了。
3、绘制这样图形的目的是什么？

ChaChing · 发表于 2011-4-11 13:45

本帖最后由 ChaChing 于 2011-4-11 14:00 编辑

真的有点受不了, 试做下 ...
不清楚LZ的PSD指的是power spectrum OR power spectrum density, 所以其他LZ自行修订

load value.txt; plot(value)
nn=length(value); sr=1/3600; DHz = sr/nn; ff=[0:nn/2-1]*DHz;
yy=fft(value); spec=abs(yy(1:nn/2)).*abs(yy(1:nn/2));
figure; loglog(ff,spec); grid on
figure; plot(log10(ff),log10(spec)); grid on

ChaChing · 发表于 2011-4-11 13:53

本帖最后由 ChaChing 于 2011-4-11 13:55 编辑

回复 17 # Apollo 的帖子

...频谱图中读不出时间信息...

週期与频率不是互为倒数!? 阁下忽略了!?

aliu · 发表于 2011-4-11 15:43

回复ChaChing
我要的PSD是power spectrum density
别生气，别生气，，呵呵，，看到你画的图了，这张图应该就是我要的效果图了，我试着改改里面的参数，看效果会是怎样的
我在这一块是个菜鸟，，所以才上咱论坛来求助高手们的....
谢谢ChaChing的热心解答...感动ing

aliu · 发表于 2011-4-11 15:49

回复Apollo
是可以读得出一些信息的，主要是周期信分布息，我贴出的样本图上面就有相应的周期提示啊
那张样本图在N多的文章中都有出现，所做的描述都是psd图，也就是功率谱密度power spectrum density图
呵呵，说实话，没想到你也那么上心了，挺有意思的，感谢....

ChaChing · 发表于 2011-4-11 23:57

回复 20 # aliu 的帖子

没生气! 我只是一时心急怕Apollo误判, so ...
若是psd图, spec=abs(yy(1:nn/2)).*abs(yy(1:nn/2))/DHz; 此时量值是有差的!

aliu · 发表于 2011-4-12 19:46

回复ChaChing
嗯，你这位朋友就是爽快，自己懂的知识和处理过程，直接就拿出来分享，呵呵，，，
谢谢....
当然，也要谢谢参与过本次求助问题的兄弟姐妹们.....
谢谢你们的建议，解答，交流....

simon_ivy · 发表于 2012-4-3 08:23

回复 1 # aliu 的帖子

你的问题解决了么

aliu · 发表于 2012-4-25 15:20

回复 24 # simon_ivy 的帖子

按照ChaChing 贴上来的代码，已经可以画出想看到的图了

irenelin · 发表于 2012-6-19 17:08

好好研究一下

mrqddp · 发表于 2012-8-15 10:39

请问大神是用什么数据做的啊，纵坐标单位是什么

ys71063947 · 发表于 2013-5-2 15:21

ChaChing 发表于 2011-4-11 13:45

登录/注册后可看大图

真的有点受不了, 试做下 ...
不清楚LZ的PSD指的是power spectrum OR power spectrum density, 所以其他LZ自 ...

我是刚接触谱密度这方面内容弱弱问句为什么ff=[0:nn/2-1],而不是ff=[0:nn-1];还有数据是一小时为单位时，横坐标怎么是负的；因为我的有些数据是以年为单位的，那么fs岂不是更小了？？

ChaChing · 发表于 2013-5-2 21:16

ys71063947 发表于 2013-5-2 15:21

登录/注册后可看大图

我是刚接触谱密度这方面内容弱弱问句为什么ff=[0:nn/2-1],而不是ff=[0:nn-1];还有数据是一小时为单位时， ...

ys71063947 · 发表于 2013-5-4 01:32

恩恩谢谢ChaChing 学习了我的上述问题解决了呵~~

账号		自动登录	找回密码
密码			我要加入

[编程技巧] (求助)：时间序列的功率谱密度PSD图怎么画？

点评