例程11-2.m

来自「小波分析Matlab实现」· M 代码 · 共 21 行

% 当前延拓模式是补零
% 装载原始图像
load sinsin;
% X 包含装载的图像
% 绘制原始图像
subplot(1,2,1);
image(X);
colormap(map);
title('原始图像');
sX = size(X);
% 使用db4对X进行单尺度分解
[cA1,cH1,cV1,cD1] = dwt2(X,'db4');
% 利用尺度为1时的系数直接重构
A0 = idwt2(cA1,cH1,cV1,cD1,'db4',sX);
% 绘制重构图像
subplot(1,2,2);
image(A0);
title('重构图像');
% 检查重构效果
max(max(abs(X-A0)))