波变换

开发平台：

Matlab

upwlev2.m：源码内容

							function [Nc,Ns,a] = upwlev2(c,s,IN3,IN4)
%UPWLEV2 Single-level reconstruction of 2-D wavelet decomposition.
%   [NC,NS,CA] = UPWLEV2(C,S,'wname') performs the single-level
%   reconstruction of wavelet decomposition structure 
%   [C,S] giving the new one [NC,NS], and extracts the last
%   approximation coefficients matrix CA.
%
%   [C,S] is a decomposition at level n = size(S,1)-2, so
%   [NC,NS] is the same decomposition at level n-1 and CA 
%   is the approximation matrix at level n.
%
%   'wname' is a string containing the wavelet name,
%   C is the original wavelet decomposition vector and
%   S the corresponding bookkeeping matrix (for 
%   detailed storage information, see WAVEDEC2).
%
%   Instead of giving the wavelet name, you can give the
%   filters.
%   For [NC,NS,CA] = UPWLEV2(C,S,Lo_R,Hi_R),
%   Lo_R is the reconstruction low-pass filter and
%   Hi_R is the reconstruction high-pass filter.
%
%   See also IDWT2, UPCOEF2, WAVEDEC2.
%   M. Misiti, Y. Misiti, G. Oppenheim, J.M. Poggi 12-Mar-96.
%   Last Revision: 17-Mar-2003.
%   Copyright 1995-2004 The MathWorks, Inc.
% $Revision: 1.10.4.2 $
% Check arguments.
if nargin == 3
    [Lo_R,Hi_R] = wfilters(IN3,'r');
else
    Lo_R = IN3; Hi_R = IN4;
end
% Extract last approximation.
nr   = s(1,1);
nc   = s(1,2);
last = nr*nc;
a    = reshape(c(1:last),nr,nc);
% One step reconstruction of the wavelet decomposition structure.
ns = size(s,1);
if ns>2
    s2 = s(2,:);
    s3 = s(3,:);
    Ns = [s3 ;s(3:ns,:)];
    nr = s2(1,1); nc = s2(1,2); nb = nr*nc;
    first = last+1; last = last+nb;
    h     = reshape(c(first:last),nr,nc);
    first = last+1; last = last+nb;
    v     = reshape(c(first:last),nr,nc);
    first = last+1; last = last+nb;
    d     = reshape(c(first:last),nr,nc);
    ra = idwt2(a,h,v,d,Lo_R,Hi_R,s3);
    Nc = [ra(:)' c(last+1:end)];
else
    Ns = [];
    Nc = [];
end

						