数字水印matlab程序.docVIP

下载本文档

120
0
约2.08万字
约 21页
2018-12-31 发布于安徽
举报
版权申诉

数字水印matlab程序.doc

1、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

WORD格式可编辑专业知识整理分享 clear all; close all; clc; V=double(imread(D:\lena\lena.jpg)); imshow(mat2gray(V)); [i u]=size(V); %计算V的规格 r=100; %设置分解矩阵的秩 W=rand(i,r) %初始化WH，为非负数 H=rand(r,u) maviter=100; %最大迭代次数 for iter=1:maviter W=W.*((V./(W*H))*H); %注意这里的三个公式和文中的是对应的 W=W./(ones(i,1)*sum(W)); H=H.*(W*(V./(W*H))); end img_V=W*H; figure; imshow(mat2gray(img_V)); 首先读入原始图象并设置参数，然后嵌入水印信息，程序代码如下：clear% size=256;block=8;blockno=size/block;LENGTH=size*size/64; Alpha1=0.02;??Alpha2=0.1; T1=3; I=zeros(size,size); D=zeros(size,size); BW=zeros(size,size); block_dct1=zeros(block,block); ? %产生水印序列并对其排序 randn(seed,10);watermark1=randn(1,LENGTH); subplot(2,2,1);plot(watermark1);title(watermarc:Gaussian noise); subplot(2,2,3);title(edge of origine image)[Y0,I0]=sort(watermark1);%%读入原图象trueImage=imread(H:\Documents and Settings\sunhw\My Documents\My Pictures\biaozhun.bmp);alfa=.1;LENGTH=2500;subplot(2,2,2);imshow(trueImage);title(origine image:I);%%对原图象进行DCT变换dctF1=dct2(H:\Documents and Settings\sunhw\My Documents\My Pictures\biaozhun.bmp);[m,n]=size(dctF1);%%找出水印嵌入位置（幅值较大的n个频域成分）A=dctF1(:);[Y1,I1]=sort(A);x=m*n;k=LENGTH;M=zeros(x,1);%%修改幅值较大的n个频域成分的幅值，嵌入水印（因为两个问题不同，所以有两个注释符）for i=1:x???if k=1????M(x)=Y1(x)*(1+alfa*Y0(k));????k=k-1;???else????M(x)=Y1(x);???end???x=x-1;endN=zeros(x,1);x=m*n;for i=1:x????N(I1(i))=M(i);enda=1;for j=1:n??for i=1:m????dctF2(i,j)=N(a);????a=a+1;??endend%%DCT反变换，得到嵌入水印的图象idctF1=idct2(dctF2);subplot(2,2,4);imshow(log(abs(idctF1)),[ ]);title(embeded image:D);%end I=imread(D:\lena\1.jpg);disp(I);I=double(I)/255;?disp(I);I=ceil(I);%%%%%%%%%%显示水印图像%%%%%%%%%%%%%figure(1);subplot(2,3,1);imshow(I),title(水印图像)dimI=size(I);rm=dimI(1);cm=dimI(2);%%%%%%%%%%%%%%%5??以下生成水印信息 %%mark=I;alpha=0.05;V=imread(D:\lena\len