基于加取模和循环左移运算的扩散算法matlab

• 基于加取模和循环左移运算的扩散算法
  实际上是对前面两种扩散算法的显著改良。
  
  LSB3表示取数据的最低三位，对于8b的灰度图像，每个像素都是8b的，所以LSB3得到的数据范围始终在0~7，是一个像素点数据循环的有效范围。
  这里注意，除了循环左移，还可以用循环右移，如果加密时用了循环左移，那么解密时应该循环右移。
  实现：

clc;clear;
P=imread('lenna.bmp');P=rgb2gray(P);
iptsetpref('imshowborder','tight');
figure(1);subplot(1,3,1);imshow(P);title('明文图像');
[M,N]=size(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);
for i=1:n+t
    K11=a*(y0-x0)+w0;K12=a*(y0-(x0+K11*h/2))+w0;
    K13=a*(y0-(x0+K12*h/2))+w0;K14=a*(y0-(x0+h*K13))+w0;
    x1=x0+(K11+K12+K13+K14)*h/6;
    
    K21=c*x1-y0-x1*z0;K22=c*x1-(y0+K21*h/2)-x1*z0<