时频分析 Gabor 变换 Matlab程序

%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（