function?[??HB?]?=?nmfdb（?X?r??）
%NMFDB?Summary?of?this?function?goes?here
%???Detailed?explanation?goes?here
dim=size（X）;????????????????????????????????????%計算x的規格
X=double（X）;
B=10*rand（dim（1）r）;????????????????????????????%初始化BH，為非負數
%B=B./（ones（dim（1）1）*sum（B））;???????????????????%歸一化B的每一列

H=10*rand（rdim（2））;
?total=sum（H2）;
for?i=1:4
????????H（i:）=H（i:）./total（i1）;
end
maxiter=2000;????????????????????????????????????%最大迭代次數
for?iter=1:maxiter
????H=H.*（B‘*X./（B‘*B*H））;
????B=B.*（X*H‘./（B*H*H‘））;
????%B=B./（ones（dim（1）1）*sum（B））;
????%total=sum（sum（sum（B）））;
????%B=B./total;
????total=sum（H2）;
????for?i=1:4
????????H（i:）=H（i:）./total（i1）;
????end
end

end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????740??2014-03-29?15:05??nmfdb.m