matlab练习程序（非负矩阵分解）

clear all;
close all;
clc;

V=double(imread('lena.jpg'));
imshow(mat2gray(V));

[i u]=size(V);                                    %计算V的规格
r=100;                                  %设置分解矩阵的秩
W=rand(i,r);                            %初始化WH，为非负数
H=rand(r,u);
maviter=100;                                    %最大迭代次数
for iter=1:maviter
W=W.*((V./(W*H))*H');           %注意这里的三个公式和文中的是对应的
W=W./(ones(i,1)*sum(W));
H=H.*(W'*(V./(W*H)));
end

img_V=W*H;
figure;
imshow(mat2gray(img_V));