program_11_06.m

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

% 当前延拓模式是补零
% 装载原始图像
load sinsin; 
% X 包含原始图像
% 绘制原始图像
subplot(3,3,1);
image(X);
colormap(map);
title('原始图像');

% 使用db1对X进行尺度为2的分解
[c,s] = wavedec2(X,2,'db1');
sizex = size(X)
sizec = size(c)
val_s = s 
% 对于各个方向，从小波分解结构[c,s]中提取尺度为2时的高频系数
[chd2,cvd2,cdd2] = detcoef2('all',c,s,2); 
sizecd2 = size(chd2)
% 绘制尺度为2时高频图像
subplot(3,3,4);
image(chd2);
title('尺度为2时高频图像的水平部分');
subplot(3,3,5);
image(cvd2);
title('尺度为2时高频图像的垂直部分');
subplot(3,3,6);
image(cdd2);
title('尺度为2时高频图像的对角部分');
% 对于各个方向，从小波分解结构[c,s]中提取尺度为1时的高频系数
[chd1,cvd1,cdd1] = detcoef2('all',c,s,1); 
sizecd1 = size(chd1)
% 绘制尺度为1时高频图像
subplot(3,3,7);
image(chd1);
title('尺度为1时高频图像的水平部分');
subplot(3,3,8);
image(cvd1);
title('尺度为1时高频图像的垂直部分');
subplot(3,3,9);
image(cdd1);
title('尺度为1时高频图像的对角部分');