function demo1_1sl_HOSVD
% 读取数据
file=load('Goldhill.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(:,:,end)))
colorbar
end