小波阈值去噪分析

数据挖掘流程

传统的去噪方法根据噪声能量一般集中于高频,而信号频谱分布于一个有限区间的特点,用傅里叶变换将含噪信号变换到频域,然后采用低通滤波器进行滤波。但基于傅里叶变换的去噪方法不能将有用信号的高频部分和由噪声引起的高频干扰加以有效地区分,存在着保护信号局部性和抑制噪声之间的矛盾。而小波变换可以很好的保存信号的尖峰和局部突出部分。小波去噪的发展可以归纳为三种方法:

(1)Mallat提出的小波系数模极大值去噪法

(2)徐长发等人提出的基于小波系数的空域相关性去噪算法

(3)Donoho与Johnstone提出的小波阈值收缩去噪法。

由于小波阈值法是其它两种方法的发展,具有一定的优势性。因而在此选择应用广泛的小波阈值去噪分析。

算法原理

小波阈值去噪原理

小波阈值去噪的基本思想是将信号通过Mallat算法处理过的小波变换后,选择产生的小波系数。因为信号经小波分解后小波系数较大,噪声的小波系数较小,并且噪声的小波系数要小于信号的小波系数。通过选取一个合适的阈值,大于阈值的小波系数被认为是有信号产生的,应予以保留,小于阀值的则认为是噪声产生的,置为零从而达到去噪的目的。

1.1.2 小波阈值收缩法

小波阈值收缩法的基本步骤为:

(1)分解:选定一种层数为N的小波对信号进行小波分解。

(2)阈值处理过程:分解后通过选取一个的阈值,用阈值函数对各层系数进行量化。

(3)重构:用处理后的系数重构信号。

其中分解和重构均为一般的小波变换,这里不具体分析,这里主要分析阈值处理过程。

小波阈值函数的选取

小波分解阈值去噪中,阈值函数体现了对超过和低于阈值的小波系数不同处理策略,是阈值去噪中关键的一步。设w表示小波系数,T为给定阈值,sign(*)为符号函数,常见的阈值函数有:

硬阈值函数:小波系数的绝对值低于阈值的置零,高于的保留不变,见下左图

软阈值函数:小波系数的绝对值低于阈值的置零,高于的系数shrinkage处理,见下右图

硬阈值法获得的重构信号具有更好的逼近性,但有附加振荡。软阈值获得的重构信号具有更好的光滑性,但误差相对较大。

小波阈值的确定

去噪选取的阈值最好刚好大于噪声的最大水平,可以证明的是噪声的最大限度以非常高的概率低于(此阈值是Donoho提出的),其中根号右边的这个参数(叫做sigma)就是估计出来的噪声标准偏差(根据第一级分解出的小波细节系数,即整个det1绝对值系数中间位置的值)。

Matlab的算法实现

Matlab调用的基本方式

在Matlab中,一维小波阈值去噪主要使用以下常用命令:

wnoisest函数 估计一维小波高频系数中的噪声偏差

wavedec函数 一维信号的多尺度分解

waverec函数 一维信号小波分解系数的重构

wdencmp函数 用于对一维或二维信号的压缩或者去噪

gbl(global的缩写)表示每层都采用同一个阈值进行处理,

lvd表示每层用不同的阈值进行处理,

N表示小波分解的层数,

THR为阈值向量,

对于格式(2)(3)每层都要求有一个阈值,因此阈值向量THR的长度为N,

SORH表示选择软阈值还是硬阈值(分别取为`s`和`h`),

参数KEEPAPP取值为1时,则低频系数不进行阈值量化处理,反之,则低频系数进行阈值量化。

XC是消噪或压缩后的信号,[CXC,LXC]是XC的小波分解结构,

PERF0和PERFL2是恢复和压缩L^2的范数百分比, 是用百分制表明降噪或压缩所保留的能量成分。

1.2.2Matlab脚本使用实例

以下是一维小波阈值去噪的具体过程和结果(完全分解图)

Matlab小波工具箱使用实例

选择App中Multiple 1-D中的Multivariate Denoising里面的Demo。


写在最后:这是我之前做的一些技术初步调研,主要目的是了解和展示一些技术,并进行了一些个人总结,其中必然有疏漏。其中引用了不少论文,百科和他人博客,但因为时间太久大都忘了出处,如果侵权提示后删除

