运用Arnold与Logistic的简单图像加密程序-AICC.m

%------------------------------------------------------
%  XOR+Arnold+XOR+G_Lorenz
%------------------------------------------------------
clear all;
clc;
%-------------------------------------------------------------
%<---DOWN
%------------------------------------------------------
%利用混沌系统生成密钥流首先生成混沌吸引子
%参数声明
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系统
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作为临时寄存变量必须在位置变换的循环之内
    end
end
Image=ImageA;%用于循环多轮加密
%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