數(shù)學建模眼科病床分配論文

1、眼科病床的合理安排模型摘要本文主要研究醫(yī)院眼科病床的合理安排問題?，F(xiàn)醫(yī)院對全體非急癥病人是按照FCFS（First come, First serve）規(guī)則安排住院，這樣雖然公平，但缺乏合理性，以致等待住院的病人隊列越來越長。我們可以運用運籌學中排隊論12的相關知識來解決此問題, 以縮短病人等待時間、提高對醫(yī)院資源的有效利用率。關于問題一 ：運用排隊理論，我們制定了平均月床位周轉次數(shù)、平均等待隊長和平均病人等待時間三個指標。運用題中所給數(shù)據(jù)，計算優(yōu)化前各個指標量：月平均床位周轉次數(shù)為3.33次/月；平均等待隊長為97人/天；病人平均等待時間為12.6天。用matlab軟件模擬出尚未住院病人數(shù)據(jù)

2、，得出各個指標量，將優(yōu)化前指標與模擬數(shù)據(jù)指標對比分析，以此評價模型好壞。關于問題二：運用6SQ統(tǒng)計軟件，對已知病人術后觀察時間進行服從均勻分布的卡方擬合優(yōu)度檢驗，結果認為各類病人的術后觀察時間均服從均勻分布。以此為依據(jù)用matlab軟件模擬尚未出院病人的術后觀察時間，得出第二天擬出院病人及出院病人數(shù)，按手術準備時間最短原則，安排各個病人入院時間、手術時間。同樣產生標準分布隨機數(shù)預測出第三段病人的出院時間，從而做出第三段數(shù)據(jù)病床的合理安排。得問題一中各個評價指標量：病人平均等待時間為11.15天，平均等待隊長72.78人/天；月平均床位周轉次數(shù)為4.76次/月?？梢钥闯瞿Ｐ偷男Ч翰∪似骄却龝r

3、間縮短1.45天，平均等待隊長縮短25人/天，月平均床位周轉次數(shù)比之前提高1.43次/月。關于問題三：依靠問題二中的優(yōu)化模型，利用matlab軟件模擬出每位病人的入院時間、第一次手術時間、出院時間，因為每次隨機模擬的結果均不相同，我們對排隊隊列中每位病人的入院時間均模擬5次，從中選擇一個最大值和一個最小值，作為等待病人的入院時間區(qū)間。關于問題四：在周六、日不做手術的情況下重新制定相應的入院規(guī)則，運用新規(guī)則做出第三段數(shù)據(jù)病床的合理安排，再分別計算其評價指標。最優(yōu)方案為：在周六、日不做手術的情況下，周三、五安排白內障手術的模型更為合理。此時，該模型中：病人平均等待時間為12.2天，平均等待隊長81

4、.18人/天；月平均床位周轉次數(shù)為3.16次/月。一、問題的重述醫(yī)院就醫(yī)排隊是大家都非常熟悉的現(xiàn)象，它以這樣或那樣的形式出現(xiàn)在我們面前，例如，患者到門診就診、到收費處劃價、到藥房取藥、到注射室打針、等待住院等，往往需要排隊等待接受某種服務。 該醫(yī)院眼科門診每天開放，住院部共有病床79張。該醫(yī)院眼科手術主要分四大類：白內障、視網膜疾病、青光眼和外傷。題目中給出了2008年7月13日至2008年9月11日這段時間里各類病人的情況。白內障手術較簡單，而且沒有急癥。目前該院是每周一、三做白內障手術，此類病人的術前準備時間只需1、2天。做兩只眼的病人比做一只眼的要多一些，大約占到60%。如果要做雙眼是周

5、一先做一只，周三再做另一只。外傷疾病通常屬于急癥，病床有空時立即安排住院，住院后第二天便會安排手術。其他眼科疾病比較復雜，有各種不同情況，但大致住院以后2-3天內就可以接受手術，主要是術后的觀察時間較長。這類疾病手術時間可根據(jù)需要安排，一般不安排在周一、周三。由于急癥數(shù)量較少，建模時這些眼科疾病可不考慮急癥。該醫(yī)院眼科手術條件比較充分，在考慮病床安排時可不考慮手術條件的限制，但考慮到手術醫(yī)生的安排問題，通常情況下白內障手術與其他眼科手術（急癥除外）不安排在同一天做。當前該住院部對全體非急癥病人是按照FCFS（First come, First serve）規(guī)則安排住院，但等待住院病人隊列卻越來

