声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 2709|回复: 8

[综合讨论] 如何用matlab做非等间隔信号的功率谱密度变换

[复制链接]
发表于 2013-4-25 11:16 | 显示全部楼层 |阅读模式

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

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

x
我用其他软件算一些量的时间历程变化,发现保存的数据时间间隔是非等距的。要对这些数据进行频谱分析,得出其功率谱密度。本人之前没接触过这些,不知怎么编程处理,请高手帮忙,先谢谢了

下面是部分数据,第一列是时间,第二列是相应信号的值:
1.00000e-04     5.87636e-05     
2.00000e-04     4.27140e-05     
4.00000e-04     2.55866e-03     
5.00000e-04     7.92671e-04     
6.00000e-04     -3.60755e-04   
7.00000e-04     -6.35785e-04   
8.00000e-04     -3.65362e-04   
9.00000e-04     2.64967e-04     
1.00000e-03     -3.01036e-04   
1.10000e-03     4.90043e-04     
1.20000e-03     -2.48811e-04   
1.30000e-03     -1.17294e-04   
1.40000e-03     -3.73790e-04   
1.50000e-03     2.59957e-04     
1.60000e-03     -9.52361e-05   
1.70000e-03     7.74584e-04     
1.80000e-03     4.34910e-04     
1.90000e-03     1.21932e-03     
2.00000e-03     7.28303e-04     
2.10000e-03     3.80541e-04     
2.20000e-03     1.10776e-03     
2.30000e-03     2.50722e-04     
2.40000e-03     4.86838e-04     
2.50000e-03     -4.34337e-04   
2.80000e-03     1.80804e-03
回复
分享到:

使用道具 举报

发表于 2013-4-25 21:06 | 显示全部楼层
个人水平有限, 只会将数据插值成等距后再fft

补充内容 (2013-5-2 23:28):
有真正的信号处理高手(impulse)赞成, 心中踏实些!

点评

赞成: 5.0
赞成: 5
  发表于 2013-5-2 22:00
 楼主| 发表于 2013-4-26 15:47 | 显示全部楼层

谢谢。但怎么插值啊?我的数据有的是等间隔的,但另一部分的间隔是不等的,象第二行和三行:
有的1.00000e-04     5.87636e-05     
2.00000e-04     4.27140e-05     
4.00000e-04     2.55866e-03     
以及最后两行:
2.50000e-03     -4.34337e-04   
2.80000e-03     1.80804e-03

我因为要对结果进行后处理才接触这方面的,可搞了很长时间,也没头绪,请高手帮忙写了程序,谢了
发表于 2013-4-26 21:23 | 显示全部楼层
怎么插值啊?
help一下interp1, 选个合适的method
 楼主| 发表于 2013-4-27 11:36 | 显示全部楼层
ChaChing 发表于 2013-4-26 21:23
help一下interp1, 选个合适的method

谢谢指点,我试一下。
另外,我想问下,做FFT变换时,时间是怎么转换成频率的?我的数据时间很小,倒数关系的话,频率是不是很大?那么,我要不不要在数据插值成等距后,先对时间(第一列)进行无量纲化后,再FFT啊?
发表于 2013-4-29 09:17 | 显示全部楼层
guohf 发表于 2013-4-27 11:36
谢谢指点,我试一下。
另外,我想问下,做FFT变换时,时间是怎么转换成频率的?我的数据时间很小,倒数关 ...

先对时间进行无量纲化, 对信号结果好像不会有影响吧! 不是吗?
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-2 17:14 | 显示全部楼层
ChaChing 发表于 2013-4-29 09:17
先对时间进行无量纲化, 对信号结果好像不会有影响吧! 不是吗?
ref: from http://home.chinavib.com/blog ...

谢谢了,我先看看,有问题再请教
发表于 2013-5-2 21:29 | 显示全部楼层

看看先学习并试下, 有问题再共同讨论!

本帖最后由 ChaChing 于 2013-5-2 21:31 编辑
guohf 发表于 2013-5-2 17:14
谢谢了,我先看看,有问题再请教

看看先学习并试下, 有问题再共同讨论!
或有结果后分享成果
发表于 2014-3-31 10:49 | 显示全部楼层
学习,支持
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

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

GMT+8, 2024-11-13 10:56 , Processed in 0.083937 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

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