排課問題分析

排課問題分析摘要：本題要求我們對多約束條件的典型組合進行分析，求解，并作最優(yōu)化處理。基于此種原因，我們先對各個元素間的沖突做預(yù)處理，進行約束條件的規(guī)劃，再通過matlab軟件將教室、教師、課程和時間間的約束條件統(tǒng)一化，構(gòu)成R-T-C表（詳見附表），再將各個元素進行優(yōu)先級的計算，從而根據(jù)排課的優(yōu)化模型，求出最優(yōu)解。經(jīng)過對所給的表格，數(shù)據(jù)的深入分析，我們可以得知，教師明顯缺少，比如課程學(xué)時要求有160個課時，然而教師能上的課時僅有116個課時，所以開始排課時，不考慮教師，向教師中安排課程。由于同類課程最好不要放在一起，同時根據(jù)老師的需求和教室的開放時間進行分配，經(jīng)過與我們實際的課表的排課情況的分析，比如隔一天排同一課，課程類別不同的課程不在同一時間上課，我們可以大致的排出一個按教室上課的表，即R-T-C表。通過對R-T-C表的分析，發(fā)現(xiàn)有很多課沒老師上和老師沒課上的情況，我們就對其進行相應(yīng)的，合理的調(diào)整。最后發(fā)現(xiàn)還是老師要外聘。將外聘14名老師去上相應(yīng)沒人上的科目，具體情況見附表。最后，我們得到了一相對優(yōu)化的，以教室為準的課表（詳見附表），從而解決問題（1）的要求。對于我們課表的安排，發(fā)現(xiàn)再沒對晚自習(xí)有其他條件約束是不會對所排的課表有所影響。關(guān)鍵詞：排課問題組合規(guī)劃多目標(biāo)函數(shù)數(shù)據(jù)量化優(yōu)先級一、問題重述對于有課程40門，教師共有25名，教室18間的條件下合理的安排課程表，而課程、教師、教室的具體屬性及要求詳見附表（表1，表2,表3）對于課表德編排，題目有如下規(guī)則：每周以5天為單位進行編排；每天最多只能編排8節(jié)課（上午4節(jié)，下午4節(jié)），特殊情況下可以編排10節(jié)課（晚上2節(jié)），每門課程以2節(jié)課為單位進行編排，同類課程盡可能不安排在同一時間。要求所要解決的問題：請你結(jié)合實際情況建立數(shù)學(xué)模型，通過編程計算，給出較為合理的課表編排方案，分析你所給出的方案的合理性。如果不準晚上排課，排課結(jié)果是否有所變化，如何變化？對教師聘用，教室配置給出合理化建議。二、問題分析隨著現(xiàn)代教學(xué)的改革及各項教育工程的實施，新的教育體制對課表的編排提出了更高的要求。但現(xiàn)實生活中，排課問題屢屢皆是，小學(xué)如此，中學(xué)如此，大學(xué)更是如此，不僅科目多樣，而且教室、老師多變，這使得排課問題往往是很令人費解的。經(jīng)過分析，排課問題就是的多資源組合問題，問題的求解就是找出各個元素間的對應(yīng)關(guān)系。進而將各個元素間的聯(lián)系進一步確定，轉(zhuǎn)化成一個可以量度其大小的值，從而確定優(yōu)先級，而我們又將如何確定各元素間的關(guān)系，目標(biāo)函數(shù)的確定？根據(jù)已有知識可以知道，本題主要分析的是建立一個排課的優(yōu)化模型。而它是一個在課程類別、教師編號、教師及時間上的一個四維空間模型，在各種約束條件下的組合規(guī)劃問題，其實質(zhì)就是解決各因素間的沖突問題。在模型建立后，我們有根據(jù)什么參量得到排課的最優(yōu)解。三、基本假設(shè)模型假設(shè)：1、學(xué)校的教師和教室資源及學(xué)生班結(jié)構(gòu)在一個學(xué)期不會有的變動2、所有的教室都在同一個校區(qū)，且1?2節(jié)課的教室到3?4節(jié)課的教室的路程不超過10min3、在一學(xué)期，任課教師身體都非常健康，不存在因病因事缺課的情況4、各種教學(xué)資源（課桌、多媒體、機房電腦）在一學(xué)期都不會發(fā)生故障，影響上課5、在上課期間，老師、學(xué)生都不遲到，不影響上課質(zhì)量6、當(dāng)有3個課時時，我們當(dāng)做2個課時處理，及3節(jié)連堂上符號說明：Rn教室編號n=1,2,3???.18Cn課程類別n=1,2,3???.40Tn教師編號n=1,2,3-..25Y—外聘外聘教師編號n=1,2,3???..14P(TRC)T老師在R教室上C課X時間段效率X教室利用率X3教師滿意度D優(yōu)化級的量度值W權(quán)衡時間段效率參數(shù)權(quán)衡教室利用率參數(shù)

