应用 信号分析中小波变换基函数选择研究 晏[摘强 周冬梅 (成都理工大学信息科学与技术学院,四川 成都 610059) 要] 不同的小波基具有不同的时频特征,用不同的小波基分析同一个问题产生的结果会有区别。因此在小波分 析应用中就存在了小波基的选取问题。本文详述了小波变换的原理以及选择小波基函数的基本方法.根据不同的应用需要 来使用相应的小波基函数。 [关键词] 小波基;信号分析;选择;小波特性 1.引言 信号分析的主要目的是寻找一种简单有效的信号变换 方法,使信号所包含的重要特征能显示出来。在信号分析 中,对信号的基本刻画往往采用两种最基本形式,即时域形 式和频域形式。把时间或空间位置作为自变量,而把信号的 某一数字化特征作为因变量描述信号是常用方式。小波变 换作为一个可以很好地完成时频分析的工具进入人们的视 野。而小波分析中小波基的选择又对分析产生举足轻重的 影响。如何能够很好地选择小波基函数,选择的要求等成为 我们研究的话题之一。 2.小波变换 小波变换是由法国科学家Morlet于1984年在进行地 震数据分析时提出的,他在研究地震波的局部特性时,发现 传统的Fourier变换和加窗Fourier变换难以达到其要求,从 而引入小波的概念。 定义:小波是一个满足条件: r L O(t)dt=-0 (1) 的函数(IlJ(t))。通过平移和伸缩可产生一函数ll, , t): .上 . Il, .b(t)=laf I1J( 二ca Jb)atbE R a≠0 (2) 其中tlJ(t)称为基小波或母小波,a为伸缩因子也称为尺 度因子,b为平移因子,式(2)称为基小波生成的连续波。 小波变换定义: ,+∞ .上,+∞ . wf(a,b)=L ift)。 (t)dt=lal J f(1)Il,( )dt (3) 如果将小波变换定义式(3)与加宽Fourier变换定义式 进行比较,不难看出,连续小波 . (t)的作用与加窗Fourier 变换中的函数gw, t)=g(t—r)e 相似,都起到窗函数的作用, 其中参数b与,r都起着平移作用,而本质不同的是a与∞, 从前面对加窗Fourier变换的讨论中,我们知道,在加窗 Fourier变换中,参数‘1)的大小与时频窗口大小有关系,即tO 的变化与时频窗口长宽有关系,但不改变窗口的大小,时频 窗口的大小只与窗函数g(t)有关,当窗函数g(t)确定,时频窗 口大小就随之确定。而在小波变换中,时频窗口的大小则与 作者简介:晏强,男,河南信1 ̄sA.,硕士,研究方向:电路与系统。 a有关。 3.小波基特性及选择条件 3.1小波函数特性 由于小波变换是将原始图像与小波基函数以及尺度函 数进行内积运算,基于小波变换的数据压缩目的就是希望经 小波分解后得到的三个方向的细节分量具有高度的局部相 关性,而整体相关性能最大限度地消除。因此对于同一幅图 像,选择不同的小波基进行分解所得到的数据压缩效果是不 同的。本文即是对小波函数的特征进行研究,分析不同小波 基对图像压缩编码的影响。 小波变换以其优异的时域和频域局部化能力、方向选择 能力和与人眼视觉特性相符的多分辨率分析能力,被广泛 应用于图像压缩领域,并取得了很好的效果。将小波变换用 于图像压缩时,并非所有的小波基都适合图像分解,小波 基帕选择直接影响到整个算法的编码能力、变换的复杂性和 重构图像的质量,因此小波基的选择是图像压缩中的一个 关键问题。对于图像压缩来说理想的小波基应该具有下列性 质: f1)正交性 用正交小波基由多尺度分解得到的各子带数据分别落 在相互正交的LZ(R2)的子空间中,使各子带数据相关性减 小。但是能准确重建的、正交的、线性相位、有限冲击响应滤 波器组是不存在的,此时一般放宽正交性条件为双正交。 (2)紧支性与衰减性 称小波 (t)是紧支的,如果它有紧支集;称小波 (t)是 急衰或急降的,当t—oo时,它快速衰减或按指数规律衰减。 紧支性与衰减性是小波的重要性质,紧支宽度越窄或衰减越 快,小波的局部化特性越好;紧支小波不需做人为的截断,应 用精度很高,但是一个函数不可能在时域和频域都是紧支 的,最多有一个是紧支的,另一个是急衰的。一般希望小波 基能够在时域上具有紧支性。 (31对称性 对称或反对称的尺度函数和小波函数是非常重要的,因 为可以构造紧支的正则小波基,而且具有线性相位。 Daubechies已经证明,除了Haar小波基,不存在对称的紧支 ・-——49--—— 学术探讨应用技术与研究 2『112率第3期 正交小波基。而对于双正交小波基,可以合成具有对称或反 对称的紧支撑小波基。 f4)正则性 (2)小波滤波器长度的选择对称的双正交小波基可以满 足图像处理中比较严格的线性相位特性要求,不仅可减少或 消除重建图像的边缘失真,而且在级联的塔型结构中无须相 正则性是函数光滑程度的一种描述,也是函数频域能量 集中的一种度量。若qs(t)具有N阶连续导数,正则性阶数r 越大,意味着llJ(t)越光滑,其频域的能量越集中,一般而言数 据压缩的结果也越好。 位补偿,同时支集较短便于快速实现和进行边界处理,因此 在图像编码领域得到广泛应用。对于非平稳图像信号,采用 非对称滤波器(AFB)比传统滤波器(如FIR)图像重建后的 视觉效果要好。小波基与待压缩图像的相似性对压缩效果 的影响实验证明,在小波函数基本图像与压缩图像的结构较 相似的情况下,可以忽略正则性,结构越相似压缩效果越 好。 f5)消失矩 小波基函数应具有消失矩性质,即: J tko(t)dt=O:k=O,1,……n一1 J・ (3)小波变换的级数由一维小波采用张量积构成的可分 离的二维小波变换是将原始图像分解成一个低频信号和三 个方向的高频分量信号。即每层分解为四个子带信号,每层 的低频信号又可以进一步分解成四个子带信号,总的子带数 为3K+1,其中K为分解层数。 D凸 其中 (t)有n个连续的零点,我们说 (t)的消失矩为 k。当k=0时有f t (t)dt=O,这表明 (t)是一个迅速衰减 且平均值为0的波。消失矩的大小决定了用小波逼近光滑 函数的收敛率。当图像光滑时,消失矩越大,将导致小波系 数越小。 4.小波变换未来发展 小波应用的范围虽广,但真正取得极佳效果的领域并 3.2小波函数选择条件 在不同的应用领域,小波基的选取标准不同,不同的小 不多,人们也正在挖掘有前景的应用领域。小波分析打开 了信号处理的大门,这个领域远比Fourier分析处理的时不 波基适应不同的具体情况。小波变换中的小波基的选择转 换为正交镜像滤波器组QMFB的选择。小波基的选取应从 一变系统复杂。在这个大领域里,小波分析是一个重要工具, 同时也需要其它的理论和工具。 参考文献: 【1】魏明果.实用小波分析[M】.北京:北京理工大学出版社 2005. 般原则和具体应用两方面考虑。根据小波函数特征,一般 (1)小波基的正则性阶数对图像压缩效果有影响 正则性是函数光滑程度的一种描述,也是函数频域能量 原则如下: 集中的一种度量。在小波变换过程中,对于正则性较弱的滤 波器,输入信号光滑部分会因变换级数的增加而发生不连 续的情况,这种非信号本身的不连续性也将反映在图像压 [2】崔锦泰,程正兴译.小波分析导论[M].西安:西安交通大学出 版社.1995. 缩时的离散小波的变换系数中,量化过程使得这些系数出 现误差,导致图像重建失真,正则性越弱,则图像变化部分 越不光滑,视觉效果越差,反之小波函数正则性越强,图像 重建效果越好。小波正交性、能量集中性和图像压缩小波基 恰当选择使得信号经小波正交变换后,其能量大部分集中 在其中k个变换系数上,若k越小,则表明该变换能量越集 中。在图像压缩中,为提高压缩比可只保留权值较大的前n 个分量。可见,能量集中特性越好,图像压缩效果越好。 [3】谭善文,秦树人,汤宝平.小波基时频特性及其在分析突变信 号中的应用[J].重庆大学学报2001,24(2):12—17. 【4]陈志奎,徐铭陶.小波基和它对一维信号的直接分解[J].重庆 大学学报,1998,21(3):131-138. [5]杨福生.小波变换的工程分析与应用[M].北京:科学出版 社.1999. [6]王建中.小波理论及其在物理和工程中的应用….数学进展, 1992,21(3):289—311. [7】邓东皋,彭立中.小渡分析[J].数学进展,1991,20(1):294-310. The Selection of Basis Functions of Wavelet Transform in Signal Analysis Yan Qiang Zhou Dongmei (Chengdu University ofTechnology,Chengdu 610059,Sichuan) 【Almr ̄】Different wavelet basis has diferent frequency characteristics.The result of a question will be diferent with different wavelet basis.Therefore,the application of wavelet analysis on the selection of wavelet basis is critica1.This paper details the principle of wavelet transform and chooses the wavelet basis function,depending on the application needs to use the corresponding wavelet basis function. 【l【‘ 喇恤】wavelet basis;signal analysis;selection;wavelet features 一50—