數(shù)學(xué)建模論文醫(yī)院眼科病床的合理安排優(yōu)化模型

1、 醫(yī)院眼科病床的合理安排優(yōu)化模型摘要：眼科病床的合理安排，既可以提高醫(yī)院的經(jīng)濟(jì)效益，又可以減小患者的排隊(duì)就醫(yī)時(shí)間，以便達(dá)到雙贏的目的。通過(guò)對(duì)醫(yī)院現(xiàn)有的綜合分析及優(yōu)化得到如下問(wèn)題的優(yōu)化結(jié)果：?jiǎn)栴}一，給出了評(píng)價(jià)模型的指標(biāo)：醫(yī)院滿意度即床位的周轉(zhuǎn)率，患者滿意度即等待隊(duì)長(zhǎng)，及逗留時(shí)間；問(wèn)題二，通過(guò)對(duì)已有數(shù)據(jù)的統(tǒng)計(jì)分析整理，得出結(jié)論：患者的到達(dá)服從泊松分布，服務(wù)時(shí)間服從負(fù)指數(shù)分布。將醫(yī)院的病床看做服務(wù)臺(tái)，建立了排隊(duì)系統(tǒng)的模型，通過(guò)對(duì)模型的求解過(guò)程，給出了患者流和服務(wù)臺(tái)的計(jì)算機(jī)模擬方法，進(jìn)而給出了根據(jù)出院患者數(shù)來(lái)安排入住患者的方案。問(wèn)題三，結(jié)合問(wèn)題二的模型，利用醫(yī)院對(duì)每類患者的服務(wù)時(shí)間，服務(wù)效率以及服務(wù)臺(tái)

2、的狀態(tài)（病床的占用和空閑）來(lái)求的他們的平均等待時(shí)間，因而可預(yù)測(cè)病人的大致入院時(shí)間。問(wèn)題四，若周六，周日不安排手術(shù)，而白內(nèi)障手術(shù)的手術(shù)時(shí)間間隔一定，通過(guò)順延的方法得出兩種優(yōu)化方案，并對(duì)兩種方案通過(guò)類比的方法得出最優(yōu)方案，即在周三，周五安排白內(nèi)障手術(shù)，而在周一，周二，周四安排其它手術(shù)。問(wèn)題五，從方便管理，醫(yī)院效益最大以及患者逗留時(shí)間最短的前提出發(fā)，對(duì)醫(yī)院現(xiàn)有的病床進(jìn)行分配建立了規(guī)劃模型，并通過(guò)對(duì)模型的求解最后給出了病床的分配方案：白內(nèi)障單眼8張床位，白內(nèi)障雙眼18張床位，青光眼11張床位，視網(wǎng)膜疾病33張床位，外傷9張床位。關(guān)鍵詞： 排隊(duì)系統(tǒng) 統(tǒng)計(jì)分析 計(jì)算機(jī)模擬 類比 規(guī)劃 一 問(wèn)題重述 醫(yī)院就

3、醫(yī)排隊(duì)是大家都非常熟悉的現(xiàn)象，我們考慮某醫(yī)院眼科病床合理安排的數(shù)學(xué)建模問(wèn)題，就要考慮醫(yī)院資源的合理有效利用及患者的等待隊(duì)長(zhǎng)等因素。已知該醫(yī)院目前情況如下：1 該醫(yī)院眼科門(mén)診每天開(kāi)放，住院部共有病床79張，眼科手術(shù)主要分四大類：白內(nèi)障、視網(wǎng)膜疾病、青光眼和外傷。附錄中給出了2008年7月13日至2008年9月11日這段時(shí)間里各類病人的情況。2白內(nèi)障手術(shù)較簡(jiǎn)單，而且沒(méi)有急癥，目前該院是每周一、三做白內(nèi)障手術(shù)，此類病人的術(shù)前準(zhǔn)備時(shí)間只需1、2天。做兩只眼的病人比做一只眼的要多一些，大約占60%。如果要做雙眼是周一先做一只，周三再做另一只。3外傷疾病通常屬于急癥，病床有空時(shí)立即安排住院，住院后第二天便

