小波模极大值除噪
模极大值与阈值决策融合的小波语音数据去噪方法
20 0 8年 1 0月
计 算 机 应 用 研 究
Ap l ai n Re e r h o o u e s p i t s a c fC mp t r c o
Vo . 5 No 1 Leabharlann 2 . O 0c . 2 0 t o8
模 极 大 值 与 阈值 决 策 融 合 的小 波 语 音 数 据 去 噪 方 法 术
a d n n i e rt rs od d n i n to n w v l t h sme h d a s r e h d a tg s a d o e c me h i d a tg s n o l a e h l e os g meh d o a ee .T i n h i t o b o b d te a v n a e n v ro d t e d s v n a e a o rd t n lmeh d .T k n e p r n ov l ae te i rv d meh d, ih s o h e t o a ea v na e o ft i o a to s a e a x e i t ai t h a i me t d mp o e t o whc h ws en w me h d h st d a tg f t h
去噪 方法 , 有效地结合 了两 者的优 点 。经 实验证 明 , 方法精度 高、 该 算法 简单 、 计算 速度快 。
关 键 词 :小 波 去 噪 ; 极 大值 ;阈值 ; 音 数 据 处 理 模 语
中图分 类号 :T 3 l P 9
文献标志 码 :A
文章 编号 :lO — 6 5 2o )0 33 —2 O 13 9 ( 0 8 l— 14 0
4 o7 1 o 3,C ia) ^n A b t ac Th sp p rprpo e e wa ee— nosn meh d, whih wa a e n wa ee rnso m du u ma i m s r t: i a e o s d a n w v ltde ii g to c s b s d o v ltta fr mo l s xmu
小波变换模极大值点在信号去噪中的应用
有I 暇 ( ( J 。 J ) l ≤ l暇 (%
。
M
口
¨
a
:
+ c一
s 。 ,
)I
,
则 称 (卟
。 ,) 为 小
竿卜 I
仔 s,
波 变换 的 模 极 人 值 点
if 2
且 重 构 后 的 信 号 非 常逼 近 原 始 信 号
。
但
小 波模 极 大 值 点 法 去噪 效 果 更好
圈■圈黧鳃翱麓骥璃刚阁骥簇黧赣 鞫 隳嘲鞫 鳓礁
小 波 变换 ;模 极 大 值 点 ;l 融 蛾
lt 随
z
指 数 ;信 号 去 噪
鞠囊嘲黪瀚豳黧嚼赣翻隧赠黧鼹罐窝赣黼鳃圈潮黼
a
j
f
rom
%∥)
:
l,
1 『 表, J÷ 为:
口 ≠
坐 嘉 f 卜6 咄 H 、 /
妒
\
“
0
2. 【 2)
’
仃崽 平 方 可 积 厢 数 / (, )
.
,
记作/ (, ) ∈ f ( 尺)
,
它 的 小波 分 析
定 义 为 以 函 数 系 眈 。( , ! 分 变换 : ) 为秋 分 核 的 移
¨ m )
1
=
(川
(2 3 )
-
‰
∽)
=
C巾
蹦
引
伽
州
。 :
[
川
d 等)
加
引言
小波 变换是 近 年 来继 傅里 叶 变换 后 的
。
叶 变换是 分 析平 稳信号的理 想 工 具
傅里 博 里 叶 变换 能提 取 出 函 数 在
信号处理之小波去噪方法介绍
本文对各种去噪方法进行了比较,总结了两大类方法的基本思想及实现流程,详细介绍了应用最广的小波阈值去噪。
一、小波去噪主要方法1、基于小波分频的去噪方法——主要用来压制面波等规则干扰;2、小波域去噪方法——主要用于压制随机干扰,目前主要有三种方法: a) 模极大值去噪方法(Mallat 和Zhang ,1992)b) 尺度相关性分析方法(Xu ,1994)c) 小波阈值收缩方法(Dohono 和Johnstone ,1994)其中,小波阈值去噪方法能在最小均方误差意义下得到信号的近似最优估计,计算速度快,适应性广,因此应用最广泛。
二、方法实现的总体流程1、基于小波分频的去噪方法小波时频分析使信号在空间域和频率域同时具有良好的局部分析性质。
小波变换可以将信号分解到各个不同的尺度或各个不同的频段上,并且通过伸缩、平移聚焦到信号的任一细节加以分析。
小波分析的这些特长,结合传统的傅立叶去噪方法,为地球物理信号去噪提供了有效途径。
对于离散序列信号,其小波变换采用 Mallat 快速算法, 信号经尺度j =1,2,…,J 层分解后,得到)(2R L 中各正交闭子空间(1W 、2W 、…、J W 、J V ), 若j j V A ∈代表尺度为j 的低频部分, j j W D ∈代表高频部分,则信号可以表示为J J D D A t f +++= 1)(,据此可重构出信号在尺度j =J 时的低频部分和j =1,2,…,J 的高频部分。
如果地震数据中的干扰波频率与有效波的频率成分是分开的,通过小波分频很容易消除干扰波;如果两种频率成分存在混叠,也可以用小波分频方法提取混叠部分,再用传统方法分离有效和干扰波。
这样可以最大限度的保留有效波能量。
2、小波域去噪方法小波域去噪方法是利用信号和噪声的小波系数在小波域不同特性来进行的。
信号和噪声的小波系数幅值随尺度变化的趋势不同,随着尺度的增加,噪声的小波系数很快衰减,而信号的小波系数基本不变。
模极大值去噪方法
3.2 基于小波变换模极大值去噪方法的研究目前利用小波变换消除噪声的方法很多,但总结起来,比较成熟的是Mallat 提出的一种多尺度小波变换模极大值的去噪方法。
3.2.1 小波变换模极大值的定义定义在尺度s 下,若0x x d " ,()()0,,Wf s x Wf s x £成立,则0x 称为模极大值点,()0,Wf s x 称为模极大值。
小波变换极大模是由信号中奇异点和噪声产生的。
根据理论分析,知道以平滑函数的一阶导数为母小波作小波变换,其小波变换在各个尺度下的模极大值对应于信号突变点的位置。
小波分析尺度越小,平滑函数的平滑区域小,小波系数模极大值点与突变点位置的对应就越准确。
但是小尺度下小波变换随噪声影响非常大,产生许多伪极值点,往往只凭一个尺度不能定位突变点的位置。
相反,在大尺度下对噪声进行了一定的平滑,极值点相对稳定,但由于平滑作用使其定位又产生了偏差。
同时,只有在适当尺度下各突变点引起的小波变换才能避免交迭干扰。
因此,在用小波变换模极大值法判断信号突变点时,需要把多尺度结合起来综合观察。
下面由小波变换模极大值在多尺度上的变化规律来表征信号突变点的性质。
在许多情况下,小波变换并不要求保留所有的连续尺度a ,为了实现快速算法,选择尺度按二进制变化,即二进制变换。
信号的突变点在不同尺度2j 上都会产生对应的模极大值。
在任意尺度2j 上模极大值对应于信号在2j 尺度上平滑后的该点一阶导数大小。
小波理论表明,模极大值的幅值随着尺度的变化规律是由信号在该突变点的局部李氏指数(Lipschitzexponent)决定的。
3.2.2 模极大值随着尺度的变化规律李氏指数的定义为,设函数在0t 附近具有下述特征:()()00,1n x t h p t h A h n n aa +-+?<+ (3-1)则称()x t 在0t 处的李氏指数为a 。
式中h 是一个充分小量,()n p t 是过()0x t 点的n 次多项式()n Z Î。
基于模极大值小波域的去噪算法研究
文章编号:1004-9037(2003)03-0315-04基于模极大值小波域的去噪算法研究张小飞1,徐大专1,齐泽锋2(1.南京航空航天大学信息科学与技术学院,南京,210016;2.武汉大学电气工程学院,武汉,430072)摘要:根据信号与噪声在小波变换下的不同特性,提出了基于模极大值小波域的去噪算法。
该算法先用Adho c 算法求出信号的模极大值,再根据模极大值小波域的定义求出信号的模极大值小波域,从而得到信号的小波系数,然后逆变换得到信号。
实例分析表明:该算法能有效消除噪声,与交替投影模极大值算法相比,该算法在原理上更简单,程序实现更容易,去噪速度更快,能满足在线监测的要求。
关键词:小波;去噪;小波变换模极大值;模极大值小波域中图分类号:T N 911.7 文献标识码:A 收稿日期:2002-08-30;修订日期:2003-02-17Denoising Algorithm Based on Modulus Maximum Wavelet FieldZ H AN G X iao -fei 1,X U Da -zhuan 1,QI Ze -f eng2(1.Colleg e of Info rma tio n Science and T echno log y ,Na nj ing U niv er sity of Aero nautics &Astro na utics ,N anjing ,210016,China ;2.Co llege o f Elect rical Engineering ,W uhan U niv e rsity ,W uha n,430072,China )Abstract:Befo re the sig nal is detected ,the noise included in signal sho uld be wiped off .A novel denoising algo rithm based on modulus m aximum w av elet field is presented.The cha racters of noise and sig nal o n w av elet transform a re discussed.The processes of the alg orith m are as fo llow s :firstly the modulus max imum is o btained acco rding to Adhoc alg orith m;seco ndly the mo dulus maximum w av elet field is com puted acco rding to its definitio n to a ttain wav elet coefficient;thirdly the singularity is reco nstructed th ro ug h inv erse w av elet transform .Com pared w ith others ,the a lg o rithm has som e adv antag es ,such as constructing efficiently ,w iping off the noise effectiv ely ,prog ramming easily .Ex amples prove that the alg orith m has better denosing perfo rmance,and ca n meet the demands o f o nline detection.The sing ula rity of sig nal is represented by Lipschitz index.Key words :w av elet ;denoising ;w av elet transform modulus maxim um ;modulus maximumw av elet field引 言在实际信号处理过程中采集到的信号包含大量噪声。
基于最优小波包基的模极大值法信号消噪研究
基于最优小波包基 的模极 大值法信号消 噪研究 *
孟庆 华 , 文 礼 赵
( 杭州电子科技 大学机 电工程学 院, 州 30 1) 杭 1 0 8
摘 要 : 给出了对信号进行基于熵准则的最优小波包基分解的基本原理, 在此基础上提出了使用模极大值法进行信号消噪
中图分类 号 :P 7 T 24
文献标 识码 :628 -3
信号处理中, 所采集到的信号总是含有不同成 分的噪声. 例如, 由于动力 电所引起的 5Hz 0 工频干 扰是 噪声 的一 个最 大来 源. 电磁 辐射 、 电子装 置 中电 子器件 的热噪声 , 对模拟信号抽样时所产生 的量化 噪声 , 有限位运算( 加减乘除) 产生的舍人误差噪声 等都是 噪声 的来 源. 目前 小 波 消噪 得 到 了广 泛 的应 用, 它充分利用 了信号与噪声 ( 或干扰) 在频谱上 的 差异 , 克服 了以往消 噪方法 的缺点 : 噪不彻底 和 当 消 含噪声 不规则 时 的过 于平 滑 , 具有 较好 的消 噪效果 . 小波消噪的原理是将含有噪声的信号在某一尺度下
分解到 不 同的频带 内 , 后再 将 噪声 所处 的频带 置 然
中” 的能力 , 它能将信号的能量集 中到少数小波系数 上 , 白噪声在任何正交基上的变换仍然是 白噪声. 而 相对来说, 信号的小波系数值必然大于那些能量分 散且幅值较小的噪声 的小波系数值. 选择一个合适 的阈值 , 对小波系数进行阈值处理 , 就可以达到消除 噪声而保留有用信号 的 目的. 该方法能得到原始信 号 的近似最优估计 , 目前应用非常广泛. 但是这种处 理 方法需 要 预先知 道 噪声 的先 验 知识 . 如果 不 清楚 的话往往得不到理想的消噪效果.
小波去噪三种方法
小波去噪常用方法目前,小波去噪的方法大概可以分为三大类:第一类方法是利用小波变换模极大值原理去噪,即根据信号和噪声在小波变换各尺度上的不同传播特性,剔除由噪声产生的模极大值点,保留信号所对应的模极大值点,然后利用所余模极大值点重构小波系数,进而恢复信号;第二类方法是对含噪信号作小波变换之后,计算相邻尺度间小波系数的相关性,根据相关性的大小区别小波系数的类型,从而进行取舍,然后直接重构信号;第三类是小波阈值去噪方法,该方法认为信号对应的小波系数包含有信号的重要信息,其幅值较大,但数目较少,而噪声对应的小波系数是一致分布的,个数较多,但幅值小。
基于这一思想,在众多小波系数中,把绝对值较小的系数置为零,而让绝对值较大的系数保留或收缩,得到估计小波系数,然后利用估计小波系数直接进行信号重构,即可达到去噪的目的。
1:小波变换模极大值去噪方法信号与噪声的模极大值在小波变换下会呈现不同的变化趋势。
小波变换模极大值去噪方法,实质上就是利用小波变换模极大值所携带的信息,具体地说就是信号小波系数的模极大值的位置和幅值来完成对信号的表征和分析。
利用信号与噪声的局部奇异性不一样,其模极大值的传播特性也不一样这些特性对信号中的随机噪声进行去噪处理。
算法的基本思想是,根据信号与噪声在不同尺度上模极大值的不同传播特性,从所有小波变换模极大值中选择信号的模极大值而去除噪声的模极大值,然后用剩余的小波变换模极大值重构原信号。
小波变换模极大值去噪方法,具有很好的理论基础,对噪声的依赖性较小,无需知道噪声的方差,非常适合于低信噪比的信号去噪。
这种去噪方法的缺点是,计算速度慢,小波分解尺度的选择是难点,小尺度下,信号受噪声影响较大,大尺度下,会使信号丢失某些重要的局部奇异性。
2:小波系数相关性去噪方法信号与噪声在不同尺度上模极大值的不同传播特性表明,信号的小波变换在各尺度相应位置上的小波系数之间有很强的相关性,而且在边缘处有很强的相关性。
基于小波变换模极大值的信号去噪方法
收 稿 日期 : 2 0 1 4 — 1 1 — 0 3
时间变化的信号 ( 平 稳随机过程 ) , 处理的理想工具仍然 是 傅立叶分 析。但是在实际应用 中绝大多数信号是非平稳的 ( 非平稳 随机过程 ) ,小波 分析 适用于非平稳 信号的处理 。 信号处理主要包括 : 信号去噪 、 特征提取 和边 缘提取[ 2 1 。信
Ke y w o r d s : s i g n a l d e n o i s i n g ; wa v e l e t t r a n s f o m; r mo d u l u s ma x i ma ; l a t e r n a t i n g p r o j e c t i o n me t h o d
( 4 ) 式给 出 了误 差 函数
) 的上界 , 不 过式 中的 指数
着尺度. 的增大而减小 ; 如果在 t o 处 的李 氏指数 a = O , 则小 波变换 的模极 大值不 随尺度 . 的变化而变化 。
n + l 是整数 , L i p s c h i t z 指数 是一个非整数 , 代替这里 的指
号去噪是信号处理 中最为 常见 的 , 本 文深入研究小 波变换
模 极 大 值 法 去 噪 的原 理 , 提 出一 种 简 单 的 基 于 小 波 变 换 模
极大值的去噪算法 。
1 1 8 l 震攘援燃 露 ・ 2 0 1 4 年 第6 期
2
度 内 则
量
数 似
有
用
鬻 剿㈣ 罄鬻 瓣激
工程 师手 记 E n g i n e e l " R e c o r d
小波去噪方法的研究
小波去噪方法的研究作者:李杨来源:《科技视界》2017年第25期【摘要】小波去噪是小波变换应用的常见用途,本文介绍了三种常用的小波去噪方式,分别是小波模极大值法、小波分解与重构法和小波阈值法。
将这几种方式分别加入白噪声并实施去噪处理,通过对比优缺点,为小波去噪方法的选取提供了参考依据。
【关键词】去噪;模极大值;小波分析0 引言在数学领域中小波分析发展迅速,它是一种时频分析的办法,小波分析具有多分辨率分析的特征,它能够对信号中任意一些细节做多分辨率的时频域分析,享有“数学显微镜”的美誉,它还有一个双重意义在于不光拥有深刻的理论还拥用广泛的应用。
1 小波去噪几种方法的基本原理1.1 小波模的极大值法去噪信号局部突变点的特征量度指数是Lip,它的定义:设有正整数n,n≤a≤n+1,如果存在正整数A>0及n次多项式Pn(x),使得|f(x)-pn(x-x0)|≤A|x-x0|a对x∈(x0-δ,x0+δ)成立,则称f(x)在x0点是Lipα。
在尺度s下的局部模的极大值点x0的定义是,若?坌x∈δx0,有|Wf(s,x)|≤|Wf(s,x0)|。
f(x)的Lip指数与小波变动模极大值条件是log2|W2jf(t)|≤log2k+jα。
对于常见的信号,由于α≥0,j的值越大小波模极大值也就越大;而对于白噪声,由于α1.2 小波阈值法小波变换阈值的方法称为“小波收缩”,这种方法是基于白噪声的小波系数的特点,大部分能量集中于较低幅度的小波系数,在任何正交基上变换,白噪声不会发生本质的变化,依然是白噪声,并具有相同的振幅。
小波系数将大于噪声信号的小波系数和相对分散的能量比较小,所以我们可以设计一个阈值,小于阈值的小波系数为零,从而有效地抑制信号中的噪声。
1.3 小波分解与重构法小波分解与重构的快速算法即Mallat算法,这是在1988年由S.Mallat在构造正交小波基的时候,第一次提出了多分辨分析的理论观点,小波的多分辨频率的特征从空间的理念上得到了解释。
强背景噪声下的小波熵模极大值消噪方法
基于小波变换模极大值的去噪方法研究
wh t o s ,wh n S i n ie e e NR o ,mo u u x mad n ii g me h d i s p r rt h e h l e n ii g o i lw s d l sma i e o sn t o u e i o t r s o d d - o sn ;f ri u s os ,wh n s o mp len ie e
要: 小波去噪 在信 号处理 中得到 广 泛 应 用。 目前 常 用 的方 法主 要 有 D n h o o o提 出 的 阈值 法和
Malt 出的模极 大值 法 。模 极 大值 法是 一种 经典 的 小波 去噪方 法 ,噪声 的模 极 大值 的 幅度 随尺 l 提 a
度 的增 大而迅 速减 小 ,而 正常信 号随尺 度的增 大而增 大 , 因此利 用合适 尺度 的 小波 变换 , 易把 噪 容
fo t e n r lsg a .Te tp o e h t o a e y g o e o sn fe to i o s n mp le n ie o u sa r m h o ma i n 1 s r v s t e me h d h s v r o d d n ii g ef c n wh t n iea d i u s o s .F rGa s in e
o d me h d p o o e y D n h n h d l sma i r p s d b l t l t o r p s d b o o o a d t emo u u x ma p o o e y Ma l .Th d l sma i t o ls ia wa e a e mo u u x ma me h d i a ca sc l v — s ltd - os n t o .Th d l sma i mp iu e o h o s e r a e a i l t h n r a e o c l ,wh l h o — e e n ii g me h d e mo u u x ma a l d ft e n ie d c e s s r p d y wih t ei c e s fs ae t i t e n r e
基于小波变换模极大值的信号去噪方法研究.
第32卷第1期2009年1月河北农业大学学报JOURNAL OF AGRICUL TURAL UNIVERSITY OF HEBEIVol. 32No. 1Jan . 2009文章编号:1000-1573(2009 01-0114-03基于小波变换模极大值的信号去噪方法研究张玉新, 滕桂法, 赵洋, 李阅历, 马建斌(河北农业大学信息科学与技术学院, 河北保定071001①摘要:信号在采集、转换和传输过程中, 由于受到设备、环境及人为因素的影响, 使信号不可避免地受到噪声干扰。
因此, 如何去除信号中的噪声, 得到感兴趣的信息是信号处理过程中的一项关键技术。
对基于小波变换模极大值的信号去噪问题进行了研究, 根据信号和噪声的小波变换模极大值在不同尺度上表现出的不同的传播特性, 给出了基于小波变换模极大值的去噪算法。
数值实验结果表明了该算法的有效性和可行性。
关键词:小波变换; 信号重构; 模极大值; 阈值; 去噪中图分类号:TP 335文献标识码:AR esearch on w avelet modulus m axim a in the de -methodZH ANG Y u 2xin , TE NG G ui 2fa , ZH ang LI ue li , (College of Information Science and , Baoding 071001, ChinaAbstract :In and transmission , signals are often corrupted by noise of ineffective equipments , environments and even human errors. De 2noising of extracting desired information has been a crucial technique in signal pro 2cessing. This paper discusses the signal de 2noising problem based on wavelet transform modulus maxima. According to the different characters of wavelet transform modulus maxima of signal and noise , a de 2noising algorithm based on wavelet transform modulus maxima is proposed. The experi 2mental results show that this method is efficient and practical.K ey w ords :wavelets transform ; signal reconstruction ; modulus maxima ; threshold ; de 2noising传统的去噪方法是将染噪信号进行傅立叶变换, 经滤波操作后再行傅立叶逆变换而获得“干净”的信号, 但由于低通平滑的作用, 使得在消除高频噪声的同时也会模糊边缘位置信息, 造成信号发生某种程度的畸变。
小波去噪三种方法
小波去噪三种方法小波去噪常用方法目前,小波去噪的方法大概可以分为三大类:第一类方法是利用小波变换模极大值原理去噪,即根据信号和噪声在小波变换各尺度上的不同传播特性,剔除由噪声产生的模极大值点,保留信号所对应的模极大值点,然后利用所余模极大值点重构小波系数,进而恢复信号;第二类方法是对含噪信号作小波变换之后,计算相邻尺度间小波系数的相关性,根据相关性的大小区别小波系数的类型,从而进行取舍,然后直接重构信号;第三类是小波阈值去噪方法,该方法认为信号对应的小波系数包含有信号的重要信息,其幅值较大,但数目较少,而噪声对应的小波系数是一致分布的,个数较多,但幅值小。
基于这一思想,在众多小波系数中,把绝对值较小的系数置为零,而让绝对值较大的系数保留或收缩,得到估计小波系数,然后利用估计小波系数直接进行信号重构,即可达到去噪的目的。
1:小波变换模极大值去噪方法信号与噪声的模极大值在小波变换下会呈现不同的变化趋势。
小波变换模极大值去噪方法,实质上就是利用小波变换模极大值所携带的信息,具体地说就是信号小波系数的模极大值的位置和幅值来完成对信号的表征和分析。
利用信号与噪声的局部奇异性不一样,其模极大值的传播特性也不一样这些特性对信号中的随机噪声进行去噪处理。
算法的基本思想是,根据信号与噪声在不同尺度上模极大值的不同传播特性,从所有小波变换模极大值中选择信号的模极大值而去除噪声的模极大值,然后用剩余的小波变换模极大值重构原信号。
小波变换模极大值去噪方法,具有很好的理论基础,对噪声的依赖性较小,无需知道噪声的方差,非常适合于低信噪比的信号去噪。
这种去噪方法的缺点是,计算速度慢,小波分解尺度的选择是难点,小尺度下,信号受噪声影响较大,大尺度下,会使信号丢失某些重要的局部奇异性。
2:小波系数相关性去噪方法信号与噪声在不同尺度上模极大值的不同传播特性表明,信号的小波变换在各尺度相应位置上的小波系数之间有很强的相关性,而且在边缘处有很强的相关性。
一种基于小波模极大值去噪的改进算法
( eatet fCv n i ei G in U i rt o e nl y G in 5 10 ) D p r n i l gn r g, ul n esyf Tc oo , ul 4 0 4 m o iE e n i v i h g i
可有效恢 复信号 。与传统的阈值 去噪法和基于平移不变量的去噪算法相 比, 该算法提高 了信噪 比。
关键 词 小波 模极大值 阈值去噪 平移不变量 形变监测
中图分类 号 : 14 0 7
文献标 识码 : A
AN M PROVED I DE. NO I I S NG ALG OI UTHM BAS ED N O 1 VELET oDULAR AXI UM l M M M
作者简介 : 何永红 , , 7 年生 , 女 1 8 9 硕士, 研究方向为 : 变形观测技术与工程信息系统 . m i hynh n20 @13 cn E— a : ogO g04 6 .o le
维普资讯
第 2 卷第2 7 期
200 4月 7年
大 地 测 量 与 地 球 动 力 学
J URNAL OF G ODE Y AND GE YNAMI S O E S OD C
V I2 . 0 . 7 No 2
Ap ., 0 r 2 07
了成 效 , 目前 国 际 上 极 为 活 跃 的 研 究 方 向 J 是 。 目前 , 常用 的小 波去 噪 方法 主要 有 基 于模 极 大值 原
理的传统的阈值法, 基于平移不变量的去噪算法 , 本 文依 据 Maa提 出 的信号 和 噪声 不 同 的奇 异指 数影 lt l 响小波变换系数模极大值随尺度传播特性的理论 , 判断噪声引起的模极大 , 并对其相邻 的模极小构成 的选通 域 进 行 直 线 平 滑 处 理 J 。考 虑 到 小 波 函数
基于三次样条插值的小波模极大值去噪算法
基于三次样条插值的小波模极大值去噪算法2014年8月第35卷第8期计算机工程与设计COMPUTERENGINEERINGANDDESIGNAug畅2014Vol畅35No畅8基于三次样条插值的小波模极大值去噪算法赵鸿图,刘云+(河南理工大学计算机科学与技术学院,河南焦作454000)摘要:为解决交替投影法计算量大、收敛速度慢的问题,提出一种基于三次样条插值的改进去噪算法。
利用Adhoc算法来搜索模极大值线,求出信号的模极大值点,通过三次样条函数对信号模极大值点进行插值,得到信号的估计小波系数进而重构信号。
Matlab仿真结果表明,该算法能在很大程度上去除噪声,提高信噪比,将仿真结果和其它算法进行比较,比较结果表明了该算法的收敛速度和信号优化质量明显优于其它算法。
关键词:交替投影法;模极大值;三次样条函数;小波系数;重构信号中图法分类号:TN911畅72文献标识号:A文章编号:1000‐7024(2014)08‐2965‐04收稿日期:2014‐01‐05;修订日期:2014‐03‐28基金项目:国家创新方法工作专项基金项目(2010IM020500)作者简介:赵鸿图(1965),男,河南长垣人,博士,副教授,研究方向为电力市场与信号处理;+通讯作者:刘云(1989),女,河南濮阳人,硕士,研究方向为现代数字信号处理。
E‐mail:lyun0304@163畅comDenoisingalgorithmofwaveletmodulusmaximumbasedoncubicsplineinterpolationZHAOHong‐tu,LIUYun+(SchoolofComputerScienceandTechnology,HenanPolytechnicUniversity,Jiaozuo454000,China)Abstract:Tosolvetheproblemsofthelargecomputationandtheslowconvergencespeedofthealternatingprojection(AP)method,animproveddenoisingalgorithmbasedoncubicsplineinterpolationwasputforward.ThemodulusmaximumlinewassearchedandthemodulusmaximumpointsofthesignalwerefoundbyusingtheAdhocalgorithm,thesignalestimatedwaveletcoefficientswereobtainedfromusingthecubicsplineinterpolationfunctionbetweenthepointsofmodulusmaximumandthesig‐nalwasreconstructed.TheMatlabsimulationresultsprovethattheimproveddenoisingalgorithmcanimprovetheSNRandre‐movethenoisetoagreatdegree.Theconvergencespeedandtheoptimalsignalqualityaresuperiortootheralgorithmsbycom‐paringwiththesimulationresults.Keywords:APalgorithm;modulusmaximum;cubicsplinefunction;waveletcoefficients;reconstructsignal0引言目前,研究者们在模极大值理论基础上对信号去噪算法提出了许多新的思想和新的改进[1‐4]。
小波变换语音降噪极大极小值算法python代码实现
小波变换语音降噪极大极小值算法概述小波变换语音降噪是一种常用的信号处理技术,可以有效地去除语音信号中的噪声,提高语音的清晰度和可懂度。
本文将详细介绍小波变换语音降噪中的一种常见算法——极大极小值算法,并给出其在Python中的实现。
小波变换简介小波变换是一种时频分析方法,能够将信号分解为不同频率的子信号。
在语音处理中,通过对语音信号进行小波变换,可以将语音信号的频谱特性和时域特性同时考虑,从而更好地去除噪声。
极大极小值算法极大极小值算法是小波变换语音降噪中常用的一种方法。
该算法的基本思想是通过寻找小波变换系数的极大值和极小值来确定信号的频段,并将低频信号与高频信号进行分离。
具体流程如下:1.选择合适的小波基函数,并将原始语音信号进行小波变换。
2.寻找小波变换后的系数序列中的极大值和极小值点。
3.根据极大值和极小值点确定信号的频段。
4.根据频段将原始信号进行分离,分别处理高频信号和低频信号。
5.对分离后的信号进行去噪处理。
6.将去噪处理后的信号进行合并,得到最终的去噪语音信号。
Python实现在Python中,可以使用pywt库来进行小波变换和相关操作。
以下为使用极大极小值算法进行小波变换语音降噪的Python代码实现:import pywtimport numpy as npdef wavelet_denoising(signal, wavelet='db4', level=4):# 进行小波变换coeffs = pywt.wavedec(signal, wavelet, level=level)# 获取小波系数的极大值和极小值索引maxima_idx = []minima_idx = []for i in range(1, len(coeffs)):maxima = np.where((coeffs[i-1][:-1] < coeffs[i][:-2]) & (coeffs[i-1] [1:] < coeffs[i][2:]))[0]minima = np.where((coeffs[i-1][:-1] > coeffs[i][:-2]) & (coeffs[i-1] [1:] > coeffs[i][2:]))[0]maxima_idx.append(maxima)minima_idx.append(minima)# 根据极大值和极小值索引确定信号的频段high_freq_mask = np.zeros_like(coeffs[0])for i in range(len(maxima_idx)):high_freq_mask[maxima_idx[i]] = 1# 分离低频信号和高频信号high_freq_coeffs = coeffs[:]high_freq_coeffs[0] = high_freq_coeffs[0] * high_freq_masklow_freq_coeffs = coeffs[:]low_freq_coeffs[0] = low_freq_coeffs[0] * (1 - high_freq_mask)# 对高频信号进行去噪处理for i in range(1, len(high_freq_coeffs)):high_freq_coeffs[i][maxima_idx[i-1]] = 0high_freq_coeffs[i][minima_idx[i-1]] = 0# 重构信号denoised_signal = pywt.waverec(low_freq_coeffs, wavelet)return denoised_signal# 测试代码input_signal = np.random.randn(1000) # 原始语音信号noisy_signal = input_signal + 0.1 * np.random.randn(1000) # 添加高斯噪声denoised_signal = wavelet_denoising(noisy_signal, wavelet='bior3.3', level=3) # 进行小波变换语音降噪# 输出结果print('原始语音信号:', input_signal)print('带噪声语音信号:', noisy_signal)print('去噪后语音信号:', denoised_signal)结论本文介绍了小波变换语音降噪中的一种常见算法——极大极小值算法,并给出了其在Python中的实现。
基于Hermite插值的小波模极大值重构滤波的肌电信号消噪方法
H
(x
)
=
⎜⎜⎜⎜⎝⎛1
+
2
x x1
− x0 −x0
⎠⎞⎟⎟⎟⎟⎝⎛⎜⎜⎜⎜
x x0
− x1 − x1
⎠⎞⎟⎟⎟⎟2
y0
+
⎛⎜⎜⎜⎜⎝1
+
2
x x1
− x1 −x0
⎞⎠⎟⎟⎟⎟⎛⎝⎜⎜⎜⎜
x x1
− x0 − x0
⎞⎠⎟⎟⎟⎟2
y1
(1)
利用式(1)对小波函数各尺度相邻模极大值之间的节点
进行插值,即能重构小波函数。显而易见,插值多项式最高
次数仅为 3,算法比交替投影重构简化了许多,节约了计算
时间。
4 实验分析
本文所述的实验系统构成为:美国奇石乐公司的肌电信 号采集仪、美国 NI 公司的数据采集卡 6024E 和 NI 的 Labview 系统。实验前,选定一组合适的肌电信号拾取电极,放于健 康男性受试者的尺侧腕伸肌和尺侧腕屈肌处,用于采集 sEMG 信号。肌电信号采集仪输出至 NI 数据采集卡 6024E, 随后用 Labview 软件获取采样的 N 点数据 xN (n) , (n=0,1, ",N − 1) , 得 到 表 面 肌 电 信 号 离 散 的 采 样 值 , 然 后 在 MATLAB7.3 平台中,完成数据的导入和小波消噪算法。
2007-12-24 收到,2008-04-22 改回 国家自然科学基金(60705010)和浙江省科技计划(2007C23088)资助 课题
2 小波系数模极大值
函数的奇异性是指函数在某处有间断点或者某阶导数
858
电子与信息学报
第 31 卷
不连续,通常用 Lipschitz 指数来度量。
基于Blackman窗函数插值的小波模极大值去噪算法
信号去噪算法的研究一直是信号处理领域的难点和热点,长期以 来人们常用傅里叶变换分析处理信号,而且取得了较好的效果。 然而傅里叶变换是一种全域的整体变换,只能显示一个域内的信 号特性,对于分析非平稳信号有很大的局限性。
小波变换是继傅里叶变换迅猛发展起来的一种分析非平稳信号 的强有力工具,具有在时频域内同时表征信号局部化特性和多分 辨分析的能力,不但在信号处理领域有极好的发展前景,而且与 科学界其他领域相互融合、相互渗透形成一种新生力量。本文 主要研究了利用小波模极大值重构小波系数实现信号重构的算 法,在分析了模极大值直接重构算法和交替投影算法之后,总结 了两种经典算法存在的优缺点,发现模极大值直接重构算法虽然 简单快捷,但重构后的信号信息不完整,有较大失真;而交替投影 算法虽然能达到很高的重构精度,但计算量太大,运算速度太慢, 实用性不强。
最后采用blocks和bump两种经典算法进行了仿真比较,结果表 明Blackman窗函数插值的去噪算法有较好的去噪效果,信号信息 完整,计算量小,且失真小,收敛速度快,实用性较强。
如何找到一种由信号的模极大值重构信号的合理有效、方便快 捷的算法,就是本文要解决的问题。考虑到Blackman窗函数有主 瓣宽且幅值大,旁瓣宽度小以及衰减速度快的特性,本文以 Blackman窗函数作为插值函数,并在此基础上,构造了一种基于 Blackman窗函数插值的小波模极大值去噪算法,该算法利用插值 技术填补了各尺度上那些不是模极大值点的值,从而得到重构信 号。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
ylabel(strcat('a ',num2str(i)));
subplot(level+1,2,2*(i)+2);
plot(swd(i,: )); axis tight;grid on;
[a,d]=swt(w0,level,Lo_D,Hi_D);
w2=d; % w2为待重建小波
for j=1:num_inter
w2=Py_Pgama(d,wpeak,wframe,1,sr); % 先进行Py投影和 Pgama投影
w0=iswt(pswa,w2,Lo_R,Hi_R); % 再进行Pv投影
snr=20*log10(norm(signal)/norm(err))
function [signal,swa,swd,ddw,wpeak]=wave_peak(points,level,Lo_D,Hi_D,Lo_R,Hi_R,offset)
% 该函数用于读取ecg信号,找到小波变换模极大序列
% signal: 原始信号; swa:小波概貌; swd:小波细节;
% ddw: 局部极大位置; wpeak:小波变换的局部极大序列]
pswa=swa(level,: ); % pswa: 为待重建的信号
wframe=(wpeak~=0);
%迭代初始化
w0=zeros(1,points);
这是用小波变换模极大值重建信号的源程序,数据是一心电信号,在matlab6。5下实现,来源于胡广书的《现代信号处理教程》附属光盘,现提供给大家供大家学习参考,滤波部分可以根据个人情况进行修改。程序包含四个部份,其中wave_peak实现信号的分解并求出模极大值
%--------------------------------------------------------------------------
ylabel(strcat('d ',num2str(i)));
end
%求小波变换的模极大值及其位置
ddw=zeros(size(swd));
pddw=ddw;
nddw=ddw;
posw=swd.*(swd>0);
pdw=((posw(:,1:points-1)-posw(:,2:points))<0);
%所处理数据的长度 分解的级数 抽样率 迭代次数 小波名称
offset=0;
[Lo_D,Hi_D,Lo_R,Hi_R]=wfilters(wf);
%计算小波分解系数和模极大序列
[signal,swa,swd,ddw,wpeak]=wave_peak(points,level,Lo_D,Hi_D,Lo_R,Hi_R,offset);
for i=1:level
subplot(level,1,i);
plot(wpeak(i,: )); axis tight;grid on;
ylabel(strcat('j= ',num2str(i)));
end
[swa,swd] = swt(signal,level,Lo_D,Hi_D);
figure;
subplot(level,1,1); plot(real(signal)); grid on;axis tight;
for i=1:level
subplot(level+1,2,2*(i)+1);
ddw=pddw|nddw;
ddw(:,1)=1;
ddw(:,points)=1;
wpeak=ddw.*swd;
wpeak(:,1)=wpeak(:,1)+1e-10;
wpeak(:,points)=wpeak(:,points)+1e-10;
ቤተ መጻሕፍቲ ባይዱ
%按级给出小波变换模极大的波形
figure;
[a,d]=swt(w0,level,Lo_D,Hi_D); % Pv
end
pswa=iswt(swa(level,: ),w2,Lo_R,Hi_R); % 计算重建信号
% 原信号和由模极大重建信号的比较
figure,
subplot(211)
plot(pswa(1:points));
subplot(212)
plot(signal(1:points),'r');
%分别计算重建小波以及原信号的信噪比
werr=w2-swd;
% 原信号的小波变换(swd)和重建后的小波变换(w2)的比较
figure,
for m=1:level
wsnr(m)=20*log10(norm(swd(m,: ))/norm(werr(m,: )))
subplot(level+1,1,m);
plot(swd(m,: )),hold on,
plot(w2(m,: ),'r');grid on;ylabel(strcat('j=',num2str(m))),axis tight;
end
err=pswa(1:points)-signal(1:points);
pddw(:,2:points-1)=((pdw(:,1:points-2)-pdw(:,2:points-1))>0);
negw=swd.*(swd<0);
ndw=((negw(:,1:points-1)-negw(:,2:points))>0);
nddw(:,2:points-1)=((ndw(:,1:points-2)-ndw(:,2:points-1))>0);
% exa130301.m 例13.3.1: 利用小波变换模极大重建原信号
%--------------------------------------------------------------------------
close all;
points=4096; level=6; sr=360; num_inter=6; wf='db3';
warning off;
ecgdata=load('x.txt' );
plot(ecgdata(1:points)),grid on,axis tight,axis([1,points,-50,300]);
signal=ecgdata(1:points)'+offset;
% 信号的小波变换,按级给出概貌和细节的波形