连续小波变换技术总结.pptVIP

　　例已知一信号f t ＝3sin 100pt ＋2sin 68pt ＋5cos 72pt ，且该信号混有白噪声，对该信号进行连续小波变换。小波函数取db3，尺度为1、1.2、1.4、1.6、…、3。其MATLAB程序如下： t＝0：0.01：1； f＝3*sin 100*pi*t ＋2*sin 68*pi*t ＋5*cos 72*pi*t ＋randn 1，length t ； coefs＝cwt f，［1：0.2：3］，￠db3￠，￠plot￠ ； title ￠对不同的尺度小波变换系数值￠ ； Ylabel ￠尺度￠ ； Xlabel ￠时间￠ ；程序输出结果如图所示。 　 图1.11 　　小波变换的系数用图所示的灰度值图表征，横坐标表示变换系数的系号，纵坐标表示尺度，灰度颜色越深，表示系数的值越大。 绘图原理 1.需要用到的小波工具箱中的三个函数? cwt（），centfrq（），scal2frq（） COEFS cwt S,SCALES,‘wname’ 说明：该函数能实现连续小波变换，其中S为输入信号，SCALES为尺度，wname为小波名称。 FREQ centfrq ‘wname’ 说明：该函数能求出以wname命名的母小波的中心频率。 F scal2frq A,‘wname’,DELTA 说明：该函数能将尺度转换为实际频率，其中A为尺度，wname为小波名称，DELTA为采样周期。? 注：这三个函数还有其它格式，具体可参阅matlab的帮助文档。 2.尺度与频率之间的关系 设a为尺度，fs为采样频率，Fc为小波中心频率，则a对应的实际频率Fa为 Fa＝Fc×fs/a 1 显然，根据采样定理，为使小波尺度图的频率范围为 0,fs/2 ，尺度范围应为 2*Fc,inf ,其中inf表示为无穷大。在实际应用中，只需取尺度足够大即可。 3.尺度序列的确定 由式 1 可以看出，为使转换后的频率序列是一等差序列，尺度序列必须取为以下形式： c/totalscal,...,c/ totalscal-1 , c 2 其中，totalscal是对信号进行小波变换时所用尺度序列的长度 通常需要预先设定好 ，c为一常数。 下面讲讲c的求法。 根据式 1 容易看出，尺度c/totalscal所对应的实际频率应为fs/2，于是可得 c 2×Fc×totalscal 3 将式 3 代入式 2 便得到了所需的尺度序列。 4.时频图的绘制 确定了小波基和尺度后，就可以用cwt求小波系数coefs（系数是复数时要取模），然后用scal2frq将尺度序列转换为实际频率序列f，最后结合时间序列t，用imagesc t,f,abs coefs 便能画出小波时频图。 注意：直接将尺度序列取为等差序列，例如1:1:64，将只能得到正确的尺度－时间－小波系数图，而无法将其转换为频率－时间－小波系数图。这是因为此时的频率间隔不为常数。。 ? 下面给出一实际例子来说明小波时频图的绘制。所取仿真信号是由频率分别为100Hz和200Hz的两个正弦分量所合成的信号。 clear; clc; fs 1024; %采样频率 f1 100; f2 200; t 0:1/fs:1; s sin 2*pi*f1*t +sin 2*pi*f2*t ;??%两个不同频率正弦信号合成的仿真信号 %%%%%%%%%%%%%%%%%小波时频图绘制 wavename cmor3-3; totalscal 256; %尺度序列的长度，即scal的长度 wcf centfrq wavename ; %小波的中心频率 cparam 2*wcf*totalscal; %为得到合适的尺度所求出的参数 a totalscal:-1:1; scal cparam./a; %得到各个尺度，以使转换得到频率序列为等差序列 coefs cwt s,scal,wavename ; %得到小波系数 f scal2frq scal,wavename,1/fs ;? ?%将尺度转换为频率 imagesc t,f,abs coefs ; %绘制色谱图 colorbar; xlabel 时间 t/s ; ylabel 频率 f/Hz ; title 小波时频图 ; 第二章 连续小波变换 2.1连续小波变换及其性质 2.1.1 连续小波基函数 小波，即小区域的波，是一种特殊的长度有限、平均值为零的波形。 小波的可容许条件： 小波特点： （一）“小”。即在时域都具有紧支集或近似紧支集。 （二）正负交替的“波动性”。即直流分量为零。 信号可分解为一系列由同一个母小波函数经平移与尺度伸缩得到的小波函数的叠加。 将小波母函数 进行伸缩和平移，就可以得到函数： 小波函数基，它们是由同一