|
马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
终于画了出来,看时频图还挺好的。不知是不是蒙对了。附上程序,希望大虾们能看看,发表一下意见!
load 3;
b=X;
b=b';
Fs=2e+6;
Ts=1/Fs;
wpt=wpdec(b,4,'sym4');
tn=leaves(wpt);
l=length(tn);
nfft=256;
TF=zeros(l,nfft/2);
for k=1:l
node=tn(k);
cfs=read(wpt,'data',node);
fftcfs=fft(cfs,nfft);
absfftcfs=abs(fftcfs);
TF(k,:)=absfftcfs(1:nfft/2);
end
TF=TF';
tnew=(1:l)*128*Ts;
fnew=(1:nfft/2)*Fs/nfft;
[T,F]=meshgrid(tnew,fnew);
contour(T,F,abs(TF)) |
|