function y=WT(level, wfil, x);
% WT(l, wfil, x) performs an l level wavelet decomposition on data x using
% wavelet filters wfil.  Acts columnwise on x.

[datalength_1,datalength_2]=size(x);

y=x;
for j=0:level-1
    
    m=datalength_1/2^j;
    y(1:m,:)=wavedec(wfil,y(1:m,:));
    
end