小波分解的Mallat快速算法

Mallat快速算法的图像分解的Matlab代码

function Y=mallatdec2（Xwnamelevel）
%输入:X????? 载入的二维图像像数值; 
%???? level? 小波分解次（级）数设定值（如果设定值超过最高可分解次数按最高分解次数分解）
%????? wname? 小波名字wavelet name
%输出：Y???? 多极小波分解后的小波系数矩阵
%wfilters：小波滤波器

[hg]=wfilters（wname‘d‘）;    %hg分别为分解的低通和高通滤波器 
X=double（X）;
t=1;
hh=size（X2）;
Y=zeros（1hh）;
while t<=level   %循环
    %先进行行小波变换
    for row=1:hh
        Y（row1:hh）=mdec1（X（row1:hh）hg） ;
    end
    %再进行列小波变换
    for col=1:hh
        temp=mdec1（ Y（1:hhcol）‘hg）; %Y（1:hhcol）为矩阵中的一列，转置后为行向量，得到的结果temp也是行向量
        Y（1:hhcol）=temp‘;             %将行向量转置为列向量存储在Y中
    end
    t=t+1;                             %循环次数加1，大于level，循环终止
    hh=hh/2;                           %分解的矩阵的大小为原来的四分之一
    X=Y;                               %将Y值赋给X等待下一次分解                      
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2012-09-12 14:06  Toolbox_Mallat\
     文件         320  2012-09-12 14:09  Toolbox_Mallat\Mallat_fusion_main.m
     文件         346  2012-06-21 09:39  Toolbox_Mallat\Untitled.m
     文件       66616  2009-04-29 10:12  Toolbox_Mallat\lena.BMP
     文件         970  2012-05-31 20:11  Toolbox_Mallat\mallatdec2.m
     文件         607  2012-05-16 14:07  Toolbox_Mallat\mallatrec2.m
     文件         891  2012-06-18 23:54  Toolbox_Mallat\mdec1.m
     文件         369  2012-05-31 20:17  Toolbox_Mallat\mrec1.m
     文件         586  2012-06-21 09:38  Toolbox_Mallat\xdec1.m