使用matlab设计小波变换程序中的若干问题

《使用matlab设计小波变换程序中的若干问题》由会员上传分享，免费在线阅读，更多相关内容在学术论文-天天文库。

1、使用MATLAB设计小波变换程序中的若干问题在使用MATLAB完成小波变换程序和通过阈值来压缩图像的过程中，我和许多同学都是边学边用，是从一个接一个的问题中逐步理解小波和MATLAB编写程序的。因此我愿意就个人遇到和解决问题的经验与大家讨论，希望能够对遇到同样问题的人有所帮助。在清华大学林福宗老师倡导的网上互动的学习方式中，老师同学的开诚布公的讨论，尤其是林老师启发大家对出现‘问题’采取的态度和做法，对我今后成长为一名合格的清华的研究生意义重大，谨以此文表示对他指导关心的敬意！1.内容简介本文分

2、为三部分:如何使用MATLAB设计小波标准与标准分解；如何完成使用小波变换压缩图像；仍需探讨的问题。每一部分主要以问题和例子的形式讲述，为了便于参照，附录部分给出部分源代码供大家参考指正。我个人在完成作业的时候，走了许多弯路，最后，反复读老师的第3章讲义的第3--5节，明白了小波变换的一些非常基础性的知识，正是如此，我主要是按照第三章的例子，做出使用MATLAB进行Haar变换的实例，这至少对完成非标准Haar小波的3级分解与合成没有任何问题，而且计算速度远远比使用一维变换快，甚至超过dwt2和

3、wavedec2。再者，是通过这些简短的例子的学习，也能从实践的角度理解小波变换的基础知识和道理，我觉得掌握小波的知识要比简单使用dwt2直接分析出结果更重要。为了让使用一维变换或打算采用卷积完成分解重构程序的同学有所参考，也给出了我以前写的代码和设想供有兴趣者参考。在图像压缩的任务中，我个人建议采用wdencmp，原因是简单可靠，能够生成老师要求小波压缩与重构的演示及PNG文件。在最后部分着重对使用小波压缩之后重构图像的PNG文件会变大问题进行了简单的猜测性的解释和讨论，希望感兴趣的同学能深入

4、研究。2.如何使用MATLAB设计小波标准与标准分解1.使用Haar和Db9小波编写图像文件的标准和非标准分解重构程序实质上，任务书中规定首先要编写４种程序之一：哈尔（haar）小波的标准分解重构程序；Daubechies9小波的标准分解重构程序；哈尔（haar）小波的非标准分解重构程序；Daubechies9小波的非标准分解重构程序。如果采用标准方法，需要生成图３－２５的系列图像。如果采用非标准方法则需要生成图３－２６的系列图像。我个人编写了非标准的分解重构程序，并总结网上同学们公开的方法，认

5、为可以有３种方法实现。但前提是必须仔细阅读林老师推荐的补充教材第三章“小波与小波变换”的3.3、3.4和3.5节（１４－２７页）。如果能使用MATLAB简单地实践一下教材的例子，至少完成哈尔小波的标准与非标准程序相当简单。为此我主要介绍如何实现老师教材给我们的例子是如何在MATLAB中实现的。然后，再引申到使用MATLAB的函数实现方法。为了描述问题简便，使用黑体表示需在MATLAB的命令窗口（commandwindow）的输入部分。如果不熟悉下面出现的函数功能和使用方法，请再命令窗口中使用he

6、lp函数名，如helpdwt2，能够得到英文的功能和用法说明。2.哈尔小波变换的MATLAB实例先介绍MATLAB的矩阵(Matrix)和向量(Vector)的赋值方法：在commandwindow的>>符号下输入：12一维向量：　Ｉ=[9735]　　　　　Ｆ=[2,5,8,9,7,4,-1,1]二维矩阵：A=[642361606757955541213515016174746202143422440262737363031333234352928383925412322444519184849

7、15145253111056858595462631]　Ｉ和Ｆ现在就是［例3.1］［例3.2］一维向量,　A就是3.5.1中的图像矩阵。在进行Haar小波的一维与二维变换（分解）之前，首先介绍一下MATLAB的矩阵运算。MATLAB是MatrixLaboratory的合写，它对矩阵运算之功能堪称一流。由于使用矩阵描述问题更象数学表达式，所以编写的程序不仅高效，更易读。所以MATLAB程序应该尽量使用矩阵直接描述。如M是一4X4的矩阵，则B=I*M就完成了使用两个for循环编写乘法的程序。分析老师

8、的［例3.1］，那么可以得到如下结果：[（9+7）*1/2（3+5）*1/2（9-7）*1/2（3-5）*1/2]如果把其看作矩阵方式的乘法，那么令M为其表述求取平均和差值的系数矩阵，则输入：M=[1/201/20121/20-1/2001/201/201/20-1/2]把上述的M输入到MATLAB中,然后使用：C=I*M得到的结果就是[841-1]。这就是非规范化Haar小波的第一级分解系数。注意C的前两项[84]就是近似系数（ApproximationCoefficients）,后两项[1_