1. 问题描述
信号处理技术中,对于非平稳信号(频域特征随着时间变化的信号),简单的傅立叶分析无法展示信号频率特征随着时间的变化。故往往采用小波分析技术。由于超声导波的频散特性,研究中往往采用小波分析,参考文献中的示例如下:
本文预期实现文献中的小波分析技术。
2. 技术背景
MATLAB小波分析函数。
3. 解决方案
已知:
fs,信号的采样频率 Hz
s,信号时程序列
wavename = ‘Morl’; % 小波类型设置
totalscal = 2048; % 尺度序列的长度
wcf = centfrq(wavename); % 小波的中心频率
cparam = 2*wcf*totalscal; % 为得到合适的尺度所求出的参数
a = totalscal:-1:0.2;
scal = cparam./a; % 得到各个尺度,以使转换得到频率序列为等差序列
coefs=cwt(s,scal,wavename); % 得到小波系数
f=scal2frq(scal,wavename,1/fs); % 将尺度转换为频率
imagesc(t2*1e6,f/100,abs(coefs)); % 绘制色谱图
3.1 分辨率控制
连续小波变换的基本函数为cwt,变换得到小波系数。我们预期的图像纵坐标为频率&#