4、會(huì)安排手術(shù)。4其他眼科疾病比較復(fù)雜，有各種不同情況，但大致住院后2-3天內(nèi)就可以接受手術(shù)，主要是術(shù)后的觀察時(shí)間較長(zhǎng)。這類疾病手術(shù)時(shí)間可根據(jù)需要安排，一般不安排在周一、周三。5該醫(yī)院眼科手術(shù)條件比較充分，通常情況下白內(nèi)障手術(shù)與其他眼科手術(shù)（急癥除外）不安排在同一天做。當(dāng)前該住院部對(duì)全體非急癥病人是按照f(shuō)cfs（first come, first serve）規(guī)則安排住院，但等待住院病人隊(duì)列卻越來(lái)越長(zhǎng)，現(xiàn)提出以下問(wèn)題：?jiǎn)栴}一：試分析確定合理的評(píng)價(jià)指標(biāo)體系，用以評(píng)價(jià)該問(wèn)題的病床安排模型的優(yōu)劣。問(wèn)題二：試就該住院部當(dāng)前的情況，建立合理的病床安排模型，以根據(jù)已知的第二天擬出院病人數(shù)來(lái)確定第二天應(yīng)該安排哪

5、些病人住院。并對(duì)你們的模型利用問(wèn)題一中的指標(biāo)體系作出評(píng)價(jià)。問(wèn)題三：作為病人，自然希望盡早知道自己大約何時(shí)能住院。能否根據(jù)當(dāng)時(shí)住院病人及等待住院病人的統(tǒng)計(jì)情況，在病人門(mén)診時(shí)即告知其大致入住時(shí)間區(qū)間。 問(wèn)題四：若該住院部周六、周日不安排手術(shù)，請(qǐng)你們重新回答問(wèn)題二，醫(yī)院的手術(shù)時(shí)間安排是否應(yīng)作出相應(yīng)調(diào)整?問(wèn)題五：有人從便于管理的角度提出建議，在一般情形下，醫(yī)院病床安排可采取使各類病人占用病床的比例大致固定的方案，試就此方案，建立使得所有病人在系統(tǒng)內(nèi)的平均逗留時(shí)間（含等待入院及住院時(shí)間）最短的病床比例分配模型。 二符號(hào)說(shuō)明 表示第類患者（=1，2，3，4，5 分別表示白內(nèi)障單眼；白內(nèi)障雙眼；青光眼；視網(wǎng)

6、膜疾病以及外傷患者） 表示服務(wù)臺(tái)個(gè)數(shù) 表示第類患者中第個(gè)患者與第個(gè)患者到達(dá)的間隔時(shí)間 表示對(duì)第類患者中第個(gè)患者服務(wù)花費(fèi)的時(shí)間 表示患者到來(lái)的時(shí)刻 表示患者等待時(shí)間 表示患者開(kāi)始接受服務(wù)（入院）時(shí)刻 表示患者結(jié)束服務(wù)（出院）時(shí)刻 表示各類患者的到達(dá)速率 表示服務(wù)臺(tái)的服務(wù)速率 表示服務(wù)強(qiáng)度 表示時(shí)間內(nèi)到來(lái)患者的總數(shù) 表示時(shí)間內(nèi)服務(wù)完（出院）的患者的總數(shù) 表示時(shí)刻系統(tǒng)中的患者數(shù) 表示等待隊(duì)長(zhǎng)（排隊(duì)等待的顧客數(shù)） 表示隊(duì)長(zhǎng) 表示服務(wù)臺(tái)的狀態(tài)三問(wèn)題的分析1背景分析本問(wèn)題的難點(diǎn)是同時(shí)考慮醫(yī)院完善醫(yī)院的管理制度，改進(jìn)患者的安置情況，提高醫(yī)院的經(jīng)濟(jì)和社會(huì)效益等諸多因素。如果僅考慮提高醫(yī)院的經(jīng)濟(jì)效益，則只要提高

