%Appl11_5
close?all???%關(guān)閉所有圖形窗口
addpath（‘tftoolbox‘）;???%添加時(shí)頻分析工具箱所在的路徑
load?hns.dat??%加載地震波數(shù)據(jù)
x=hns;
clear?hns;??
dt=0.02;???%初始采樣間隔
y=decimate（x2510）;????%對信號按原采樣率的1/25轉(zhuǎn)換轉(zhuǎn)換前采用10階Chebyshev低通濾波器進(jìn)行濾波
clear?x?hns;???%清除x和hns變量
sampinterval=dt*25;
N=length（y）;
Nf=300;???%頻率軸的數(shù)據(jù)個(gè)數(shù)
[tfr1]=tfrgabor（yNfNfgausswin（301））;
T1=[0:Nf-1]/Nf*N*sampinterval;F1=[0:N-1]./N/sampinterval;
subplot（221）pcolor（T1F1（1:N/2/10）abs（[tfr1（1:N/2/10:）]））
shading?interp?%將圖像進(jìn)行平滑
title（‘Gabor變換‘）
xlabel（‘時(shí)間/s‘）ylabel（‘頻率/Hz‘）;
colorbar??%加上色標(biāo)
clear?tfr1?T1?F1;??%清除tfr1T1F1變量
[tfr2T2F2]=tfrwv（y1:length（y）Nf）;??%Wigner-Ville分布
T2=T2*sampinterval;F2=[0:Nf-1]/Nf;??%縱橫軸標(biāo)記
subplot（222）pcolor（T2F2（1:Nf/10）abs（[tfr2（