混沌保密通信传输图像数据.txt

通过细胞神经网络混沌系统对图像数据进行加密和解密

a=imread('in.jpg');
aa=double(a);
k=0;
for i=1:99
    for j=1:99
        k=k+1;
        in(k)=aa(i,j)/255;
    end
end
dt=0.01;
x1(1)=0;x2(1)=0;x3(1)=0;x4(1)=0;
for i=1:9801
    y4(i)=(abs(x4(i)+1)-abs(x4(i)-1))/2;
    s(i)=80.1900*x1(i)+20.1546*x2(i)+11.9363*x3(i)-89.8*x4(i)+200*y4(i)+in(i);
    x1(i+1)=x1(i)+(-x3(i)-x4(i))*dt;
    x2(i+1)=x2(i)+(2*x2(i)+x3(i))*dt;
    x3(i+1)=x3(i)+(14*x1(i)-14*x2(i))*dt;
    x4(i+1)=x4(i)+(100*x1(i)-100*x4(i)+200*y4(i)+in(i))*dt;    
end

x11(1)=0;x12(1)=0;x13(1)=0;x14(1)=0;
for i=1:9801
    y14(i)=(abs(x14(i)+1)-abs(x14(i)-1))/2; 
    in1(i)=s(i)-(80.1900*x11(i)+20.1546*x12(i)+11.9363*x13(i)-89.8*x14(i)+200*y14(i));   
    x11(i+1)=x11(i)+(-x13(i)-x14(i))*dt;
    x12(i+1)=x12(i)+(2*x12(i)+x13(i))*dt;
    x13(i+1)=x13(i)+(14*x11(i)-14*x12(i))*dt;
    x14(i+1)=x14(i)+(100*x11(i)-100*x14(i)+200*y14(i)+in1(i))*dt;    
end
e=in-in1;
k=0;
for i=1:99
    for j=1:99
        k=k+1;
        bb(i,j)=in1(k)*255;
    end
end
b=uint8(bb);

subplot(3,5,1)
plot(in);xlabel('n');ylabel('in(n)');title('输入');
subplot(3,5,6)
plot(s);xlabel('n');ylabel('s(n)');title('混沌');
subplot(3,5,11)
plot(in1);xlabel('n');ylabel('in1(t)');title('输出');
subplot(3,5,2)
plot(x1,x2);xlabel('x1');ylabel('x2');title('x1-x2相图');
subplot(3,5,7)
plot(x1,x3);xlabel('x1');ylabel('x3');title('x1-x3相图');
subplot(3,5,12)
plot(x1,x4);xlabel('x1');ylabel('x4');title('x1-x4相图');
subplot(3,5,3)
plot(x11,x12);xlabel('x11');ylabel('x12');title('x11-x12相图');
subplot(3,5,8)
plot(x11,x13);xlabel('x11');ylabel('x13');title('x11-x13相图');
subplot(3,5,13)
plot(x11,x14);xlabel('x11');ylabel('x14');title('x11-x14相图');
subplot(3,5,4)
imshow(a);xlabel('x');ylabel('y');title('原灰度图像');
subplot(3,5,9)
plot(e);xlabel('n');ylabel('e(n)');title('误差');
subplot(3,5,14)
imshow(b);xlabel('x');ylabel('y');title('还原灰度图像');

for i=1:9801
    s(i)=1*s(i)+0.01*2*(rand(1)-0.5);
end

x11(1)=0;x12(1)=0;x13(1)=0;x14(1)=0;
for i=1:9801
    y14(i)=(abs(x14(i)+1)-abs(x14(i)-1))/2; 
    in1(i)=s(i)-(80.1900*x11(i)+20.1546*x12(i)+11.9363*x13(i)-89.8*x14(i)+200*y14(i));   
    x11(i+1)=x11(i)+(-x13(i)-x14(i))*dt;
    x12(i+1)=x12(i)+(2*x12(i)+x13(i))*dt;
    x13(i+1)=x13(i)+(14*x11(i)-14*x12(i))*dt;
    x14(i+1)=x14(i)+(100*x11(i)-100*x14(i)+200*y14(i)+in1(i))*dt;    
end
e=in-in1;
k=0;
for i=1:99
    for j=1:99
        k=k+1;
        bb(i,j)=in1(k)*255;
    end
end
b=uint8(bb);

subplot(3,5,5)
imshow(a);xlabel('x');ylabel('y');title('原灰度图像');
subplot(3,5,10)
plot(e);xlabel('n');ylabel('e(n)');title('误差');
subplot(3,5,15)
imshow(b);xlabel('x');ylabel('y');title('还原灰度图像（0.01噪声）');