7、醫(yī)院病床的使用率，運(yùn)用數(shù)據(jù)分析方法可以給出患者的安排情況；如果僅考慮患者的安置情況，需要增加病床，顯然，這兩種方案都不是最佳的。于是既要考慮醫(yī)院的效益（即醫(yī)院的滿意度），又要考慮患者的等待時(shí)間最短（即患者的滿意度）。 醫(yī)院的滿意度取決于病床的使用率及其服務(wù)強(qiáng)度，使用率和服務(wù)強(qiáng)度越高，醫(yī)院的經(jīng)濟(jì)效益越好，醫(yī)院的滿意度越高；患者的滿意度取決于排隊(duì)的隊(duì)長(zhǎng)和等待時(shí)間的長(zhǎng)短，等待時(shí)間越短患者的滿意度越高。所以我們需要在兩個(gè)因素之間找出一個(gè)合理的匹配關(guān)系使雙方的滿意度達(dá)到最高。2數(shù)據(jù)分析通過(guò)數(shù)據(jù)分析，形成了如圖所示的柱狀圖及其以下數(shù)據(jù)表：通過(guò)柱狀圖可以比較清楚的看出，每天各類疾病的患者到達(dá)數(shù)量及所占比重：

8、視網(wǎng)膜疾病患者最多，相當(dāng)于青光眼和外傷的總和，其次是白內(nèi)障（雙眼）和白內(nèi)障。從上圖中還可以看出患者人數(shù)基本以每?jī)芍転橐粋€(gè)循環(huán)周期。表 1：病種白內(nèi)障白內(nèi)障（雙眼）青光眼視網(wǎng)膜疾病外傷人數(shù) 1001336616764平均患者到達(dá)率2 2 1 3 1 表 2：病種白內(nèi)障白內(nèi)障（雙眼）青光眼視網(wǎng)膜疾病外傷病種平均服務(wù)時(shí)間59101379 平均等待時(shí)平均逗留時(shí)利用spss軟件驗(yàn)證門(mén)診到達(dá)人數(shù)服從泊松分布。observed nexpected nresidual1.0012.9-1.92.0023.8-1.83.0031.81.24.0044.8-.8

9、5.0051.83.2total15test statisticsxchi-square(a)8.640df4asymp. sig.071a 5 cells (100.0%) have expected frequencies less than 5. the minimum expected cell frequency is 1.8. 通過(guò)對(duì)數(shù)據(jù)的分析整理，利用spss得到上面的兩表。由于顯著性水平為0.05，而上表中所得顯著性水平0.0710.05且到達(dá)門(mén)診的患者是隨機(jī)的，所以到達(dá)門(mén)診的患者人數(shù)服從泊松分布，認(rèn)為理論分布是符合實(shí)際情況的。四模型的假設(shè)1. 患者的到達(dá)是隨機(jī)的，且患者源是無(wú)

10、限的。2. 忽略門(mén)診時(shí)間及手術(shù)時(shí)間。3. 將病床看為服務(wù)臺(tái)，總病床數(shù)即為服務(wù)臺(tái)數(shù)。4. 同一天看病的患者到達(dá)間隔時(shí)間為零。五模型的建立與求解問(wèn)題一 評(píng)價(jià)性指標(biāo)大致可以從患者的滿意度和醫(yī)院的滿意度兩個(gè)方面出發(fā)：患者滿意度方面：要求醫(yī)院充分調(diào)動(dòng)和利用現(xiàn)有的一切資源，使得到醫(yī)院的每一位患者滿意而歸，這就要求醫(yī)院的資源合理分配，使患者在醫(yī)院停留的時(shí)間最短，因?yàn)榉?wù)時(shí)間一定，故只要求等待時(shí)間最短。醫(yī)院滿意度方面：達(dá)到醫(yī)院現(xiàn)有資源的充分利用，即效益最大，以服務(wù)更多的患者為宗旨，即要求患者的等待隊(duì)長(zhǎng)盡可能小。資源的利用率和服務(wù)強(qiáng)度盡可能高。 因此，該模型的評(píng)價(jià)指標(biāo)體系是等待隊(duì)長(zhǎng)，等待時(shí)間，病床使用率，服務(wù)強(qiáng)

