newton_NEWTONODE_newton_DECRYPT_chaotic_Encrypt_源码

clc; close all; clear all; % 读取图像'lena.bmp' I = imread('lena.jpg'); figure; imshow(I); title('original image'); drawnow; figure; imhist(I) drawnow;% 产生实数混沌序列x,y,z n =size(I,1)*size(I,2)*size(I,3); level =0; s=10; r=3.55; b=0.25; x0=1.1840000001; y0=1.362700000; z0=1.2519000001; h=0.01; [x,y,z]=newton(n,level,s,r,b,x0,y0,z0,h); drawnow; % 对实数混沌序列进行预处理 X=10*x-round(10*x); Y=10*y-round(10*y); Z=10*z-round(10*z); drawnow; % % 产生索引序列lx,ly,lz [X1,lx]=sort(X); [Y1,ly]=sort(Y); [Z1,lz]=sort(Z); drawnow; % % % 对图像进行加密 A=I(:,:,1); B=I(:,:,2); C=I(:,:,3); A=repmat(A,3); B=repmat(B,3); C=repmat(C,3); A1=reshape(A,1,[]); B1=reshape(B,1,[]); C1=reshape(C,1,[]); A2(1:n)=A1(lx); B2(1:n)=B1(ly); C2(1:n)=C1(lz); A2=repmat(A2,1,3); B2=repmat(B2,1,3); C2=repmat(C2,1,3); % A3=reshape(A2,1,[]); % B3=reshape(B2,1,[]); % C3=reshape(C2,1,[]); A3=reshape(A2,size(A,1),size(A,2)); B3=reshape(B2,size(B,1),size(B,2)); C3=reshape(C2,size(C,1),size(C,2)); drawnow; % 显示加密后的图像 figure; A3=A3(1:128,1:128); B3=B3(1:128,1:128); C3=C3(1:128,1:128); I5(:,:,1)=A3; I5(:,:,2)=B3; I5(:,:,3)=C3; zz=randi([0,255],128,128); %zz=repmat(zz,1,3) zz=uint8(zz); zzz=bitxor(zz,I5); imhist(zzz) figure; imshow(zzz) title('encrypted image'); imwrite(zzz,'lennaencrypt.jpg'); drawnow; % 对图像进行解密 I11=imread('lennaencrypt.jpg'); I1=bitxor(zz,I11); A4=I1(:,:,1); B4=I1(:,:,2); C4=I1(:,:,3); A4=repmat(A4,3); B4=repmat(B4,3); C4=repmat(C4,3); A5=reshape(A4,1,[]); B5=reshape(B4,1,[]); C5=reshape(C4,1,[]); % A5=repmat(A5,1,3); % B5=repmat(B5,1,3); % C5=repmat(C5,1,3); llx=double(lx);lly=double(ly);llz=double(lz); % A6=A6(llx); % A6=A6(llx); % A6=A6(llx); A6(1:n)=A5(lx); B6(1:n)=B5(ly); C6(1:n)=C5(lz); % A6=double(A6); % B6=double(B6); % C6=double(C6); % A6=A6(llx); % B6=B6(lly); % C6=C6(llz); A6=repmat(A6,1,3); B6=repmat(B6,1,3); C6=repmat(C6,1,3); A7=reshape(A6,size(A4,1),size(A4,2)); B7=reshape(B6,size(B4,1),size(B4,2)); C7=reshape(C6,size(C4,1),size(C4,2)); A7=A7(1:128,1:128); B7=B7(1:128,1:128); C7=C7(1:128,1:128); I1(:,:,1)=A7; I1(:,:,2)=B7; I1(:,:,3)=C7; drawnow; % 显示解密后的图像 figure; imshow(I1); title('decrypted image'); imwrite(I1,'lennadecrypt.png'); drawnow; figure; imhist(I1); % 错误密码解密后的图像 x1=1.1840000001 y1=1.3627000001 z1=1.2519000001 I2=imread('lennaencrypt.jpg'); n =size(I2,1)*size(I2,2)*size(I2,3); level =0; s=10; r=3.55; b=0.25; x1=1; y1=1; z1=0; h=0.01 [x2,y2,z2]=newton(n,level,s,r,b,x1,y1,z1,h); %[x2,y2,z2]=Arneodo(n,level,s,r,b,x1,y1,z1,h); X2=10*x2-round(10*x2); Y2=10*y2-round(10*y2); Z2=10*z2-round(10*z2); [X3,lx1]=sort(X2); [Y3,ly1]=sort(Y2); [Z3,lz1]=sort(Z2); A8=I2(:,:,1); B8=I2(:,:,2); C8=I2(:,:,3); A8=repmat(A8,1,3); B8=repmat(B8,1,3); C8=repmat(C8,1,3); A9=reshape(A8,1,[]); B9=reshape(B8,1,[]); C9=reshape(C8,1,[]); A10(lx1)=A9(1:n); B10(ly1)=B9(1:n); C10(lz1)=C9(1:n); A11=reshape(A10,size(A8,1),size(A8,2)); B11=reshape(B10,size(B8,1),size(B8,2)); C11=reshape(C10,size(C8,1),size(C8,2)); A11=A11(1:128,1:128); B11=B11(1:128,1:128); C11=C11(1:128,1:128); I2(:,:,1)=A11; I2(:,:,2)=B11; I2(:,:,3)=C11; figure; imshow(I2); title('error decrypted image'); imwrite(I2,'errordecryptlenna.jpg'); drawnow; figure; imhist(I2) drawnow;