6、越長，醫(yī)院方面希望你們能通過數(shù)學建模來幫助解決該住院部的病床合理安排問題，以提高對醫(yī)院資源的有效利用。需要解決的問題如下：問題一：試分析確定合理的評價指標體系，用以評價該問題的病床安排模型的優(yōu)劣。問題二：試就該住院部當前的情況，建立合理的病床安排模型，以根據(jù)已知的第二天擬出院病人數(shù)來確定第二天應該安排哪些病人住院。并對你們的模型利用問題一中的指標體系作出評價。問題三：作為病人，自然希望盡早知道自己大約何時能住院。能否根據(jù)當時住院病人及等待住院病人的統(tǒng)計情況，在病人門診時即告知其大致入住時間區(qū)間。 問題四: 周六和周日不安排手術時間（外傷除外），那么在白內障手術時間不變的情況下，重新確定各類病人

7、入院規(guī)則。在第二問中已經隨機產生門診人數(shù)的基礎上，按照新規(guī)則確定病人入院時間，隨后確定手術時間，并用Matlab軟件隨機模擬產生各個病人的出院時間。問題五： 問題的基本假設和符號系統(tǒng)2.1 問題的基本假設1) 已經門診的患者不會因等待時間過長而改奔其他醫(yī)院；2) 2008-7-13開始統(tǒng)計數(shù)據(jù)時醫(yī)院病滿床；3）第四問中周六、周日可以安排外傷急癥手術，但不安排其他疾病手術；4）不考慮重大突發(fā)事件引起的急病情況。2.2 符號系統(tǒng) ：白內障、視網膜、青光眼、外傷四類病人（=1，2，3，4）；：第類病人的平均等待時間，即病人門診時間到入院時間均值；：第類病人的平均逗留時間，即病人門診時間到出院時間均值

8、；：第類病人的平均住院時間，即病人入院時間到出院時間均值；：第類病人的住院等待手術平均時間，即病人入院時間到手術時間均值；：第類病人的病情平均恢復時間，即病人第一次手術時間到出院時間均值；：病人平均住院天數(shù)。：病床的平均月周轉次數(shù)；：第天第類病人的門診人數(shù)；：第天第類病人的入院人數(shù)；：第天第類病人的等待人數(shù)；：第天第類病人的排隊人數(shù)，即隊長；：第類病人的平均排隊人數(shù)，即平均隊長；三、問題一的解答3.1 模型評價指標的選取以及對模型一（FCFS）的評價一、平均月床位周轉次數(shù)在一定時間內，每張病床住過幾次病人即反映了醫(yī)院治療病人的服務效率。 平均床位周轉次數(shù)越大說明醫(yī)院的醫(yī)病效率越高，則模型越優(yōu)。

9、本文定義： ； (1)其中（病人平均住院天數(shù)）=每人住院天數(shù)求和/總人數(shù)。通過統(tǒng)計篩選后的數(shù)據(jù)得出模型一=9，則月平均床位周轉次數(shù)為3.33次/月。 二、平均等待隊長第天的等待隊長為： （2） 平均排隊長為系統(tǒng)內各天等待服務的病人數(shù)的數(shù)學期望。默認=0；其他具體結果見附錄一。則平均等待隊長為： (3)其中為選定的某段時間（盡量長），該時間長度可根據(jù)等待人數(shù)穩(wěn)定時間段來確定。經過數(shù)據(jù)處理，可得出模型一（FCFS）等待隊長時間的折線圖（如圖3.1所示） 圖3.1 等待隊長隨時間分布折線圖注：圖3.1中，時間軸上0表示2008-7-12.圖3.1中顯示，前10天左右等待人數(shù)呈上升趨勢，并不穩(wěn)定，這是

