program_11_13.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

% 载入原始图像
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');
% 第二种用法
a0bis = 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)))