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

dsatttt

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

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

ChaChing

搜过了吗? 已经很多了!

dsatttt

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

ChaChing

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

xiaoice_cui

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

905lili

二种基本的方法，1，自相关函数法，程序如下，

1. clear;
   
2. Fs=1000; %采样频率
   
3. n=0:1/Fs:1;
   
4. 
   
5. %产生含有噪声的序列
   
6. xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
   
7. 
   
8. nfft=1024;
   
9. cxn=xcorr(xn,'unbiased'); %计算序列的自相关函数
   
10. CXk=fft(cxn,nfft);
    
11. Pxx=abs(CXk);
    
12. index=0:round(nfft/2-1);
    
13. k=index*Fs/nfft;
    
14. plot_Pxx=10*log10(Pxx(index+1));
    
15. plot(k,plot_Pxx);
    

复制代码

2.周期图法，程序如下，

1. clear;
   
2. Fs=1000; %采样频率
   
3. n=0:1/Fs:1;
   
4. 
   
5. %产生含有噪声的序列
   
6. xn=cos(2*pi*40*n)+3*cos(2*pi*100*n)+randn(size(n));
   
7. 
   
8. window=boxcar(length(xn)); %矩形窗
   
9. nfft=1024;
   
10. [Pxx,f]=periodogram(xn,window,nfft,Fs); %直接法
    
11. plot(f,10*log10(Pxx));
    

复制代码

xiaoice_cui

回复 6 # 905lili 的帖子

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

图1

图2

图3

无一涯

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

ChaChing

回复 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-s ... -blog-id-18250.html

jifeiyun

junyanmm2008

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

mrqddp

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

pengzc

学习一下

清晨露珠

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

                               
登录/注册后可看大图

回复 6 # 905lili 的帖子

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

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

gaogaola

遇到了跟楼主一样的问题，求解答啊！！