图像NMF分解MATLAB代码

function NMF

clear all;
close all;
clc;

X=double（imread（‘lena.bmp‘））;
imshow（mat2gray（X））;

[n m]=size（X）;                                    %计算X的规格
r = （n*m）/（m+n）                                  %设置分解矩阵的秩
U=rand（nr）;                            %初始化UV，为非负数