meyer小波的分解与重构

i1=imread（‘G:\tupian\hongwai.jpg‘）;
i=rgb2gray（i1）;
figure（1）;
imshow（i）;
[cs]=wavedec2（i3‘meyer‘）;
sizei=size（i）;
sizec=size（c）;
val_s=s;
%提取1，2，3层的低频系数
%提取尺度为3时的低频系数
ca3=appcoef2（cs‘meyer‘3）;
sizeca3=size（ca3）;
%绘制尺度为3时的低频系数
figure（2）;
subplot（221）;
image（ca3）;
colormap（gray（256））;
%imshow（ca3）;
title（‘尺度3时低频‘）;
%提取尺度为2时的低频系数
ca2=appcoef2（cs‘meyer‘2）;
sizeca2=size（ca2）;
%绘制尺度为2时的低频图像
subplot（222）;
image（ca2）;
colormap（gray（256））;
%imshow（ca2）;
title（‘尺度2时低频‘）;
%提取尺度为1时的低频图像
ca1=appcoef2（cs‘meyer‘1）;
sizeca1=size（ca1）;
%绘制尺度为1时的低频图像
subplot（223）;
image（ca1）;
colormap（gray（256））;
%imshow（ca1）;
title（‘尺度为1时低频‘）;
%提取1，2，3层的各个方向的高频系数
%从小波分解结构[cs]中提取尺度为3时的各个方向的高频系数
[chd3cvd3cdd3]=detcoef2（‘all‘cs3）;
sizechd3=size（chd3）;
%绘制尺度为3时的水平方向高频图像
figure（3）;
subplot（331）;
imshow（chd3）;
title（