905lili 发表于 2012-4-26 09:59

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

这 是求功率谱密度函数的程序，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));

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

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

问题已解决，{:4_63:}

kyu16866 发表于 2012-11-8 19:57

905lili 发表于 2012-4-26 14:35 static/image/common/back.gif
问题已解决，

怎么解决的？给点提示？

粤语残片 发表于 2012-11-11 15:15

同求提示啊

wy145223 发表于 2012-12-2 21:48

是怎么解决的呢？求教。

happy 发表于 2012-12-3 09:50

trapz可以实现积分

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

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

hcharlie 发表于 2013-3-31 15:47

沐雨柠檬 发表于 2013-3-28 13:47 static/image/common/back.gif
其功率谱密度为什么要给信号加窗啊？功率谱密度是解决什么问题的？求教！

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