島權(quán)衡教師滿意度參數(shù)Q安排課程的最優(yōu)值相關(guān)名詞解釋:上課總?cè)藬?shù)教室最大容納量時間段效率：經(jīng)上網(wǎng)查詢及對相關(guān)資料的查閱，我們得知一天聽課效率最高的是上午8?10,下午1?3,故我們定義上午1?2的聽課效率為3,其上課總?cè)藬?shù)教室最大容納量教室利用率：為充分利用教室資源，我們定義：教室利用率二四、問題的分析及模型的建立問題分析(1)從數(shù)學(xué)角度上講，本題主要分析建立一個排課模型，而它是一個在課程類別、教師編號、教師及時間上的一個四維空間模型，在各種約束條件下的組合規(guī)劃問題，其實質(zhì)就是解決各因素間的沖突問題。在此為了簡化處理，先從課程類別、教室編號入手，建立一個關(guān)于C-R的關(guān)系表，再采用化零為整的思路建立我們的目標(biāo)函數(shù)一一優(yōu)化模型，最后，我們根據(jù)各因素對排課模型的優(yōu)先度，求解出排課模型的最優(yōu)解。在對問題初始化分析時，我們發(fā)現(xiàn)課程類別、教室編號、教師、上課時間存在這么一個對關(guān)系：1—1的對應(yīng)關(guān)系)1—n的對應(yīng)關(guān)系)n—n的對應(yīng)關(guān)系進而，我們再對它們之間的屬性分析，根據(jù)它們間的聯(lián)系求出一種相對合理的排課方案，最后，對方案的合理性進行分析。模型的建立經(jīng)過分析，我們需將所有課程盡量合理的安排在一個星期。首先我們將一個星期劃分為五天，記作1、2、3、4、5,將一天分為四個部分，記作1、2、3、4,進而，我們將得到一個5X4的矩陣。其中，j=1、2、3、4、5分別表示星期一、星期二、星期三、星期四、星期五；i二1、2、3、4分別表示1?2節(jié)課、3?4節(jié)課、5?6節(jié)課、7?8節(jié)課。即有：力ISP(TRC)t=1r=1c=1我們記作P(TRC)是一個TXRXC維的數(shù)列矩陣，表示T老師在R教室上C課，我們定義？(TRC)=1時，即老師、教室、課程三者都相互符合是記作1而P(TRC)=0時，即老師、教室、課程三者中有一項不符合記作0規(guī)定：A(TR)表示T老師到R教室上了一次課是，即2個節(jié)B(TR)表示T老師到R教室上總課時約束條件：每一個時間段都不能多于一個老師在一個教室上課，此時應(yīng)滿足的條件是：IP(TRC)<1ngN(N=1、2、3-25)t=1每位老師在每一天不能同時對同一個班上上兩次課，IP(TRC)<1i=1某位老師在某一間教室上課時，安排的課就該在這間教室排完，此時應(yīng)滿足的條件是：x室4P(TRC)=2xA(TR)xB(TR)J=1i=1當(dāng)一位老師連續(xù)兩天對同一個班上兩次課以上的次數(shù)越少、課程安排在聽課效率高的時間段次數(shù)越多和老師與老師之間的沖突(滿意度)次數(shù)越少，此時定義：Q=IIIP(TRC)t=1r=1c=1同時有滿足以上的約束條件，Q將取到最優(yōu)值，即此時安排的課表最優(yōu)化排課的預(yù)處理同一教師在同一時間不能安排兩門課同一教室在同一時間不能安排兩門課同一時間安排的課時總數(shù)不能大于教室的課時總數(shù)所提供的教室屬性安排課程的所需教室屬性一致優(yōu)化級的計算：考慮到課表的安排是為了按時保量的完成教學(xué)任務(wù)，而其影響因素有眾多，這里我們主要考慮一下三種：1)時間段效率2）教室利用率3）教師滿意度所以我們對課表的安排原則將依據(jù)優(yōu)先級的大小進行排課。先對以上三個因素進行量化處理可得到：時間段效率：上午1?2上午3?4下午1?2下午3?43021教室利用率:1?0.80.6?0.80.4?0.60.2?0.43210教師滿意度:很滿意滿意一般不滿意3210我們定義優(yōu)先級的函數(shù)表達式為：D=W1XX1+XXUW3XX3其中W1W2W3表示相應(yīng)參數(shù)，可根據(jù)實際情況進行調(diào)整的通過上式，我們可以很容易得知上課的時間段效率、教室利用率、教師滿意度越大，D值就越大即優(yōu)先級越大，就優(yōu)先排課；反之，D值越小，優(yōu)先級就越小，就相對后排課。通過模型得到的T-R-C表（見附表），再根據(jù)優(yōu)先級的排課，我們可得到一個相對優(yōu)化的課表：R3教室一周課程安排星期一星期二星期三星期四星期五上午1-2C12-T9上午3-4C18-T24C18-T24下午1-2C6-T6下午3-4這里就不一一羅列，其余詳見附表。模型的評價與改進優(yōu)點：到了一種排課的方法，即將課程表按課時數(shù)撤分，再按沒門課程的優(yōu)先順序依次放入課表中。在考慮其他因素（某課程只能安排在固定的某時間段上課，教師只在固定時間段安排課程）的情況下，本模型便于進一步發(fā)展、完善（只需改變最優(yōu)解的參數(shù)加相應(yīng)限制就可解決）。比起目前最相對合理的遺產(chǎn)法簡單易懂很多。缺點：有很多待改進之處，該算法的優(yōu)先級只是根據(jù)經(jīng)驗常識來進行優(yōu)先級安排課程，具有一定的非科學(xué)性，如果能根據(jù)一定的科學(xué)知識，對課程的優(yōu)先級進行科學(xué)的評定，科學(xué)的地算出其優(yōu)先級系數(shù)，再確定其優(yōu)先級，會具有更好的合理性。再解題過程中限制條件多，數(shù)據(jù)量大，過程相對較繁瑣。模型推廣本模型有編程優(yōu)化的模塊，但排課表是人為因素較多，也是為了滿足課程的安排過程中更人性化。因此本模型適用于中小規(guī)模的排課，其主要限制因素是程序完善度不高，如要進行推廣，則必須增加程序的模塊，使其更完善，相信改進完善后，即可更科學(xué)的實現(xiàn)大規(guī)模大排參考文獻數(shù)學(xué)建模（上冊），電子機械高等?？茖W(xué)校，二零一零年三月數(shù)學(xué)建模資料——最優(yōu)化模型.pdf2010.7回溯法baike.baidu./view/45.htm?fr=ala012010.7貪婪法——baike.baidu./view/112297.html?tp=0102010.7一天的最有聽課效率時間段zhidao.baidu.2010.7附表：表1:教師屬性:教師編號能勝任課程類別周最大課時數(shù)對教室類別要求上課時間要求盡可能不同天上課的教師T011，84多媒體教室或機房上午T4T0214普通教室下午T0316普通教室上午T0424多媒體教室上午T1T0524普通教室下午T0626普通教室無T0734普通教室上午T083，83普通教室或機房下午T0934普通教室上午T11,T18T1036多媒體教室上午T1148普通教室下午T9,T1244普通教室無T1346多媒體教室下午T1452普通教室上午T155，83普通教室或機房下午T23T1654普通教室上午T1764普通教室下午T1866普通教室無T9T1964多媒體教室下午T2074普通教室上午T2176普通教室下午T2276多媒體教室上午T233，84普通教室或機房無T15T244,86普通教室或機房上午T256，84普通教室或機房下午