11、度。問(wèn)題二1模型分析與建立從所要解決的問(wèn)題和對(duì)問(wèn)題所做的假設(shè)出發(fā)，在前面的數(shù)據(jù)處理中我們確定了患者的到達(dá)率服從泊松分布，有個(gè)服務(wù)臺(tái)，患者源的數(shù)量是無(wú)限的。因?yàn)榛颊叩竭_(dá)的速率為，每個(gè)服務(wù)臺(tái)的服務(wù)速率為，則整個(gè)系統(tǒng)的最大速率為服務(wù)強(qiáng)度，顯然，要使系統(tǒng)穩(wěn)定運(yùn)行，必須有。該模型的特點(diǎn)在于整個(gè)系統(tǒng)的服務(wù)速率與系統(tǒng)中患者到達(dá)特點(diǎn)有關(guān)。如果系統(tǒng)中只有一個(gè)患者，則系統(tǒng)的服務(wù)率等于，因?yàn)槠渌姆?wù)臺(tái)處于閑置狀態(tài)，如果體統(tǒng)中有2個(gè)患者，則系統(tǒng)的服務(wù)速率等于，依次類推，如果系統(tǒng)中的患者有個(gè)，則系統(tǒng)的最大服務(wù)速率為，所有的服務(wù)臺(tái)均投入服務(wù)；而當(dāng)系統(tǒng)中的患者數(shù)目超過(guò)個(gè)時(shí)，多余的患者只能進(jìn)入派對(duì)系統(tǒng)等待服務(wù)，此時(shí)系統(tǒng)的服

12、務(wù)速率為。該模型的狀態(tài)轉(zhuǎn)移圖如下： （排隊(duì)服務(wù)系統(tǒng)狀態(tài)轉(zhuǎn)移圖）通過(guò)狀態(tài)轉(zhuǎn)移圖可得如下方程： （1） （2） （3）而 聯(lián)立可解得： （4） （5） （6） （7）設(shè)排隊(duì)隊(duì)長(zhǎng)為，當(dāng)系統(tǒng)有個(gè)患者到達(dá)時(shí)，若，則必有個(gè)患者在排隊(duì)。這樣的話，能求出每種病的等待時(shí)長(zhǎng)，服務(wù)時(shí)長(zhǎng)及停留的時(shí)間。而當(dāng)時(shí)，服務(wù)臺(tái)空閑個(gè)數(shù)為，既有個(gè)空病床，然后可以安排其他的患者入院。 從而建立如下的模型： （8） （9） （10） （11） （12）2模型的模擬求解2.1 病人流的模擬由于已認(rèn)定病人達(dá)到流服從泊松分布，在模擬的時(shí)候可以通過(guò)隨機(jī)數(shù)生成函數(shù)來(lái)產(chǎn)生模擬病人的到達(dá)，基本算法如下：步驟一：在泊松分布中，求出 x取何值時(shí)，p(x

13、=k)取最大值，設(shè)最大值為，這相當(dāng)于求解在k取何值時(shí)有最大值，可以通過(guò)一個(gè)循環(huán)來(lái)得到取最大值時(shí)的整數(shù)自變量。步驟二：通過(guò)迭代，不斷生成0-1區(qū)間上的隨機(jī)數(shù)。當(dāng)隨機(jī)數(shù)小于時(shí)，則終止迭代，否則重復(fù)步驟二。步驟三：記錄迭代過(guò)程的次數(shù)，即為所需要得到的符合泊松分布的隨機(jī)量。2.2 服務(wù)臺(tái)的模擬服務(wù)臺(tái)可由隨機(jī)數(shù)生成函數(shù)來(lái)模擬，其算法為：令滿足均勻分布的隨機(jī)變量 ，y的概率是，因?yàn)?x是服從均勻分布的，所以它的概率分布函數(shù)是當(dāng)x的范圍在0-1之間的時(shí)候，概率分布函數(shù)即是x，因此上式的結(jié)果是，即為負(fù)指數(shù)分布。由此滿足負(fù)指數(shù)分布的隨機(jī)變量可以用得出。（具體的計(jì)算機(jī)編程見(jiàn)附錄1）方案的給出通過(guò)數(shù)據(jù)的統(tǒng)計(jì)和計(jì)算機(jī)

