cwq3.m

来自「matlab程序源码,这个是我和同学一起做的,大家可以参考看看怎么样的,谢谢大家」· M 代码 · 共 31 行

%基于小波的图象消噪
load wifs;
subplot(221);
image(X);
colormap(map);
title('原图 ');
axis square; %画出原图象
%生成含噪图象并显示
init=2055615866;
randn('seed',init)
XX=X+randn(size(X));
subplot(222);
image(XX);
colormap(map);
title('含噪声图象 ');
axis square; %画出含噪声图象
%下面对图象消噪
%首先用sym4小波函数对x进行两层分解
[c,s]=wavedec2(XX,2,'sym4');
a1=wrcoef2('a',c,s,'sym4',1);%第一次低通滤波消噪
subplot(223);
image(a1);
title('第一次消噪后图象 ');
axis square; %画出第一次低通滤波消噪后图象
%再次实现低通滤波消噪
a2=wrcoef2('a',c,s,'sym4',2); 
subplot(224);
image(a2);
title('第二次消噪后图象 ');
axis square; %画出第二次低通滤波消噪后图象