遗传算法matlab源代码(遗传算法的matlab代码)

关于遗传算法matlab源代码，遗传算法的matlab代码这个很多人还不知道，今天小爱来为大家解答以上的问题，现在让我们一起来看看吧！

1、一、嵌入水印信息的MATLAB程序首先读入原始图象并设置参数，然后嵌入水印信息，程序代码如下：clear%%读入原图象trueImage=imread('C:Documents and Settingsks001My DocumentsMy Pictureslean.tif');alfa=.1;LENGTH=0;subplot(2,2,1);imshow(trueImage);title('原始图象');%%对原图象进行DCT变换dctF1=dct2('C:Documents and Settingsks001My DocumentsMy Pictureslean.tif');subplot(2,2,2);imshow(log(abs(dctF1)),[ ]);title('DCT cofficient matrix');[m,n]=size(dctF1);%%产生水印序列并对其排序radon('copyright',10);watermark1=radon(LENGTH,1);subplot(2,2,3);title('watermark seqence')[Y0,I0]=sort(watermark1);%%找出水印嵌入位置（幅值较大的n个频域成分）A=dctF1(:);[Y1,I1]=sort(A);x=mn;k=LENGTH;M=zeros(x,1);%%修改幅值较大的n个频域成分的幅值，嵌入水印（因为两个问题不同，所以有两个注释符）for i=1:xif k>=1M(x)=Y1(x)(1+alfaY0(k));k=k-1;elseM(x)=Y1(x);endx=x-1;endN=zeros(x,1);x=mn;for i=1:xN(I1(i))=M(i);enda=1;for j=1:nfor i=1:mdctF2(i,j)=N(a);a=a+1;endend%%DCT反变换，得到嵌入水印的图象idctF1=idct2(dctF2);subplot(2,2,4);imshow(idctF1,[ ]);title('嵌入水印后的图象');end二、提取恢复水印信息的MATLAB程序水印提取过程是水印嵌入过程的逆过程，相对嵌入过程来说比较复杂，难度较大，下面是水印提取检测的MATLAB程序代码：function watermark_detect(image,Y1,I0,waterMark1)%image:嵌入水印的图象%Y1:原始图象的序列排序%I0:原始水印的序列排序%waterMark1:原始水印序列%%对嵌入水印图象进行DCT变化dctW1=dct2(image);%%找出幅值较大的系数B=dtW1(:);[Y1,I2]=sort(B);[m1,n1]=size(dctW1);y=m1n1;k=length(waterMark1);N0=zeros(k,1);%%提取水印序列while k>=1N0(k)=(Y2(y)-Y1(y))/alfa/Y1(y);k=k-1;y=y-1;endk=length(waterMark1);waterMark2=zeros(k,1);for i=1:kwaterMark2(I0(i))=N0(i);end%%选取50个测试序列，其中第10个为提取出的水印figure;for i=1;50if i==10;waterMark=waterMark2;elsewaterMark=rand(k,1);end%计算各个序列与原来水印序列的相关值c=waterMark'waterMark1/sqrt(waterMark'waterMark);stem(i,c);hold on;end%三、接下来对嵌入水印的图象进行不同的攻击，用以测试水印的鲁棒性。

2、程序的目的和程序代码如下：%%攻击实验disp('input you cho according to the followingimage processing operation:');disp('0--exit');disp('1--oothing patterns');%添加噪音disp('2--adding uniorm noise 添加噪音');%滤波disp('3--adding filter [10 10] 滤波');%剪切disp('4--cutting part of the image 剪切');%压缩disp('5--10 quality JPEG compressing 压缩');%旋转45度disp('6--rotate 45 旋转');%d=input('please input you cho(请输入您的选择):');while d~=0switch dcase 1watermark_detect(idctF1,Y1,I0,waterMark1);case 2WImage2=idctF1;noise0=10rand(size(WImage2));WImage2=WImage2+noise0;figure;imshow(WImage2,[ ]);title('adding uniform noise 添加噪音');watemark_detect(WImage2,Y1,I0,waterMark1);case 3WImage3=idctF1;H=fspcial('gaussian高斯',[10,10],5);WImage3=imfilter(WImage3,H);figure;imshow(WImage3,[ ]);title(through filter [10,10] 滤波');watemark_detect(WImage3,Y1,I0,waterMark1);case 4WImage4=idctF1; WImage4(1:128,1;128)=256;figure;imshow(WImage4);title('cutting part of the image 剪切');watemark_detect(WImage4,Y1,I0,waterMark1);case 5WImage5=idctF1;WImage5=im2double(WImage5);cnum=10;dctm=dctmtx(8);p1=dctm;p2=dctm.';imageDCT=blkproc(WImage5,[8,8],'p1p2x',dctm,dctm.');DCTvar=im2col(imageDCT,[8,8],'distinct').';n=size(DCTvar,1);DCTvar=(sum(DCTvar.DCTvar)-(sum(DCTvar)/n).^2)/n;[dum,order]=sort(DCTvar);cnum=64-cnum;mask=ones(8,8);mask(order(1:cnum))=zeros(1,cnum);im88=zeros(9,9);im88(1:8,1:8)=mask;im128128=kron(im88(1:8,1:8),ones(16));dctm=dctmtx(8);p1=dctm.';p2=mask(1;8,1:8);p3=dctm;Wimage5=bikproc(imageDCT,[8,8],'p1(x.8p2)p3',dctm.',mask(1:8,1:8),dctm);figure;imshow(Wimage5);title('JPEG Image 压缩');watemark_detect(WImage5,Y1,I0,waterMark1);case 6 WImage6=idctF1;WImage6=imrotate(WImage6,45,'bilinear','corp');figure;imshow(Wimage6);title('rotate 45 旋转');watemark_detect(WImage6,Y1,I0,waterMark1);case 0break;otherwiseerror('you he a valid value(您的输入错误)');endd=input('please input you cho(请输入您的选择):');end%结束。

本文到这结束，希望上面文章对大家有所帮助。