14、的模擬可以給出病人出院后床位的空閑安置情況按照下面的原則：服務(wù)時(shí)間最短的優(yōu)先考慮，即先安排白內(nèi)障單眼患者入住，其次安排服務(wù)時(shí)間短（即住院時(shí)間短）的病人入住。該模型對(duì)于醫(yī)院而言，極大地提高了醫(yī)院病床的周轉(zhuǎn)率，使得醫(yī)院的效益達(dá)到最大。而對(duì)于患而言，醫(yī)院的病床周轉(zhuǎn)率提高了，服務(wù)的人多了，等待的隊(duì)列就減短了，能滿足大部分的患者滿意，當(dāng)然這樣就犧牲了部分需要服務(wù)時(shí)間長(zhǎng)的患者的利益。問(wèn)題三根據(jù)題目所給的數(shù)據(jù)整理和統(tǒng)計(jì)，結(jié)合問(wèn)題二的模型得出如下表表示的結(jié)果。平均等待時(shí)間在床等待時(shí)間平均服務(wù)時(shí)間白內(nèi)障11.240513.9493675.236111111白內(nèi)障(雙眼)11.384622.7788468.560

15、97561青光眼11.354172.08333310.48717949視網(wǎng)膜疾病9.9849621.48120312.54455446外傷117.036363636由上表可以大致確定每種患者入院前的平均等待時(shí)間，具體如下表：病類白內(nèi)障白內(nèi)障（雙眼）青光眼視網(wǎng)膜疾病外傷平均等待時(shí)間11.2411.3811.359.981根據(jù)每種病在入院前的平均等待時(shí)間，只要確定了是那類患者，便可以估計(jì)其入院時(shí)間。白內(nèi)障，青光眼患者大概在門(mén)診后11天入院，視網(wǎng)膜疾病患者大概在門(mén)診后10天入院，而外傷患者在門(mén)診后的第二天就可以入院。問(wèn)題四周六、周日不安排手術(shù)，故按照原有的方案，白內(nèi)障病人手術(shù)在周一、周三進(jìn)行，而其它

16、的手術(shù)只能在周二、周四、周五進(jìn)行，這樣，以“1”表示白內(nèi)障手術(shù)時(shí)間，“0”表示其它手術(shù)時(shí)間，“”表示不進(jìn)行手術(shù)，得下表：一二三四五六七10100從該表可以看出若在周六、周日不進(jìn)行手術(shù)，則其它手術(shù)只能在周二、周四、周五進(jìn)行，這樣，本在周六、周日進(jìn)行的手術(shù)只能在周二、周四、周五進(jìn)行，這樣，使得患者等待期間增加病床利用率較低。而從題目中給出的數(shù)據(jù)可以統(tǒng)計(jì)出周一到周五的各類患者的人數(shù)情況如下表：（具體數(shù)據(jù)見(jiàn)附錄2）星期一星期二星期三星期四星期五星期六星期日白內(nèi)障（單）1621221181012白內(nèi)障（雙）23122124251315青光眼疾病99515889視網(wǎng)膜疾病25193131221824外傷1

17、06121012410現(xiàn)考慮白內(nèi)障手術(shù)安排相隔時(shí)間一定，采取將白內(nèi)障手術(shù)時(shí)間順延的方法，得到如下兩種方案，如圖所示：方案1一二三四五六七01010方案2一二三四五六七00101對(duì)于方案一，可以看出白內(nèi)障病人在周二、周四進(jìn)行，其它手術(shù)安排在周一、周三、周五，這樣，周六、周日累積下來(lái)的病人要求在周一內(nèi)做完，這樣對(duì)醫(yī)院、醫(yī)生的負(fù)荷太大，否則要等到周三、周五來(lái)做，使患者的等待時(shí)間延長(zhǎng)。而對(duì)于方案二，白內(nèi)障病人安排在周三、周五，其它病人安排在一、三、四三天進(jìn)行。這樣，周六、周日的病人可在周一、周二內(nèi)做，較方案一而言，方案二的患者的等待時(shí)間較短，故醫(yī)院應(yīng)采用方案二。問(wèn)題五1分析：通過(guò)數(shù)據(jù)的整理統(tǒng)計(jì)，病床在

