dsatttt 发表于 2009-2-26 17:22

怎么用MATLAB绘制离散时间序列的功率谱密度函数？

已知一个离散时间序列，画成以时间为横坐标的曲线如图一：

怎么通过这个序列画出它的功率谱密度曲线（横坐标为频率），如图二所示？谢谢！用MATLAB

ChaChing 发表于 2009-2-26 18:40

回复 楼主 dsatttt 的帖子

搜过了吗? 已经很多了!

dsatttt 发表于 2009-2-26 19:26

我看了一些，我先给序列做了自相关分析xcorr()
然后做傅里叶变换，然后再求幅值，绘制曲线，对吗？

ChaChing 发表于 2009-2-26 22:19

个人印象中有现成的函数!

xiaoice_cui 发表于 2012-4-27 19:40

楼主，这个应该怎么做呀？我现在在做这个，不懂呢

905lili 发表于 2012-4-27 20:15

二种基本的方法，1，自相关函数法，程序如下，clear;
Fs=1000; %采样频率
n=0:1/Fs:1;

%产生含有噪声的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));

nfft=1024;
cxn=xcorr(xn,'unbiased'); %计算序列的自相关函数
CXk=fft(cxn,nfft);
Pxx=abs(CXk);
index=0:round(nfft/2-1);
k=index*Fs/nfft;
plot_Pxx=10*log10(Pxx(index+1));
plot(k,plot_Pxx);
2.周期图法，程序如下，clear;
Fs=1000; %采样频率
n=0:1/Fs:1;

%产生含有噪声的序列
xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));

window=boxcar(length(xn)); %矩形窗
nfft=1024;
=periodogram(xn,window,nfft,Fs); %直接法
plot(f,10*log10(Pxx));

xiaoice_cui 发表于 2012-4-28 01:56

回复 6 # 905lili 的帖子

谢谢解答！不过我还是不大明白。就是我的数据是离散时间数据的，1年365天的数据，1天1个，那这样的话，我的参数应该怎样设置呀？ Fs=1/3600？nfft=365 ?
如果是这样设置，并且数据没有做标准化处理的话，出来的图是图1；
如果用zscore函数对数据先做标准化处理，就是图2；
可是，这都不太对。我想要的效果是大概图3这样的。为什么图3的纵坐标都是小于1的？横坐标怎么设置在一定的范围内的？

无一涯 发表于 2012-4-28 16:09

学习一下，增加体能+知识

ChaChing 发表于 2012-4-29 14:07

回复 5 # xiaoice_cui 的帖子

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

jifeiyun 发表于 2012-5-3 14:58

{:{13}:}

junyanmm2008 发表于 2012-6-11 21:24

我也在做这个，求30个离散点的功率谱，谁能给个详细点的解释啊

mrqddp 发表于 2012-9-3 17:04

我也有这方面的困扰，楼主的问题解决了嘛

pengzc 发表于 2013-1-4 21:32

{:{39}:}学习一下

清晨露珠 发表于 2013-1-24 15:44

xiaoice_cui 发表于 2012-4-28 01:56 static/image/common/back.gif
回复 6 # 905lili 的帖子

谢谢解答！不过我还是不大明白。就是我的数据是离散时间数据的，1年365天的数据， ...

您好！怎么画功率谱模型函数啊？是否可以指点下！或者推荐下着方面的书籍！谢谢您

gaogaola 发表于 2013-11-21 14:56

遇到了跟楼主一样的问题，求解答啊！！
页: [1]
查看完整版本: 怎么用MATLAB绘制离散时间序列的功率谱密度函数？