非線性整數(shù)規(guī)劃的遺傳算法matlab程序

《非線性整數(shù)規(guī)劃的遺傳算法matlab程序》由會(huì)員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在應(yīng)用文檔-天天文庫(kù)。

1、羇肈蒆薄蚆芃莂薃蝿肆羋螞袁芁膄蟻羃肄蒃蝕蚃袇葿蠆裊膂蒞蠆羈羅芁蚈蚇膁膇蚇螀羄蒅蚆袂腿莁螅羄羂芇螄蚄膇膃螃螆羀薂螃羈芆蒈螂肁肈莄螁螀芄芀莇袃肇膆莇羅節(jié)蒅蒆蚅肅莁蒅螇芀芇蒄衿肅節(jié)蒃肂袆薁蒂螁膂蕆蒁襖羄莃蒁羆膀艿蒀蚅羃膅蕿螈膈蒄薈袀羈莀薇肂膆莆薆螂聿節(jié)薅襖芅膈薅羇肈蒆薄蚆芃莂薃蝿肆羋螞袁芁膄蟻羃肄蒃蝕蚃袇葿蠆裊膂蒞蠆羈羅芁蚈蚇膁膇蚇螀羄蒅蚆袂腿莁螅羄羂芇螄蚄膇膃螃螆羀薂螃羈芆蒈螂肁肈莄螁螀芄芀莇袃肇膆莇羅節(jié)蒅蒆蚅肅莁蒅螇芀芇蒄衿肅節(jié)蒃肂袆薁蒂螁膂蕆蒁襖羄莃蒁羆膀艿蒀蚅羃膅蕿螈膈蒄薈袀羈莀薇肂膆莆薆螂聿節(jié)薅襖芅膈薅羇肈蒆薄蚆芃莂薃蝿肆羋螞袁芁膄蟻羃肄蒃蝕蚃袇葿蠆

2、裊膂蒞蠆羈羅芁蚈蚇膁膇蚇螀羄蒅蚆袂腿莁螅羄羂芇螄蚄膇膃螃螆羀薂螃羈芆蒈螂肁肈莄螁螀芄芀莇袃肇膆莇羅節(jié)蒅蒆蚅肅莁蒅螇芀芇蒄衿肅節(jié)蒃肂袆薁蒂螁膂蕆蒁襖羄莃蒁羆膀艿蒀蚅羃膅蕿螈膈蒄薈袀羈莀薇肂膆莆薆螂聿節(jié)薅襖芅膈薅羇肈蒆薄蚆芃莂薃蝿肆羋螞袁芁膄蟻羃肄蒃蝕蚃袇葿蠆裊膂蒞蠆羈羅芁蚈蚇膁膇非線性整數(shù)規(guī)劃的遺傳算法Matlab程序（引自人工智能及數(shù)據(jù)挖掘論壇）這是一個(gè)具有200個(gè)01決策變量的多目標(biāo)非線性整數(shù)規(guī)劃，編寫優(yōu)化的目標(biāo)函數(shù)如下，其中將多目標(biāo)轉(zhuǎn)化為單目標(biāo)采用簡(jiǎn)單的加權(quán)處理。functionFitness=FITNESS(x,FARM,e,q,w)%%適應(yīng)度函數(shù)

3、%輸入?yún)?shù)列表%x決策變量構(gòu)成的4×50的0-1矩陣%FARM細(xì)胞結(jié)構(gòu)存儲(chǔ)的當(dāng)前種群，它包含了個(gè)體x%e4×50的系數(shù)矩陣%q4×50的系數(shù)矩陣%w1×50的系數(shù)矩陣%%gamma=0.98;N=length(FARM);%種群規(guī)模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)]));針對(duì)問(wèn)題設(shè)計(jì)的遺傳算法如下，其中對(duì)模型約束的處理是重點(diǎn)考慮的地方function[Xp,LC1,LC2,LC3,LC4]=MYGA(M,N,Pm)%%求解01整數(shù)規(guī)劃的遺傳算法%%輸入?yún)?shù)列表%M遺傳進(jìn)化迭代次數(shù)%N種群規(guī)模%Pm變異概率%%輸出參數(shù)列表%Xp最優(yōu)個(gè)體%LC1子目標(biāo)1的收斂曲線%LC2子目標(biāo)2的收斂曲線%LC3平均適應(yīng)度函數(shù)的收斂曲

5、線%LC4最優(yōu)適應(yīng)度函數(shù)的收斂曲線%%參考調(diào)用格式[Xp,LC1,LC2,LC3,LC4]=MYGA(50,40,0.3)%%第一步：載入數(shù)據(jù)和變量初始化loadeqw;%載入三個(gè)系數(shù)矩陣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;%%第二步：隨機(jī)產(chǎn)生初始種群farm=cell(1,N);%用于存儲(chǔ)種群的細(xì)胞結(jié)構(gòu)k=0;whilek%以下是一個(gè)合法個(gè)體的產(chǎn)生過(guò)程x=zeros(4,50);%x每一列的1的個(gè)

6、數(shù)隨機(jī)決定fori=1:50R=rand;Col=zeros(4,1);ifR<0.7RP=randperm(4);%1的位置也是隨機(jī)的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%下面是檢查行和是否滿足約束的過(guò)程，對(duì)于不滿足約束的予以拋棄Temp1=sum(x,2);Temp2=find(Temp1>20);iflength(Temp2)==0k=k+1;farm{k}=x;endend%

7、%以下是進(jìn)化迭代過(guò)程counter=0;%設(shè)置迭代計(jì)數(shù)器whilecounter%%第三步：交叉%交叉采用雙親雙子單點(diǎn)交叉newfarm=cell(1,2*N);%用于存儲(chǔ)子代的細(xì)胞結(jié)構(gòu)Ser=randperm(N);%兩兩隨機(jī)配對(duì)的配對(duì)表A=farm{Ser(1)};%取出父代AB=farm{Ser(2)};%取出父代BP0=unidrnd(49);%隨機(jī)選擇交叉點(diǎn)a=[A(:,10),B(:,(P0+1):end)];%產(chǎn)生子代ab=[B(:,10),A(:,(P0+1):end)];%產(chǎn)生子代bnewfarm{2*N-1}=a;%加入子代種群new

8、farm{2*N}=b;%以下循環(huán)是重復(fù)上述過(guò)程fori=1N-1)A=farm