10、因為題中之前數(shù)據(jù)未給出。在計算平均等待隊長時，選取為12-61天，代入式（3）計算得出為97人/天，即=97人/天。從圖中觀察得等待人數(shù)呈波動上升。故醫(yī)院原系統(tǒng)不可取，即模型一（FCFS）不合理，需要建立新規(guī)則。三、平均病人等待時間對不同病人分別求解，其中病人等待時間過長說明醫(yī)院系統(tǒng)不完善。病人平均逗留時間=平均等待時間+平均住院時間即： （4）統(tǒng)計可得模型一（FCFS）各類指標（如下表3.1所示）表3.1 各類病人的各項時間指標病類指標白內障（單眼）白內障（雙眼）視網膜疾病青光眼外傷平均等待時間（）12.712.712.712.31平均住院時間（）5.28.612.510.57.0平均逗留時

11、間（）17.921.125.122.78.0入院至手術時間（）2.383.632.372.421 表3.1中數(shù)據(jù)結果顯示：除外傷以外，其他四類病人的平均等待時間均為12天以上，而住院時間最多的（視網膜病人）只有12.5天，等待時間占去了各類病人住院時間的50%以上，因此可判定模型一（FCFS）存在較大的缺陷。3.2對醫(yī)院的原系統(tǒng)分析： 平均月床位周轉次數(shù)低。 平均等待隊長隨時間變化波動上升，均值大概為97人。病人的等待時間較長，各類病人住院后平均等待手術時間較長，如：白內障病人入院后1-2天即可做手術，但是統(tǒng)計該指標得白內障單眼病人需要等待2.38天，白內障雙眼病人需要等待3.63天，遠大于1

12、-2天，延長了病人平均逗留時間，造成病床資源浪費。 四、 問題二的解答4.1 問題的分析 以平均月床位周轉次數(shù)最多，平均等待隊長最短，病人消耗時間最少為目標，對原模型“全體非急癥病人是按照FCFS（First come, First serve）規(guī)則安排住院”的排隊模型進行優(yōu)化，建立模型二。運用6SQ統(tǒng)計軟件，對已知病人術后觀察時間進行服從均勻分布的卡方擬合優(yōu)度檢驗，結果認為各類病人的術后觀察時間均服從均勻分布。以此為依據(jù)用matlab軟件模擬尚未出院病人的術后觀察時間，得出第二天擬出院病人及出院病人數(shù)，按手術準備時間最短原則，安排各個病人入院時間、手術時間。同樣產生標準分布隨機數(shù)預測出第三段

13、病人的出院時間，從而做出第三段數(shù)據(jù)病床的合理安排。根據(jù)題中所說手術時間的限制，仍將病人分為白內障（單眼）、白內障（雙眼）、視網膜、青光眼、外傷五類病人。題中所述白內障病人的手術時間為每周一、周三，外傷病人手術時間無限制，其他情況病人不能在周一、周三做手術。結合此規(guī)定需要重新設定入院規(guī)則，并以此安排各類病人入院時間。由于醫(yī)院的病人出入情況是動態(tài)問題，所以在安排第三段病人的入院時間時應該先預測出9月11號以后門診的人數(shù)，這就需要對不同類病人到達門診人數(shù)概率分布進行擬合。利用MATLAB軟件隨機模擬出門診人數(shù)的分布函數(shù)，嚴格按照此函數(shù)預測題中所給第三段病人的入院時間。依照第一問中的評標指標對該模型進

14、行評價，并與模型一作比較，評價模型二的優(yōu)劣。4.2 模型的建立和求解 模型二（病人入院規(guī)則） 根據(jù)“準備時間最短原則”制定模型二具體方案：規(guī)則一：外傷病人當天入院，第二天手術；規(guī)則二：白內障單眼患者一般安排在周一、周二、周六、周日入院；規(guī)則三：白內障雙眼患者一般安排在周六、周日入院；規(guī)則四：視網膜和青光眼患者安排在周二，周三、周四、周五，周日入院；規(guī)則五：同一天同一類病患入院時，采用FCFS服務制；規(guī)則六：當病人等待時間達到20天時，只要有空病床立即安排入院。 根據(jù)以上規(guī)則，通過各類病人術后恢復時間的概率分布函數(shù)（正態(tài)分布），利用Matlab隨機模擬產生其出院時間，填補滿第三段數(shù)據(jù)。 各概率分

15、布函數(shù)的擬合檢驗對已知數(shù)據(jù)，利用6SQ統(tǒng)計軟件進行卡方擬合優(yōu)度檢驗得到如下結果。.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結果接受零假設結果接受零假設結果接受零假設青光眼患者到達時間：外傷患者到達時間：假設檢驗假設檢驗零假設服從泊松分布零假設服從泊松分布