18、7月13號(hào)到9月11號(hào)安排的每類患者人數(shù)的如下統(tǒng)計(jì)表：病類白內(nèi)障（單）白內(nèi)障（雙）青光眼視網(wǎng)膜疾病外傷服務(wù)人數(shù)791035113163服務(wù)一個(gè)患者的時(shí)間5.238.5610.4912.547.04服務(wù)總時(shí)間413.17881.68534.991642.72443.52根據(jù)上表做直方圖如下（圖中橫坐標(biāo)1，2，3，4，5分別代表上表中5種?。河芍狈綀D可以直觀的看出每種病在醫(yī)院接受服務(wù)的人數(shù)，服務(wù)的平均時(shí)間和服務(wù)的總時(shí)間的分布情況。根據(jù)每類病人占用時(shí)間的比例，以及題目中所給的醫(yī)院的目前現(xiàn)有的床位數(shù)，建立規(guī)劃模型。設(shè)白內(nèi)障（單眼），白內(nèi)障（雙眼），青光眼，視網(wǎng)膜疾病及外傷分配的床位數(shù)分別為。床位的服

19、務(wù)效率分別為。2.模型的建立建立模型如下： 且為整數(shù)模型的求解：利lingo軟件求解得：，故從便于管理的角度出發(fā)，將病床采取病人占用病床的比例可將醫(yī)院現(xiàn)有病床按照各類病人的情況分為：白內(nèi)障單眼8張，白內(nèi)障雙眼18張，青光眼11張，視網(wǎng)膜疾病33張，外傷9張。六 模型的改進(jìn)與評(píng)價(jià) 對(duì)醫(yī)院的排隊(duì)系統(tǒng)而言，要使醫(yī)院在經(jīng)濟(jì)上獲得最大效益，可以等價(jià)的轉(zhuǎn)化為用時(shí)間費(fèi)用來(lái)衡量經(jīng)濟(jì)效益，將等待時(shí)間轉(zhuǎn)化為等待費(fèi)用。根據(jù)統(tǒng)計(jì)的資料估計(jì)，由于隊(duì)列過(guò)長(zhǎng)導(dǎo)致等待費(fèi)用增加。因而該問(wèn)題費(fèi)用函數(shù)期望取最小的問(wèn)題屬于非線性規(guī)劃問(wèn)題。對(duì)于眼科這種類型的離散型隨機(jī)變量，也可以用邊際分析法進(jìn)行求解。該模型綜合考慮了醫(yī)院的滿意度和患者

20、的滿意度，從雙方互利的目的出發(fā)給出了病床的合理安排方案；還可以利用計(jì)算機(jī)模擬的方法由患者的就診時(shí)間預(yù)測(cè)住院時(shí)間和出院時(shí)間。該模型的不足在于：1.沒(méi)能用量化的方法準(zhǔn)確給出指標(biāo)體系；2.忽略了一些偶然因素，導(dǎo)致結(jié)果與實(shí)際情況有一定的偏差；但都不影響模型的實(shí)用性。七 參考文獻(xiàn)【1】 鄭阿奇，matlab實(shí)用教程，北京：電子工業(yè)出版社，2004。【2】 李大潛，中國(guó)大學(xué)生數(shù)學(xué)建模競(jìng)賽，北京：高等教育出版社，1998。【3】 盧向南，應(yīng)用運(yùn)籌學(xué)，杭州：浙江大學(xué)出版社，2005?！?】 陳杰，matlab寶典，北京：電子工業(yè)出版社，2007?！?】 吳建國(guó)，數(shù)學(xué)建模案例精編，北京：中國(guó)水利水電出版社，20

