声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

12
返回列表 发新帖
楼主: aliu

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

  [复制链接]
发表于 2011-4-10 23:25 | 显示全部楼层

并非我画的图!? 是Apollo
可能是log scale的差异吧!?
回复 支持 反对
分享到:

使用道具 举报

发表于 2011-4-11 09:15 | 显示全部楼层
我有一下疑问:
1、分析频率是多少?楼主提供的是一年中每小时每点的数据,难道实际采样间隔为1小时?肯定不是
2、频率能与时间在频谱中同时出现吗?像楼主要的效果图也只能如我贴的波谱图那样分别显示,在频谱图中读不出时间信息。从图形看,楼主画的并非是功率谱密度图,而是时域图,只是坐标单位改了。
3、绘制这样图形的目的是什么?
发表于 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
zzz.jpg
发表于 2011-4-11 13:53 | 显示全部楼层
本帖最后由 ChaChing 于 2011-4-11 13:55 编辑

回复 17 # Apollo 的帖子
...频谱图中读不出时间信息...
週期与频率不是互为倒数!? 阁下忽略了!?
 楼主| 发表于 2011-4-11 15:43 | 显示全部楼层
回复ChaChing
我要的PSD是power spectrum density
别生气,别生气,,呵呵,,看到你画的图了,这张图应该就是我要的效果图了,我试着改改里面的参数,看效果会是怎样的
我在这一块是个菜鸟,,所以才上咱论坛来求助高手们的....
谢谢ChaChing的热心解答...感动ing
 楼主| 发表于 2011-4-11 15:49 | 显示全部楼层
回复Apollo
是可以读得出一些信息的,主要是周期信分布息,我贴出的样本图上面就有相应的周期提示啊
那张样本图在N多的文章中都有出现,所做的描述都是psd图,也就是功率谱密度power spectrum density图
呵呵,说实话,没想到你也那么上心了,挺有意思的,感谢....
发表于 2011-4-11 23:57 | 显示全部楼层
回复 20 # aliu 的帖子

没生气! 我只是一时心急怕Apollo误判, so ...
若是psd图, spec=abs(yy(1:nn/2)).*abs(yy(1:nn/2))/DHz; 此时量值是有差的!
 楼主| 发表于 2011-4-12 19:46 | 显示全部楼层
回复ChaChing
嗯,你这位朋友就是爽快,自己懂的知识和处理过程,直接就拿出来分享,呵呵,,,
谢谢....
当然,也要谢谢参与过本次求助问题的兄弟姐妹们.....
谢谢你们的建议,解答,交流....
发表于 2012-4-3 08:23 | 显示全部楼层
回复 1 # aliu 的帖子

你的问题解决了么

点评

18F?  发表于 2012-4-5 08:39
 楼主| 发表于 2012-4-25 15:20 | 显示全部楼层
回复 24 # simon_ivy 的帖子

按照ChaChing  贴上来的代码,已经可以画出想看到的图了
发表于 2012-6-19 17:08 | 显示全部楼层
好好研究一下
发表于 2012-8-15 10:39 | 显示全部楼层
请问大神是用什么数据做的啊,纵坐标单位是什么
发表于 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岂不是更小了??
发表于 2013-5-2 21:16 | 显示全部楼层
ys71063947 发表于 2013-5-2 15:21
我是刚接触谱密度这方面内容 弱弱问句为什么ff=[0:nn/2-1],而不是ff=[0:nn-1];还有数据是一小时为单位时, ...

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
发表于 2013-5-4 01:32 | 显示全部楼层
恩恩 谢谢ChaChing 学习了 我的上述问题解决了呵~~
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-14 14:02 , Processed in 0.086140 second(s), 24 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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