小波包分析后的功率谱

etet-1234

function xbb2()
N=60;
t=0:0.1:N;
x1=sin(2*pi*t)+randn(1,601)+sin(6*pi*t);
fs=10;
nfft=1024;
xd=(x1-mean(x1))/std(x1,1);
figure(1)
subplot(2,2,1)
plot(t,xd)
Y=abs(fft(xd,nfft));
subplot(2,2,2)
plot((0:nfft/2-1)/nfft*fs,Y(1:nfft/2));
T=wpdec(xd,5,'db4');
y1=wprcoef(T,[5,1]);
subplot(2,2,3)
plot(t,y1);
Y1=abs(fft(y1,nfft));
subplot(2,2,4)
plot((0:nfft/2-1)/nfft*fs,Y1(1:nfft/2));

请大家 看一下我感觉我的小波包分析后的功率谱 感觉有些问题 希望大家指导 我是个新手 最好能给讲解一下
小波包分析后 的功率谱与原来的差别太大 谢谢

etet-1234

自己先回答一下 呵呵 我想了想感觉 是我重够后的图象效果不好的原因 造成的 不知道对不对 请指点

taoboy

Y1=abs(fft(y1,nfft))/1024;

qqvirile

xd=(x1-mean(x1))/std(x1,1);是什么意思？

luzq1982

楼主还在不 ？路过交流下