全國一等獎數(shù)學建模論文-眼科病床安排的數(shù)學模型

眼科病床安排的數(shù)學型本文解決的是醫(yī)院眼科病床的安排問題，現(xiàn)醫(yī)院安排病人入院的原則是先來先服務，這樣雖然公平，但缺乏合理性以致等待住院的病人隊列越來越長，為解決此問題，我們建立了三個最優(yōu)化模型。對于問題一：我們確定了三個評價指標：手術前的平均逗留時T，平均每天出院人數(shù)NO，病人手術前的準備時間。然后計算出在原來先來先服務的原則下各指標值g為T7.8605。q對于問題二：我們采用優(yōu)先級原則動態(tài)地對病床進行安排。首先，統(tǒng)計初始數(shù)據，通過6SQ軟件進行分布的卡方擬合檢驗得：每類病人的到來均服從泊松分布、術后觀察時間服從均勻分布。然后，我們發(fā)現(xiàn)合理的調度方案必須使得病人的術前準備時間盡量短。因此，重新制定入院規(guī)則：外傷優(yōu)先級始終最高；其它病的優(yōu)先級隨時間的變化而變化。接著，再以三個指標為目標函數(shù)，病人入院規(guī)則為約束建立了多目標的最優(yōu)化模型，最后，根據入隊與服務時間服從的分布，用計算機隨機模擬，得到在隊列穩(wěn)定時，此規(guī)則下三個指標值為：T10.311，9.633T1.6526；這樣手術前的q平均逗留時間減少21.6%，平均每天出院人數(shù)增加22.55%，平均術前準備時間減少了32.31%。對于問題三：在問題二的計算機隨機模擬的基礎上，已經可以求得對應的等待隊列中病人的入院時間的模擬結果，因為存在一定隨機性，我們模擬次，取出每次所得結果中的模擬入院時間，作為病人的一個大致入院時間。對于問題四于星期六與星期日不安排除了外傷手術的其它手術安排在周四，五住院的視網膜和青光眼病人的手術要推遲到下周二、四，以此我們同樣建立了多目標的最優(yōu)化模型，得出在隊列穩(wěn)定時，三個指標值分別為：T10.436，NOqT；g對于問題五：為便于醫(yī)院的管理，可根據各類病人服從的分布按照比例給各類病人安排固定的病床數(shù)，但要先單獨分配外傷類的病床，因為醫(yī)院要保證有足夠的床鋪滿足外傷類病人，據統(tǒng)計結果知外傷病人到達和外傷病人被服務的時間都是服從泊松分布，則先建立排隊論中的M/M/C模型求出分配給外傷病人的病床數(shù)余下的病床按照一定的比例分配給其它類的病人。為得到平均逗留時間最短，我們建立了單目標最優(yōu)化模型。關鍵詞：優(yōu)先級

調度排隊論

計算機模擬

最優(yōu)化1

醫(yī)院就醫(yī)排隊是大家都非常熟悉的現(xiàn)象，它以這樣或那樣的形式出現(xiàn)在我們面前，例如，患者到門診就診、到收費處劃價、到藥房取藥、到注射室打針、等待住院等，往往需要排隊等待接受某種服務。在本文中，我們考慮某醫(yī)院眼科病床的合理安排的數(shù)學建模問題。該醫(yī)院眼科門診每天開放，住院部共有病79張。該醫(yī)院眼科手術主要分四大類：白內障、視網膜疾病、青光眼和外傷。附錄一中給出了2008年月13日至20089月11日這段時間里各類病人的情況。各類眼病手術的安排情況：白內障手術：較簡單，而且沒有急癥。目前該院是每周一、三做白內障手術，此類病人的術前準備時間只需1、2天。做兩只眼的病人比做一只眼的要多一些，大約占到60%。如果要做雙眼是周一先做一只，周三再做另一只。外傷手術通常屬于急癥床有空時立即安排住院住院后第二天便會安排手術。視網膜、青光眼手術：比較復雜，有各種不同情況，但大致住院以2-3天內就可以接受手術，主要是術后的觀察時間較長。這類疾病手術時間可根據需要安排，一般不安排在周一、周三。由于急癥數(shù)量較少，建模時這些眼科疾病可不考慮急癥。該醫(yī)院眼科手術條件比較充分，在考慮病床安排時可不考慮手術條件的限制，但考慮到手術醫(yī)生的安排問題，通常情況下白內障手術與其他眼科手術（急癥除外）不安排在同一天做。當前該住院部對全體非急癥病人是按照（Firstcome,Firstserve規(guī)則安排住院，但等待住院病人隊列卻越來越長，醫(yī)院方面希望你們能通過數(shù)學建模來幫助解決該住院部的病床合理安排問題，以提高對醫(yī)院資源的有效利用。本文需解決的問題有：問題一分析確定合理的評價指標體系以評價該問題的病床安排模型的優(yōu)劣。問題二：試就該住院部當前的情況，建立合理的病床安排模型，以根據已知的第二天擬出院病人數(shù)來確定第二天應該安排哪些病人住院并對你們的模型利用問題一中的指標體系作出評價。問題三：作為病人，自然希望盡早知道自己大約何時能住院。能否根據當時住院病人及等待住院病人的統(tǒng)計情況，在病人門診時即告知其大致入住時間區(qū)間。問題四：若該住院部周六、周日不安排手術，請你們重新回答問題二，醫(yī)院的手術時間安排是否應作出相應調整?問題五：有人從便于管理的角度提出建議，在一般情形下，醫(yī)院病床安排可采取使各類病人占用病床的比例大致固定的方案，試就此方案，建立使得所有病人在系統(tǒng)內的平均逗留時間（含等待入院及住院時間）最短的病床比例分配模型。模型的假2.1模型的假設假設1：題目所給數(shù)據是合理、正確的假設2：視網膜與青光眼兩類病不考慮急癥假設3：白內障手術與其他眼科手術（急癥除外）不安排在同一天假設4：該醫(yī)院眼科手術條件比較充分，在安排病床時不考慮手術條件的限制假設5：對于問題四，假定周六和周日可安排外傷手術，但不能安排其他手術2.2符號說明2

