资源描述：
《非线性整数规划的遗传算法matlab程序》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。
1、羇肈蒆薄蚆芃莂薃蝿肆芈蚂袁芁膄蚁羃肄蒃蚀蚃袇葿虿袅膂莅虿羈羅芁蚈蚇膁膇蚇螀羄蒅蚆袂腿莁螅羄羂芇螄蚄膇膃螃螆羀薂螃羈芆蒈螂肁肈莄螁螀芄芀莇袃肇膆莇羅节蒅蒆蚅肅莁蒅螇芀芇蒄衿肃节蒃肂袆薁蒂螁膂蒇蒁袄羄莃蒁羆膀艿蒀蚅羃膅蕿螈膈蒄薈袀羁莀薇肂膆莆薆螂聿节薅袄芅膈薅羇肈蒆薄蚆芃莂薃蝿肆芈蚂袁芁膄蚁羃肄蒃蚀蚃袇葿虿袅膂莅虿羈羅芁蚈蚇膁膇蚇螀羄蒅蚆袂腿莁螅羄羂芇螄蚄膇膃螃螆羀薂螃羈芆蒈螂肁肈莄螁螀芄芀莇袃肇膆莇羅节蒅蒆蚅肅莁蒅螇芀芇蒄衿肃节蒃肂袆薁蒂螁膂蒇蒁袄羄莃蒁羆膀艿蒀蚅羃膅蕿螈膈蒄薈袀羁莀薇肂膆莆薆螂聿节薅袄芅膈薅羇肈蒆薄蚆芃莂薃蝿肆芈蚂袁芁膄蚁羃肄蒃蚀蚃袇葿虿
2、袅膂莅虿羈羅芁蚈蚇膁膇蚇螀羄蒅蚆袂腿莁螅羄羂芇螄蚄膇膃螃螆羀薂螃羈芆蒈螂肁肈莄螁螀芄芀莇袃肇膆莇羅节蒅蒆蚅肅莁蒅螇芀芇蒄衿肃节蒃肂袆薁蒂螁膂蒇蒁袄羄莃蒁羆膀艿蒀蚅羃膅蕿螈膈蒄薈袀羁莀薇肂膆莆薆螂聿节薅袄芅膈薅羇肈蒆薄蚆芃莂薃蝿肆芈蚂袁芁膄蚁羃肄蒃蚀蚃袇葿虿袅膂莅虿羈羅芁蚈蚇膁膇非线性整数规划的遗传算法Matlab程序（引自人工智能及数据挖掘论坛）这是一个具有200个01决策变量的多目标非线性整数规划，编写优化的目标函数如下，其中将多目标转化为单目标采用简单的加权处理。functionFitness=FITNESS(x,FARM,e,q,w)%%适应度函数
3、%输入参数列表%x决策变量构成的4×50的0-1矩阵%FARM细胞结构存储的当前种群，它包含了个体x%e4×50的系数矩阵%q4×50的系数矩阵%w1×50的系数矩阵%%gamma=0.98;N=length(FARM);%种群规模F1=zeros(1,N);F2=zeros(1,N);fori=1:Nxx=FARM{i};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(ppp));f2=sum(
4、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变异概率%%输出参数列表%Xp最优个体%LC1子目标1的收敛曲线%LC2子目标2的收敛曲线%LC3平均适应度函数的收敛曲
5、线%LC4最优适应度函数的收敛曲线%%参考调用格式[Xp,LC1,LC2,LC3,LC4]=MYGA(50,40,0.3)%%第一步：载入数据和变量初始化loadeqw;%载入三个系数矩阵e,q,w%输出变量初始化Xp=zeros(4,50);LC1=zeros(1,M);LC2=zeros(1,M);LC3=zeros(1,M);LC4=zeros(1,M);Best=inf;%%第二步：随机产生初始种群farm=cell(1,N);%用于存储种群的细胞结构k=0;whilek%以下是一个合法个体的产生过程x=zeros(4,50);%x每一列的1的个
6、数随机决定fori=1:50R=rand;Col=zeros(4,1);ifR<0.7RP=randperm(4);%1的位置也是随机的Col(RP(1))=1;elseifR>0.9RP=randperm(4);Col(RP(1:2))=1;elseRP=randperm(4);Col(RP(1:3))=1;endx(:,i)=Col;end%下面是检查行和是否满足约束的过程，对于不满足约束的予以抛弃Temp1=sum(x,2);Temp2=find(Temp1>20);iflength(Temp2)==0k=k+1;farm{k}=x;endend%
7、%以下是进化迭代过程counter=0;%设置迭代计数器whilecounter%%第三步：交叉%交叉采用双亲双子单点交叉newfarm=cell(1,2*N);%用于存储子代的细胞结构Ser=randperm(N);%两两随机配对的配对表A=farm{Ser(1)};%取出父代AB=farm{Ser(2)};%取出父代BP0=unidrnd(49);%随机选择交叉点a=[A(:,10),B(:,(P0+1):end)];%产生子代ab=[B(:,10),A(:,(P0+1):end)];%产生子代bnewfarm{2*N-1}=a;%加入子代种群new
8、farm{2*N}=b;%以下循环是重复上述过程fori=1N-1)A=farm