16、自由度3自由度2卡方統(tǒng)計量4卡方統(tǒng)計量1p值0p值1顯著性水平0顯著性水平0結果接受零假設結果接受零假設根據(jù)以上假設檢驗的結果知：。根據(jù)指數(shù)分布與泊松分布的關系：如相繼兩個事件出現(xiàn)的間隔時間服從參數(shù)為的指數(shù)分布，則在單位時間間隔內事件出現(xiàn)的次數(shù)服從參數(shù)為的泊松分布，即單位時間內該事件出現(xiàn)次的概率為 （5）將各類病人平均每天到達的人數(shù)，即值代入式（5），可得相應的概率密度函數(shù)，計算結果見表4.2。表4.2：各類病人達到時間的概率密度函數(shù)和相應的值類別白內障單眼白內障雙眼視網膜青光眼外傷1.62.22.811.2 各類患者的術后觀察時間服從均勻分布首先統(tǒng)計出各類病人的術后觀察時間，根據(jù)統(tǒng)計結果，我

17、們假定各類病人的術后觀察時間服從均勻分布，然后通過6SQ統(tǒng)計軟件進行卡方擬合優(yōu)度檢驗，檢驗結果見表4.3。表4.3：各類患者的術后觀察時間服從均勻分布的卡方擬合優(yōu)度檢驗白內障單眼術后觀察時間白雙術后觀察時間視網膜術后觀察時間假設檢驗假設檢驗假設檢驗零假設服從均勻分布零假設服從均勻分布零假設服從均勻分布自由度34自由度39自由度99卡方統(tǒng)計量7.858300287卡方統(tǒng)計量2.093007186卡方統(tǒng)計量55.04926108p值0.999999107p值1p值0.999897347顯著性水平0.05顯著性水平0.05顯著性水平0.05結果接受零假設結果接受零假設結果接受零假設青光眼術后觀察時間

18、外傷術后觀察時間假設檢驗假設檢驗零假設服從均勻分布零假設服從均勻分布自由度37自由度53卡方統(tǒng)計量11.66013072卡方統(tǒng)計量29.25609756p值0.999978078p值0.996696665顯著性水平0.05顯著性水平0.05結果接受零假設結果接受零假設根據(jù)以上假設檢驗的結果知：各類病人的術后觀察時間均服從均勻分布。 指標計算和模型評價選擇模型一中所用指標，對模型二進行評價?？傻媚Ｐ投牡谝活愔笜巳缦卤恚Ｐ鸵粩?shù)據(jù)作為比較）所示：1） 平均月床位周轉次數(shù)通過matlab擬合住院天數(shù)關于時間的結果，如下圖（圖4.1）： 圖4.1可知在37天以后，系統(tǒng)開始穩(wěn)定，此時的平均住院天數(shù)約為

19、6.3天 選擇模型一中所用指標，對模型二進行評價?？傻媚Ｐ投牡谝活愔笜巳缦卤?表4.4 平均月床位周轉次數(shù)比較標指型模模型一模型二病床月平均周轉率3.33次/月4.76次/月 2）平均等待隊長 經過數(shù)據(jù)處理得出等待人數(shù)時間的折線圖（如圖4.2所示）： 圖4.2 模型二-隊長-時間 上圖顯示模型二的等待人數(shù)在10天后趨于穩(wěn)定，并在平均隊長72上下波動。3）平均病人等待時間表4.5 平均等待時間的比較型模類病白內障（單）白內障（雙）視網膜疾病青光眼外傷模型一 12.712.712.712.31模型二 11.910.611.210.914.3 問題結果分析及討論 根據(jù)改進后的模型二，通過數(shù)據(jù)統(tǒng)計可

20、得到以上各項指標值，并與模型一作比較，可得出結論：1.病床的月平均周轉次數(shù)增至4.45次/月，較模型一增多了病床使用次數(shù)，縮小了病床使用周期，顯然提高了病床的工作效率，醫(yī)院的病床配置數(shù)量得到優(yōu)化。2.平均等待隊長為72.4，較模型一中的平均隊長97大大降低，且可觀察得出隨著時間延長，等待人數(shù)在穩(wěn)定后有下降的趨勢，且下降點對應模型二計算隊長的起始時間，明顯得出模型二相比模型一較優(yōu)。3.各類病人的平均等待時間 11.15天。 綜合以上各項指標分析可得模型二較模型一有很大改進。五、問題三的解答5.1問題分析依靠問題二中的優(yōu)化模型，利用matlab軟件模擬出每位病人的入院時間、第一次手術時間、出院時間