符號(i)(jq(i)(jg(i)(j)fNO

符號說明患i類病的第j個病人的等待入院時間患i類病的第j個病人的術前的準備時間患i類病的第j個病人的住院時間平均每天出院人數(shù)T

手術前的平均逗留時間，即從門診到第一次手術的平均時間Tg

病人的平均術前的準備時間T

f

病人的平均住院時間T(i)(j)hT(i)(j)oi)N(i)(j)NIi)(j(i)(j)

患i類病的第j個病人的入院時T(ij)1,2,3,(ij)示2008hh年7月13日(i)(j)表示年7月14日，依此類推?）h患第i類病的第j個病人的手術時間（(i)(j)的數(shù)方式與(i)(j同，即ohT(ij)示日T(i)(j)表示714，依oo此類推??）一段時間內到門診看病的i病人的人數(shù)第ji病人的在院人數(shù)（不包括當天新入院的人數(shù)）第ji病人的新入院人數(shù)第ji病人的出院人數(shù)

i

i類病人平均每天到門診看病的人數(shù)Ci

分配給i類病人的病床數(shù)，單位為張%求余符號，等價于modi=1類病表示白障（單眼）疾病，i=2類病表示白內障（雙眼）疾病，i=3類病表示視膜疾病，i=4類病表示青光眼疾病，i=5類病表示外傷疾病此題研究的是某醫(yī)院眼科病床合理安排的數(shù)學建模問題。要對病床進行合理的安排，就要有合理的安排規(guī)則，尤其是在醫(yī)院病床不夠的時候。當前該住院部對全體非急癥病人是按照FCFS（FirstFirstserve）規(guī)則安排住院，這樣雖然對病人很公平，但缺乏合理性，例如根據FCFS原則，白內障雙眼的病人可能會在星期二入院，但醫(yī)院規(guī)定“白內障雙眼的患者在星期一做一只眼，在星期三做另一只眼”，所以該患者的術前準備時間就變?yōu)?天，而一般情況下，白內障患者的術前準備時間只12天；3

這樣的情況會延遲其它類病人的入院時間，進而使得病人隊列越來越長。在病床不夠的情況下從醫(yī)院的角度講院自然希望在多做手術的同時減小病人占用病床的時間。為了得到合理的安排規(guī)則，首先要確定合理的評價指標體系，用此評價按該規(guī)則建立的病床安排模型的優(yōu)劣。針對問題一：從病人的角度看，病人到醫(yī)院看病分為以下幾個階段:掛號看病時間即門診時間，入院，手術前的準備，手術，手術后的觀察，出院。合理的安排就是讓病人從掛號看病到出院的時間盡量的短但根據實際情況知病人的術后觀察時間是由病情決定的故所建立的模型只能縮短門診看病到接受手術的時間間隔即病人手術前的逗留時間，所以模型的評價指標可以是病人手術前的平均逗留時間，平均術前準備時間。從醫(yī)院的角度看，我們可以將病床的周轉次數(shù)作為評價指標。由于病床的周轉次數(shù)與醫(yī)院每天出院人數(shù)是密切相關的——在病床不夠的情況下，醫(yī)院每天出院的人數(shù)越多，能夠入院的病人就越多，病床周轉次數(shù)就越多，醫(yī)院的效益就越好。所以，綜合考慮病人和醫(yī)院的利益，我們把病人手術前的平均逗留時間，平均術前準備時間，平均每天出院人數(shù)作為評價指標，當前兩個指標值越小，最后一個指標值越大的時候，病床安排模型越好。針對問題二：在確定病人入院規(guī)則時，要考慮以下幾點:白內障病人只能安排在周一與周三做手術，而其它病人除外傷病人外不能在周一和周三做手術，還有不同的病人的術前準備時間不一樣；使得建立的入院規(guī)則能夠讓病人的等待時間盡量短，這就可以縮短病人在醫(yī)院的時間，亦可縮短病人從門診到入院的時間。再對題給數(shù)據進行分析，得出各類病人到門診看病的統(tǒng)計規(guī)律，然后就可以以問題一的評價指標作目標函數(shù)，以安排病人入院規(guī)則為約束建立一個病床安排的多目標最優(yōu)化模型。針對問題三：根據各類病人的統(tǒng)計結果，可得出每天有多少人患病以及患什么病，找出其分布規(guī)律，然后根據病人的入院規(guī)則，可以得出各類病人大致在星期幾入院，再根據術后觀察時間的統(tǒng)計規(guī)律，便可以得到病人的出院時間，從而可安排病人入院，這樣就可在病人門診時告知其大致的入院時間。針對問題四：同問題二一樣，以問題一的評價標準作為目標函數(shù)，建立一個病床安排的多目標最優(yōu)化模型。但由于周六、日不安排手術，會使得約束條件發(fā)生改變。針對問題五：從便于管理的角度醫(yī)院可以根據各類病人的到達規(guī)律安排病床，故先統(tǒng)計出各類病人的到達服從什么樣的分布再建立模型求出平均逗留時間最短時的病床分配方案，但在分配時要把外傷類除外，因外傷類病人不允許等待，故分派給外傷病人的病床必須保證每天都能滿足需入院的外傷病人，因此先分配外傷類的病床，再統(tǒng)一分配余下各類病的病床。定義1術后觀察時間指病人出院與第一次手術的時間間隔定義2手術前的平均逗留時間指門診到第一次手術的平均時間根據對題給數(shù)據的統(tǒng)計結果知：4.1各類病的術前準備時間時間為1-7（包括1天和7天）4.2等待住院病人隊列越來越長的原因因為08年713日-08年9月11平均每天到門診看病的人數(shù)為8.6885人，平均每天出院人數(shù)為8.1163人，所以近似認為平均每天到門診看病的人數(shù)為人，平均每天出院人數(shù)為8.1163人。正因為每天到門診看病的人數(shù)大于每天出院的人數(shù)，所以才導致了等待住院病人隊列越來越長。4.3各類病人的到達（病人到達時間指病人的門診時間）服從泊松分布4

由于醫(yī)院就醫(yī)排隊是典型的排隊論問題，而一般的排隊論模型都是泊松輸入，所以我們先假定病人的到達服從泊松分布然后根據附錄一給出的數(shù)據求出每天到門診看的各類病人的人數(shù)（統(tǒng)計結果見附錄二）再利用6SQ統(tǒng)計軟件進行卡方擬合優(yōu)度檢驗得到如下結果（結果見表4.1）：表4.1：各類病人的到達服從泊松分布的卡方擬合優(yōu)度檢白內障單眼患者到達時間：白內障雙眼患者到達時間：視網膜患者到達時間：假設檢驗

假設檢驗

假設檢驗零假設服從松分布零假設服從泊松分布零假設服從泊分布自由度3自由度4自由度4卡方統(tǒng)計量1.8卡方統(tǒng)計量1.9卡方統(tǒng)計量3.7p值0.6p值0.8p值0.5顯著性水平0.1顯著性水平0.1顯著性水平0.1結果接受假設結果接受零假設結果接受零設青光眼患者到達時間：外傷患者到達時間：假設檢驗

假設檢驗零假設服從松分布零假設服從泊松分布自由度3自由度2卡方統(tǒng)計量4卡方統(tǒng)計量1p值0p值1顯著性水平0顯著性水平0結果接受假設結果接受零假設根據以上假設檢驗的結果知：各類病人的到達時間均符合泊松分布。根據指數(shù)分布與泊松分布的關[1]：如相繼兩事件出現(xiàn)的間隔時間服從參數(shù)的指數(shù)分布，則在單位時間間隔內事件出現(xiàn)的次數(shù)服從參數(shù)泊松分布，即單位時間內該事件出k次的概率為e()(k,n)k!

⑴將各類病人平均每天到達的人數(shù)，值代入式⑴，可得相應的概率密度函數(shù)，計算結果見表4.2。表4.2：各類病人達到時間的概率密度函數(shù)和相應類別白內障單眼白內障雙眼視網膜青光眼外傷1.62.22.811()

1.6

kk!

2.2

k!

2.8

kkee!k!k!k!!4.4各類患者的術后觀察時間服從均勻分布首先統(tǒng)計出各類病人的術后觀察時間（統(tǒng)計結果見附錄三），根據統(tǒng)計結果，我們假定各類病人的術后觀察時間服從均勻分布然后通過6SQ統(tǒng)計軟件進行卡方擬合優(yōu)度檢驗，檢驗結果見表4.3。表4.3：各類患者的術后觀察時間服從均勻分布的方擬合度檢驗白內障單眼術后觀察時間白雙術后觀察時間視網膜術后觀察時間假設檢驗

假設檢驗5

假設檢驗

n(in(i)零假設自由度

服從均分布

零假設34自由度

服從均分布零假設39自由度

服從均分布99卡方統(tǒng)計量

7.858300287卡方統(tǒng)計量2.093007186卡方統(tǒng)計量55.04926108p值顯著性水平

0.999999107p值0.05顯著性水平

1p值0.05顯著性水平

0.9998973470.05結果

接受零設

結果

接受零設

結果

接受零設青光眼術后觀察時間假設檢驗

外傷術后觀察時間假設檢驗零假設自由度

服從均分布

零假設37自由度

服從均分布53卡方統(tǒng)計量

11.66013072卡方統(tǒng)計量29.25609756p值顯著性水平

0.999978078p值0.05顯著性水平

0.9966966650.05結果

接受零設

結果

接受零設根據以上假設檢驗的結果知：各類病人的術后觀察時間均服從均勻分布。4.5外傷病人住院時間服從泊松分布首先統(tǒng)計出外傷病人的住院時間，根據統(tǒng)計結果，我們假定其服從泊松分布，然后通過6SQ統(tǒng)計軟件進行卡方擬合優(yōu)度檢驗，檢驗結果見表表4.4:外傷病人住院時間服從泊松分布的卡方擬合優(yōu)度檢驗結果假設檢驗零假設服從松分布自由度6卡方統(tǒng)計量11.01007449p值0.088065567顯著性水平0.05結果接受假設5問本文研究的是某醫(yī)院眼科病床的合理安排的數(shù)學建模問題對于病床安排模型的優(yōu)劣，不能憑人們的主觀感受進行判斷，而要確定合理的評價指標體系進行判斷，為此我們確定了如下的評價指標體系：指標1：手術前的平均逗留時間手術前的平均逗留時間指門診到第一次手術的平均時間，其數(shù)學表達式為T

ij

i)(j(i)(j))g(i)i該指標值越小，表示病床安排模型越好。6

5n(i)55n(i)5指標二：平均每天出院人數(shù)平均每天出院人數(shù)的數(shù)學表達式為：NO

5ij

(i)(j)n該指標值越大，表示病床安排模型越優(yōu)。指標三：病人平均術前的準備時間病人平均術前準備時間的數(shù)學表達式為Tg

(i)(j)gij(i)i該指標值越小，表示病床分配模型越優(yōu)。對于題目中給出的以FCFS為原則（急癥除外）的病床安排模型，我們通過求解得到其三個評價指標分別為：T天q338NO7.8605人（從題給數(shù)據，在7月30日之前，只有外傷病人出院，但44在7月30日之后，各類病人均出院達到穩(wěn)定，所以7月30日之后的統(tǒng)計數(shù)據求解平均每天出院人數(shù)）T天g問題二的針對問題二我們建立了模型一。6.1模型一的建立6.1.1確定目標函數(shù)該模型是為了解決醫(yī)院的病床安排問題，為了使病床的安排更加合理，我們只需使三個評價指標——手術前的平均逗留時T最小，平均每天出院人NO最大，病人平均術前的準備時T最小即可，所以我們建立了如下的目標函數(shù)：g

TqNOT6.1.2確定約束條件ⅰ由于白內障手術比較簡單，此類病人的術前準備時間只需、2天，而且根據附表一知，各類病人的術前準備時間均在天之內（包括1天和7天），所以白內障單眼和雙眼病人術前準備時間為1-7天，即(i)(jiⅱ由于視網膜和青光眼疾病比較復雜，大致住院以2-3天內就可以接受手術，而7

且根據附表一知，各類病人的術前準備時間均在天之內（包括1天和7天），所以這兩類疾病的術前準備時間為2-7天，即(i)(j)igⅲ外傷疾病有空床時立即安排住院，且住院后第二天便會安排手術，所以此類病的術前準備時間為1天，即(5)(jgⅳ根據數(shù)據分析結果，為了縮短等待入院病人的隊列長度，我們制定了如下的病床安排原則：原則1白內障單眼患者一般安排在周一、周二、周六、周日入院原則2白內障雙眼患者一般安排在周六、周日入院原則3視網膜和青光眼患者安排在周三、周四、周五入院原則4外傷病人當天入院，第二天手術原則5當病人等待時間達到25天時，只要有空病床立即安排入院說明：外傷病人在任意一天優(yōu)先級是最高的，在某一天病人可安排住院此時此類?。ǔ鈧┤说膬?yōu)先級第二高，且若有兩類或兩類以上的?。ǔ鈧┤硕伎稍谕惶烊朐?，則這些病人優(yōu)先級第二高且相等，在這天不安排住院的病人的優(yōu)先級最低。例如白內障單眼患者在周一周二周六周日的優(yōu)先級第二高在其它的時間優(yōu)先級最低。故患者的優(yōu)先級隨時間動態(tài)的發(fā)生改變。因為T(i)(j示第i第j病（i)(j)，hhT(i)(j)表示7月13日T(i)(j)示2008年7月14日此類推??h又由于2008年7月13日是星期日，所以示該病人的入院時間為周一i)(j

2,表示該病人的入院時間為周表示該病人的入院時為周三%7表示該病人的入院間為周四表示該病人的入院時間為周五表示該病人的入院時為周六0,表示該病人的入院時間為周日

（2）根據式（2）和病床安排原則，我們確定了如下的約束條件：原1(1)(j)h原則j)h原則ijih原則(5)(j)g原則(i)(j)25q整約：(i)(j，(ij)1,2,3q6.1.3綜上所述，得到問題二的多目標最優(yōu)化模型

運8

5(i)55n5n(i)55(i)55n5n(i)5ghq

minTq

ij

(i)(j)(i)(j))qg(ii

NO

ij

NO(i)(j)

minTg

(i)(j)gij(iii)(ji2(ij)igTj)gj)s.t.(2)(j)hT(i)(j%73,4(ij)j，T(i)(j)1,2,3

運算

6.2模型一的求解首先我們由數(shù)據分析可知，對于每一類病人，從第一次接受手術到出院之間經歷的時間服從均勻分布，每天到來的各類病人的數(shù)目分別服從各自參數(shù)possion分布，我們首先可以模擬得79個已住院的人的出院時間（肯定在9.11號之后），就可以將排隊隊列的一定數(shù)量的病人送入服務隊列并且可以按照每天到來的各類病人的流來補充每天的排隊隊列人數(shù)。接著，由于本題涉及到大量的隨機現(xiàn)象，故用一般的規(guī)劃方法難以求得最優(yōu)解，所以我們采用一種計算機隨機模擬（相關程序見附錄五）的算法求得本題的最優(yōu)解。算法思想：從912開始，根據服務隊列中正在接受服務的對象的服務時間分布，模擬得到正在接受服務的對象的出隊時間（即病人的出院時間），然后按照約束條件中給定的的原則，將等待隊列中的對象分配進入服務隊列，然后改變等待隊列與服務隊列的狀態(tài)，進行下一步的模擬，直到我們得到一定數(shù)量滿足我們需要的對象（即：滿足退出條件），得出此時的天數(shù)即可。模擬的流程如下所示（注：退出條件即為我們需要用模擬分配病人的天數(shù)20或60）9

開始賦f為9月號其初f足出件

yes

輸結

結

No統(tǒng)f天時出院數(shù)按原進服務即：束條）改服隊與等隊記狀與出圖：計算機模擬的簡化流程圖模擬后的數(shù)據的存儲格式為：表6-1：模擬后數(shù)據的示例存儲格式編號

生病類別

診斷時間

入院時間

手術（一）時間

手術（二）時間

出院時間12345?10

白內障(眼視網膜疾病青光眼視網膜疾病視網膜疾病????視網膜疾病

4949494949?61

6362626262??74

6564646464??76

'/''/''/''/''/'??'/'

7171727779??90210

白內障627274[]76310

白內障(眼627779[]83410

10

白內障(眼627779[]85510

白內障(眼627779[]84610

白內障(眼627879[]84710

青光眼627476[]808?

?????????????注：第六列在102號以前是’’表示第二次手術時間，模擬時，我們沒有考慮第二次手術的時間。102以后是[]，表明該行數(shù)據是我們自己通過possion流對隊列追加顧客后的效果，是新加入排隊系統(tǒng)的，方便進行連續(xù)模擬。統(tǒng)計從2008年的9月12日始后60的隊列的長度，我們即可認為在這段時間內的接受服務的對象達到了穩(wěn)定狀態(tài)。然后作出隊列長度的隨時間的變化圖示：如下所示：隊長隨天數(shù)的變化趨勢1009080

7060504030

0102030天數(shù)

405060圖：隊長隨時間變化趨勢此時，在我們的分配原則下，可以發(fā)現(xiàn)隊列的長度不斷減小，并且有趨于穩(wěn)定的趨勢，我們可以認為此時的穩(wěn)態(tài)50，因此我們認為這種方案是比較合理的?，F(xiàn)在，我們取得從9月12日后的天一個模擬的入院安排情況，因為20已經足夠確定我們需要填滿102個排隊對象的的數(shù)據計算得我們定義的三個目標函數(shù)的最優(yōu)解為：病人手術前的平均逗留時間：12.1（越小越好）病人平均術前準備時間：1.6722天；（越小越好）平均每天出院人數(shù)：9人（越大越好）6.3模型一結果分析在這種分配方案下，我們發(fā)現(xiàn)與原模型的手術前的平均逗留時間：13.1519平均術前準備時間2.4413，平均每天出院人數(shù)7.8605相比，并沒有太大的差別，我們認為這是在隊列沒有進入穩(wěn)定狀態(tài)時的統(tǒng)計數(shù)據造成的，我們用以下的方法處理：我們可以認為當系統(tǒng)服務了100個顧客后，它已經進入了穩(wěn)定狀態(tài)，又由于我們要11

去足夠的數(shù)據才能具有說服力，因此我們定義一個評價區(qū)間：從100個排隊等待手術的病人到第30天結束時最后一個出院的病人。在此區(qū)間上，我們再用同樣的方法進行評價，會發(fā)現(xiàn)我們的三個評價指標值為：病人手術前的平均逗留時間：10.311病人平均術前準備時間：1.6526天；平均每天出院人數(shù)：9.633人此時我們就可以發(fā)現(xiàn)，當這個排隊系統(tǒng)在盡量趨于穩(wěn)定狀態(tài)時，它的手術前的平均逗留時間、術前準備時間、平均每天出院人數(shù)均比前邊的結果有了一定的優(yōu)化，這是由于9月12日后的20天的排隊系統(tǒng)受醫(yī)院最初的先來先服務的影響較大而當系統(tǒng)服務了100個病人后，此時的排隊系統(tǒng)趨于穩(wěn)定，所以求得的結果較優(yōu)。從而進一步證明我們的排隊系統(tǒng)比原有的效率更高。問題三的根據問題二的模型，我們已經完全模擬出來了每位病人的入院時間、第一次手術時間、出院時間（如上表所示），所以我們可以求得排隊隊列102人的入院時間，但是每次隨機模擬的結果均不相同，所以，我們可以通過模擬若干次，求出每次每一個病人的出院時間，從中選擇一個最大值和一個最小值，將它作為病人的一個大概的入院時間的區(qū)間。我們取模擬的次數(shù)為10，下面是我們所得的一個近似的結果：表7-1：10次模擬后產生的10個模擬的出院時間編號生病類別一二三四五六七八九十1

白內障(雙眼)63636363

6363636363

632視網膜疾病6262626262626262623青光眼6262626262626262624視網膜疾病6262626262626262625視網膜疾病626262626262626262??????????????99視網膜疾病767474697474747574100白內障667372727278727171101視網膜疾病767574697474747574102視網膜疾病767574747474747574（注：具體數(shù)據見附錄一）那么我們就可以根據以上表格中的數(shù)據確定出病人的大致入院區(qū)間：表7-2：病人的大致入院區(qū)間

62626262?69726974編號12345?99

生病類別白內障雙)視網膜疾病青光眼視網膜疾視網膜疾病????視網膜疾病

最佳入院時間6362626262?,75

對應的日期9-139-129-129-129-12?月日月24日月25日月日100

白內障

,789月日,月21,月日,月23日月28日101

視網膜疾病

,75

月日月24日月25日月日12

5(i)55n5(i)55n5n(i)5102

視網膜疾病74---769-24---9-26日（注：完整數(shù)據見附錄四）問題四的針對問題四我們建立了模型二。問題四與問題二的區(qū)別在于：在問題二中，醫(yī)院每天都可以安排手術，而在問題四中，只能在周一至周五安排手術（外傷每天均可安排手術）。8.1模型二的建立8.1.1確定目標函數(shù)（同模型一的目標函數(shù)）以三個評價指標最優(yōu)為目標函數(shù)：

TqNOT8.1.2確定約束條件該模型的約束條件除了包含模型一的約束條件外，還有以下幾個：由于白內障手術之后安排在周一和周三，外傷手術每天都可以安排，所以周六和周日不安排手術只會影響視網膜和青光眼的手術安排。在模型一給出的病床安排原則下，對于視網膜病人和青光眼病人，將其中周三入院的手術安排在同一周的周五，在周四和周五入院的手術安排在下周周二。據此，我們又建立了如下的約束條件：T(ij)i)(j)()(j)hhT(ij)i)(j)(i)(j)hhT(ij)i)(j)()(j)hh8.1.3綜上所述，得到問題四的多目標優(yōu)化模型

minTq

ij

(i)(j)(i)(j))qg(ii

NO

ij

NO(i)(j)

minTg

(i)(j)gij(ii13

gqqhohhgqqhohh(i)(j)i(i)(j)ig(5)(j)gi)(j)i)(j，(i)(j)1,2,3j)hj)h(i)(j)ij)(i)(j)(i)(j)ij)(T(i)(j)oh(i)(j)(i)(j)(i)(j)oh8.2.模型二的求解

%表示求余運算）我們按照與問題二相同的思想，按照同樣的原理進行計算機模擬（相關程序見附錄五），在這里只是改變了對于青光眼和視網膜疾病的分配方案：若這兩類患者在周四、周五分配入院，則他們均到下周二進行手術。統(tǒng)計從2008年的9月12號始后60的隊列的長度，我們即可認為在這段時間內的接受服務的對象達到了穩(wěn)定狀態(tài)。然后作出隊列長度的隨時間的變化圖示：如下所示：隊長隨天數(shù)的變化趨勢120110100度長列隊的數(shù)天應對

9080706050

0102030天數(shù)

405060圖8-1:隊長隨時間的變化趨勢此時，我們同時發(fā):當我們忽略外傷病人的等待時間時，隨著時間的推移，病人的等待時間的規(guī)律如下圖所示：14

0kC!0kC!!K等待時間的變化趨勢

2520151050

050100150200250300病例序號

350400450500圖8-2:病人的等待時間隨時間的變化趨勢此時，在我們的分配原則下，隊列的長度不斷減小，并且有趨于穩(wěn)定的趨勢，我們可以認為此時的穩(wěn)態(tài)等待人數(shù)是60而且人從就診到第一次手術的時間也有一個逐步下降的趨勢，所以因此我們認為這種分配的方案是比較合理的：據此我們可以在這種條件下求得的三個指標值：現(xiàn)在，我們同樣取得9月12號后的20天的一個模擬的入院安排情況，然后計算得我們定義的三個目標函數(shù)的最優(yōu)解為：手術前的平均逗留時間：12.691天病人平均術前準備時間：2.2191天；平均每天出院人數(shù)：8.9人8.3模型二的結果分析在這種分配方案下，我們發(fā)現(xiàn)與原模型的手術前的平均逗留時間：13.1519平均術前準備時間2.4413，平均每天出院人數(shù)7.8605相比，并沒有太大的差別，所以可以同樣按照問題二中的解決辦法，我們會發(fā)現(xiàn)我們的三個評價指標值為：手術前的平均逗留時間：10.432天病人平均術前準備時間：2.017天；平均每天出院人數(shù)：9.1667人此時我們就可以發(fā)現(xiàn)，當這個排隊系統(tǒng)在盡量趨于穩(wěn)定狀態(tài)時，它的各個指標均比前邊的結果有了一定的優(yōu)化，從而進一步證明我們的排隊系統(tǒng)比原有的效率更高。問題五的針對問題五我們建立了模型三。相關知識引入排隊中有一種多服務臺多顧客的模型機制模型.顧客到達服務臺具有隨機性，服從泊松分布，服務臺對顧客服務服從指數(shù)分布，則：P

11

15

!!(i)5n(i)5n1nCC!其中：

：單位時間內顧客到達服務臺的平均人數(shù)：單位時間內服務臺服務的人數(shù)

，C:服臺的個數(shù)P():P():顧在問題五中，有人從便于管理的角度建議將各類病人占用的病床數(shù)大致固定”，就此方案建立使得所有病人在系統(tǒng)內的平均逗留時間最短的病床比例分配模型。在分配病床時因外傷病人特殊，故應先分配病床，從數(shù)據統(tǒng)計知患外傷的病人到醫(yī)院的人數(shù)服泊松分布，醫(yī)院對其服務的時間同樣股從

17

的泊松分布，故可把病床看作排隊系統(tǒng)中的服務臺，把外傷病人看作顧客，這樣就可以建立一個M/M/C的模型，則只要C(5)使得P(5)我們設定(n(5))，就能求出C，即為安排給外傷病人的病床數(shù)。那Ci)9.1模型三的建立9.1.1確定目標函數(shù)

)4)i

