NMF matlab代码

 [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）;
 
%----------将矩阵转化为非负矩阵-------------%
 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）;