小波变换与小波滤波

1、 汇报人：汇报人：*指导老师：指导老师：* 汇报日期：汇报日期：2014.4.212基本思想： 将信号分解成一系列不同频率的连续正弦波的叠加。缺陷：丢掉了时间信息，无法根据变换结果判断一个特定的信号是在什么时候发生的。傅立叶变换傅立叶变换FT变换适于分析平稳信号。实际中大多数信号含有大量的非平稳信号，例如：突变，奇异，事件的起始与终止等情况。这些情况反映了信号的重要特征，是分析的对象。例如下图：典型的地震信号。3典型的地震记录 4实际采集的地震信号实际采集的地震信号它们的频域特性都随时间而变化。分析它需要提取某一时间段的频域信息或某一频率段所对应的时间信息如何完成只分析数据中的一小部分?5基本

2、思想： 给信号加一个小窗，主要集中在对小窗内的信号进行变换，因此反映了信号的局部特征。6缺陷： 其窗函数的大小和形状均与时间和频率无关，保持固定不变，对于分析时变信号不利。（高频信号持续时间短，低频长。我们希望对于高频采用小的时间窗，低频使用大时间窗进行分析。）STFT无能为力了！ 不能构成不能构成正交基正交基，给数值计算带来不便，给数值计算带来不便。 7登场原因： （1）继承和发展了STFT的局部化思想。 （2）克服了窗口大小不随频率变化、缺乏离散正交基的缺点。8若一物体可用颜色和大小表示，我们称颜色和大小为特征基，构成此物体特征描述空间。大小和颜色是互不相干的2种描述，我们称其为正交。 同

3、时若这些基能够完全表示所有物体，我们称其为完备特征基。 因为特征基表现了物体特征，因而可以用更简洁的描述表示物体。 91984年法国的年轻的地球物理学家Jean Morlet在进行石油勘探的地震数据处理分析时与法国理论物理学家A.Grossman一起提出了小波变换(wavelet transform, WT)的概念 10小波（小波（WaveletWavelet），），即小区域的波，是一种特殊的即小区域的波，是一种特殊的长度有限、平均值为长度有限、平均值为0 0的波形。的波形。特点：特点：（1 1）“小小”，即在时域都具有紧支集或近似紧支集，即在时域都具有紧支集或近似紧支集 （2 2）正负交替的

4、）正负交替的“波动性波动性”，也即直流分量为零，也即直流分量为零1112傅立叶分析所用的正弦波在时间上没有限制，从负无穷傅立叶分析所用的正弦波在时间上没有限制，从负无穷到正无穷，但小波倾向于不规则与不对称。到正无穷，但小波倾向于不规则与不对称。FTFT将信号分解成一系列不同频率正弦波的叠加将信号分解成一系列不同频率正弦波的叠加，小波分小波分析是将信号分解成一系列小波函数的叠加析是将信号分解成一系列小波函数的叠加。而这些小波。而这些小波函数都是由一个母小波函数经过平移与尺度伸缩得来的。函数都是由一个母小波函数经过平移与尺度伸缩得来的。用不规则的小波函数来逼近尖锐变化的信号显然要比光用不规则的小波

5、函数来逼近尖锐变化的信号显然要比光滑的正弦曲线要好滑的正弦曲线要好，同样，信号局部的特性用小波函数，同样，信号局部的特性用小波函数来逼近显然要比光滑的正弦函数来逼近要好。来逼近显然要比光滑的正弦函数来逼近要好。1314 连续小波变换（连续小波变换（Continuous Wavelet TransformContinuous Wavelet Transform， CWTCWT）用下式表示）用下式表示： (,)( ) (, )scale positionf tscale positionCt dt表示小波变换是表示小波变换是信号信号f f( (x x) )与与被缩放和平移被缩放和平移的的小波函数小

6、波函数()()之积在信号存在的整个期间里求和的结果。之积在信号存在的整个期间里求和的结果。1.4 1.4 连续小波变换连续小波变换CWTCWT的的变换结果变换结果是许多是许多小波系数小波系数C C，这些系数是缩放因，这些系数是缩放因子（子（scalescale）和平移（）和平移（positionposition）的函数。）的函数。 15 基本小波函数()的缩放和平移操作 (1) (1) 缩放缩放。就是压缩或伸展基本小波，缩放系数越小，则小波越窄小波的缩放操作 OOOf (t)f (t)f (t)tttf (t)(t)；scale1f (t)(2t)；scale0.5f (t)(4t)；scal

