randcs.m

matlab语言写的细胞自动机数字水印算法

%%%产生一个矩阵
function randmtx=randCS(row,col,seed)
A=4;
Xb=2.5;
randmtx(1,1)=seed;
for i=2:row*col
    randmtx(1,i)=A*sin(randmtx(1,i-1)-Xb)*sin(randmtx(1,i-1)-Xb);
    
end
a=randmtx;
a=reshape(a,row,col);
%%%转变为二值矩阵%%%%
for i=1:row*col
 if(randmtx(1,i)>=2*A/3)
     randmtx(1,i)=1;
 else
     randmtx(1,i)=0;
 end
end
randmtx=reshape(randmtx,row,col);
b=randmtx;

%%%subplot(133);imshow(k1);title('展示矩阵效果');
subplot(231);imshow(a);title('展示随机矩阵效果');
subplot(232);imshow(b);title('展示二值矩阵效果');