%該函數針對其中的近似系數和高頻系數進行加密，以達到加密的效果（采用函數調用的形式，即構建一個hundungen函數）
clear?all;
t0?=?clock;%測試程序運行時間
im=imread（‘e:\lena.jpg‘）;
im1=rgb2gray（im）;%圖像灰度化
im1=medfilt2（im1[3?3]）;%圖像平滑處理
%subplot（331）;
%figure（1）;
%imshow（im）;%X為原始圖像
%title（‘原始圖像‘）;
figure（2）;
imshow（im1）;
title（‘灰度化處理‘）;
im1=double（im1）;
[ca1ch1cv1cd1]=dwt2（im1‘bior3.7‘）;%小波變換，獲取圖像的低頻高頻系數
figure（3）;
subplot（231）;
imshow（ca1[]）;
title（‘圖像近似‘）;
subplot（232）;
imshow（ch1）;
title（‘低頻水平分量‘）;
subplot（233）;
imshow（cv1）;
title（‘低頻垂直分量‘）;
subplot（234）;
imshow（cd1）;
title（‘高頻分量‘）;
%以下為混沌加密算法
%以下為混沌加密算法
%以下為混沌加密算法
%以下為混沌加密算法
%以下為混沌加密算法
%以下為混沌加密算法
[MN]=size（ca1）;
e=hundungen（MN0.1）;
tt=0.1;
fca1=mod（tt*ca1+（1-tt）*e256）;
subplot（235）;
imshow（fca1[]）;
title（‘加密‘）;
im2=idwt2（ca1ch1cv1cd1‘bior3.7‘）;
figure（4）;
imshow（uint8（im2）[]）;
title（‘灰度圖像小波重構‘）;
im3=idwt2（fca1ch1cv1cd1‘bior3.7‘）;
figure（5）;
imshow（uint8（im3）[]）;
title（‘加密圖像小波重構‘）;
%以下為混沌解密算法
%以下為混沌解密算法
%以下為混沌解密算法
%以下為混沌解密算法
%以下為混沌解密算法
%以下為混沌解密算法
e=hundungen（MN0.1）;
[fca1ch1cv1cd1]=dwt2（im3‘bior3.7‘）;
fca2=（fca1-（1-tt）*e）/tt;
im4=idwt2（fca2ch1cv1cd1‘bior3.7‘）;
figure（6）;
imshow（uint8（im4）[]）;
title（‘解密圖像小波重構‘）;
%置亂后圖像的均值
figure（7）;
subplot（221）
imhist（uint8（im1））;
title（‘初始圖像的直方圖‘）;
subplot（222）
imhist（uint8（fca1））;
title（‘ca1系數加密之后的直方圖‘）;
subplot（223）
imhist（uint8（im3））;
title（‘加密之后的直方圖‘）;
subplot（224）
imhist（uint8（im4））;
title（‘解密之后的直方圖‘）;
ssy=sum（sum（im3））;
uy=ssy/（M*N）;%置亂后圖像的均值
vy=sum（sum（（im3-uy）^2））;
ssx=sum（sum（im1））;
ux=ssx/（M*N）;%原圖像的均值
vx=sum（sum（（im1-ux）^2））;
Variancey=vy/uy;%置亂后圖像的方差
Variancex=vx/ux;%原圖像的方差
DDD=Variancey/Variancex;%置亂度
etime（clockt0）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2008-08-07?16:22??tuxiangjiami8105\
?????文件????????2024??2008-04-24?12:25??tuxiangjiami8105\diaoyong.m
?????文件?????????686??2008-04-27?20:38??tuxiangjiami8105\hundungen.m
?????文件????????1223??2007-11-29?11:22??tuxiangjiami8105\說明.txt