7、e0.2516(2) (2) 平移平移。小波的延迟或超前。在数学上小波的延迟或超前。在数学上, ,函数函数f f( (t t) )延延迟迟k k的表达式为的表达式为f f( (t-kt-k) )，小波的平移操作(a) 小波函数(t)； (b) 位移后的小波函数(t-k) 基本小波函数()的缩放和平移操作小波变换的步骤: 一 取一个小波与信号的最前面部分比较; 二 计算相关因子C,C代表小波和这段数据的相关性即:C越大,两者越相似;17 三 移动小波,重复步骤一和二,一直遍历整个数据;1819小波尺度和信号频率的关系小波尺度和信号频率的关系小尺度小尺度 信号的高频信号的高频大尺度大尺度 信号的低

8、频信号的低频20 在在每个可能的缩放因子和平移参数下计算小波系每个可能的缩放因子和平移参数下计算小波系数数，其计算量相当大，将产生，其计算量相当大，将产生惊人的数据量惊人的数据量，而且有，而且有许多数据是无用许多数据是无用的。的。1.6 1.6 离散小波变换（离散小波变换（DWTDWT）如果如果缩放因子和平移参数都选择为缩放因子和平移参数都选择为2 2j j（j j00且为且为整数）的倍数整数）的倍数， 即只选择部分缩放因子和平移参数即只选择部分缩放因子和平移参数来进行计算，来进行计算， 就会使分析的就会使分析的数据量大大减少数据量大大减少。21使用这样的缩放因子和平移参数的小波变换称为使用这

9、样的缩放因子和平移参数的小波变换称为双尺度小波变换（双尺度小波变换（Dyadic Wavelet TransformDyadic Wavelet Transform），它，它是是离散小波变换离散小波变换（Discrete Wavelet TransformDiscrete Wavelet Transform， DWTDWT）的一种形式。）的一种形式。1.6 1.6 离散小波变换（离散小波变换（DWTDWT）通常离散小波变换就是指双尺度小波变换。通常离散小波变换就是指双尺度小波变换。 22 执行离散小波变换的执行离散小波变换的有效方法是使用滤波器有效方法是使用滤波器， 该方法是该方法是Malla

10、tMallat于于19881988年提出的，称为年提出的，称为MallatMallat算法算法( (马马拉拉) )。这种方法实际上是一种信号分解的方法，。这种方法实际上是一种信号分解的方法， 在数在数字信号处理中常称为字信号处理中常称为双通道子带编码双通道子带编码。1.6 1.6 离散小波变换（离散小波变换（DWTDWT）23一个滤波器为低通滤波器，通过该滤波器可得到信号的近似值A（Approximations）另一个为高通滤波器， 通过该滤波器可得到信号的细节值D（Detail）。1.6 1.6 离散小波变换（离散小波变换（DWTDWT）24实际应用中，实际应用中，信号的低频分量往往是最重要

11、的信号的低频分量往往是最重要的，而高频分量只，而高频分量只起一个修饰的作用。如同一个人的声音一样，起一个修饰的作用。如同一个人的声音一样， 把高频分量去掉把高频分量去掉后，听起来声音会发生改变，但还能听出说的是什么内容，但后，听起来声音会发生改变，但还能听出说的是什么内容，但如果把低频分量删除后，就会什么内容也听不出来了。如果把低频分量删除后，就会什么内容也听不出来了。 1.6 1.6 离散小波变换（离散小波变换（DWTDWT）25图图 多级信号分解示意图多级信号分解示意图（a a） 信号分解；信号分解； (b) (b) 小波分树；小波分树； （c c）小波分解树）小波分解树 26在使用滤波器

12、对真实的数字信号进行变换时，在使用滤波器对真实的数字信号进行变换时，得到的数据将是得到的数据将是原始数据的两倍原始数据的两倍。 根据根据尼奎斯特尼奎斯特(Nyquist)(Nyquist)采样定理就采样定理就提出了降采样的方提出了降采样的方法，即在每个通道中每两个样本数据取一个，得到的法，即在每个通道中每两个样本数据取一个，得到的离散小波变换的系数离散小波变换的系数(coefficient)(coefficient)分别用分别用cDcD和和cAcA表示表示 1.6 1.6 离散小波变换（离散小波变换（DWTDWT）27 将信号的小波分解的分量进行处理后，一般还要将信号的小波分解的分量进行处理后

