function?[SlopeScheduleNewTimeArrayStartTimeFinishTimeMakespan]?=?Ch7_Palmer（TimeArray）
%%?程序功能說明
%第7章：帕爾瑪算法計(jì)算工件的近優(yōu)排序和最長流程時(shí)間，并繪制甘特圖
%====輸出參數(shù)====
???%Slope???????????：工件的斜度
???%Schedule????????：近優(yōu)排序
???%NewTimeArray????：近優(yōu)排序下的加工時(shí)間矩陣
???%StartTime???????：開始時(shí)間
???%FinishTime??????：完工時(shí)間
???%Makespan????????：最長流程時(shí)間
%====輸入?yún)?shù)====
???%TimeArray???????：加工時(shí)間矩陣
%編寫時(shí)間：2012年04月；完善時(shí)間：2012年12月。

%%?程序主體部分
%Get?the?near?optimal?schedule?by?the?Palmer?rule
%（1）根據(jù)帕爾瑪法則確定工件的近優(yōu)加工排序
[mn]?=?size（TimeArray）;??????????????%加工時(shí)間矩陣的行數(shù)和列數(shù)，行數(shù)m為機(jī)器數(shù)，列數(shù)n為工件數(shù)
if?m?<=?1
????error（‘機(jī)器數(shù)必須大于等于2‘）
end

for?i?=?1:n
????SlopeIndex（i）?=?0;????????????????????????????????????????????????%工件i的斜度指標(biāo)的初值
????for?j=1:m
????????SlopeIndex（i）?=?SlopeIndex（i）+0.5*（2*j-m-1）*TimeArray（ji）;???%計(jì)算工件i的斜度指標(biāo)
????end
end
TimeArray（m+1:）?=?1:n;???????????????????????????????????????????%第m+1行為工件的代號(hào)
[SlopeBestIndex]?=?sort（SlopeIndex‘descend‘）;???????????????????%給出斜度及斜度的非遞增順序——近優(yōu)加工順序
Schedule?=?BestIndex;?????????????????????????????????????????????%工件的近優(yōu)加工排序
TimeArray?=?TimeArray（:BestIndex）;???????????????????????????????%根據(jù)近優(yōu)加工順序重新排列加工時(shí)間矩陣
NewTimeArray?=?TimeArray（1:m:）;??????????????????????????????????%重排后的加工時(shí)間矩陣

%（2）計(jì)算開工時(shí)間、完工時(shí)間、最長流程加工時(shí)間
WorkTime（1:m:）?=?TimeArray（1:m:）;???????????????????????????????%重新定義的加工時(shí)間矩陣
StartTime（11）?=?0;??????????????????????????????????????????