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（