论文部分程序仿真代码

designed by lilizong
lilizong[at]gmail
实验仿真部分代码
row=1;
col=1;
%控制最终显示窗口中显示的图像行列数
myinforarray=linspace(0,0,5);
mypsnrarray=linspace(0,0,5);
w=1;
formyouter=1:1%此处设置迭代次数，多次迭代可以计算优化效果
A=imread('e:oat.bmp');
%subplot(row,col,[124578]),imshow(A),title('原始载体图像');
%subplot(row,col,1),imshow(A),title('原始载体图像');
[m,n]=size(A);
B=cell(8,8);%这里的k为你要分成大小为k*k的块
AA=A;
myinfor=0;
formyinter=1:1%此处设置插值次数

fori=1:512/8
forj=1:512/8
B{i,j}=A((i-1)*8+1:i*8,(j-1)*8+1:j*8);
end
end


%将B化为一维
c=cell(8,8);
k=1;
fori=1:512/8
forj=1:512/8
c{k}=B{i,j};
k=k+1;
end
end

%下面的代码用于生成混沌序列，生成4096个数字，然后排序
%得到一个位于4096之间的无重复值的混沌序列

x=linspace(0,0,4096);

x(1)=0.75+0.001+myinter/10+myouter/1000;
%0.73得值获得好结果,其他结果不佳
t=linspace(0,0,4096);
m=0;
n=0;
fori=1:4095
x(i+1)=1-2*x(i)*x(i);
k=k+1;
end

[pm,pn]=sort(x);
%pn为排序后结果，pm为原始数列

tod=cell(8,8);
%定义tod为需要处理的图像块
fori=1:4096
tod{i}=c{pn(i)};

test=tod;%添加一个测试块，看看问题出现在什么地方
end
%对图像进行拼接,显示一下置乱后的图像
toad=cell(8,8);
fori=1:4096
toad{i}=tod{i};
end
k=1;
toyuan=cell(8,8);
fori=1:64
forj=1:64
toyuan{i,j}=toad{k};
k=k+1;
end
end


toyuan=cell2mat(toyuan);

%subplot(row,col,2),imshow(toyuan),title('分块置乱后载体图像');

%temp=cell(8,8);
%得到像素值数组xiang
xiangsu=linspace(0,0,256);
%对取得的图像快进行分析，找到最大像素值
form=1:4096
temp=tod{i};
%temp=cell2mat(temp);
fori=1:8
forj=1:8
forn=1:256
iftemp(i,j)==n
xiangsu(n)=xiangsu(n)+1;
end
end
end
end
end
%求出最大像素的值
maxxiangsu=xiangsu(1);
fori=1:256
ifmaxxiangsu < xiangsu (i)
maxxiangsu =xiangsu(i);
maxvalueofxiangsu =i;
end
end

maxvalueofxiangsu;
%开始对需要嵌入信息的图像快进行移位处理
form =1:4096 %表示需要处理的块数
temp =tod{m};
fori =1:8
forj =1:8
iftemp(i,j)<maxvalueofxiangsu
%temp(i,j) =maxvalueofxiangsu-1;
temp(i,j) =temp(i,j)-1;
end
end

end
%tod{m} =mat2cell(temp); 更改4.28
tod{m} =temp;
end
%处理完毕，开始嵌入水印

%读入水印信息

b =imread('e:logo1.bmp');
b =im2bw(b);
[p,q] =size(b);
oo =p*q; %用来标识水印的总像素数

%watermark =reshape(b,1,oo);

watermark =zeros(1,oo);

n =1;
fori =1:p
forj =1:q
watermark(n) =b(i,j);
n =n+1;
end
end
%subplot(row,col,6),imshow(b),title('水印图像');
watermark;
waterlogistic =zeros(1,1600);
waterlogistic(1) =0.8;
fori =2:1600
waterlogistic(i) =1-2*waterlogistic(i-1)*waterlogistic(i-1);
end

fori =1:1600
ifwaterlogistic(i) > 0
waterlogistic(i)=1;
else
waterlogistic(i)=0;
end
end

waterlogistic;

watermark=xor(watermark,waterlogistic);
n=1;
fori=1:40
forj=1:40
towatermark(i,j)=watermark(n);
n=n+1;
end
end


%subplot(row,col,8),imshow(towatermark),title('加密后水印图像');



yuanwatermark=xor(watermark,waterlogistic);
n=1;
fori=1:40
forj=1:40
yuanwatermark1(i,j)=yuanwatermark(n);
n=n+1;
end
end

