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

aliu 发表于 2011-4-10 19:49 static/image/common/back.gif
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=*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 static/image/common/back.gif
真的有点受不了, 试做下 ...
不清楚LZ的PSD指的是power spectrum OR power spectrum density, 所以其他LZ自 ...

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

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

ys71063947 发表于 2013-5-2 15:21 static/image/common/back.gif
我是刚接触谱密度这方面内容 弱弱问句为什么ff=,而不是ff=;还有数据是一小时为单位时， ...

ref: from http://home.chinavib.com/blog-63979-18250.html
5.[原创]关于FFT的频谱对应关系 http://www.chinavib.com/thread-23614-1-1.html
Matlab编程实现FFT实践及频谱分析 http://www.chinavib.com/thread-51332-1-1.html
FFT方法和FFT的应用注意点 http://www.chinavib.com/thread-77716-1-1.html
[原创]关于FFT的结果为什么要除以N http://www.chinavib.com/thread-23665-1-1.html
频谱与功率谱有什么区别? http://www.chinavib.com/thread-8098-1-1.html
信号经过FFT后，该怎样计算幅值和相位 http://www.chinavib.com/thread-53683-1-1.html
matlab中关于FFT的使用(理解频率分辨率、补零问题) http://www.chinavib.com/thread-98535-1-1.html
如何用matlab求PSD曲线 http://www.chinavib.com/thread-72741-1-1.html

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

恩恩 谢谢ChaChing 学习了 我的上述问题解决了呵~~
页: 1 [2]
查看完整版本: (求助)：时间序列的功率谱密度PSD图怎么画？