13、，一般还要根据需要把根据需要把信号恢复出来信号恢复出来，也就是利用信号的小波分，也就是利用信号的小波分解的系数还原出原始信号，这一过程称为解的系数还原出原始信号，这一过程称为小波重构小波重构（Wavelet ReconstructionWavelet Reconstruction）或叫做小波合成或叫做小波合成（Wavelet SynthesisWavelet Synthesis）。）。1.7 1.7 小波重构小波重构 这一合成过程的数学运算叫做这一合成过程的数学运算叫做逆离散小波变换逆离散小波变换（Inverse Discrete Wavelet TransformInverse Discre

14、te Wavelet Transform，IDWTIDWT）。）。 28小波重构算法示意图 SHLHL1.7 1.7 小波重构小波重构 29 (1) (1) 重构近似信号与细节信号重构近似信号与细节信号由小波分解的近似系数和细节系数可以重构出原由小波分解的近似系数和细节系数可以重构出原始信号。始信号。同样，可由近似系数和细节系数分别重构出信号同样，可由近似系数和细节系数分别重构出信号的的近似值近似值或或细节值细节值，这时只要近似系数或细节系数置，这时只要近似系数或细节系数置为零即可。为零即可。 1.7 1.7 小波重构小波重构 30重构近似和细节信号示意（a） 重构近似信号； (b) 重构细节

15、信号 A1HL1000个 样 点0约 500个 0cA1约 500个 近 似 分 量(a)D1HL1000个 样 点(b)约 500个 0约 500个 近 似 分 量0cD11.7 1.7 小波重构小波重构 31 (2)多层重构重构出信号的近似值A1与细节值D1之后，则原信号可用A1D1S重构出来。对应于信号的多层小波分解，小波的多层重构图: 1.7 1.7 小波重构小波重构 32重构过程为：A3D3A2；A2D2A1；A1+D1S。A3D3A2D2SA1D11.7 1.7 小波重构小波重构 心电信号的噪声特点小波分析与传统信号处理方法的比较小波去噪的基本原理小波去噪的基本步骤小波去噪中的阈值

16、函数和阈值的选取小波去噪中小波函数的选择去噪效果的评价程序说明总结心电信号（ECG）是典型的强噪声的非平稳的随机信号。正常心电信号的频率范围在0.01Hz-100Hz之间，而90%的ECG频谱能量又集中在0.25Hz-35Hz之间。在心电信号的采集和A/D转换过程中,心电信号不可避免地受到各种类型的噪声干扰,概括起来主要包括以下三类噪声:一、心电信号的噪声特点一、心电信号的噪声特点（1）由于电源磁场作用于心电图机与人体之间的环形电路所致的50Hz/60Hz工频干扰;（2）由于病人肌肉紧张产生的肌电干扰;（3）由于病人呼吸运动或者由电极电极皮肤之间界面阻抗所致的频响,一般小于1Hz的基线漂移；