表2：課程屬性及要求:課程編號課程類別周課時數(shù)對教室座位最大要求數(shù)對教室類別要求時間要求C011450多媒體教室上午C021430普通教室下午C031640普通教室下午C041425多媒體教室上午C051360普通教室下午C0624100普通教室下午C072450多媒體教室上午C082230普通教室上午C092440普通教室下午C102325多媒體教室上午C113660普通教室上午C123480普通教室上午C133650多媒體教室下午C143230普通教室下午C153340普通教室下午C164425多媒體教室上午C174460普通教室下午C184690普通教室上午C194450多媒體教室上午C204250普通教室上午C215430普通教室上午C225440多媒體教室下午C235325普通教室上午C245660普通教室下午C2554100多媒體教室上午C266350多媒體教室下午C276430普通教室下午C286440普通教室下午C296450多媒體教室上午C306430普通教室下午C317640普通教室上午C327425多媒體教室下午C337360普通教室無C347480多媒體教室上午C357450普通教室下午C368630機房無C378440機房上午C388225機房上午C398460機房下午C408450機房無表3：教室屬性:教室編號最大座位數(shù)教室類別R01100多媒體教室R02100普通教室R03100普通教室R0450多媒體教室R0550普通教室R0650普通教室R0750普通教室R0860普通教室R0960普通教室R1060多媒體教室R1160普通教室R1260機房R1340機房R1440機房R1540機房R1640多媒體教室R1750普通教室R1840普通教室表4：T-R-C表課程老師教室容量C1T1R(10)50/60C2T2R(18)30/40C3T3R(18)40/40C4C5C6T6R(2)R(3)100/100C7T4R(4)50/50C8T6R(18)30/40C9T5R(18)40/40C10C11C12T7R(2)R(3)80/100C13T10R(4)R(16)50/50C14T8R(18)40/40C15T9R(18)40/40C16T13R(16)30/40C17T11R(8)R(9)R(11)60/60C18T24R(2)R(3)90/100C19T11R(5)R(6)R(7)R(17)50/50C20T12R(5)R(6)R(7)R(17)50/50C21T16R(18)30/40C22

