地震波fft变换、小波变换对比分析matlab运行文件

利用小波变换、fft实现地震波的时域、频域特征的对比分析

fs=10000;
%采样频率
t=（DUT20180413215506040（:1））;
s=（DUT20180413215506040（:2））;
%两不同频率正弦信号合成的信号
%%%%小波时频图绘制%%%%
wavename=‘cmor3-3‘;
totalscal=256;
%尺度序列的长度，即scal的长度
wcf=centfrq（wavename）;
%小波的中心频率
cparam=2*wcf*totalscal;
%为得到合适的尺度所求出的参数
a=totalscal:-1:1;
scal=cparam./a;
%得到各个尺度，以使转换得到频率序列为等差序列
coefs=cwt（sscalwavename）;
%得到小波系数
f=scal2frq（scalwavename1/fs）;
%将尺度转换为频率
figure（1）plot（ts）;
xlabel（‘时间/ms‘）;
ylabel（‘幅值‘）;
title（‘时域波形图‘）;
%绘制原始信号图
figure（2）;
%绘制信号的傅里叶变换的频谱图
N=length（t）;
y=fft（sN）;
amg=abs（y）;
plot（（1:N/2）*fs/Namg（1:N/2）*2/N）;
xlabel（‘频率 f/Hz‘）;
ylabel（‘幅值‘）;
title（‘fft频域图‘）;
figure（3）;
c=cwt（sawavename‘plot‘）;
%直接使用尺度
figure（4）imagesc（tfabs（coefs））;
%绘制色谱图
colorbar;
xlabel（‘时间 t/ms‘）;
ylabel（‘频率 f/Hz‘）;
title（‘小波时频图‘）;
figure（5）;
mesh（tfabs（coefs））; 
axis tight;    
colorbar;
xlabel（‘时间 t/ms‘）;
ylabel（‘频率 f/Hz‘）;     
title（[‘小波时频幅值图‘‘（‘num2str（wavename）‘）‘]）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件         431  2018-03-29 11:35  cwt2.m
     文件        1109  2018-06-23 10:49  cwt1.m