量子图像加密仿真实现

量子图像加密算法仿真实现，包含经典图像仿真与量子图像仿真实现代码

function N=Arnold（ImageabFrequencycrypt）
%图像数值矩阵 Arnold转换函数
%输入参数
%Image：待加密（待解密）图像（注意书写格式，后缀），只能为二维
%Frequency：图像需要变换迭的次数
%crypt    0～加密；1～解密
%输出参数
%M：转换后图像数据矩阵
%   输出由M对应的图像矩阵
if nargin<3
    disp（‘请按程序的输入参数格式输入参数！！！‘）;
    return;
end
if crypt~=0&crypt~=1
    disp（‘encrypt必须为0或1！‘）;
end
%将Q赋值给M，计算Q的大小
Q=Image;
M=Q;
Size_Q=size（Q）;
%如果不是二维或三维数组，则不处理，返回
if（length（Size_Q）<=2）
    if Size_Q（1）~=Size_Q（2）
        disp（‘不是方阵，不能Arnold变换‘）;
        return;
    end
else
    disp（‘不是二维数组，不能Arnold变换‘）;
    return;
end
%...................
%Arnold转换
n=0;
k=Size_Q（1）;
M1_t=Q;
M2_t=Q;
if crypt==1                   %解密
    Frequency=ArnoldPeriod（Size_Q（1））-Frequency;
end
for s=1:Frequency
    n=n+1;
    if mod（n2）==0;
        for i=1:k;
            for j=1:k;
                c=M2_t（ij）;
                M1_t（mod（i+a*j-1-ak）+1mod（b*i+（a*b+1）*j-b-（a*b+1）k）+1）=c;
            end
        end
    else
        for i=1:k;
            for j=1:k;
                c=M1_t（ij）;
                M2_t（mod（i+a*j-1-ak）+1mod（b*i+（a*b+1）*j-b-（a*b+1）k）+1）=c;
            end
        end
    end
end
if mod（Frequency2）==0;
    M=M1_t;
else
    M=M2_t;
end
%...................
%imwrite（double（M）/255strcat（‘Arnold_‘num2str（Frequency）‘_‘Image）‘bmp‘）;
N=M;
function Period=ArnoldPeriod（N）;
%求周期
if （N<2）
    Period=0;
    return;
end
n=1;
x=1;
y=1;
while n~=0
    xn=x+y;
    yn=x+2*y;
    if（mod（xnN）==1&&mod（ynN）==1）
        Period=n;
        return;
    end
    x=mod（xnN）;
    y=mod（ynN）;
    n=n+1;
end

    

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1764  2014-05-22 10:30  方案二\Arnold.m

     文件     292638  2014-03-31 15:20  方案二\baboon.tiff

     文件     262278  2014-04-06 14:46  方案二\boat.tiff

     文件     176914  2014-03-31 15:20  方案二\bridge.tiff

     文件        918  2014-05-26 16:25  方案二\cuojiemi.m

     文件       1685  2014-05-27 21:41  方案二\cuowujiemi.m

     文件       1260  2014-05-27 17:04  方案二\cuowujiemi1.m

     文件        677  2014-05-08 11:13  方案二\imgentropy.m

     文件        591  2014-05-08 15:22  方案二\imgentropy1.m

     文件       1303  2014-07-03 10:14  方案二\jiemi.m

     文件        656  2014-05-26 16:15  方案二\jiemi1.m

     文件     249438  2013-10-08 10:04  方案二\lena.tiff

     文件       4166  2015-03-31 16:55  方案二\maohesui.m

     文件     269806  2014-03-31 15:20  方案二\peppers.tiff

     文件     214598  2014-03-31 15:20  方案二\plane.tiff

     文件       1124  2013-10-08 10:04  方案二\relativity1.m

     文件       1124  2013-11-20 15:25  方案二\relativity2.m

     目录          0  2015-03-09 16:08  方案二

----------- ---------  ---------- -----  ----

              1480940                    18