利用insfreq、tfrstft、winger－ville方法进行时频分析出现的问题讨论
下面是利用光学多普勒效应测量的喇叭振动信号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=;%1对应奈奎斯特频率
a=;
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);
=cart2pol(x,y);
figure(6)
polar(TH,R);
=min(R);
=instfreq(sign2);
figure(7)
plot(f)
页:
[1]