%% dwtmatrix.m
% usage: Y = dwtmatrix(N,lev,la,ha) 
% Compute the matrix for the multilevel DWT linear map 
% N = length of finite signal
% lev number of levels of the DWT
% la = low pass analysis filter
% ha = high pass analysis filter
% la and ha must have the same length
% the periodic transform is used 
% N must be divisible by 2^lev
% uses the wavelet tool box

% S. Allen Broughton - 12 Nov 2010

function W=dwtmatrix(N,lev,la,ha)
dwtmode('per')
E = eye(N); 
W = zeros(N,N);
 for  j = 1:N 
  W(:,j) = wavedec(E(:,j),lev,la,ha);
 end