matlab 功率谱mse,功率谱和功率谱密度

宋老师，您好：

1、关于语音信号，求功率谱和功率谱密度，是不是一样的，没有区别的？

2、看到您书中求取功率谱密度(p28)用的是pwech方法。看到还有用spectrogram函数，用[S,F,T,P]=spectrogram(x,window,noverlap,nfft,fs)或者

[S,F,T,P]=spectrogram(x,window,noverlap,F,fs)。还看到用STFT中的[s, f, t] = stft(x, wlen, h, nfft, fs)。

这些是不是都是求功率谱的？

3、最近看到您在论坛里的一个psd的帖子，有段代码

clc;clear all;

%读入语音信号

[x,fs]=wavread('F:female1.wav');

t=(0:length(x)-1)/16000;

%采样频率

fs=16000;

%绘制时域波形

figure(1);

plot(t,x,'LineWidth',2);

xlabel('时间 t');

ylabel('幅值 A');

%计算功率谱

nfft=2048;

psd=abs(fft(x,nfft));

%绘制功率谱图

figure(2)

plot((0:nfft/2-1)/nfft*fs,psd(1:nfft/2),'LineWidth',2);

xlabel('频率 f');

ylabel('功率谱 PSD');

%计算短时傅里叶变换

%窗函數

h=window('hamming',85);

[tfr,t,f]=tfrstft(x,1:length(x),256,h);

%绘制等高线图

figure (3);

contour(t/fs,(0:127)/256*fs,abs(tfr(1:128,

).^2);

xlabel('时间 t');

ylabel('频率 f');

%三维图

figure (4);

mesh(t/fs,(0:127)/256*fs,abs(tfr(1:128,

).^2);

xlabel('时间 t');

ylabel('频率 f');

zlabel('幅值 A');

试着运行后出现[img=0,1]file:///C:UsersAdministratorAppDataRoamingTencentUsers123273043QQWinTempRichOleUQ2`683%_[O[~839R@XIPLU.jpg[/img]

UQ2`683%_[O[~839R@XIPLU.jpg (67.35 KB, 下载次数: 55)

2014-10-31 10:35 上传

错误，这改如何修改的？

还有，这个里面用psd=abs(fft(x,nfft))求功率谱，怎么要加绝对值的？

谢谢宋老师