% 使用Haar小波，得到相應的提升方案
lshaar = liftwave('haar');
% 添加ELS 到提升方案
els = {'p',[-0.125 0.125],0};
lsnew = addlift(lshaar,els);
% 對于簡單圖像，尺度為2進行LWT
x = reshape(1:16,4,4);
xDec = lwt2(x,lsnew,2)
% 提取第1層的低頻系數
ca1 = lwtcoef2('ca',xDec,lsnew,2,1)
% 重構低頻和高頻
a1 = lwtcoef2('a',xDec,lsnew,2,1)
a2 = lwtcoef2('a',xDec,lsnew,2,2)
h1 = lwtcoef2('h',xDec,lsnew,2,1)
v1 = lwtcoef2('v',xDec,lsnew,2,1)
d1 = lwtcoef2('d',xDec,lsnew,2,1)
h2 = lwtcoef2('h',xDec,lsnew,2,2)
v2 = lwtcoef2('v',xDec,lsnew,2,2)
d2 = lwtcoef2('d',xDec,lsnew,2,2)
% 檢查重構效果
err = max(max(abs(x-a2-h2-v2-d2-h1-v1-d1)))