matlab中功率谱密度函数的讨论

905lili

这 是求功率谱密度函数的程序，

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));
    

复制代码

Pxx就是功率谱，f为频率。我现在要求 。怎么实现呢，求高手指导。（只知道Pxx关于f的图形，没有函数关系式）

905lili

问题已解决，

kyu16866

905lili 发表于 2012-4-26 14:35

                               
登录/注册后可看大图

问题已解决，

怎么解决的？给点提示？

粤语残片

同求提示啊

wy145223

是怎么解决的呢？求教。

happy

trapz可以实现积分

沐雨柠檬

其功率谱密度为什么要给信号加窗啊？功率谱密度是解决什么问题的？求教！

hcharlie

沐雨柠檬 发表于 2013-3-28 13:47

                               
登录/注册后可看大图

其功率谱密度为什么要给信号加窗啊？功率谱密度是解决什么问题的？求教！

你是研究振动信号处理的，这个问题对你很重要。
功率谱密度是为解决随机振动问题而建立的方法，在分析随机信号时要用功率谱密度来表示。
一个正弦信号用一个频率和振幅表示就够了；
一个（无穷长的）周期信号由一个周期信号的无限重复构成，它包含有限个正弦分量，取一个或几个整周期信号的离散采样做DFT分析，就可以得出这几个频率分量值。
而一个宽带随机信号，它没有可重复的周期，它的频率分量是连续分布的，所以必须用功率谱密度来描述，意即在各频率上单位频率内含有的功率（例如一个Hz的频率范围内含有多少g平方）。所以功率谱密度是定义在随机信号里，功率谱密度方法适用与随机信号的分析。
随机信号求功率谱密度时也用DFT或FFT方法，先加窗是为了减小功率泄漏，这不是一两句话说得清楚，找点资料看看。