tu3.m

来自「该程序基于小波变换」· M 代码 · 共 30 行

%装载并显示原始图象
load flujet;
subplot(1,3,1);
image(X);
colormap(map);
title('原始图象');
axis square;
%=============================================
%加噪并显示图象
init=2055615866;
rand('seed',init);
XX=X+50*randn(size(X));
subplot(1,3,2);
image(XX);
colormap(map);
title('含噪图象');
axis square;
%===============================================
%对图象消噪
%用sym4小波函数对XX进行3层分解
[c,l]=wavedec2(XX,3,'sym4');
%实现低通滤波消噪
a2=wrcoef2('a',c,l,'sym4');
%==============================================
%显示消噪后的结果
subplot(1,3,3);
image(a2);
colormap(map);
title('消噪图象');
axis square;