《流水線車間生產(chǎn)調(diào)度的遺傳算法》由會員分享,可在線閱讀,更多相關《流水線車間生產(chǎn)調(diào)度的遺傳算法(6頁珍藏版)》請在裝配圖網(wǎng)上搜索。
.
流水線車間生產(chǎn)調(diào)度的遺傳算法MATLAB源代碼
n個任務在流水線上進行m個階段的加工,每一階段至少有一臺機器且至少有一個階段存在多臺機器,并且同一階段上各機器的處理性能相同,在每一階段各任務均要完成一道工序,各任務的每道工序可以在相應階段上的任意一臺機器上加工,已知任務各道工序的處理時間,要求確定所有任務的排序以及每一階段上機器的分配情況,使得調(diào)度指標(一般求Makespan)最小。
function [Zp,Y1p,Y2p,Y3p,Xp,LC1,LC2]=JSPGA(M,N,Pm,T,P)
%--------------------------------------------------------------------------
% JSPGA.m
% 流水線型車間作業(yè)調(diào)度遺傳算法
% GreenSim團隊——專業(yè)級算法設計&代寫程序
% 歡迎訪問GreenSim團隊主頁→http://blog.sina.com.cn/greensim
%--------------------------------------------------------------------------
% 輸入?yún)?shù)列表
% M 遺傳進化迭代次數(shù)
% N 種群規(guī)模(取偶數(shù))
% Pm 變異概率
% T mn的矩陣,存儲m個工件n個工序的加工時間
% P 1n的向量,n個工序中,每一個工序所具有的機床數(shù)目
% 輸出參數(shù)列表
% Zp 最優(yōu)的Makespan值
% Y1p 最優(yōu)方案中,各工件各工序的開始時刻,可根據(jù)它繪出甘特圖
% Y2p 最優(yōu)方案中,各工件各工序的結束時刻,可根據(jù)它繪出甘特圖
% Y3p 最優(yōu)方案中,各工件各工序使用的機器編號
% Xp 最優(yōu)決策變量的值,決策變量是一個實數(shù)編碼的mn矩陣
% LC1 收斂曲線1,各代最優(yōu)個體適應值的記錄
% LC2 收斂曲線2,各代群體平均適應值的記錄
% 最后,程序還將繪出三副圖片:兩條收斂曲線圖和甘特圖(各工件的調(diào)度時序圖)
%第一步:變量初始化
[m,n]=size(T);%m是總工件數(shù),n是總工序數(shù)
Xp=zeros(m,n);%最優(yōu)決策變量
LC1=zeros(1,M);%收斂曲線1
LC2=zeros(1,N);%收斂曲線2
%第二步:隨機產(chǎn)生初始種群
farm=cell(1,N);%采用細胞結構存儲種群
for k=1:N
X=zeros(m,n);
for j=1:n
for i=1:m
X(i,j)=1+(P(j)-eps)*rand;
end
end
farm{k}=X;
end
counter=0;%設置迭代計數(shù)器
while counter
rand;%變異概率為Pm
X=farm{i};
I=unidrnd(m);
J=unidrnd(n);
X(I,J)=1+(P(J)-eps)*rand;
farm{i}=X;
end
end
farm{pos(1)}=Xp;
counter=counter+1
end
%輸出結果并繪圖
figure(1);
plotif=1;
X=Xp;
[Zp,Y1p,Y2p,Y3p]=COST(X,T,P,plotif);
figure(2);
plot(LC1);
figure(3);
plot(LC2);
function [Zp,Y1p,Y2p,Y3p]=COST(X,T,P,plotif)
% JSPGA的內(nèi)聯(lián)子函數(shù),用于求調(diào)度方案的Makespan值
% 輸入?yún)?shù)列表
% X 調(diào)度方案的編碼矩陣,是一個實數(shù)編碼的mn矩陣
% T mn的矩陣,存儲m個工件n個工序的加工時間
% P 1n的向量,n個工序中,每一個工序所具有的機床數(shù)目
% plotif 是否繪甘特圖的控制參數(shù)
% 輸出參數(shù)列表
% Zp 最優(yōu)的Makespan值
% Y1p 最優(yōu)方案中,各工件各工序的開始時刻
% Y2p 最優(yōu)方案中,各工件各工序的結束時刻
% Y3p 最優(yōu)方案中,各工件各工序使用的機器編號
%第一步:變量初始化
[m,n]=size(X);
Y1p=zeros(m,n);
Y2p=zeros(m,n);
Y3p=zeros(m,n);
%第二步:計算第一道工序的安排
Q1=zeros(m,1);
Q2=zeros(m,1);
R=X(:,1);%取出第一道工序
Q3=floor(R);%向下取整即得到各工件在第一道工序使用的機器的編號
%下面計算各工件第一道工序的開始時刻和結束時刻
for i=1:P(1)%取出機器編號
pos=find(Q3==i);%取出使用編號為i的機器為其加工的工件的編號
lenpos=length(pos);
if lenpos>=1
Q1(pos(1))=0;
if lenpos>=2
for j=2:lenpos
Q1(pos(j))=Q2(pos(j-1));
Q2(pos(j))=Q2(pos(j-1))+T(pos(j),1);
end
end
end
end
Y1p(:,1)=Q1;
Y3p(:,1)=Q3;
%第三步:計算剩余工序的安排
for k=2:n
R=X(:,k);%取出第k道工序
Q3=floor(R);%向下取整即得到各工件在第k道工序使用的機器的編號
%下面計算各工件第k道工序的開始時刻和結束時刻
for i=1:P(k)%取出機器編號
pos=find(Q3==i);%取出使用編號為i的機器為其加工的工件的編號
lenpos=length(pos);
if lenpos>=1
EndTime=Y2p(pos,k-1);%取出這些機器在上一個工序中的結束時刻
POS=zeros(1,lenpos);%上一個工序完成時間由早到晚的排序
for jj=1:lenpos
POS(jj)=ppp(1);
EndTime(ppp(1))=Inf;
end
%根據(jù)上一個工序完成時刻的早晚,計算各工件第k道工序的開始時刻和結束時刻
Q1(pos(POS(1)))=Y2p(pos(POS(1)),k-1);
Q2(pos(POS(1)))=Q1(pos(POS(1)))+T(pos(POS(1)),k);%前一個工件的結束時刻
if lenpos>=2
for j=2:lenpos
Q1(pos(POS(j)))=Y2p(pos(POS(j)),k-1);%預定的開始時刻為上一個工序的結束時刻
if Q1(pos(POS(j)))
下載提示(請認真閱讀)
- 1.請仔細閱讀文檔,確保文檔完整性,對于不預覽、不比對內(nèi)容而直接下載帶來的問題本站不予受理。
- 2.下載的文檔,不會出現(xiàn)我們的網(wǎng)址水印。
- 3、該文檔所得收入(下載+內(nèi)容+預覽)歸上傳者、原創(chuàng)作者;如果您是本文檔原作者,請點此認領!既往收益都歸您。
文檔包含非法信息?點此舉報后獲取現(xiàn)金獎勵!
下載文檔到電腦,查找使用更方便
10
積分
- 配套講稿:
如PPT文件的首頁顯示word圖標,表示該PPT已包含配套word講稿。雙擊word圖標可打開word文檔。
- 特殊限制:
部分文檔作品中含有的國旗、國徽等圖片,僅作為作品整體效果示例展示,禁止商用。設計者僅對作品中獨創(chuàng)性部分享有著作權。
- 關 鍵 詞:
-
流水線
車間
生產(chǎn)
調(diào)度
遺傳
算法
- 溫馨提示:
1: 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
2: 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權益歸上傳用戶所有。
3.本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預覽,若沒有圖紙預覽就沒有圖紙。
4. 未經(jīng)權益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
5. 裝配圖網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負責。
6. 下載文件中如有侵權或不適當內(nèi)容,請與我們聯(lián)系,我們立即糾正。
7. 本站不保證下載資源的準確性、安全性和完整性, 同時也不承擔用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
裝配圖網(wǎng)所有資源均是用戶自行上傳分享,僅供網(wǎng)友學習交流,未經(jīng)上傳用戶書面授權,請勿作他用。
鏈接地址:http://www.3dchina-expo.com/p-12932940.html