%------------------------------------------------------
%??XOR+Arnold+XOR+G_Lorenz
%------------------------------------------------------
clear?all;
clc;
%-------------------------------------------------------------
%<---DOWN
%------------------------------------------------------
%利用混沌系統(tǒng)生成密鑰流首先生成混沌吸引子
%參數聲明
a1=-2;a2=6.7;a13=0;a23=-1;
b1=1;b2=-0;b13=0;b23=0;d2=0;
c3=-1;c12=0;c11=0;c22=1;c33=0;d3=0;
detat=0.001;
%--------------------------------------------------------
%初始值
x（1）?=?0.1;
y（1）?=?1.5;
z（1）?=?3.1;
%--------------------------------------------------------
%離散化形式的廣義Lorenz系統(tǒng)
for?i?=?1:100000
????x（i+1）?=?（1+detat*a1）*x（i）...
????????+detat*a2*y（i）...
????????+detat*a13*x（i）*z（i）...
????????+detat*a23*y（i）*z（i）;
????y（i+1）?=?detat*b1*x（i）...
????????+（1+detat*b2）*y（i）...
????????+detat*b13*x（i）*z（i）...
????????+detat*b23*y（i）*z（i）...
????????+d2;
????z（i+1）?=?（1+detat*c3）*z（i）...
????????+detat*c12*x（i）*y（i）...
????????+detat*c11*（x（i）^2）...
????????+detat*c22*（y（i）^2）...
????????+detat*c33*（z（i）^2）...
????????+d3;
end
%------------------------------------------------------
%讀入原圖
Image=double（imread（‘cameraman3.tif‘））;
[mn]=size（Image）;
N=m*n;
e=mod（N2）;%判斷奇偶數，E為1則奇數，反之則偶數
%------------------------------------------------------
disp（‘加密中，請稍候.........‘）
tic
%------------------------------------------------------
%XOR1，從上到下從左到右依次進行
S1=119;
Image1=reshape（Image1N）;
for?i=1:N
????if?i==1
????????Image1（i）=bitxor（Image1（i）S1）;
????%elseif?i==n+2
????????%Image1（i）=bitxor（Image1（i）Image1（N））;
????else
????????Image1（i）=bitxor（Image1（i）Image（i-1））;
????end
end
ImageX1=reshape（Image1mn）;
figure（1）
imshow（uint8（ImageX1））;
%for?k=1:m
????%for?l=1:n
????????%Image（11）=bitxor（Image（11）S1）;%第一個元素與固定密鑰XOR
????????%if?k==m?&?l==n
????????????%Image（22）=bitxor（Image（22）Image（mn））;
????????%else
????????????%Image（kl）=bitxor（Image（kl）Image（kl-1））;
????????%end
????%end
%end
%------------------------------------------------------
%Arnold
for?k=1:m
????for?l=1:n
????????k1=k+l;
????????l1=k+2*l;%Arnold
????????if?k1>m
????????????k1=mod（k1m）;
????????end
????????if?l1>n
????????????l1=mod（l1m）;
????????end
????????if?k1==0
????????????k1=m;
????????end
????????if?l1==0
????????????l1=m;
????????end
????????ImageA（k1l1）=ImageX1（kl）;
????????%Image1作為臨時寄存變量必須在位置變換的循環(huán)之內
????end
end
Image=ImageA;%用于循環(huán)多輪加密
%Image2=Image1;
figure（2）
imshow（uint8（ImageA））;
%------------------------------------------------------
%XOR2，從左到右從上到下依次進行先行后列
S2=207;
ImageA=ImageA‘;
Image2=reshape（ImageA1N）;
for?i=1:N
????if?i==1
????????Image2（i）=bitxor（Image2（i）S2）;
????%elseif?i==n+2
????????%Image2（i）=bitxor（Image2（i）Image2（N））;
????else
????????Image2（i）=bitxor（Image2（i）Image2（i-1））;
????end
end
ImageX2=reshape（Image2mn）;
figure（3）
imshow（uint8（ImageX2））;
%Image2=reshpmn21N（Image）;%用自定義的函數
%for?i=1:N
????%if?i==1
????????%Image2（i）=b