很好的一个小波变换图像压缩程序

可以压缩达到12-13倍以上，相信对你会很有帮助的

function [x1phix2psi] = cascade（hgdx）
% Wavelet cascade algorithm

c = h{1}*2/sum（h{1}）;
x = 0:dx:length（c） - 1;
x1 = x - h{2};
phi0 = 1 - abs（linspace（-11length（x）））.‘;

ii = []; jj = []; s = [];

for k = 1:length（c）
   xk = 2*x - （k-1）;
   i = find（xk >= 0 & xk <= length（c） - 1）;
   ii = [iii];
   jj = [jjfloor（xk（i）/dx）+1];
   s = [sc（k）+zeros（size（i））];
end

% Construct a sparse linear operator that iterates the dilation equation
Dilation = sparse（iijjslength（x）length（x））;

for N = 1:30
   phi = Dilation*phi0;
   if norm（phi - phi0inf） < 1e-5 break; end
   phi0 = phi;
end

if norm（phi） == 0
   phi = ones（size（phi））*sqrt（2）;   % Special case for Haar scaling function
else
   phi = phi/（norm（phi）*sqrt（dx））;  % Rescale result
end

if nargout > 2
   phi2 = phi（1:2:end）;  % phi2 is approximately phi（2x）

   if length（c） == 2
      L = length（phi2）;
   else
      L = ceil（0.5/dx）;
   end

   % Construct psi from translates of phi2
   c = g{1};
   psi = zeros（length（phi2）+L*（length（c）-1）1）;
   x2 = （0:length（psi）-1）*dx - g{2} - 0*h{2}/2;

   for k = 1:length（c）
      i = （1:length（phi2）） + L*（k-1）;
      psi（i） = psi（i） + c（k）*phi2;
   end
end
return;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1235  2008-12-22 10:11  wavelets\cascade.m

     文件       1077  2008-12-28 08:45  wavelets\dec.asv

     文件       1159  2008-12-28 08:50  wavelets\dec.m

     文件        859  2008-12-27 20:09  wavelets\decode.asv

     文件        218  2008-12-27 10:43  wavelets\double2uint.asv

     文件        251  2008-12-27 11:02  wavelets\double2uint.m

     文件        137  2008-12-27 15:36  wavelets\Dvalue.asv

     文件        305  2008-12-27 16:24  wavelets\Dvalue.m

     文件       1211  2008-12-22 10:11  wavelets\filterstr.m

     文件       9779  2008-12-22 10:08  wavelets\getwavelet.m

     文件     262874  2008-12-22 09:01  wavelets\lena.bmp

     文件        446  2008-12-22 10:09  wavelets\lp_lift.m

     文件        985  2008-12-28 09:05  wavelets\main.asv

     文件        873  2008-12-28 09:17  wavelets\main.m

     文件        243  2008-12-22 10:11  wavelets\numvanish.m

     文件       3498  2008-12-27 13:52  wavelets\quantize.asv

     文件       3846  2008-12-27 15:29  wavelets\quantize.m

     文件        206  2008-12-27 16:06  wavelets\Recover_Dvalue.asv

     文件        214  2008-12-27 16:21  wavelets\Recover_Dvalue.m

     文件       1521  2008-12-27 21:47  wavelets\RLC.asv

     文件       1699  2008-12-28 00:23  wavelets\RLC.m

     文件         31  2008-12-27 16:42  wavelets\select.asv

     文件        230  2008-12-29 08:50  wavelets\select.m

     文件        979  2008-12-22 10:08  wavelets\seqdual.m

     文件         16  2008-12-27 16:50  wavelets\shiyan.m

     文件         82  2008-12-27 10:58  wavelets\uint2double.asv

     文件         82  2008-12-27 10:58  wavelets\uint2double.m

     文件         74  2008-12-27 21:49  wavelets\Untitled.m

     文件      12764  2008-12-24 12:35  wavelets\wavelet.asv

     文件      12764  2008-12-23 07:52  wavelets\wavelet.m

............此处省略5个文件信息