基于混沌的一般彩色图像加密算法的Matlab实现

基于混沌的一般彩色图像加密算法的Matlab实现

说明：

1. 混沌系统使用Loren连续混沌系统，使用四阶龙格库塔法将其离散化
2. 加密方法使用基于空域的像素置乱-扩散结构
3. 密文只能以无损压缩方式保存或分发，例如.bmp

Matlab代码（密文无损压缩图像保存）

clc;clear;close all;
P=imread('房屋.tiff'); %载入图片
figure;
subplot(1,3,1);
imshow(uint8(P));
title('房屋明文');

P_R=P(:,:,1);P_G=P(:,:,2);P_B=P(:,:,3);


[M,N]=size(P_R);  %计算图片P尺寸大小
P=double(P);    %将图像转化为双精度类型

%% 第一部分
%产生混沌随机序列
n=2*M*N;      %计算所需要的序列长度
h=0.002; t=800;    
a=10; b=8/3; c=28; r =-1;   
x0=1.1; y0=2.2; z0=3.3; w0=4.4;
s=zeros(1,n);

tic;
%利用四阶龙格库塔法将超Lorenz混沌系统离散化的计算过程
for i=1:n+t
    K11=a*(y0-x0);               %x0表示xi
    K12=a*(y0-(x0+K11*h/2))+w0;
    K13=a*(y0-(x0+K12*h/2))+w0;
    K14=a*(y0-(x0+K13*h))+w0;
    x1=x0+(K11+2*K12+2*K13+K14)*h/6;    %x1表示xi+1
    
    K21=c*x1-y0-x1*z0;               
    K22=c*x1-(y0+K21*h/2)-x1*z0;
    K23=c*x1-(y0+K22*h/2)-x1*z0;
    K24=c*x1-(y0+K23*h)-x1*z0;
    y1=y0+(K21+2*K22+2*K23+K24)*h/6; 
    
    K31=x1*y1-b*z0;               
    K32