利用insfreq、tfrstft、winger－ville方法进行时频分析出现的问题讨论

wdtabcd

      下面是利用光学多普勒效应测量的喇叭振动信号

clear all
cd E:\data\6.21;
load w6212c000.dat;
sign=w6212c000;
figure(1)
plot(sign);
fsign=fft(sign);
figure(2)
plot(abs(fsign));
Fs=2.5e6;
f=[0 0.05 0.05 1];%1对应奈奎斯特频率
a=[1 1 0 0];
b=firls(100,f,a);
wvtool(b);
sign1=filtfilt(b,1,sign);
meansign=mean(sign1);
sign1=sign1-meansign;
figure(4)
plot(sign1);
cd D:\MATLAB6p5p1\toolbox\时频分析工具箱;
sign2=hilbert(sign1); %解析信号的产生方法
figure(5)
plot(sign2);
x=real(sign2);
y=imag(sign2);
[TH,R]=cart2pol(x,y);
figure(6)
polar(TH,R);
[Y,I]=min(R);
[f,t]=instfreq(sign2);
figure(7)
plot(f)