C23C24C25C26T19R(4)50/50C27T18R(18)40/40C28T17R(18)40/40C29C30C31C32C33T20R(8)R(9)R(11)60/60C34T22R(1)80/100C35T21R(5)R(6)R(7)R(17)50/50C36C37C38T15R(13)R(14)R(15)25/40C39T25R(12)60/60C40T23R(12)60/60表5：各教室課程安排:R1教室一周課程安排星期一星期二星期三星期四星期五上午1-2C34-T22C34-T22C25上午3-4C25下午1-2下午3-4R2教室一周課程安排星期一星期二星期三星期四星期五上午1-2C18-T24上午3-4C12-T9下午1-2C6-T6下午3-4R3教室一周課程安排星期一星期二星期三星期四星期五上午1-2C12-T9上午3-4C18-T24C18-T24下午1-2C6-T6下午3-4R4教室一周課程安排

星期一星期二星期三星期四星期五上午1-2C13-T10C19-T11C13-T10上午3-4C19-T11下午1-2C26-T19下午3-4R5教室一周課程安排星期一星期二星期三星期四星期五上午1-2C21-T16C35-T21上午3-4C35-T21C21-T16下午1-2C30C17-T12C30下午3-4R6教室一周課程安排星期一星期二星期三星期四星期五上午1-2C31-T21上午3-4下午1-2C28-T17C9-T5C9-T5下午3-4C28-T17R7教室一周課程安排星期一星期二星期三星期四星期五上午1-2C8-T6上午3-4下午1-2C20-T12C27-T17下午3-4C27-T17C3-T3R8教室一周課程安排星期一星期二星期三星期四星期五上午1-2C31上午3-4C31下午1-2C5C33-T20下午3-4R9教室一周課程安排星期一星期二星期三星期四星期五上午1-2C11C11C11上午3-4下午1-2C24下午3-4C24C24R10教室一周課程安排星期一星期二星期三星期四星期五

