遗传算法程序定义.pdfVIP

遗传算法 （一）目的： 通过遗传算法实现函数 F(x1,x2)=100*(x1^2-x2)^2+(1-x2)^2 的最优解， 其中 x1,x2 在 -2.048 和 2.048 之间。 （二）程序代码： clear clc G=100; % 循环次数 %/////////////////// 随机初始化 ////////////////////////////////////////////////////// E1=rand(80,20);E=round(E1); for p=1:G % 迭代次数 q=size(E); row=q(1); for i=1:row for j=1:10 x1(i,j)=E(i,j); end end for i=1:row for j=11:20 x2(i,(j-10))=E(i,j); end end %***** 主要是编程实现二进制换十进制 ****** for i=1:row c=0; for j=1:10 a=x1(i,j)*2^(10-j); c=c+a; y1(i)=c;% 十进制数据存 y1 end end for i=1:row d=0; for j=1:10 b=x2(i,j)*2^(10-j); d=d+b; y2(i)=d;% 十进制数据存 y2 end end %***** 编码（转换区间） ***** for i=1:row x11(i)=y1(i)/1023*4.096-2.048; x22(i)=y2(i)/1023*4.096-2.048; end %/////////////////////////////////////// 选择 ///////////////////////////////////////////////////// for i=1:row f(i)=100*(x11(i)^2-x22(i))^2+(1-x11(i))^2;% 求 f(i) 值 end [Or,idx1]=sort(f);% 对 f 排序 f1(p)=Or(row);%f 最大值 r1(p)=x11(idx1(row));%f 最大时对应 x11 的与 x22 值 r2(p)=x22(idx1(row)); %***** 求 f 平均值 ***** s=0; for i=1:row s=s+f(i); end avr=s/row;% 平均适应度 %///////////////////////////////////////// 繁殖 /////////////////////////////////////////////////// for i=1:row k(i)=f(i)/avr;% 选中概率 g(i)=round(k(i));% 四舍五入 end %///////////////////////////////////////// 新建矩阵 Et/////////////////////////////////////// k=1; for i=1:row n=g(i); if n0 for j=1:n for l=1:20 Et(k,l)=E(i,l); end k=k+1; end end end %////////////////////////////////////////// 交叉 ///////////////////////////////////////// o1=size(Et); for i=1:o1 r3=rand(1); if r3=0.6