function demo2_1sl_HOSVD
% 读取数据
file=load('Peppers.mat');
X=file.XData;

% 截断值60，HO-SVD分解
[G,U]=HO_SVD(X,60);
% 使用模-n积还原拟合张量
% Gx_1U_1 x_2U_2 x_3U_3
apprX=modeN_prod(modeN_prod(modeN_prod(G,U{1},1),U{2},2),U{3},3);

% 展示拟合结果
imshow(uint8(apprX(:,:,1)))
% 展示与原图像差值
figure()
imagesc(abs(X(:,:,1)-apprX(:,:,1)))

colormap(PYCM().inferno())% 使用python colormap渲染，可删去
colorbar
end