时频分析作业.docVIP

时频分析作业 学 院：电气工程与信息工程学院 专 业： 检测技术与自动化装置 姓 名： 陈 彦 峰 学 号： 082081102005 指导教师： 周 红 莉 小波（wavelet)是定义在有限间隔且平均值为0的函数，小波函数多以开发者名字命名，如图1所示： 图1 部分小波 众所周知，傅立叶分析是把一个信号分解成各种不同频率的正弦波，因此正弦波是傅立叶变换的基函数。同样，小波分析是把一个信号分解成由原始小波经过移位和缩放后的一系列小波，因此小波是小波变换的基函数，即小波可用作表示一些函数的基函数。 小波是近十几年才发展并迅速应用到图像和语音分析等众多领域的数学工具，是继110多年前建立傅立叶（Joseph Fourier）分析之后的一个重大突破。经过十几年的努力，小波理论基础已经基本建立并成为应用数学的一个新领域，引起了众多数学家和工程技术人员的极大关注，是国际上科技学术界高度关注的前沿领域。本文试图从工程和实验角度出发，较为直观地探讨小波变换在图像压缩中的应用。 2、小波变换和重构 小波变换的基本思想是用一组小波或基函数表示一个函数或信号，例如图像信号。以哈尔（Haar）小波基函数为例，基本哈尔小波函数（Haar wavelet function）定义如下： Ψ(x) = 1, 当0≤x1/2 Ψ(x) = -1, 当1/2≤x1 Ψ(x) = 0, 其他 设有一幅分辨率只有4个像素的一维图像，对应像素值为：[9 7 3 5]。用哈尔小波变换的过程是：计算相邻像素对的平均值（averaging，亦可称之为近似值approximation），得到一幅分辨率为原图像1/2的新图像：[8 4]。这时图像信息已部分丢失，为了能从2个像素组成的图像重构出4个像素的原图像，必须把每个像素对的第一个像素值减这个像素的平均值作为图像的细节系数（detail coefficient）保存。因此，原图像可用下面的两个平均值和两个细节系数表示：[8 4 1 -1]。可以把第一步变换得到的图像进一步变换，原图像两级变换的过程如表1所示： 表1 哈尔小波变换过程 分辨率 平均值 细节系数 4 [9 7 3 5] ? 2 [8 4] [1 -1] 1 [6] [2] ? 哈尔变换过程事实上是用求均值和差值的方法对函数或图像进行分解，对于f(x)=[9 7 3 5]，我们可作最多2层的分解。对于2维图像，同样可以用依次对行列进行小波变换得到2维图像的分解。这时经过一次小波变换得到是2维图像的近似值(CA)以及水平(CH)、垂直(CV)和对角(CD)细节分量值。显然，从2维图像的CA、CH、CV和CD值可以重构出原来的2维图像。 3、图像压缩 事实上，去掉某些经过小波变换得到的细节分量值对重构图像的质量影响不大。具体的做法是设置一个阈值δ，例如把≤δ的经小波变换得到的水平(CH)、垂直(CV)和对角(CD)细节分量值细节分量值置为0，进行行变换和列变换，这样就实现了图像压缩。 例1. 压缩前图像X的大小： Name Size Bytes Class X 256x256 524288 double array Grand total is 65536 elements using 524288 bytes 第一次压缩图像的大小为： Name Size Bytes Class ca1 135x135 145800 double array Grand total is 18225 elements using 145800 bytes 第二次压缩图像的大小为： Name Size Bytes Class ca2 75x75 45000 double array Grand total is 5625 elements using 45000 bytes 程序如下： %装入图像 load wbarb; %显示图像 subplot(221);image(X);colormap(map) title(原始图像); axis square disp(压缩前图像X的大小：); whos(X) %对图