非负矩阵分解 matlab

本程序为非负矩阵分解，适用于高光谱解混。图形图像处理方面的 matlab程序

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