MATLAB在小波理论学习中应用.docVIP

MATLAB在小波理论学习中应用 　　摘要 在信号分析领域,与一些经典技术相比,小波变换提供了关于信号数据的更加准确的信息。小波工具箱是在MATLAB技术计算环境基础上建立的一个与小波有关的函数的集合,它提供了在MATALB框架内使用小波和小波包进行信号和图像的分析与综合的工具。为了有助于充分利用该工具箱,也为了介绍商业科学软件的开发理念,在深入分析其目录结构和源代码的基础上,研究该工具箱的体系结构和实现技术。 　　关键词 小波工具箱;小波滤波器; 面向对象编程 　　中图分类号:TP319 文献标识码:B 文章编号:1671-489X(2010)33-0150-02 　　Application of MATLAB in Learning of Wavelet Theory//Wang Xin 　　Abstract In the field of signal analysis, wavelet transform provides more precise information about the signal data than some classic techniques. Wavelet toolbox is a set of functions related to wavelet, which are based on the MATLAB technical computing environment. It provides tools for the analysis and synthesis of signals and images in the MATLAB framework. In order to take full advantage of the toolbox and to introduce the development philosophy of commercial scientific software, the architecture and implementation of wavelet toolbox are studied by analyzing its directory structures and source codes. 　　Key words wavelet toolbox; wavelet filter; object-oriented programming 　　Author’s address School of Internet of Things Engineering, Jiangnan University, Wuxi, Jiangsu, China 214122 　　 　　1 引言 　　小波工具箱使用图形化工具和命令行函数扩展MATLAB技术的计算环境,可以用于开发基于小波的算法,进行信号和图像的分析、综合、去噪和压缩等运算[1]。 　　小波工具箱支持对小波属性和小波应用的交互式研究。这对于语音和音频处理、图像和视频处理、生物成像以及通信和地球物理学中的一维和二维应用是非常有益的[2]。 　　通过对小波工具箱的目录结构和源代码的深入分析,研究小波工具箱的体系结构和实现技术。这样的实践有助于学生对小波理论的深入理解。 　　2 小波滤波器的计算 　　从小波理论可以知道,如果尺度函数φ是紧支的,那么对应的尺度滤波器W的长度是有限的,就可以被看做一个FIR滤波器[3]。小波工具箱将这个FIR滤波器保存为系统预定义的数据文件。载入这个数据文件,即可获得尺度滤波器W。从尺度滤波器W出发,定义4个FIR滤波器,如表1所示。 　　从计算小波滤波器的wfilters函数的源代码可知,上述4个滤波器的计算方案如图1所示。首先,对W进行归一化,即可求得Lo_R。从Lo_R出发,使用qmf函数和wrev函数,即可求得Hi_R、Lo_D和Hi_D。其中,qmf函数使得Hi_R和Lo_R是正交镜像滤波器,wrev的作用是翻转Hi_R和Lo_R的系数。因此,Hi_D和Lo_D也是正交镜像滤波器。 　　在完成小波函数ψ(t)、尺度函数φ(t)或者小波滤波器的定义的基础上,即可进行各种小波分析的任务:连续小波分析、离散小波分析和小波包分析。 　　3 离散小波分析 　　将小波理论应用于实际的工程问题中,使用较多的是离散小波变换。在小波工具箱中,离散小波变换的分析-分解函数主要有dwt函数(单级离散小波变换)和wavedec函数(多级小波分解),重构-综合函数主要有idwt函数(单级逆离散小波变换)、waverec函数(多级小波重构)和wrcoef函数(从小波系数重构单个分解分支)。这些函数之间的计算关系如图2和图3所示。 　　图2