?[Amap]=imread（‘E:\Hydrangeas.jpg‘）;
?A=rgb2gray（A）;
?wpt=wpdec2（A3‘Haar‘‘shannon‘）;
?s00=wprcoef（wpt[00]）;
?s00=s00./255;
?M?=?size（s001）;
?N?=?size（s002）;
?
%----------將矩陣轉(zhuǎn)化為非負矩陣-------------%
?X?=?min（s00）;
?m?=?min（X）?
?F?=?ones（MN）?*（-m）;
?s00?=?s00?+?F;
?X?=?min（s00）;
?m?=?min（X）
%------------------------------------------%

D?=?s00;
r?=?3;
W?=?abs（rand（Mr））;
H?=?abs（rand（rN））;
%WH?=?W?*?H;
%V?=（D-WH）.^2;
%e=sum（（sum（V））‘）;

%while?（e>=?0.001）
for??iterator?=?1?:1000
????WD?=?W‘*D;
????WWH?=?W‘*W*H;
????
????for?a?=?1?:?r
?????????for?u?=?1:N
?????????????H（au）?=?H（au）?*WD（au）/WWH（au）;