小波分解矩阵Matlab.doc

文档介绍：%----------------------------------------------------------%
小波图像分解 Matlab 程序-
http://blog./chenyusiyuan/archive/2008/06/05/
 
小波图像重构 Matlab 程序-
http://blog./chenyusiyuan/archive/2008/06/05/
%----------------------------------------------------------%
 
 
%----------------------------------------------------------% 
     小波分解矩阵 Matlab 程序-
%----------------------------------------------------------%
 
function [coef,scf]=mywavedec2(x,N,wname)
%----------------------------------------------------------%
% 函数 MYWAVEDEC2() 对输入矩阵 x 进行 dim 层分解,得到相应的分解系数矩阵 y
% 输入参数:x ——输入矩阵
%          N ——分解级数
%          wname ——分解所用的小波函数
% 输出参数:scf ——存储各级分解系数矩阵的大小以及原始(图像)矩阵的大小
%           coef ——分解系数矩阵,其结构如下:
% coef = {cA_N;cV_N;cH_N;cD_N;cV_N-1;cH_N-1;cD_N-1;……;cV_1;cH_1;cD_1}
%
% Copyright by Zou Yuhua ( chenyusiyuan )
% Version: , Date: 2008-07-08
%----------------------------------------------------------%
 
 
% 求出小波函数的滤波器组系数向量
[Lo_D,Hi_D] = wfilters(wname,'d');
lf=length(Lo_D);
% 画出原始图像
imshow(x);title('Original Image');
% 标明图像大小
[r,c]=size(x);
xlabel(['Size : ',num2str(r),'*',num2str(c)]);
 
 
% 将矩阵x的数据格式转换为适合数值处理的double格式
xd=double(x);              
[rx,cx]=size(x);
[o2sa,f1sa,rsx]=sizcoef([rx,cx],lf,N);  a=[o2sa,f1sa,rsx]
 
coef=[];
scf=[rx,cx];
for i=1:N
    [cA,cV,cH,cD]=mydwt2(xd,Lo_D,Hi_D);     %