%subplot(row,col,10),imshow(yuanwatermark1),title('测试：恢复后水印图像');
%嵌入方式为，找到最高位，如果水印为0，不发生变化，如果水印为1，最高位减1
%令tt为嵌入水印循环量，如果t的值超过水印的信息量，则开始进行新的循环。
t=1;%用来标识嵌入位
infor=0;%用来标识总共可以嵌入的信息
form=1:4096
%temp=cell2mat(toa{m});
temp=tod{m};
fori=1:8
forj=1:8
iftemp(i,j)==maxvalueofxiangsu
ifwatermark(t)==1%水印的信息如果是1，则进行原信息的改变
temp(i,j)=temp(i,j)-1;%此处测试一下
%temp(i,j)=256;
end
t=t+1;
ift==1600
t=1;
end
infor=infor+1;
end
end
end
%tod{m}=mat2cell(temp);
tod{m}=temp;
end
myinfor=myinfor+infor;

%对图像进行拼接,显示一下添加水印以后的图像
toad=cell(8,8);
fori=1:4096
toad{i}=tod{i};
end
k=1;
toyuan=cell(8,8);
fori=1:64
forj=1:64
toyuan{i,j}=toad{k};
k=k+1;
end
end
toyuan=cell2mat(toyuan);
%subplot(row,col,1),imshow(toyuan),title('添加水印后的置乱图像');
%对图像进行拼接，显示添加水印后原始图像效果
toad=cell(8,8);
fori=1:4096
toad{pn(i)}=tod{i};
end
k=1;
toyuan=cell(8,8);
fori=1:64
forj=1:64
toyuan{i,j}=toad{k};
k=k+1;
end
end


toyuan=cell2mat(toyuan);


%subplot(row,col,2),imshow(toyuan),title('测试：加水印后原始图像');
%figure,imshow(toyuan);
%此处请屏蔽外层循环后开启！！！！！
A=toyuan;
end
%%%%%以下为提取水印过程
%开始读取水印
%1首先将嵌入水印后图像按照加密过程进行置乱，并从中提取嵌入到其中的加密后的水印信息
%2对提取出来的加密后的水印与混沌徐序列进行相与操作，得到原始水印图像
B=cell(8,8);%这里的k为你要分成大小为k*k的块
AA=A;
myinfor=0;
formyinter=1:1%此处设置插值次数

fori=1:512/8
forj=1:512/8
B{i,j}=A((i-1)*8+1:i*8,(j-1)*8+1:j*8);
end
end
%将B化为一维
c=cell(8,8);
k=1;
fori=1:512/8
forj=1:512/8
c{k}=B{i,j};
k=k+1;
end
end
x=linspace(0,0,4096);

x(1)=0.75+0.001+myinter/10+myouter/1000;
%0.73得值获得好结果,其他结果不佳
t=linspace(0,0,4096);
m=0;
n=0;
fori=1:4095
x(i+1)=1-2*x(i)*x(i);
k=k+1;
end

[pm,pn]=sort(x);
%pn为排序后结果，pm为原始数列

tod=cell(8,8);
%定义tod为需要处理的图像块
fori=1:4096
tod{i}=c{pn(i)};

test=tod;%添加一个测试块，看看问题出现在什么地方
end
watermark=zeros(1,1600);
n=1;
form=1:4096
%temp=cell2mat(toa{m});
temp=tod{m};
fori=1:8
forj=1:8
iftemp(i,j)==maxvalueofxiangsu
watermark(n)=0;
n=n+1;
end
iftemp(i,j)==maxvalueofxiangsu-1
watermark(n)=1;
n=n+1;
end

infor=infor+1;
end
end
end
%tod{m}=mat2cell(temp);
tod{m}=temp;
end
towatermark=zeros(40,40);
n=1;
fori=1:40
forj=1:40
towatermark(i,j)=watermark(n);
n=n+1;
end
end

%subplot(row,col,2),imshow(towatermark),title('测试：加水印后原始图像');

watermark2=zeros(1,1600);
fori=1:1600
watermark2(i)=watermark(i);
end


towatermark=xor(watermark2,waterlogistic);

towatermark2=zeros(40,40);

n=1;
fori=1:40
forj=1:40
towatermark2(i,j)=towatermark(n);
n=n+1;
end
end
subplot(row,col,1),imshow(towatermark2),title('提取水印图像');

%%%%%以上为提取水印过程
%figure,imshow(toyuan),title('加载了隐藏信息的图片');
%toyuan
%上述语句用于显示每次加密后的载体图片
x=AA;
x1=toyuan;
x=double(x);
x1=double(x1);
fori=1:256;
forj=1:256;
MYa(i,j)=x(i,j)^2;%afterfiltering'ssingle
MYb(i,j)=x(i,j)-x1(i,j);
end
end
d=0;
e=0;
fori=1:256;
forj=1:256;
d=d+MYa(i,j);
e=e+MYb(i,j)*MYb(i,j);
end
end
f=log10(d/e);
g=10*f;
%下面开始同时他的信息值和psnr值
myinforarray(w)=myinfor;
mypsnrarray(w)=uint8(g);
w=w+1;


end
%myinter
%fori=1:myouter


%myinforarray
%mypsnrarray
%figure,plot(myinforarray);
%figure,plot(mypsnrarray);


%myinforarray
%mypsnrarray