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