21、，因為每次隨機模擬的結果均不相同，我們對排隊隊列中每位病人的入院時間均模擬5次，從中選擇一個最大值和一個最小值，作為等待病人的入院時間區(qū)間。表7-1：5次模擬后產生的出院時間編號一二三四五大概入院時間（上下限）163636464632008-9-132008-9-14262626262622008-9-122008-9-12381767683762008-9-262008-10-3462626262622008-9-122008-9-12562626262622008-9-122008-9-129983838382832008-10-22008-10-310072727273722008-9-

22、222008-9-2310183838383832008-10-32008-10-310283838383832008-10-32008-10-3（注：具體數(shù)據(jù)見附錄二 ）那么我們就可以根據(jù)以上表格中的數(shù)據(jù)確定出病人的大致入院區(qū)間：即為模擬出入院天數(shù)的最小最大值。六、問題四的解答6.1 問題的分析及模型三的建立若周六和周日不安排手術時間（外傷除外），那么在白內障手術時間不變的情況下，需重新確定各類病人入院規(guī)則。模型三具體方案：規(guī)則一：外傷病人當天入院，第二天手術；規(guī)則二：白內障單眼患者在周一、周二、周六、周日優(yōu)先安排入院；規(guī)則三：白內障雙眼患者在周六、周日優(yōu)先安排入院；規(guī)則四：視網膜和青光眼患

23、者在周三、周四、周五優(yōu)先安排入院；規(guī)則五：同一天同一類病患入院時，采用FCFS服務制；規(guī)則六：當病人等待時間超過20天時，優(yōu)先其安排入院。 規(guī)則七：視網膜和青光眼患者在周六周日不進行手術。在第二問中已經隨機產生門診人數(shù)的基礎上，按照新規(guī)則確定病人入院時間，隨后確定手術時間，并用Matlab軟件隨機模擬產生各個病人的出院時間。6.2 指標計算及模型評價表6.1 第一類指標比較指標模型模型一模型二模型三病床月平均周轉率3.33次/月4.76次/月3.16次/月 與問題二類似，根據(jù)新規(guī)則用matlab模擬出每天的隊長，得折線圖（詳見圖6.1）：圖6.1 模型三隊長關于時間分布 能看出來隊長隨著天數(shù)的

24、增加趨于穩(wěn)定，得其平均81.18人/天 表6.2 平均等待時間的比較 病類指標白內障（單）白內障（雙）視網膜疾病青光眼外傷模型一 平均等待時間（）12.712.712.712.31模型二 平均等待時間（）11.910.611.210.91模型三 平均等待時間（）12.311.911.513.41以上數(shù)據(jù)分析可得，模型三相較與模型一，平均等待時間和每日平均隊長都有了很大的改善，因此此模型可行。七模型的評價8.1 模型的優(yōu)點：1) 從本文所選指標看，該文模型大大降低了醫(yī)院資源浪費現(xiàn)象，同時減短了病人的等待時間，緩解了病人擁擠現(xiàn)狀。從醫(yī)院的長遠發(fā)展角度看是比較樂觀的；2) 較為合理的提出病人入院規(guī)則

25、,該模型有較強的可行性；3）給出了系統(tǒng)穩(wěn)定性的定義，對文中指標的存在合理性給出了證明；4）動態(tài)隨機產生符合模型分布的數(shù)據(jù)，較好地滿足了問題原數(shù)據(jù)的趨勢； 5）采用非線性規(guī)劃模型給出模型的最優(yōu)解，較為合理。8.2 模型的缺點：1）該模型未減少白內障病人的等待時間，仍存在一定的資源浪費現(xiàn)象；2）模型的平均等待人數(shù)仍很多，分析其主要原因為醫(yī)院本身的病床資源欠缺。如果可以增加病床個數(shù)，模型可以進一步優(yōu)化。3）由于用matlab產生模擬隨機數(shù)容易產生奇異數(shù)，導致系統(tǒng)的穩(wěn)定性下降。 參考文獻：1 醫(yī)院排隊論模型，2 排隊論模型，3運籌學教材編寫組，運籌學第三版，清華大學出版社，2005

