matlab三相功率测量不对,matlab中welch法出处理信号的功率谱怎么看着不对呢？

我用welch法对一个震动信号进行分析，我看到其他人做的功率谱分析时候，功率谱(图像的y轴)都是负的，而且波动很大。而我做的结果却不是这样的。不知道原因出到哪里。跪求众大神，帮忙。谢谢了。求大家帮我找找问题啊。我已经上传了结果和原始数据还有程序。

下边是我的程序：

C=textread('zhendong.txt','','headerlines',200);

Fs=1024;

nfft=512;

window=boxcar(100); %矩形窗

window1=hamming(100); %海明窗

window2=blackman(100); %blackman窗

noverlap=60; %数据无重叠

range='half'; %频率间隔为[0 Fs/2]，只计算一半的频率

[Pxx,f]=pwelch(C(:,2),window,noverlap,nfft,Fs,range);

[Pxx1,f]=pwelch(C(:,2),window1,noverlap,nfft,Fs,range);

[Pxx2,f]=pwelch(C(:,2),window2,noverlap,nfft,Fs,range);

plot_Pxx=10*log10(Pxx);

plot_Pxx1=10*log10(Pxx1);

plot_Pxx2=10*log10(Pxx2);

subplot(2,2,1);

plot(f,plot_Pxx);

title('矩形窗');

xlabel('频率/Hz');ylabel('功率谱/dB');

subplot(2,2,2);

plot(f,plot_Pxx1);

title('海明窗')

xlabel('频率/Hz');ylabel('功率谱/dB');

subplot(2,2,3);

plot(f,plot_Pxx2);

title('blackman窗')

xlabel('频率/Hz');ylabel('功率谱/dB');

subplot(2,2,4);

plot(C(:,1),C(:,2));

title('原始信号')

xlabel('时间t');ylabel('复制');

we.png

(38.46 KB, 下载次数: 0)

2014-9-16 22:34 上传

结果图

2014-9-16 22:36 上传

点击文件名下载附件

124.19 KB, 下载次数: 10

原始数据