`
- 浏览:
1816172 次
- 性别:
- 来自:
北京
-
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
分享到:
Global site tag (gtag.js) - Google Analytics
相关推荐
按照 认知无线电系统频谱感知技术研究 硕士论文中的描述, 自己做的论文第三章图的matlab仿真程序
无线传感器网络所有定位算法matlab仿真代码 很全的无线传感器网络定位算法仿真代码 并且包含各个算法的论文
我发表小论文用的程序代码,永磁直线电机的速度位置跟踪,用的无模型自适应控制
matlab图像处理边缘算子程序仿真加PPT加论文
数控云台proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资料 数控云台proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资料 数控云台proteus仿真+程序...
电子秤proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资料 电子秤proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资料 电子秤proteus仿真+程序资料.zip...
量程自动切换数字电压表proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资料 量程自动切换数字电压表proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB...
基于STM32流水灯八个灯仿真代码以及相关报告,能够仿真成功。还包括相关软件下载,和学习资料,包括仿真文件,hen文件以及源代码可以自我修改方便和学习
单片机C语言程序设计实训100例——基于8051+Proteus仿真.zip电子设计项目毕业设计及产品设计资料论文,源代码及PC单片机C语言程序设计实训100例——基于8051+Proteus仿真.zip电子设计项目毕业设计及产品设计资料论文,...
电压检测系统(含VB上位机)proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资料 电压检测系统(含VB上位机)proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及...
一种智能频率计的设计与制作(AVR)proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文,源代码及PCB资一种智能频率计的设计与制作(AVR)proteus仿真+程序资料.zip电子设计项目毕业设计及产品设计资料论文...
在解决基本算法问题的基础上进一步优化程序,编写出具备泛用性质而且友好的图形用户界面程序。具有实际意义,可为管理者进行决策提供辅助。 创新使用MATLAB数学工具代码编写和GUI设计来解决二次指派问题。能够借助...
通信论文仿真代码CL-Paper code:his 是与以下论文相关的代码包: Sanjeev Sharma、Anubha Gupta 和 Vimal Bhatia,“使用信号簇稀疏性的 IR-UWB 通信中的脉冲噪声缓解”,IEEE 通信快报,第一卷。 22,没有。 3, pp....
基于《基于Matlab的曲柄摇杆机构运动学仿真》郭飞,韩冰这篇论文完成的simulink程序。新增根据四个摇杆距离可以得知初试速度。
本代码实现了基于机器学习方法(SVM和NB)的MIMO天线选择仿真,针对Transmit Antenna Selection in MIMO Wiretap Channels: A Machine Learning Approach论文进行了复现,绘制出信道容量与信噪比SNR关系图。
北京航空航天大学大学出本系统辨识与自适应控制一书的MATLAB源程序
论文的matlab仿真程序代码
2 程序及仿真图形的建立 2 2.1 ASK调制解调 2 2.1.1 ASK程序 2 2.1.2 ASK图形 2 2.2 PSK调制解调 4 2.2.1 PSK程序 4 2.2.2 PSK图形 5 2.3 FSK调制解调 7 2.3.1 FSK程序 7 2.3.2 FSK图形 7 2.4 误码率分析 9 2.4.1 ...
论文的matlab仿真程序代码多天线用户的大规模 MIMO 这是一个与以下科学论文相关的代码包: Xueru Li、Emil Björnson、Shidong Zhou、Jing Wang,“”,国际电信会议 (ICT) 论文集,希腊塞萨洛尼基,2016 年 5 月。 ...