版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
1、精选优质文档-倾情为你奉上非线性整数规划的遗传算法Matlab程序(附图)通常,非线性整数规划是一个具有指数复杂度的NP问题,如果约束较为复杂,Matlab优化工具箱和一些优化软件比如lingo等,常常无法应用,即使能应用也不能给出一个较为令人满意的解。这时就需要针对问题设计专门的优化算法。下面举一个遗传算法应用于非线性整数规划的编程实例,供大家参考!模型的形式和适应度函数定义如下:这是一个具有200个01决策变量的多目标非线性整数规划,编写优化的目标函数如下,其中将多目标转化为单目标采用简单的加权处理。function Fitness=FITNESS(x,FARM,e,q,w)% 适应度函数
2、% 输入参数列表% x 决策变量构成的4×50的0-1矩阵% FARM 细胞结构存储的当前种群,它包含了个体x% e 4×50的系数矩阵% q 4×50的系数矩阵% w
3、60; 1×50的系数矩阵%gamma=0.98;N=length(FARM);%种群规模F1=zeros(1,N);F2=zeros(1,N);for i=1:N xx=FARMi; ppp=(1-xx)+(1-q).*xx; F1(i)=sum(w.*prod(ppp); F2(i)=sum(sum(e.*xx);endppp=(1-x)+(1-q).*x;f1=sum(w.*prod(pp
4、p);f2=sum(sum(e.*x);Fitness=gamma*sum(min(sign(f1-F1);zeros(1,N)+(1-gamma)*sum(min(sign(f2-F2);zeros(1,N); 针对问题设计的遗传算法如下,其中对模型约束的处理是重点考虑的地方function Xp,LC1,LC2,LC3,LC4=MYGA(M,N,Pm)% 求解01整数规划的遗传算法% 输入参数列表% M 遗传进化迭代次数% N 种群规模% Pm
5、60; 变异概率% 输出参数列表% Xp 最优个体% LC1 子目标1的收敛曲线% LC2 子目标2的收敛曲线% LC3 平均适应度函数的收敛曲线% LC4 最优适应度函数的收敛曲线% 参考调用格式Xp,LC1,LC2,LC3,LC4=MYGA(50,40,0.3)% 第一步:载入数据和变量初始化load eqw;%载入三个系数矩阵e,q,w%输出变量初始化Xp=zeros(4,5
6、0);LC1=zeros(1,M);LC2=zeros(1,M);LC3=zeros(1,M);LC4=zeros(1,M);Best=inf;% 第二步:随机产生初始种群farm=cell(1,N);%用于存储种群的细胞结构k=0;while k %以下是一个合法个体的产生过程 x=zeros(4,50);%x每一列的1的个数随机决定 for i=1:50 R=rand;
7、0; Col=zeros(4,1); if R<0.7 RP=randperm(4);%1的位置也是随机的 Col(RP(1)=1;
8、60; elseif R>0.9 RP=randperm(4); Col(RP(1:2)=1; else
9、; RP=randperm(4); Col(RP(1:3)=1; end x(:,i)=Col; end %下面是检查行和是否满足约束的过程,对于不满足约束的予以抛弃 Temp1=sum(x,
10、2); Temp2=find(Temp1>20); if length(Temp2)=0 k=k+1; farmk=x; endend % 以下是进化迭代过程counter=0;%设置迭代计数器while counter % 第三步:交叉
11、 %交叉采用双亲双子单点交叉 newfarm=cell(1,2*N);%用于存储子代的细胞结构 Ser=randperm(N);%两两随机配对的配对表 A=farmSer(1);%取出父代A B=farmSer(2);%取出父代B P0=unidrnd(49);%随机选择交叉点 a=A(:,1:P0),B(:,(P0+1):end);%产生子代a
12、 b=B(:,1:P0),A(:,(P0+1):end);%产生子代b newfarm2*N-1=a;%加入子代种群 newfarm2*N=b; %以下循环是重复上述过程 for i=1:(N-1) A=farmSer(i); B=farmSer(i+1);
13、60; P0=unidrnd(49); a=A(:,1:P0),B(:,(P0+1):end); b=B(:,1:P0),A(:,(P0+1):end); newfarm2*i-1=a; ne
14、wfarm2*i=b; end FARM=farm,newfarm;%新旧种群合并 % 第四步:选择复制 FLAG=ones(1,3*N);%标志向量,对是否满足约束进行标记 %以下过程是检测新个体是否满足约束 for i=1:(3*N) x=FARMi;
15、; sum1=sum(x,1); sum2=sum(x,2); flag1=find(sum1=0); flag2=find(sum1=4); flag3=find(sum2>20); &
16、#160; if length(flag1)+length(flag2)+length(flag3)>0 FLAG(i)=0;%如果不满足约束,用0加以标记 end end NN=length(find(FLAG)=1);%满足约束的个体数目,它一定大于等于N
17、; NEWFARM=cell(1,NN); %以下过程是剔除不满主约束的个体 kk=0; for i=1:(3*N) if FLAG(i)=1 kk=kk+1;
18、0; NEWFARMkk=FARMi; end end %以下过程是计算并存储当前种群每个个体的适应值 SYZ=zeros(1,NN); syz=zeros(1,N); for i=1:NN x=NEWFARMi;&
19、#160; SYZ(i)=FITNESS2(x,NEWFARM,e,q,w);%调用适应值子函数 end k=0; %下面是选择复制,选择较优的N个个体复制到下一代 while k minSYZ=min(SYZ);
20、posSYZ=find(SYZ=minSYZ); POS=posSYZ(1); k=k+1; farmk=NEWFARMPOS; syz(k)=SYZ(POS); SYZ(POS)
21、=inf; end %记录和更新,更新最优个体,记录收敛曲线的数据 minsyz=min(syz); meansyz=mean(syz); pos=find(syz=minsyz); LC3(counter+1)=meansyz; if minsyz Best=m
22、insyz; Xp=farmpos(1); end LC4(counter+1)=Best; ppp=(1-Xp)+(1-q).*Xp; LC1(counter+1)=sum(w.*prod(ppp); LC2(counter+1)=sum(sum(e.*Xp); % 第五步:变异
23、0; for i=1:N if Pm>rand%是否变异由变异概率Pm控制 AA=farmi;%取出一个个体 POS=unidrnd(50);%随机选择变异位
24、160; R=rand; Col=zeros(4,1); if R<0.7 RP=randperm(4
25、); Col(RP(1)=1; elseif R>0.9 RP=randperm(4); &
26、#160; Col(RP(1:2)=1; else RP=randperm(4); &
27、#160; Col(RP(1:3)=1; end %下面是判断变异产生的新个体是否满足约束,如果不满足,此次变异无效 AA(:,POS)=Col; T
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 云南省昆明八中2025年高考物理模拟试卷(含答案)
- 安徽省联考2024-2025学年高一上学期1月期末英语试题(无答案)
- 生物-内蒙古鄂尔多斯市西四旗2024-2025学年2025届高三第一(上)学期期末联考试卷试题和答案
- 2024重庆离婚协议书格式范文
- 2024自动售货机租赁合同
- 2025年度5G通信基站建设与维护合同范本3篇
- 2025年度宠物猫寄养与宠物用品租赁一体化协议3篇
- 2024颐和园职工制服采购环保责任及回收利用合同3篇
- 2024甲乙双方关于购买房产的合同
- 2024版公私合营养老院服务协议示范文本版B版
- 2024年萍乡卫生职业学院单招职业适应性测试题库参考答案
- 中国农业银行信用借款合同
- ISO 56001-2024《创新管理体系-要求》专业解读与应用实践指导材料之9:“5领导作用-5.3创新战略”(雷泽佳编制-2025B0)
- 江苏省连云港市2023-2024学年八年级上学期期末数学试题(原卷版)
- 初中英语听力高频词
- 2024智能变电站新一代集控站设备监控系统技术规范部分
- 2024年建筑业10项新技术
- 语文七年级下字帖打印版
- 二年级上册音乐课件---选唱歌曲-我们和祖国最亲亲-西师大版(共8张PPT)
- 设备租赁服务方案
- 最新中石油带压作业技术规程
评论
0/150
提交评论