26、 附錄一：總體病人的到達人數(shù)（沒有分類）時間 門診人數(shù)入院人數(shù)等待人數(shù)2008-7-137072008-7-1491152008-77-1670322008-777-19102622008-7-2091702008-7-2192772008-7-2262812008-7-23161962008-7-24721012008-7-25541022008-7-2648982008-7-27911962008-7-281211972008-7-2959932008-7-3069902008-7-31118932008-8-1

27、117972008-8-2612912008-8-3116962008-8-466962008-8-555962008-8-61491012008-8-71581082008-8-88151012008-8-9419862008-8-1089852008-8-1166852008-8-1272902008-88-1469942008-88-16713902008-8-1785932008-8-181251002008-8-191471072008-8-209101062008-8-21341052008-8-22581022008-8-238

28、18922008-8-24711882008-8-2596912008-8-2673952008-8-2788952008-8-28912922008-8-291110932008-8-30914882008-8-3182942008-9-1116992008-9-2621032008-9-3951072008-9-41591132008-9-513121142008-9-6522972008-9-736942008-9-894992008-9-9751012008-9-10913972008-9-11137103附錄二：前來門診病人大致入院的時間區(qū)間編號第一次第二次第三次第四次第五次入院下限

29、163636464632008-9-132008-9-14262626262622008-9-122008-9-12381767683762008-9-262008-10-3462626262622008-9-122008-9-12562626262622008-9-122008-9-12664636464632008-9-132008-9-14765656565652008-9-152008-9-15881767683742008-9-242008-10-3964636464632008-9-132008-9-141062626262622008-9-122008-9-12116262626

30、2622008-9-122008-9-121262626262622008-9-122008-9-121381767483742008-9-242008-10-31465656565652008-9-152008-9-151562626762622008-9-122008-9-171662626762622008-9-122008-9-171781766983742008-9-192008-10-31864646464642008-9-142008-9-141964647064642008-9-142008-9-202064647064642008-9-142008-9-20216762696

31、9622008-9-122008-9-192265656565652008-9-152008-9-152367686969672008-9-172008-9-192467686969672008-9-172008-9-192565656565652008-9-152008-9-152665656565652008-9-152008-9-152764647064642008-9-142008-9-202865656565652008-9-152008-9-152967696969672008-9-172008-9-193069696969672008-9-172008-9-19316969696

32、9692008-9-192008-9-193264647064642008-9-142008-9-203365666566652008-9-152008-9-163469696969692008-9-192008-9-193565666566652008-9-152008-9-163669696969692008-9-192008-9-193769696969692008-9-192008-9-193864647070702008-9-142008-9-203965666566652008-9-152008-9-164081766983692008-9-192008-10-3416969757

33、5692008-9-192008-9-254269757575692008-9-192008-9-254369757575692008-9-192008-9-254481767683752008-9-252008-10-34564647070702008-9-142008-9-204670647070702008-9-142008-9-204781767676752008-9-252008-10-14875767676762008-9-252008-9-264975767676762008-9-252008-9-265075767676762008-9-252008-9-26517064707

34、0702008-9-142008-9-205270647170702008-9-142008-9-215370647170702008-9-142008-9-215476767676762008-9-262008-9-265570707170702008-9-202008-9-215676767676762008-9-262008-9-265770707170702008-9-202008-9-215865666566652008-9-152008-9-165970707170712008-9-202008-9-216071707170712008-9-202008-9-21617170717

35、0712008-9-202008-9-216277767777772008-9-262008-9-276377767777772008-9-262008-9-276471707170712008-9-202008-9-216578787878782008-9-282008-9-286671707170712008-9-202008-9-216778787878782008-9-282008-9-286865666566662008-9-152008-9-166979797979792008-9-292008-9-297079797979792008-9-292008-9-29716566666

36、6662008-9-152008-9-167271717171712008-9-212008-9-217366666666662008-9-162008-9-167479797979792008-9-292008-9-297566666666662008-9-162008-9-167680808080802008-9-302008-9-307780808080802008-9-302008-9-307880808080802008-9-302008-9-307966666666662008-9-162008-9-168066726666662008-9-162008-9-22818181828