17、这些噪声干扰与心电信号混杂,引起心电信号的畸变,使整个心电信号波形模糊不清,对随后的信号分析处理,尤其是计算机自动识别诊断造成误判和漏判,因此,心电信号的消噪有重要的意义。00.511.522.533.544.55-0.8-0.6-0.4-0.200.20.40.60.8128 1111111811111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111111Time / sVoltage / mVECG signal 100.dat)(*)()(tetfts1_, 1

18、, 0nt在实际工程应用中，通常所分析的信号具有非线性，非平稳,并且奇异点较多的特点。含噪的一维信号模型可表示为：其中，f(t)为真实信号，s(t)为含噪信号，e(t)为噪声， 为噪声标准偏差。 有用信号通常表现为低频信号或是相对比较平稳。而噪声信号通常表现为高频分解后，含噪部分主要集中在高频小波系数中，并且，包含有用信号的小波系数幅值较大，但数目少；而噪声对应的小波系数幅值小，数目较多。 基于上述特点，可以应用门限阈值法对小波系数进行处理。（即对较小的小波系数置为，较大的保留或削弱），然后对信号重构即可达到消噪的目的。CD3CD2CA3CA2sCA1CD1小波分解的结构示意图小波分解系数示意

19、图三、小波分解示意图：1.小波变换去噪的流程示意图：预处理小波变换多尺度分解各尺度小波系数除噪小波逆变换重构信号除噪后的信号含噪信号小波除噪的具体步骤：(1) 对含噪信号进行预处理,并进行小波分解。选择小波确定分解的层数N,然后对信号s进行N层分解。(2) 小波分解的高频系数的阈值量化。对第一层到第N层高频系数，选择软阈值或硬阈值量化处理。() 一维小波重构。根据小波分解的第N层低频系数和第一层到第N层的高频系数，进行一维重构。 在上面的步骤中，最为关键的就是如何选取阈值和如何阈值量化，从某种意义上讲，它直接影响信号去噪的质量。为阈值。阈值后的小波系数，为小波系数，设ww阈值函数阈值函数分为软

20、阈值和硬阈值两种。(1).硬阈值(hard threshol ding)当小波系数的绝对值大于等于给定阈值时，保持不变，而小于时，令其为。即：(2).软阈值(soft threshol ding)当小波系数的绝对值大于等于给定的阈值时，令其值为减去阈值；而小于时，令其为即：采用这种阈值方法去噪在实际应用中，已取得了较好的效果，但也存在着一些潜在的缺点，如硬阈值在阈值点不连续，重构可能产生一些震荡；软阈值连续，但估计的小波系数和分解的小波系数有恒定的偏差，直接影响重构信号对真实信号的逼近程度阈值的选取 阈值的选择是小波去噪和收缩最关键的一步，在去噪过程中阈值起着决定性的作用：如果太小，施加阈值后

21、小波系数包含太多的噪声分量，达不到去噪效果；反之，则去除了有用部分，使信号失真。阈值选择方案及对应的MATLAB命令(1) 固定阈值(sqtwolog)选取的算法是： (2) Stein无偏似然估计阈值(rigrsure)对于给定一个阈值t，得到它的似然估计，再将非似然的t最小化，就得到了所选的阈值。(3) 启发式阈值(heursure)它是前两种阈值的综合，是最优预测变量阈值选择，如果信噪比很小时，无偏似然估计的误差较大，此时，采用固定阈值。令：进行比较，如果时采用固定阈值，反之，选择无偏似然估计。(4) 极大极小阈值(minimaxi)它的原理是令估计的最大风险最小化，其阈值选取的算法是：

22、小波变换不像傅里叶变换是由正弦函数唯一决定的，小波基可以有很多种，不同的小波适合不同的信号去噪，对于确定的信号，如果小波选择不当，去噪结果可能相差很远，还有可能丢失有用的信息。面对各种小波，到底选择哪一种来处理心电信号才能满足医疗上的需要，必须经过大量的仿真研究结果来进行筛选 。根据大量文献记录B样条函数适合心电去噪：样条函数是一种非紧支撑正交的对称小波，有较高的光滑性，频率特性好，分频能力强，频带相干小的特性。在信号处理中小波的作用是带通滤波器，且对称和反对称性分别等价为线性相位和广义线性相位。当一个带通滤波器不是线性相位或广义线性相位时，它将使通过的信号产生畸变。从理论和实际应用的观点出发

23、，具有紧支集的小波是最富吸引力的。 B样条是一类基本的样条函数，而它的支撑区是最小的所以，B样条小波是一种合适的选择。式中yi表示标准原始信号, xi 表示经处理后的估计信号。其中，SNR越大越好, MSE 越小越好。输出评价结果及去噪后的心电信号去噪效果的评价根据选取的阈值去噪及重构求各层的阈值提取各尺度小波系数对被选的心电信号进行小波分解选取其中一个导联的数据调用心电数据开始去噪程序流程图八、小波去噪程序%应 用 db5作 为 小 波 函 数 进 行 3层 分 解%利 用 无 偏 似 然 估 计 阈 值%对 100.dat from MIT-BIH-DB的 单 导 联 数 据 进 行 去 噪 处 理clear;clcload(D:/matlab/matlab7.2/work/M.mat);E=M(:,2);E=E;n=size(E);s=E(1:2000);%小 波 分 解C L=wavedec(E,3,db5);% 从 c中 提 取 尺 度 3下 的 近 似 小 波 系 数cA3=appcoef(C,L,db5,3);%从 信 号 c中 提 取 尺 度 1,2,3下 的 细 节 小 波 系 数cD1=detcoef(C,