ploterror.m

非负矩阵分解的matlab代码

function ok = plotError(error,indx,dim)
%08.08.2006, Bjarni Bodvarsson ([e-mail])

%dim = [x y slice frame];
x = dim(1);
y = dim(2);
slice = dim(3);
frame = dim(4);

rowsColumns = ceil(sqrt(slice));
imgTmp = zeros(x*rowsColumns,y*rowsColumns);

n=1;
sourceImg = zeros(x,y,slice);
sourceImg(indx) = error;

for j=1:rowsColumns
    for i=1:rowsColumns
        if n>slice;
            imgTmp((j*y-y+1):j*y,(i*x-x+1):i*x) = zeros(y,x);
        else
            imgTmp((j*y-y+1):j*y,(i*x-x+1):i*x) = rot90(sourceImg(:,:,n));
        end
        n = n+1;
    end
end
figure,
imagesc(imgTmp), axis equal, ylim([0 rowsColumns*y]), xlim([0 rowsColumns*x]);
title('Average absolute error')