 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1990 | 回復: 5 | ||
dellluoluo金蟲 (小有名氣)
|
[求助]
求助：資源優(yōu)化問題的遺傳算法matlab編程
|
|
最近在研究項目資源優(yōu)化問題，其中一個問題即工期固定，資源均衡問題，用遺傳算法來做，使用matlab的遺傳算法工具箱自己編程來實現(xiàn)，有詳細的文獻和例子，只做單資源的，需要將目標函數(shù)寫成.m 文件， 遇到困難：目標函數(shù)是使得資源均衡，滿足工序的合適開始時間即基因值 ，困難是基因值不體現(xiàn)在目標函數(shù)的顯式表達式中，而是在約束條件中，不知道怎么處理？ 如果有哪位蟲友感興趣并愿意提供幫助，將不勝感激！如果需要，我可以把例子附上。 圖1.gif 圖2.gif 圖片不清晰，附上文獻的鏈接，文中的測試問題1，圖5 http://www.sciencedirect.com/sci ... i/S0926580599000114 [ Last edited by dellluoluo on 2013-4-2 at 17:45 ] |
回復此樓» 猜你喜歡
281求調(diào)劑（0805）
已經(jīng)有8人回復
-
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名
已經(jīng)有3人回復
-
材料專碩306英一數(shù)二
已經(jīng)有10人回復
-
301求調(diào)劑
已經(jīng)有6人回復
-
一志愿天津大學化學工藝專業(yè)（081702）315分求調(diào)劑
已經(jīng)有7人回復
-
302求調(diào)劑
已經(jīng)有6人回復
-
26博士申請
已經(jīng)有3人回復
-
268求調(diào)劑
已經(jīng)有3人回復
-
311求調(diào)劑
已經(jīng)有10人回復
-
被我言中：新模板不強調(diào)格式了，假專家開始管格式了
已經(jīng)有4人回復
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- matlab論文分享（包括神經(jīng)網(wǎng)絡(luò)，遺傳算法，粒子群算法）
已經(jīng)有442人回復
- matlab中是否可以用simulink實現(xiàn)遺傳算法
已經(jīng)有6人回復
- 如何MATLAB程序?qū)崿F(xiàn)遺傳算法
已經(jīng)有11人回復
- 一個關(guān)于matlab的遺傳算法的問題
已經(jīng)有3人回復
- 求助有關(guān)Matlab有約束非線性最優(yōu)化問題
已經(jīng)有9人回復
- 【分享】matlab遺傳算法工具箱與應(yīng)用（雷英杰）【已搜索無重復】
已經(jīng)有74人回復
- 【求助】matlab 二次規(guī)劃的優(yōu)化的問題
已經(jīng)有4人回復
- 【求助】matlab 遺傳算法工具箱繪圖命令是什么？
已經(jīng)有6人回復
libralibra
至尊木蟲 (著名寫手)
驃騎將軍
- 程序強帖: 40
- 應(yīng)助: 817 (博后)
- 金幣: 12914.1
- 紅花: 64
- 帖子: 2238
- 在線: 287.3小時
- 蟲號: 696514
- 注冊: 2009-02-05
- 專業(yè): 計算機軟件
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
dellluoluo(ben_ladeng代發(fā)): 金幣+8, 應(yīng)樓主要求代發(fā)金幣 2013-06-01 13:31:46
ben_ladeng: 金幣+2, 謝謝專家 2013-06-01 13:32:36
ben_ladeng: 專家考核, 代發(fā)金幣 2013-06-01 13:33:33
ben_ladeng: 2013-06-01 13:33:59
dellluoluo(ben_ladeng代發(fā)): 金幣+2 2013-06-01 13:34:11
感謝參與，應(yīng)助指數(shù) +1
dellluoluo(ben_ladeng代發(fā)): 金幣+8, 應(yīng)樓主要求代發(fā)金幣 2013-06-01 13:31:46
ben_ladeng: 金幣+2, 謝謝專家 2013-06-01 13:32:36
ben_ladeng: 專家考核, 代發(fā)金幣 2013-06-01 13:33:33
ben_ladeng: 2013-06-01 13:33:59
dellluoluo(ben_ladeng代發(fā)): 金幣+2 2013-06-01 13:34:11
| GA都一個流程,約束一般都跟penalty function有關(guān)系,如果不滿足約束條件,增加penalty值就會讓這些'非優(yōu)良'個體避免遺傳給下一代. |
2樓2013-04-02 18:37:53
dellluoluo
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1398.8
- 散金: 21
- 帖子: 146
- 在線: 238.1小時
- 蟲號: 550976
- 注冊: 2008-04-25
- 專業(yè): 項目管理
3樓2013-04-03 10:46:51
dellluoluo
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1398.8
- 散金: 21
- 帖子: 146
- 在線: 238.1小時
- 蟲號: 550976
- 注冊: 2008-04-25
- 專業(yè): 項目管理
|
想得到您的相助！ %programm: %計算工期固定，資源均衡的程序，為資源均衡之前的結(jié)果，得到RLI=76 clear; T0=16;%固定工期，由project軟件求得關(guān)鍵路徑 n=11;%活動數(shù) for j=1:T0 for i=1:n ri=[2 1 4 4 2 4 6 0 2 1 2];%每個活動日資源 % ri=[2 1 4 4] di=[2 4 1 4 3 6 6 1 4 5 1];%每個活動持續(xù)工期 %di=[2 4 1 4] esi=[0 0 0 2 4 4 1 7 7 10 15];%每個活動最早開始時間 %esi=[0 0 0 2] efi=esi+di;%每個活動最早完成時間 xij=zeros(i,j);%系數(shù)矩陣 r_=sum(ri*di')/T0;%每個活動平均日資源需求 ti=esi;%每個活動實際開始時間 rij=0; end; end; for j=1:T0 for i=1:n if (j>= ti(i)+1)&(j<=ti(i)+di(i)) xij(i,j)=1; rij=ri*xij; tii(i)=ti(i); RLI=sum(abs(rij-r_));%目標函數(shù),絕對誤差和 end end end x=1:T0; bar(x,rij);%直方圖 xlabel('Time'); ylabel('Daily Resource Demand Before Leveling'); 上面這個已經(jīng)得到結(jié)果！下面這個是想要得到的結(jié)果： %programm:想要得到的結(jié)果，即想要通過遺傳算法的計算，找到一組合適的基因值，也即是 %每個活動合適的實際開始時間ti屬于[esi(i);esi(i)+tfi(i)],得到如下的結(jié)果，使RLI=25。 clear; T0=16;%固定工期 n=11;%活動數(shù) for j=1:T0 for i=1:n ri=[2 1 4 4 2 4 6 0 2 1 2];%每個活動日資源 % ri=[2 1 4 4] di=[2 4 1 4 3 6 6 1 4 5 1];%每個活動持續(xù)工期 %di=[2 4 1 4] esi=[0 0 0 10 7 4 1 14 11 10 15];%每個活動最早開始時間 %esi=[0 0 0 2] efi=esi+di;%每個活動最早完成時間 xij=zeros(i,j);%系數(shù)矩陣 r_=sum(ri*di')/T0;%每個活動平均日資源需求 ti=esi;%每個活動實際開始時間%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %（需要經(jīng)過遺傳算法計算后得到的結(jié)果)% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% rij=0; end; end; for j=1:T0 for i=1:n if (j>= ti(i)+1)&(j<=ti(i)+di(i)) xij(i,j)=1; rij=ri*xij; tii(i)=ti(i); RLI=sum(abs(rij-r_));%目標函數(shù) end end end x=1:T0; bar(x,rij);%直方圖 xlabel('Time'); ylabel('Daily Resource Demand After Leveling'); 下面想用GA求：程序沒有寫明白？？？？？？？？？？？？？？？？？ %GA programm T0=16;%固定工期 %T0=6; n=11;%活動數(shù) for j=1:T0 for i=1:n ri=[2 1 4 4 2 4 6 0 2 1 2];%每個活動日資源 % ri=[2 1 4 4] di=[2 4 1 4 3 6 6 1 4 5 1];%每個活動持續(xù)工期 %di=[2 4 1 4] esi=[0 0 0 2 4 4 1 7 7 10 15];%每個活動最早開始時間(遺傳算法應(yīng)用之前) %esi=[0 0 0 2] efi=esi+di;%每個活動最早完成時間 lsi=[6 0 3 10 8 4 4 14 11 10 15];%最遲開始時間 tfi=lsi-esi;%總時差 xij=zeros(i,j);%系數(shù)矩陣 r_=sum(ri*di')/T0;%每個活動平均日資源需求 ti=esi;%每個活動實際開始時間 rij=0; end; end; for j=1:T0 for i=1:n if (j>= ti(i)+1)&(j<=ti(i)+di(i)) xij(i,j)=1; rij=ri*xij; tii(i)=ti(i); RLI=sum(abs(rij-r_));%目標函數(shù)，絕對誤差和 end end end %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %想把上面這一段寫成function的形式來調(diào)用！不清楚怎么寫？% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %定義 遺傳算法參數(shù) NVAR=11;%變量維數(shù) NIND=50;%基因規(guī)模 PRECI=20;%基因長度 GGAP=0.9;%代溝 Pc=0.5;%交叉率 Pm=0.4;%變異率 MAXGEN=100;%演化代數(shù) %建立區(qū)域掃描器%%%%%%%%%%%% %想在下面這段遺傳算法的程序基礎(chǔ)上修改，這個程序是雷英杰那個書上的的 %7.3求多元多峰函數(shù)優(yōu)化的實例%%%%%%%%%%%%%，不會寫目標函數(shù)？？？？ %這個目標函數(shù)在遺傳算法中不知道怎么寫 RLI=sum(abs(rij-r_)) %關(guān)鍵是基因值沒有體現(xiàn)在目標函數(shù)中，只是隱含在xij中，當ti在活動的 %【esi，esi+tfi】之間時，xij=1，才計算RLI. else xij=0不計算 %最終得到的結(jié)果是找到一組合適的ti，使得計算的RLI=25 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% trace = zeros(MAXGEN,NVAR); for i=1:NVAR FieldD = [rep([PRECI],[1,NVAR]);rep([esi(i);esi(i)+tfi(i)],[1,NVAR]);rep([1;0;1;1],[1,NVAR])]; end; Chrom=crtbp(NIND,NVAR*PRECI);%創(chuàng)建初始種群 gen = 0; trace=zeros(MAXGEN,2); x=bs2rv(Chrom,FieldD) ObjV= RLI;%%%%%%%%%%%%%%%%%%這個函數(shù)不會寫？ while gen SelCh=select('sus',Chrom,FitnV,GGAP); SelCh=recombin('xovsp',SelCh,0.5); SelCh=mut(SelCh,0.4); x=bs2rv(Chrom,FieldD) ObjVSel=RLI;%%%%%%%%%%%%%%%%%%%%%%%%？？？？？？？？ [Chrom ObjV]=reins(Chrom,SelCh,1,1,ObjV,ObjVSel); gen=gen+1; %輸出最優(yōu)解 [Y,I] = min(ObjVSel); Y,bs2rv(Chrom(I,  ,FieldD);% plot(variable(I),Y,'bo'); trace(gen,1)=min(ObjV); trace(gen,2)=sum(ObjV)/length(ObjV); if (gen==100) figure(2); plot(ObjV);hold on; plot(ObjV,'b*');grid; end end figure(3);clf; plot(trace(:,1));hold on; plot(trace(:,2)','-.');grid; %legend('aa','bb').  |
4樓2013-04-08 13:57:39
dellluoluo
金蟲 (小有名氣)
- 應(yīng)助: 1 (幼兒園)
- 金幣: 1398.8
- 散金: 21
- 帖子: 146
- 在線: 238.1小時
- 蟲號: 550976
- 注冊: 2008-04-25
- 專業(yè): 項目管理
5樓2013-04-10 21:14:10
6樓2015-11-21 22:10:49
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 dellluoluo 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|---|---|---|---|---|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[考研] 085601材料工程專碩求調(diào)劑 +5 | 慕寒mio 2026-03-16 | 5/250 |
|
|
[考研] 299求調(diào)劑 +4 | △小透明* 2026-03-17 | 4/200 |
|
|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿，福州大學材料專碩339分求調(diào)劑 +3 | 木子momo青爭 2026-03-15 | 3/150 |
|
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
|
[考研] 藥學383 求調(diào)劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +6 | 大可digkids 2026-03-16 | 6/300 |
|
|
|
[考研] 326求調(diào)劑 +4 | 諾貝爾化學獎覬?/a> 2026-03-15 | 7/350 |
|
|
|
[考研] 070300化學學碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 085600材料與化工 求調(diào)劑 +13 | enenenhui 2026-03-13 | 14/700 |
|
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
|
[考研] 材料工程調(diào)劑 +9 | 咪咪空空 2026-03-12 | 9/450 |
|
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考研] 0817化學工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|