function?ZZP_Arnold（Imagekey）
%輸入?yún)?shù)說(shuō)明：
%		Image:		????????待加密的圖像文件名（注意寫格式后綴），只能為二維
%		??key:	?????????	圖像需要變換迭代的次數(shù)（密鑰）
%輸出參數(shù)說(shuō)明：
%????????	M:				置亂后圖像數(shù)據(jù)矩陣
%?????zzp.bmp：?????????????輸出由M對(duì)應(yīng)的圖像文件
%中間變量說(shuō)明
%???????????N:??????????????圖像的階數(shù)
%???????????T:??????????????圖像的周期??
%*********************圖像數(shù)值矩陣Arnold轉(zhuǎn)換函數(shù)****************************
%輸入格式控制
if?nargin<2????????????????????????????%如果輸入?yún)?shù)小于2
????disp（‘請(qǐng)按程序的輸入?yún)?shù)格式輸入?yún)?shù)！！！‘）;
????return;
end
%----------------------------Arnold?變換-----------------------------------
message=double（imread（Image））;??????????%讀入文件
Size_Q?=?size（message）;?????????????????%計(jì)算原始圖像的大小
Mm=size（message1）;?????????????????????%計(jì)算原始圖像矩陣的行數(shù)???????????
Nm=size（message2）;?????????????????????%計(jì)算原始圖像矩陣的列數(shù)???
%如果不是二維或三維數(shù)組，則不處理，返回
if?（length（Size_Q）?==?2）?
???if??Mm~=Nm??
??????disp（‘不是方陣，不能Arnold轉(zhuǎn)換‘）;
??????return
???end
else
???disp（‘不是二維數(shù)組，不進(jìn)行Arnold變換‘）;
???return?
end
tempImg=message;???????????????????????%圖像矩陣賦給tempImg
for?n=1:key????????????????????????????%置亂次數(shù)
??for?u=1:Mm
????for?v=1:Nm
??????temp=tempImg（uv）;
??????ax=mod（（u-1）+（v-1）Mm）+1;????????%新像素行位置
??????ay=mod（（u-1）+2*（v-1）Nm）+1;??????%新像素列位置
??????outImg（axay）=temp;
????end
??end
tempImg=outImg;
end
%?顯示原始圖像
figure（1）
subplot（121）;
imshow（message[]）;
title（‘原始圖像‘）;
%顯示置亂后的圖像
subpl