图像处理之像素级环形扩散

刚入门学习matlab，自己根据《位级同步置乱扩散和像素级环形扩散图像加密算法》这个文章实现了一下像素级环形扩散这一小节的代码，和大家分享一下，能力有限，有不足之处还望见谅。代码如下：

clear all
clc
X=imread('lena.jpg');
X=double(X);
figure()
subplot(221)
imshow(X,[]);
title("原图");
S=size(X,2);

%密钥生成
% % % % % % % % 
x2=0.4;
L2=0.3;
x2 = mod(x2+sum(sum(X))/256,1);

for i=1:300
    x2=F(x2,L2);
end

y3=zeros(1,S*S);
y4=zeros(1,S*S);
y3(1)=x2;
for i=1:S*S-1
  y3(i+1)=F(y3(i),L2);
end
y4(1)=y3(end);
for i=1:S*S-1
    y4(i+1)=F(y4(i),L2);
end
%混沌序列标准化
z3=floor(mod(y3.*10^14,256));
z4=floor(mod(y4.*10^14,256));
% % % % % % % % 

H=reshape(flipud(rot90(X)),1,S^2);
H1=[];
H1(1)=bitxor(bitxor(sum(H),H(1)),z3(1));
for i=1:S^2-1
    H1(i+1)=bitxor(bitxor(H(i+1),H1(i)),z3(i+1));
end

H1=reshape(H1,S,S);

H1=flipud(fliplr(H1));
V=reshape(H1,1,S^2);
E=[];
E(1)=bitxor(bitxor(V(1),sum(V