小波变换后功率谱问题
我是做轴承故障诊断的，数据是正确的，但是功率谱却是这样的，几乎没有。是程序错了吗？load 279
fs=5000;
N=2048;
t=(0:1/fs:(N-1)/fs);
x=X279_FE_time ;
data=x(1:N,1);
=wavedec(data,5,'db10');
ca5=appcoef(c,l,'db10',5);
cd5=appcoef(c,l,'db10',5);
cd4=appcoef(c,l,'db10',4);
cd3=appcoef(c,l,'db10',3);
cd2=appcoef(c,l,'db10',2);
cd1=appcoef(c,l,'db10',1);
d1=wrcoef('d',c,l,'db10',1);
d2=wrcoef('d',c,l,'db10',2);
d3=wrcoef('d',c,l,'db10',3);
d4=wrcoef('d',c,l,'db10',4);
d5=wrcoef('d',c,l,'db10',5);
a5=wrcoef('a',c,l,'db10',5);
figure(1);
subplot(321);plot(t,a5);title('a5');
subplot(322);plot(t,d5);title('d5');
subplot(323);plot(t,d4);title('d4');
subplot(324);plot(t,d3);title('d3');
subplot(325);plot(t,d2);title('d2');
subplot(326);plot(t,d1);title('d1');
y=hilbert(cd1);
y1=imag(y).^2;
w=sqrt(cd1.^2+y1);
y2=abs(w);
w1=abs(w-mean(w));
psd=abs(fft(y2,N)).^2./N;
figure(2);
plot((0:N/2-1)/N*fs,w1(1:N/2));
figure(3);
plot((0:N/2-1)/N*fs,psd(1:N/2));
包络图如下；
我做的功率谱如下：坐标我不想换成对数坐标，为什么出现的功率谱图是这样的呢？
页:
[1]