*(79(5))i1,2,3,設病人在在系統(tǒng)內的平均逗留時間，ww

ij

)(j)(i)

f

(i)(j)

i其(i)(j)表示患i類病的第j個病人的等待入院時間T(i)(j)示患i病qf的第j個病人的住院時間n(i表示一段時間內患i病的人數(shù)）所以，目標函數(shù)為T

ij

j)(i)(j)qf(i

9.1.2確定約束條件

i因各類病人都分配了相應的病床故每一類病人在可在任意一天住院且是先來先服務的，這一點和模型一，二不一樣，故約束條件只有：(1)各類病人的術前準備時間(2)白內障的手術只能安排在周一三。（3）由于分配給各類病人的病床數(shù)的比例大致固定，所以每天各類病人分別占用的病床數(shù)不能超過醫(yī)院分配的數(shù)目，由此得到如下的約束：NI(i)(j)(i)(j)NO(i)(ji),ij1,2,…16

n(i)ggn(i)gg9.1.3綜上所述，得到問題五的單目標非線性優(yōu)化模型T

ij

)(j)(i)(j)f(i)

ii)(ji1,22(i)(j)i3,4g(5)(j)（i4,5;j1,2,）T(ij20qT(ij)1,2,3h(i)(j(i)(ji)(j)(i9.2模型三的求解與結果分析問題五：P(0.05，即外傷病人到醫(yī)院不需等待的概率小0.05，通過計算得當分配12張床給外傷病人時則外傷病人等待入院的概率小于0.035（基本為不發(fā)生）。按照比例求得其它類病分配的病床如下病型白內障白內障視青（雙）（單）網膜光眼分配1914259的床至于此條件下的安排病人入院規(guī)則是各類病人都各自服從先來先到的原則視網膜和青光眼病人的逗留時間會減少，因為這兩種病的術前準備時間可以控制在兩天內，這樣就可使病人的住院時間縮短，使排隊的隊長變短?？墒前變日喜∪说钠骄毫魰r間會增多因是這兩類病只能在周一三做手術，這就使得在周三——周日入院的病人要等到周一，三才能做手術，這樣使得病人的住院時間增長，使排隊的隊長變長。10.1模型評價優(yōu)點：（）根據我們定下的安排病人入院規(guī)則，建立的模型在一定程度上縮短了病人排隊的隊長，因為原模型的三個指，7.8605；我們建的qg模型三個指標：T10.3119.6331.6526這樣手術前的平均逗留時間減qg少21.6%，均每天出院人數(shù)增加了22.55%，平均術前準備時間減少了32.31%；（2）根據模型可推算出當前病人的出院時間，故我們把表二出院時間的填充了；（）利用我們建立的模型二，可根據第二天擬出院的病人確定病人入院的最佳時間；缺點：由于所給數(shù)據太少以致在統(tǒng)計數(shù)據時不是很準確，又由于計算機模擬帶有一定的隨機性，以致得到模型的三個指標不是很讓人滿意。10.2模型改進（1）查詢更多的數(shù)據，以使得統(tǒng)計結果更正確，也可使計算機模更少的數(shù)據或不模擬以減少不確定性因我們建模時沒有考慮到經濟性若考慮到不同的手術經費不17

一樣，則在制定安排病人入院規(guī)則時要考慮一定的優(yōu)先級，即手術費用高的優(yōu)先級高。（2）所建模型是針對當前所給數(shù)據的，對長遠病人入院和出院的測并不能很準確，故建模時應把時間加上去建立動態(tài)規(guī)劃模型。使所建模型能準確的預測出病人的入院和出院時間。10.3模型推廣我們建的模型不僅可用于醫(yī)院病床安排也可用于其它資源的安排還可用于諸如像試卷評價模型的其它類型的問題。[1]宋來忠，王志明,數(shù)學建模與實驗北京：科學出版社,2005。[2]《運籌學》教材編寫組編，運籌學（版），北京：清華大學出版社，[3]王玉升，排隊論模型及其在醫(yī)院管理中的作用中國醫(yī)院管理,，附錄一：2008-07-13到2008-09-11的該醫(yī)院的病人信息（略）附錄二：08年7月13日-08年9月11，各類病人每天到門診看病的人數(shù)7137147157167177187197207217227237247257267277287297307318月日8月日8月日8月日8月日

第1天第2天第3天第4天第5天第6天第7天第8天第9天第10天第11天第12天第13天第14天第15天第16天第17天第18天第19天第20天第21天第22天第23天第24天

白內障單眼113112332352122012001112

白內障雙眼132264115020202111141330

青光眼112111200121021420410201

視網膜343233330154202511434512

外傷100112122120002202230010

合計79107121210996167549125611116116518

8月日8月日8月日8月日8108118128138148158168178188198208218228238248258268278288298308319月日9月日9月日9月日9月日9月日9月日9月日9月日910911

第25天第26天第27天第28天第29天第30天第31天第32天第33天第34天第35天第36天第37天第38天第39天第40天第41天第42天第43天第44天第45天第46天第47天第48天第49天第50天第51天第52天第53天第54天第55天第56天第57天第58天第59天第60天第61天合計

4401203101114221110321121213211004241100

5411213532324410232211332131137211034133

011011011200110002210012121004110021263

4452141615221741113223425351551123115170

121020001013202111012302001112310120164

141584867136107812149358797891198116915135397913530附錄三：各類病人術后觀察時間白內障單眼術后觀察時間

白內障(眼)術后觀察時間

視網膜術后觀察時間

青光眼術后觀察時間

外傷術后觀察時間35

12

9

419

2222442343443322343433244333332342433323223

4555555445464654555555556445455554456555556

81381011814111010111081071011125889111010129111410131071210101113131212812

6787778119988879878710106479878612889799128

496547664475477610557765777348686946371044410720

2323323343333343433322222333

55565555655545545556565555656645454565

812141113158713911129912139810513101111111479914136511910156

46151215101021

88788118118910797附錄四：未入院的102名病人的最佳入院時間編號生病類別

最佳入院時間1白內障(雙眼)9月13日2視網膜疾病9月12日3

青光眼9月12日4視網膜疾病9月12日5視網膜疾病9月12日6白內障(雙眼)9月13日78

白內障9月13日青光眼9月12日9白內障(雙眼)9月13日10視網膜疾病9月12日11視網膜疾病9月12日12視網膜疾病9月12日1314

青光眼9月12日，9月17日白內障9月13日，9月14日15視網膜疾病9月12日，9月17日16視網膜疾病9月12日，9月17日17

青光眼9月12日，9月17日18白內障(雙眼)9月13日，9月14日19白內障(雙眼)9月13日，9月14日20白內障(雙眼)9月13日，9月14日21視網膜疾病9月12日，9月17日，9月18日22

白內障9月13日，9月14日23視網膜疾病9月12日，9月17日，9月18日24視網膜疾病9月12日，9月17日，9月18日2526

白內障9月13日，9月14日白內障9月13日，9月14日27白內障(雙眼)9月13日，9月14日，9月20日22

2829303132333435363738394041424344454647484950515253545556575859606162636465666768697071

白內障視網膜疾病視網膜疾病視網膜疾病白內障(雙眼)白內障視網膜疾病白內障視網膜疾病視網膜疾病白內障(雙眼)白內障青光眼視網膜疾病視網膜疾病視網膜疾病青光眼白內障(雙眼)白內障(雙眼)青光眼青光眼視網膜疾病視網膜疾病白內障(雙眼)白內障(雙眼)白內障(雙眼)視網膜疾病白內障(雙眼)青光眼白內障(雙眼)白內障白內障(雙眼)白內障(雙眼)白內障(雙眼)視網膜疾病青光眼白內障(雙眼)視網膜疾病白內障(雙眼)視網膜疾病白內障視網膜疾病視網膜疾病白內障

9月14日，9月15日9月12日，9月17日，9月18日9月12日，9月17日，9月18日9月12日，9月17日，月18日，9月19日9月14日，9月20日9月12日，9月17日，月18日，9月19日9月12日，9月17日，月18日，9月19日9月14日，9月15日9月17日，9月18日，9月19日9月17日，9月18日，9月19日9月14日，9月20日9月14日，9月15日，9月16日9月17日，9月18日，9月19日9月17日，9月18日，9月19日9月17日，9月18日，9月19日9月17日，9月18日，9月19日9月18日，9月19日9月14日，9月20日9月14日，9月20日9月18日，9月19日9月18日，9月19日9月18日，9月19日9月18日，9月19日9月14日，9月20日9月14日，9月20日，9月21日9月14日，9月20日，9月21日9月18日，9月19日9月20日，9月21日9月18日，9月19日9月20日，9月21日9月15日，9月16日9月20日，9月21日9月20日，9月21日9月20日，9月21日9月18日，9月19日9月19日9月20日，9月21日9月19日9月20日，9月21日9月19日9月15日，9月16日，9月21日9月19日9月19日，9月24日9月15日，9月16日，9月21日23

72白內障(雙眼)9月20日，9月21日73

白內障9月15日，9月16日，9月21日74視網膜疾病9月19日，9月24日757677

白內障9月15日，9月16日，月20日，9月21日，922日青光眼9月19日，9月24日青光眼9月19日，9月24日78視網膜疾病9月19日，9月24日7980

白內障9月15日，9月16日，月20日，9月21日，922日白內障9月15日，9月16日，月20日，9月21日，922日81視網膜疾病9月19日，9月24日82

白內障9月16日，9月20日，9月21日，9月22日83白內障(雙眼)9月20日，9月21日，9月22日8485

白內障9月16日，9月20日，9月21日，9月22日白內障9月16日，9月21日，9月22日，9月23日86白內障(雙眼)9月20日，9月21日，9月27日8788

白內障9月16日，9月21日，9月22日，9月23日青光眼9月19日，9月24日89白內障(雙眼)9月20日，9月21日，9月27日90視網膜疾病9月19日，9月24日，9月25日91視網膜疾病9月19日，9月24日，9月25日92

青光眼9月19日，9月24日，9月25日93白內障(雙眼)9月20日，9月21日，9月27日94白內障(雙眼)9月20日，9月21日，9月27日95

青光眼9月19日，9月24日，9月25日96白內障(雙眼)9月20日，9月21日，9月27日97

外傷9月12日98白內障(雙眼)9月21日，9月27日99視網膜疾病9月19日，9月24日，月25日，9月26日100白內障9月16日，9月21日，9月日，9月23日，928日101視網膜疾病9月19日，9月24日，月25日，9月26日102視網膜疾病9月24日，9月25日，9月26日附錄五：模擬所用的程序（1）函數(shù)M文件jsdata.m

用來計算連個日期之間的天數(shù)（2）模擬的程序：clear;clc24

loadt1,loadt2,loadt3;%t1t2t3分別是給出的三部分的數(shù)據%時間統(tǒng)一化處理fori=1:size(t1,1)t1{i,2}=jsdata('2008-7-12',t1{i,2});t1{i,3}=jsdata('2008-7-12',t1{i,3});t1{i,4}=jsdata('2008-7-12',t1{i,4});t1{i,6}=jsdata('2008-7-12',t1{i,6});if~strcmp(t1{i,5},'/')t1{i,5}=jsdata('2008-7-12',t1{i,5});endendfori=1:size(t2,1)t2{i,2}=jsdata('2008-7-12',t2{i,2});t2{i,3}=jsdata('2008-7-12',t2{i,3});t2{i,4}=jsdata('2008-7-12',t2{i,4});if~strcmp(t2{i,5},'/')t2{i,5}=jsdata('2008-7-12',t2{i,5});endendfori=1:size(t3,1)t3{i,2}=jsdata('2008-7-12',t3{i,2});end%計算機模擬產生對于t2中的79個人的出院情況%已知到第61天都還沒有出院，所以最早的出院時間是第天fori=1:size(t2,1)ifstrcmp(t2{i,1},'白內障')time=unidrnd(3)+1;t2{i,6}=t2{i,4}+time;elseifstrcmp(t2{i,1},'白內障雙眼)')time=unidrnd(3)+3;t2{i,6}=t2{i,4}+time;elseifstrcmp(t2{i,1},'青光眼)time=unidrnd(9)+3;t2{i,6}=t2{i,4}+time;elseifstrcmp(t2{i,1},'視網膜疾病)time=unidrnd(11)+4;t2{i,6}=t2{i,4}+time;elseifstrcmp(t2{i,1},'外傷')time=unidrnd(8)+2;t2{i,6}=t2{i,4}+time;endift2{i,6}<62t2{i,6}=62;25

endendfori=1:size(t1,1)t1{i,7}=mod(t1{i,2}-1,7);t1{i,8}=mod(t1{i,3}-1,7);t1{i,9}=mod(t1{i,4}-1,7);t1{i,10}=mod(t1{i,6}-1,7);t1{i,11}=11;endfori=1:size(t2,1)t2{i,7}=mod(t2{i,2}-1,7);t2{i,8}=mod(t2{i,3}-1,7);t2{i,9}=mod(t2{i,4}-1,7);t2{i,10}=mod(t2{i,6}-1,7);t2{i,11}=11;endfori=1:size(t3,1)t3{i,7}=mod(t3{i,2}-1,7);t3{i,8}=-1;t3{i,9}=-1;t3{i,10}=-1;t3{i,11}=11;endtt=t2;%初始的序列；k=1;in={};inn=[];long=102;lon=[];f=62;%9月12號的序號，以后以此遞增；whilek<=30index=[];%存儲需要記錄的下標號。n=0;%統(tǒng)計每天調度時的可用的病床數(shù)fori=1:size(tt,1);iftt{i,6}==fn=n+1;endendlong=long-n;fori=1:size(t3,1)%首先安排傷的病人ifstrcmp(t3{i,1},'外傷')&t3{i,11}~=0&n>0index=[index,i];t3{i,11}=0;t3{i,3}=f;t3{i,4}=f+1;t3{i,6}=f+1+unidrnd(8)+2;t3{i,8}=mod(t3{i,3}-1,7);26

t3{i,9}=mod(t3{i,4}-1,7);t3{i,10}=mod(t3{i,6}-1,7);n=n-1;endendifmod(f-1,7)==2%這天是周二時：fori=1:size(t3,1)%此時確定是安排白內障還是（青光眼、視網膜）ifstrcmp(t3{i,1},'白內障)&t3{i,11}~=0&n>0index=[index,i];t3{i,11}=0;t3{i,3}=f;t3{i,4}=f+1;t3{i,6}=f+1+unidrnd(3)+1;t3{i,8}=mod(t3{i,3}-1,7);t3{i,9}=mod(t3{i,4}-1,7);t3{i,10}=mod(t3{i,6}-1,