上午1-2C1-T1C29C29C7-T4上午3-4C7-T4C1-T1下午1-2下午3-4R11教室一周課程安排星期一星期二星期三星期四星期五上午1-2上午3-4下午1-2C17-T11下午3-4C17-T11R12教室一周課程安排星期一星期二星期三星期四星期五上午1-2上午3-4下午1-2C40-T23C39-T25下午3-4C39-T25C40-T23R13教室一周課程安排星期一星期二星期三星期四星期五上午1-2上午3-4下午1-2C28-T18下午3-4C36C28-T18R14教室一周課程安排星期一星期二星期三星期四星期五上午1-2C37上午3-4C37下午1-2C36C36下午3-4R15教室一周課程安排星期一星期二星期三星期四星期五上午1-2C38-T24上午3-4下午1-2下午3-4R16教室一周課程安排星期一星期二星期三星期四星期五

上午1-2C10C16-T13C4上午3-4C4C16-T13下午1-2C22下午3-4C32C22C32R17教室一周課程安排星期一星期二星期三星期四星期五上午1-2上午3-4下午1-2C15-T9C14-T08C2-T2下午3-4R18教室一周課程安排星期一星期二星期三星期四星期五上午1-2C23上午3-4C8-T6下午1-2C3-T3下午3-4C2-T2C3-T3C14-T8表6：外聘教師課程安排表:外聘教師課程課時教室時間YT1C254R1二34，五12YT2C304R5二56，五56YT3C314R8一12，三34YT4C53R8二123YT5C116R9一12，三12，五12YT6C244R9一56，三78，五78YT7C294R10二12，四12YT8C366R14一56，三78，五56YT9C103R16一123YT10C44R16二34，五12YT11C224R16一56，四78YT12C324R16二78，五78YT13C233R18三123YT14C374R14二12，四34clear；C=[1221121213121211123222422221211222122211333232423321311232124211423243424221412252125211521253325241622162126212程序附錄:-1-11-1-111-1111-1-1-11-11111-1-11-1-1-1

TOC\o"1-5"\h\z1-11-101-1011-101402-1002100110-10020002100-100211118100-1902000-1001005-123021002-1002090-1001002-1001100501505100622621731721723724722831821811823822];T=[184104106204204206304383304306408404406502583504604606604704706706384486684];R=[414242212222225-100KRc-t-*#網(wǎng)>出身建>部Rc-t-nzeros(25“18)'riHr一8ormHr25for3Hr40ifR(r-1)〉HC(F3)ilHCp1)'ifTGIhlYHil一一T(nh2)Ril*T(m-4)\switch(i4)case1ifRF2YH1Rr+(m“iYrendcase2ifR(r-2)H2Rr+(m“iYrendcase3ifR(r-2)R3Rr+(m3)Hl'endcase4ifR(r-2)H1--R(r-2)R3Rr+(m3)Hl'endcase5ifR(r-2)R2--R(r-2)H3Rr+(m“iYrendendendendendendendRt；%Rc為教室與課程的關(guān)系Rc二zeros(18,40)；fori=1:18forn=1:40ifR(i,1)>=C(n,3)&&R(i,2)==C(n,4)Rc(i,n)=1;endendendRc；%Tcform=1:25forn=1:40ifT(m,1)=C(n,1)||T(m,2)=C(n,1)ifT(m,3)>=2*C(n,2)q=T(m,4);t=T(m,5);switchqcase1ifC(n,4)==1switchtcase1ifC(n,5)==1Tc(m,n)=1;endcase-1ifC(n,5)==-1Tc(m,n)=1;endcase0ifC(n,5)==1||C(n,5)=T||C(n,5)=0Tc(m,n)=1;endendendcase2ifC(n,4)==2switchtcase1ifC(n,5)==1Tc(m,n)=1；endcase-1ifC(n,5)==-1Tc(m,n)=1;endcase0ifC(n,5)==1||C(n,5)=T||C(n,5)=0Tc(m,n)=1;endendendcase3ifC(n,4)==3switchtcase1ifC(n,5)==1Tc(m,n)=1;endcase-1ifC(n,5)==-1Tc(m,n)=1;endcase0ifC(n,5)==1||C(n,5)=T||C(n,5)=0Tc(m,n)=1;endendendcase4ifC(n,4)==1||C(n,4)==3switchtcase1ifC(n,5)==1Tc(m,n)=1;endcase-1ifC(n,5)==-1Tc(m,n)=1;endcase0ifC(n,5)==1||C(n,5)=T||C(n,5)=0endendendcase5ifC(n,4)==2||C(n,4)==3switchtcase1ifC(n,5)==1Tc(m,n)=1;endcase-1ifC(n,5)==-1Tc(m,n)=1;endcase0ifC(n,5)==1||C(n,5)=T||C(n,5)=0Tc(m,n)=1;endendendendendendendendTc;%在完全滿足的情況下課程選擇最優(yōu)化的老師%Ct說明行代表課程，列代表能上該課的老師編號%if有相同編號的老師那針對該門課程進行求最優(yōu)值，將值最大的老師放入相應(yīng)課程中Ct=[100;00;000;TOC\o"1-5"\h\z00;00;60;00;00;60;400;000;923;000;230;00;000;TOC\o"1-5"\h\z120；00；000；240；00；000；00；000；000；00；1825；171825；000；171825；000；000；000；00；00；000；12324；12324；250；2300；］；%具體實現(xiàn)Ct1二zeros(40,1)；%用于存放最終值%先判斷只有一個老師上課的課程forn=1:40fornn=n+1:40iflength(Ct(Ct(n,:)~=0))==1ifnn~=nforu=1:3ifCt(n,1)==Ct(nn,u)%判斷老師編號相同u1二Ct(n,1)；%獲取老師編號a1=n；%獲取課程號a2二nn；%獲取課程號b1二T(u1,3)-2*C(a1,2)；%計算老師課時與課程課時的差b2二T(u1,3)-2*C(a2,2)；%計算老師課時與課程課時的差ifb1<=b2ifCt1(a1)==0Ct1(a1)二u1；endelseifCt(a2)==0Ct(a2)=u1；endendifCt(n,2)==Ct(nn,u)%判斷老師編號相同u1二Ct(n,1)；%獲取老師編號a1=n；%獲取課程號a2=nn;%獲取課程號b1二T(u1,3)-2*C(a1,2);%計算老師課時與課程課時的差b2二T(u1,3)-2*C(a2,2);%計算老師課時與課程課時的差ifb1<=b2ifCt1(a1)==0Ct1(a1)二u1;endelseifCt(a2)==0Ct(a2)=u1;endendendifCt(n,3)==Ct(nn,u)%判斷老師編號相同u1二Ct(n,1);%獲取老師編號a1=n;%獲取課程號a2=nn;%獲取課程號b1二T(u1,3)-2*C(a1,2);%計算老師課時與課程課時的差b2二T(u1,3)-2*C(a2,2);%計算老師課時與課程課時的差ifb1<=b2ifCt1(a1)==0Ct1(a1)二u1;endelseifCt(a2)==0Ct(a2)=u1;endendendendendiflength(Ct(Ct(n,:廣二0)廣二1ifnn~=nforu=1:3ifCt(n,1)==Ct(nn,u)%判斷老師編號相同u1二Ct(n,1);%獲取老師編號a1=n;%獲取課程號a2=nn;%獲取課程號b1二T(u1,3)-2*C(a1,2);%計算老師課時與課程課時的差b2二T(u1,3)-2*C(a2,2);%計算老師課時與課程課時的差ifb1<=b2ifCt1(a1)==0endelseifCt(a2)==0Ct(a2)=u1；endendendifCt(n,2)==Ct(nn,u)%判斷老師編號相同u1二Ct(n,1)；%獲取老師編號a1=n;%獲取課程號a2=nn;%獲取課程號b1二T(u1,3)-2*C(a1,2);%計算老師課時與課程課時的差b2二T(u1,3)-2*C(a2,2);%計算老師課時與課程課時的差ifb1<=b2ifCt1(a1)==0Ct1(a1)二u1;endelseifCt(a2)==0Ct(a2)=u1;endendendifCt(n,3)==Ct(nn,u)%判斷老師編號相同u1二Ct(n,1);%獲取老師編號a1=n;%獲取課程號a2=nn;%獲取課程號b1二T(u1,3)-2*C(a1,2);%計算老師課時與課程課時的差b2二T(u1,3)-2*C(a2,2);%計算老師課時與課程課時的差ifb1<=b2ifCt1(a1)==0Ct1(a1)二u1;endelseifCt(a2)==0Ct(a2)=u1;endendendendendendendendend%T-R-C優(yōu)化［老師教室］TRC=［110；%1-3218；TOC\o"1-5"\h\z18；2；%6-94；18；18；2；%12-214；18；18；16；8；2；5；5；18；4；%26-2818；18；8；%33-351；5；1513；%38-4012；12；］;%定義時間片%B［周幾上課時間段老師一周幾節(jié)］B二zeros(26,4)；%課程1%課程1B二zeros(10,4,20)；n=1；fora=1:5a1二a；forb=1:2b1二b；forn=n+1:20B(1,:,n)=［a1b112］；breakendendend%課程2n=1；fora=1：5al二a；forb=3：4b1=b；forn=n+1：20B(2,：,n)=[alb122]；breakendendend%課程3n=1;fora=1:5a1二a;forb=3:4b1二b;forn=n+1:20B(3,:,n)=[a1b133];breakendendend%課程6n=1;fora=1:5a1二a;forb=3:4b1二b;forn=n+1:20B(4,:,n)=[a1b162];breakendendend%課程7n=1;fora=1:5a1二a;forb=1:2b1=b；forn=n+1：20B(5,：,n)=[albl42]；breakendendend%課程8n=1;fora=1:5a1二a;forb=1:2b1二b;forn=n+1:20B(6,:,n)=[a1b161];breakendendend%課程9n=1;fora=1:5a1二a;forb=3:4b1二b;forn=n+1:20B(7,:,n)=[a1b152];breakendendend%課程12n=1;fora=1:5a1二a;forb=1:2b1二b;forn=n+1:20B(8,:,n)=[a1b172];breakendendendn=1；fora=1：5al二a；forb=3：4b1=b；forn=n+1：20B(9,：,n)=[alb1103]；breakendendend%課程14n=1;fora=1:5a1二a;forb=1:2b1二b;forn=n+1:20B(10,:,n)=[a1b181];breakendendend%課程15n=1;fora=1:5a1二a;forb=3:4b1二b;forn=n+1:20B(11,:,n)=[a1b192];breakendendend%課程16n=1;fora=1:5a1二a;forb=1:2b1二b;forn=n+1:20B(12,:,n)=[a1b1132];breakendend%課程17n=1；fora=1：5al二a；forb=3：4b1=b；forn=n+1：20B(13,：,n)=[a1b1112];breakendendend%課程18n=1;fora=1:5a1二a;forb=1:2b1二b;forn=n+1:20B