例程11-13.m

来自「matlab7小波分析书中的例子」· M 代码 · 共 30 行

% 载入原始图像
load sinsin;
% X 包含装载的图像
% 绘制原始图像
subplot(2,2,1);
image(X);
colormap(map);
title('原始图像');
% 使用sym4对X在第3层进行SWT分解
swc = swt2(X,3,'sym4');
% 第二种用法
[ca,chd,cvd,cdd] = swt2(X,3,'sym4');
% 从小波分解结构swc，重构s
a0 = iswt2(swc,'sym4');
% 绘制重构图像1
subplot(2,2,3);
image(a0);
colormap(map);
title('重构图像1');
% 第二种用法
a0 = iswt2(ca,chd,cvd,cdd,'sym4');
% 绘制重构图像2
subplot(2,2,4);
image(a0);
colormap(map);
title('重构图像2');
% 检查重构的效果
err = max(max(abs(X-a0)))
errbis = max(max(abs(X-a0bis)))