37、1812008-10-12008-10-28266727266662008-9-162008-9-228371717171712008-9-212008-9-218466727272662008-9-162008-9-228566727272722008-9-162008-9-228671717171712008-9-212008-9-218766727272722008-9-162008-9-228881818281812008-10-12008-10-28971717771712008-9-212008-9-279082828282822008-10-22008-10-2918382838

38、2832008-10-22008-10-39283828382832008-10-22008-10-39371717771712008-9-212008-9-279471717771712008-9-212008-9-279583828382832008-10-22008-10-39671717771712008-9-212008-9-279762626262622008-9-122008-9-129871717771712008-9-212008-9-279983838382832008-10-22008-10-310072727273722008-9-222008-9-2310183838

39、383832008-10-32008-10-310283838383832008-10-32008-10-3附錄三：函數(shù)：function f=zhuanhua(data1,data2)if str2num(data1(6)-str2num(data2(6)<0 f1=31*(str2num(data2(6)-str2num(data1(6);else f1=0;endn1=str2num(data1(8:end);n2=str2num(data2(8:end);f=n2-n1+f1; 模型二主程序：clear;clc;%a 在預測新來的病人是作中間變量%b 白內障預測到達人數(shù)%s 視網

40、膜疾病預測到達人數(shù)%w 外傷病人預測到達人數(shù)%q 青光眼病人預測到達人數(shù)%bs 白內障（雙眼）預測到達人數(shù)%chuyuan 每一天的出院人數(shù)%cishu 仿真模擬循環(huán)次數(shù)%da2 入住人數(shù)據(jù)%da3 等待人數(shù)據(jù)%day 第幾天%deng 等待人數(shù)據(jù)錄入后的調用位置% ruzhu 分別表示已入住人數(shù)據(jù)錄入后的原始數(shù)據(jù)%f da3中已經有多少人入住%l 表示新增加病人前da3中的人數(shù)%lo 隊長%m da3 中數(shù)據(jù)轉向da2中時的中轉數(shù)據(jù)%n 當天出院人數(shù)%nn 當天未安排的病床數(shù)%t 從原始數(shù)據(jù)中提取數(shù)據(jù)時的中轉數(shù)據(jù)%time 計算出院時間時的中轉數(shù)據(jù)%week 當天是周幾%dadj 每個人的等

41、待時間%zysj 每個人的住院天數(shù)nn=0;load('F:ruzhu.mat')load('F:deng.mat')%總數(shù)據(jù)轉換for i=1:79 t=char(ruzhu(i,2); da2(i,1)=zhuanhua('2008-7-12',t);%第幾天門診 t=char(ruzhu(i,3); da2(i,2)=zhuanhua('2008-7-12',t);%第幾天入院 t=char(ruzhu(i,4); da2(i,3)=zhuanhua('2008-7-12',t);%第幾天第一次手術 t=ch

42、ar(ruzhu(i,5); if strcmp(t,'/') da2(i,4)=zhuanhua('2008-7-12',t);%若需要，第幾天第二次手術 end t=char(ruzhu(i,1); if strcmp(t,'白內障') da2(i,11)=1; elseif strcmp(t,'白內障(雙眼)') da2(i,11)=2; elseif strcmp(t,'青光眼') da2(i,11)=3; elseif strcmp(t,'視網膜疾病') da2(i,11)=4; else

43、 da2(i,11)=5; endendfor i=1:102 t=char(deng(i,2); da3(i,1)=zhuanhua('2008-7-12',t); da3(i,10)=0; t=char(deng(i,1); if strcmp(t,'白內障') da3(i,11)=1; elseif strcmp(t,'白內障(雙眼)') da3(i,11)=2; elseif strcmp(t,'青光眼') da3(i,11)=3; elseif strcmp(t,'視網膜疾病') da3(i,11)=4; else da3(i,11)=5; endend%預測已入住者的出院時間%1 白內障 2 白內障（雙眼） 3 青光眼 4 視網膜疾病 5 外傷 x=1;day=62;cishu=0;lo(1)=102;f=0;while f<102for i=1:size(da2,1) if da2(i,11)=1 time=unidrnd(3)+1; da2(i,5)=da2(i,3)+time; elseif da2(i,11)=2 time=unidrnd(3)+3; da2(i,5)=da2(i,4)+time; elseif da2(i,11)=3