21、05?！?】 葉其孝，大學(xué)生數(shù)學(xué)建模競(jìng)賽輔導(dǎo)教材，湖南：湖南教育出版社，2000。附錄1（算法模擬）第一周期病床使用情況：x = 117l = 38程序：r=input(r=);x=13 1 1 0 0 1 1 1 0 1 0 0 1 0 0 1 0 1 0 0 1 0 0 0 0 1 0 1 1 0;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=18 16 15 15 15 18 17 19 20 20 19 20 19 18 18

22、 16 16 18 19 18 18 17 17 21 23 22 21 20 19 19;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x1=n-m;x=13 1 0 0 1 0 1 0 1 0 0 0 0 0 1 0 0 0 1 1 1 0 0 0 0 2 0 2 0;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=20 18 19 19 18 25 24 2

23、4 22 22 23 22 24 22 23 22 21 22 21 20 19 19 19 19 26 25 23 21 22;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x2=n-m;x=13 1 1 0 1 1 1 1 0 3 1 0 1 2 0;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=23 21 21 22 22 21 21 23 28 25 2

24、4 24 24 23 22;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x3=n-m;x=13 0 0 1 0 0 0 1 0 0 1 0 1 0 0 1 0 0 1 0 0 1 0 0 2 1 0 0 0 0 1 0 0 0 1 0;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=26 22 29 23 25 26 23 28 25 26 25 26 24 2

25、2 24 21 24 26 26 20 23 24 25 27 26 26 28 25 27 31 26 29 26 24 28 26;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x4=n-m;x=13 3 1 1 0 1 1 0 1 0 1 1 0;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=6 6 11 8 7 6 9 8 8 6 6 9 7;e=;fo

26、r i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x5=n-m;x=x1+x2+x3+x4+x5l=max(x-79,0)輸入 r=20,第二周期病床使用情況：x = 100l =21程序：r=input(r=);x=1 0 2 1 0 3 1 1 1 0 1 0 0 0 1 0 0 0 2;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=18 17 16 16 14 20 1

27、8 17 16 15 15 14 14 17 20 21 20 21 19;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x1=n-m;x=27 0 1 1 1 1 1 0 0 0 1 1 0 0 1 0 0 2 0 0 0 0 1 0 0 0 1 1;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=19 20 19 18 17 15 21 22 23 22 21

28、 20 20 20 20 19 19 17 18 24 24 24 23 24 23 23 22 20;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x2=n-m;x=27 1 0 0 0 1 0 2 0 0 0 1 2 0 2 2 1;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=23 23 22 23 22 24 24 18 16 19 23 22 20 2

29、1 20 29 26;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x3=n-m;x=27 0 1 0 0 0 0 1 1 1 0 0 0 1 0 0 1 0 0 0 1 0 0 0 0 1 1 0 1 0 0 0 1 0 0 0 1 0 0 0 0 1 0;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=i-1;d=26 27 28 28 27 27 23 26 21 24

30、 28 25 27 28 22 21 27 23 25 26 23 23 26 27 25 23 24 19 26 23 26 26 26 29 22 24 25 30 29 22 21 28 26;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x4=n-m;x=27 0 1 0 2 0 1 0 1 0 0 3 2 1 0 1;tt=size(x,2);y=;for i=1:1:tt if i=1 y(i)=x(i); else y(i)=x(i-1)+x(i); end if y(i)r break; endendn=

31、i-1;d=6 9 9 8 12 7 7 9 9 8 7 9 9 9 5 6;e=;for i=1:1:tt e(i)=y(i)+d(i); if e(i)r break; endendm=i-1;x5=n-m;x=x1+x2+x3+x4+x5l=max(x-79,0)輸入 r=31附錄2門(mén)診匯總2008-7-28 星期一11452122008-8-18 星期一144112122008-9-1 星期一113151112008-7-14 星期一1131492008-7-21 星期一125292008-8-25 星期9-8 星期一1413192008-8-4 星期一1131162008-8-11 星期一111461 匯總162392510832008-8-19 星期二22417142008-7-15 星期二23223102008-8-12 星期二2331