www.pudn.com > spiht-0.3.rar > test97.m

fid=fopen('lena512.raw','rb'); 
a=fread(fid,[512,512]); 
a=a'; 
dwt_img=d2fldwt9_7(a,5); 
imshow(dwt_img,[]) 
img=d2bldwt9_7(dwt_img,5); 
figure 
imshow(img,[]) 
figure 
d=reshape(a-img,1,512*512); 
hist(d,100)