信号去噪方法综述

信号去噪方法综述

【摘要】在信号传输过程中往往会因为噪声的干扰而影响信号的质量，为了改善这种情况，往往需要对信号进行噪声处理。本文对空域相关法，阈值法等与小波相关的典型算法进行了论述，并将其和传统的滤波器法进行对比，总结出了这些方法在信号去噪方面的优缺点。

【关键词】小波；阈值；空域；信号去噪

The Summarization of signal denoising methods

A bstract:In the process of signal transmission because of the interference of noise ,the quality of the signal often be affected. in order to improve the situation, We need to dispose the noise that mixed in the signal.In this paper ,several typical methods are introduced ,including the spatial filtering method,the threshold method and so on.Those methods were compared with the traditional filter method and the advantages and disadvantages in these methods are summarized in this paper.

Key words: The wavelet ;The threshold value; Airspace;Signal denoising

引言

如何获得一个高质量的信号是信号处理领域一个孜孜不倦的研究方向，而人们在这一领域也取得了巨大的成就。长久以来，人们用傅里叶变换对信号进行相关的处理，并且也取得了一系列的成就。但是，一种方法并不能在任何情况下都适用，傅里叶变换在信号去噪方面也有很多的局限性。其中傅里叶变换在处理这类问题时的一个缺陷就是，用傅里叶进行分析时，它的构造函数是周期性的正弦波和余弦波[1]。鉴于其局限性，它只适合对那些具有周期性或者是具有近似周期性的信号进行滤波或压缩，而在对那些具有非周期或者局部特征很明显的信号的处理上效果就不是很好。

虽然傅里叶变换在信号去噪方面存在局限性，但是由其发展来的小波变换则能很好地解决上述问题。作为在信号处理领域中的一种新的分析方法，它不仅保留了傅里叶变换的许多优点，而且在原来的基础上进行了改进和发展，使其能够在时频域对信号进行处理。小波变换的显著特优点是通过变换可以将信号进行更细微的处理，并且能够将信号的某些特征较好的表现出来，实现了在时频域对信号进行局部化、多尺度的分析的要求。在小波基础上发展来的信号去噪方法表现出了良好的去噪效果，是Fourier变换在信号处理领域的完善和发展。1小波基础知识

1.1小波变换原理

定义1：)

(

)(2R

L

t

f∈

∀平方平方可积空间，连续小波变换为：

dt

a

t

t

f

a

a

W R

R

f)

(

)

(

1

)

,

(⎰--

=

τ

ψ

τ（1）

其中：)

,

(τ

a

w f是小波变换系数；)

(t

ψ是小波函数。

离散小波变换式定义为：

)

2(

)

(

2

)

,

(

1

2k

n

n

f

k

j

W j

N

n

j

f-

=-

-

=

-

∑ψ（2）其中，)

,

(k

j

W f表示小波系数，N是采样点数，j为分解层数。

在使用小波对信号进行处理的过程中，任何一个信有效信号都可以用下式来表示：

∑∑

∑

=∈

∈

+

=

j

m z

k

k

m

f

z

k

k

j

f k

m

W

t

k

j

A

t

f

1

,

,)

,

(

)(

)

,

(

)(ψ

φ

(3)

其中，f(t)是原信号，)

,

(k

j

A f表示尺度系数，)

,

(k

j

W f表示小波系数。

1.2多分辨率分析

定义2：令j V，j=…，-2，-1,0,1,2，…为

)(2R L 中的一函数子空间序列。如果能够满足以下条件，

（嵌套性）1+⊂j j V V （稠密性）)(2R L V j =

（分立性）}

0{=j V

（尺度性）j V x f ∈)(，当且仅当

0)2(V x f j

∈-（Vj 是近似空间）

那么},{Z j V j ∈称为依φ的多分辨率分析，φ是尺度空间。

多分辨率分析的原理是：把全空间

)(2

R L 依分辨率j

2分解成一系列嵌套的闭的子空间序列},{Z j V j ∈，再根据正交补分解，

将)(2

R L 分解为一系列正交的小波子空间},{Z j W j ∈，最后将信号)()(2R L t f ∈投影

分解到不同分辨率的小波子空间上来分析[2]

。

2信号去噪的方法

2.1 滤波器滤波法

在信号传输过程中，如果信号和噪声的频谱不在相同的频段，我们经常采用滤波的方法将噪声去除。这是一种在频域对信号进行处理的方法，在使用这种方法进行信号去噪时，我们要根据信号与噪声的频带特点来选取合理的滤波器来进行去噪。该方法结构简单、易于操作，只需选择合适的滤波器让含噪信号通过即可，滤波效果也很好。但是该方法的使用是有条件的，因为该方法是在频域进行的，因此只有在有效信号和噪声的频谱没有重叠时才可以把信号和噪声分开。但是在实际中我们遇到的往往是信号和噪声的频谱相差不大或重叠，甚至是噪声频谱分布在整个频域中。在这种情况下我们很难使用滤波的方法将它们彼此分开，亦或是以损失有效信号为代价，所以这种方法有很大的局限性。

2.2 基于小波的信号去噪方法

2.2.1 小波的分解与重构

根据Mallat 多分辨率分析的思想和小波分解与重构的算法[12]可以得到：

若信号f(t)的离散采样数据是fk,fk=c0,k,对信号f(t)进行正交小波变换之后得到：

⎪⎩⎪⎨

⎧==∑∑----n k

n n j k i n k n n j k j g d d h c c 2,1,2,1,

其中,k i c ,为尺度系数;k j d ,为小波系

数;h 、g 表示正交镜像滤波器组;j 为分解层数;N 为离散采样点数。

因为小波的分解和重构过程是互逆的，则重构的数学运算可以表示为：

∑∑---+=n

n

k n j n

n k n j n j g d h c c 2,2,,1

（5）

在对信号进行小波分解与重构的过程中首先要对其进行初始化，小波的初始化包含两个方面。第一是找出小波的近似空间j V ,使它可以最好的反应分解信号f 的各种信息。然后再选择一个fj ，并且 j j V f ∈这样就可以最好的向信号f 逼近。初始化之后就可以使

j f f ≈。然后可以把j f 再分解为两部分，分

别是较低级的近似部分

11--∈j j V f 和小波部

分11--∈j j W ω，即

11--+=j j j f f ω。然后再

对其他较低级别采用相同的步骤进行分解，直到分解到第0级，算法流程图如下：

图1小波分解流程图

在用小波对信号进行分解时，待处理信号可以用父函数和母函数来表示。这样我们可以从不同的尺度和频带上对信号进行处理。从而可以将含有噪声的频带置零而将有效信号的频带保留，最后再把处理过的信号进行重构，恢复出原始信号。

重构的过程和上图相反再此不详细赘述。

2.2.2小波分析的空域相关法去噪

在实际中，信号的不连续点往往具有较好的局部特性，如果一个信号的边缘是正常的，则它的李氏指数不小于零。当用小波的相关方法对信号进行处理之后，原信号和噪声的幅值与分解尺度会呈现出相反的变化趋势，具体是原信号与之成正相关，噪声则是负相关。这样以来，与噪声相关点的模极值点可能会随着分解尺度的增加而消失。小波变换之后，有效信号的边缘特征会呈现出很强的相关性，但是噪

K=1,2,…，N-1 （4）