小波调研(三):小波阈值去噪分析相关推荐

  1. 傅里叶变换 matlab代码_基于机器学习的心律失常分类(三)——小波阈值去噪[MATLAB]...

    在对心电信号进行分析前,需要对其进行去噪预处理.因为在心电信号采集或者数模转换的过程中,会产生噪声,包括工频干扰.肌电干扰.基线漂移这三种.本文将介绍目前十分常用且有效的一种去噪方法--小波阈值去噪. ...

  2. 自适应小波阈值去噪python_基于python的小波阈值去噪算法

    小波图像去噪原理 图像和噪声在经小波变换后具有不同的统计特性:图像本身的能量对应着幅值较大的小波系数,主要集中在低频(LL)部分:噪声能量则对应着幅值较小的小波系数,并分散在小波变换后的所有系数中.基 ...

  3. 小波阈值去噪原理及实现

    小波阈值去噪原理及实现 1.小波阈值去噪的基本原理 2.小波阈值去噪需要考虑的问题 3.小波阈值的实现 1.小波阈值去噪的基本原理 将信号通过小波变换(采用Mallat算法)后,信号产生的小波系数含有 ...

  4. matlab 小波启发式阈值滤波,小波阈值去噪

    小波阈值去噪 小波阈值去噪 目录 1.概念 2.原理 3.影响降噪效果的因素 3.1小波基的选择 3.2分解层数的选择 3.3阈值的选择 3.4阈值函数的选择 4.MATLAB代码 参考文献 小波分析 ...

  5. 一维信号小波阈值去噪

    转载:http://blog.csdn.net/ebowtang/article/details/40481393 一,小波阈值去噪基本理论       本博文根据小波的分解与重构原理,实现了基于硬阈 ...

  6. 基于MATLAB的小波阈值去噪

    非常感谢: http://blog.sina.com.cn/u/1861445474 http://blog.chinaaet.com/detail/3083.html MATLAB中实现了信号的阈值 ...

  7. 小波阈值去噪的原理及程序

    小波阈值去噪的原理及程序 一.小波阈值去噪基本思想 二.小波变换理论基础 三.小波函数和阈值的选择 四.程序 4.1 MATLAB程序 4.2 Python程序 4.3 效果图 一.小波阈值去噪基本思 ...

  8. 小波变换和小波阈值去噪

    1. 小波变换 小波变换是一种信号的时间--尺度(时间--频率)分析方法,它具有多分辨分析的特点,而且在时频两域都具有表征信号局部特征的能力,是一种窗口大小固定不变但其形状可改变,时间窗和频率窗都可以 ...

  9. sEMG信号小波阈值去噪

    1.sEMG信号简介 表面肌电信号是从皮肤表面获取的关于肌肉活动的微弱生物电信息,能够为肌肉活动状态分析.神经肌肉功能分析及病理分析等研究提供有效的帮助.由于受到检测仪器.环境和肌肉自身生理特性等多种 ...

最新文章

  1. 2011-2012年《软件工程》本科期末试卷
  2. 红曲面怎么做_新年就要“红”运当头 买款高颜值红色旗舰让年味更浓
  3. 兵家必争之地!爱尔兰是如何成为物联网的先锋?
  4. java list 转 map_高并发下的Java数据结构(List、Set、Map、Queue)
  5. android wifi驱动_OTT盒子WiFi方案首选:博通2T2R WiFi模块
  6. hbase 学习(十四)Facebook针对hbase的优化方案分析
  7. mysql online ddl和pt_MySQL变更之:Online DDL 和 PT-OSC 该选谁?
  8. 洛谷P2016 战略游戏【树形dp】
  9. Mr.J--JavaScript-恶搞小代码
  10. vue中 v-show和v-if的区别?
  11. 【js高级程序设计】迭代器
  12. Drools7中文教程 文档 指南
  13. 记录手机连接笔记本热点无法上网问题
  14. 他靠写程序封神,买不起操作系统,就动手写一个!曾拒绝乔布斯的offer
  15. [Typecho]更换头像源及显示QQ显示头像
  16. android极光推送 消息穿透广告弹窗,极光推送弹出两次消息,请大佬们帮忙指点迷津...
  17. 【STM32F429开发板用户手册】第33章 STM32F429的SPI总线应用之驱动DAC8563(双通道,16bit分辨率,正负10V)
  18. IPFS -- 节点搭建
  19. 万字长文详解目标检测算法,超大超全(2022最新)
  20. REW声学测试(三):生成测试信号

热门文章

  1. php new object delete,php – S3 DeleteObject – DeleteMarker始终返回空
  2. if [ $# -ne 1 ] 作用
  3. 系统迁移必知会(多年总结)
  4. GNS模拟linux,GNS3模拟器官方版最新版及安装说明
  5. 遥感影像条带噪声去除
  6. IMDB影评的情感分析任务:双层bi-LSTM模型(双向LSTM)
  7. golang实现稀疏数组(Sparse array)
  8. ABP官方文档(十三)【对象之间的映射】
  9. 利用集合进行数组的排序
  10. implicitly has an ‘any‘ type...以及suppressImplicitAnyIndexErrors版本错误TypeScript 5.5问题