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，為非負數