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

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

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

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

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

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

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

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

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

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

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

12、答4.1 問題的分析 以平均月床位周轉次數(shù)最多，平均等待隊長最短，病人消耗時間最少為目標，對原模型“全體非急癥病人是按照FCFSFirst come, First serve規(guī)那么安排住院的排隊模型進行優(yōu)化，建立模型二。運用6SQ統(tǒng)計軟件，對病人術后觀察時間進行服從均勻分布的卡方擬合優(yōu)度檢驗，結果認為各類病人的術后觀察時間均服從均勻分布。以此為依據(jù)用matlab軟件模擬尚未出院病人的術后觀察時間，得出第二天擬出院病人及出院病人數(shù)，按手術準備時間最短原那么，安排各個病人入院時間、手術時間。同樣產(chǎn)生標準分布隨機數(shù)預測出第三段病人的出院時間，從而做出第三段數(shù)據(jù)病床的合理安排。根據(jù)題中所說手術時間的限

13、制，仍將病人分為白內(nèi)障單眼、白內(nèi)障雙眼、視網(wǎng)膜、青光眼、外傷五類病人。題中所述白內(nèi)障病人的手術時間為每周一、周三，外傷病人手術時間無限制，其他情況病人不能在周一、周三做手術。結合此規(guī)定需要重新設定入院規(guī)那么，并以此安排各類病人入院時間。由于醫(yī)院的病人出入情況是動態(tài)問題，所以在安排第三段病人的入院時間時應該先預測出9月11號以后門診的人數(shù)，這就需要對不同類病人到達門診人數(shù)概率分布進行擬合。利用MATLAB軟件隨機模擬出門診人數(shù)的分布函數(shù)，嚴格按照此函數(shù)預測題中所給第三段病人的入院時間。依照第一問中的評標指標對該模型進行評價，并與模型一作比較，評價模型二的優(yōu)劣。4.2 模型的建立和求解4.2.1

14、模型二病人入院規(guī)那么 根據(jù)“準備時間最短原那么制定模型二具體方案：規(guī)那么一：外傷病人當天入院，第二天手術；規(guī)那么二：白內(nèi)障單眼患者一般安排在周一、周二、周六、周日入院；規(guī)那么三：白內(nèi)障雙眼患者一般安排在周六、周日入院；規(guī)那么四：視網(wǎng)膜和青光眼患者安排在周二，周三、周四、周五，周日入院；規(guī)那么五：同一天同一類病患入院時，采用FCFS效勞制；規(guī)那么六：當病人等待時間到達20天時，只要有空病床立即安排入院。 根據(jù)以上規(guī)那么，通過各類病人術后恢復時間的概率分布函數(shù)正態(tài)分布，利用Matlab隨機模擬產(chǎn)生其出院時間，填補滿第三段數(shù)據(jù)。4.2.2 各概率分布函數(shù)的擬合檢驗對數(shù)據(jù)，利用6SQ統(tǒng)計軟件進行卡方擬

15、合優(yōu)度檢驗得到如下結果。 各類病人的到達時間均符合泊松分布表4.1：各類病人的到達服從泊松分布的卡方擬合優(yōu)度檢驗白內(nèi)障單眼患者到達時間：白內(nèi)障雙眼患者到達時間：視網(wǎng)膜患者到達時間：假設檢驗假設檢驗假設檢驗零假設服從泊松分布零假設服從泊松分布零假設服從泊松分布自由度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值0

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

17、然后通過6SQ統(tǒng)計軟件進行卡方擬合優(yōu)度檢驗，檢驗結果見表4.3。表4.3：各類患者的術后觀察時間服從均勻分布的卡方擬合優(yōu)度檢驗白內(nèi)障單眼術后觀察時間白雙術后觀察時間視網(wǎng)膜術后觀察時間假設檢驗假設檢驗假設檢驗零假設服從均勻分布零假設服從均勻分布零假設服從均勻分布自由度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ù)以上假設檢驗的結果知：各類病人的術后觀察時間均服從均勻分布。4.2.3 指標計算和模型評價選擇模型一中所用指標，對模型二進行評價。可得模型二的第一類指標如下表模型一數(shù)據(jù)作為比較所示：平均月床位周轉次數(shù)通過matlab擬合住院天數(shù)關于時間的結果，如以下圖圖4.1： 圖4.1可知在37天以后，系統(tǒng)開始穩(wěn)定，此時的平均住院天數(shù)約為6.3天 選擇模型一中所用指標，對模型二進行

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

20、月平均周轉次數(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次模擬后產(chǎn)生的出院時間編號一二三四五大概入院時間上下限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-222008-9-2310183838383832008-10-

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

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

24、等待時間的比較 病類指標白內(nèi)障單白內(nèi)障雙視網(wǎng)膜疾病青光眼外傷模型一 平均等待時間12.712.712.712.31模型二 平均等待時間11.910.611.210.91模型三 平均等待時間12.311.911.513.41以上數(shù)據(jù)分析可得，模型三相較與模型一，平均等待時間和每日平均隊長都有了很大的改善，因此此模型可行。七模型的評價8.1 模型的優(yōu)點：從本文所選指標看，該文模型大大降低了醫(yī)院資源浪費現(xiàn)象，同時減短了病人的等待時間，緩解了病人擁擠現(xiàn)狀。從醫(yī)院的長遠開展角度看是比較樂觀的；較為合理的提出病人入院規(guī)那么,該模型有較強的可行性；3給出了系統(tǒng)穩(wěn)定性的定義，對文中指標的存在合理性給出了證明；

25、4動態(tài)隨機產(chǎn)生符合模型分布的數(shù)據(jù)，較好地滿足了問題原數(shù)據(jù)的趨勢； 5采用非線性規(guī)劃模型給出模型的最優(yōu)解，較為合理。8.2 模型的缺點：1該模型未減少白內(nèi)障病人的等待時間，仍存在一定的資源浪費現(xiàn)象；2模型的平均等待人數(shù)仍很多，分析其主要原因為醫(yī)院本身的病床資源欠缺。如果可以增加病床個數(shù)，模型可以進一步優(yōu)化。3由于用matlab產(chǎn)生模擬隨機數(shù)容易產(chǎn)生奇異數(shù)，導致系統(tǒng)的穩(wěn)定性下降。 參考文獻：1醫(yī)院排隊論模型， :/ docin /p-249477.html&pageType=previewUser，2011-8-202排隊論模型， HYPERLINK :/ /lxy

26、/shuxuejm/jiao%20cai/6%B5%DA%C1%F9%D5%C2%20%C5%C5%B6%D3%C2%DB.pdf :/ /lxy/shuxuejm/jiao%20cai/6%B5%DA%C1%F9%D5%C2%20%C5%C5%B6%D3%C2%DB.pdf，3?運籌學?教材編寫組，運籌學第三版，清華大學出版社，2005 附錄一：總體病人的到達人數(shù)沒有分類時間 門診人數(shù)入院人數(shù)等待人數(shù)2008-7-137072008-7-1491152008-77-1670322008-77-18121542

27、008-7-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-1117972008-8-2612912008-8-3116962008-8-466962008-8-555962008-8-61491012008-8-71581082008-8-88151012008-8-9419862008-8

28、-1089852008-8-1166852008-8-1272902008-88-1469942008-88-16713902008-8-1785932008-8-181251002008-8-191471072008-8-209101062008-8-21341052008-8-22581022008-8-23818922008-8-24711882008-8-2596912008-8-2673952008-8-2788952008-8-28912922008-8-291110932008-8-30914882008-8-318294200

29、8-9-1116992008-9-2621032008-9-3951072008-9-41591132008-9-513121142008-9-6522972008-9-736942008-9-894992008-9-9751012008-9-10913972008-9-11137103附錄二：前來門診病人大致入院的時間區(qū)間編號第一次第二次第三次第四次第五次入院下限163636464632008-9-132008-9-14262626262622008-9-122008-9-12381767683762008-9-262008-10-3462626262622008-9-122008-9-12

30、562626262622008-9-122008-9-12664636464632008-9-132008-9-14765656565652008-9-152008-9-15881767683742008-9-242008-10-3964636464632008-9-132008-9-141062626262622008-9-122008-9-121162626262622008-9-122008-9-121262626262622008-9-122008-9-121381767483742008-9-242008-10-31465656565652008-9-152008-9-1515626

31、26762622008-9-122008-9-171662626762622008-9-122008-9-171781766983742008-9-192008-10-31864646464642008-9-142008-9-141964647064642008-9-142008-9-202064647064642008-9-142008-9-202167626969622008-9-122008-9-192265656565652008-9-152008-9-152367686969672008-9-172008-9-192467686969672008-9-172008-9-1925656

32、56565652008-9-152008-9-152665656565652008-9-152008-9-152764647064642008-9-142008-9-202865656565652008-9-152008-9-152967696969672008-9-172008-9-193069696969672008-9-172008-9-193169696969692008-9-192008-9-193264647064642008-9-142008-9-203365666566652008-9-152008-9-163469696969692008-9-192008-9-1935656

33、66566652008-9-152008-9-163669696969692008-9-192008-9-193769696969692008-9-192008-9-193864647070702008-9-142008-9-203965666566652008-9-152008-9-164081766983692008-9-192008-10-34169697575692008-9-192008-9-254269757575692008-9-192008-9-254369757575692008-9-192008-9-254481767683752008-9-252008-10-345646

34、47070702008-9-142008-9-204670647070702008-9-142008-9-204781767676752008-9-252008-10-14875767676762008-9-252008-9-264975767676762008-9-252008-9-265075767676762008-9-252008-9-265170647070702008-9-142008-9-205270647170702008-9-142008-9-215370647170702008-9-142008-9-215476767676762008-9-262008-9-2655707

35、07170702008-9-202008-9-215676767676762008-9-262008-9-265770707170702008-9-202008-9-215865666566652008-9-152008-9-165970707170712008-9-202008-9-216071707170712008-9-202008-9-216171707170712008-9-202008-9-216277767777772008-9-262008-9-276377767777772008-9-262008-9-276471707170712008-9-202008-9-2165787

36、87878782008-9-282008-9-286671707170712008-9-202008-9-216778787878782008-9-282008-9-286865666566662008-9-152008-9-166979797979792008-9-292008-9-297079797979792008-9-292008-9-297165666666662008-9-152008-9-167271717171712008-9-212008-9-217366666666662008-9-162008-9-167479797979792008-9-292008-9-2975666

37、66666662008-9-162008-9-167680808080802008-9-302008-9-307780808080802008-9-302008-9-307880808080802008-9-302008-9-307966666666662008-9-162008-9-168066726666662008-9-162008-9-228181818281812008-10-12008-10-28266727266662008-9-162008-9-228371717171712008-9-212008-9-218466727272662008-9-162008-9-2285667

38、27272722008-9-162008-9-228671717171712008-9-212008-9-218766727272722008-9-162008-9-228881818281812008-10-12008-10-28971717771712008-9-212008-9-279082828282822008-10-22008-10-29183828382832008-10-22008-10-39283828382832008-10-22008-10-39371717771712008-9-212008-9-279471717771712008-9-212008-9-2795838

39、28382832008-10-22008-10-39671717771712008-9-212008-9-279762626262622008-9-122008-9-129871717771712008-9-212008-9-279983838382832008-10-22008-10-310072727273722008-9-222008-9-2310183838383832008-10-32008-10-310283838383832008-10-32008-10-3附錄三：函數(shù)：function f=zhuanhua(data1,data2)if str2num(data1(6)-str

40、2num(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 白內(nèi)障預測到達人數(shù)%s 視網(wǎng)膜疾病預測到達人數(shù)%w 外傷病人預測到達人數(shù)%q 青光眼病人預測到達人數(shù)%bs 白內(nèi)障雙眼預測到達人數(shù)%chuyuan 每一天的出院人數(shù)%cishu 仿真模擬循環(huán)次數(shù)%da2 入住人數(shù)據(jù)%da3 等待人數(shù)據(jù)%day 第幾天%deng 等待人

41、數(shù)據(jù)錄入后的調(diào)用位置% ruzhu 分別表示已入住人數(shù)據(jù)錄入后的原始數(shù)據(jù)%f da3中已經(jīng)有多少人入住%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 每個人的等待時間%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);

42、%第幾天門診 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=char(ruzhu(i,5); if strcmp(t,/) da2(i,4)=zhuanhua(2008-7-12,t);%假設需要，第幾天第二次手術 end t=char(ruzhu(i,1); if strcmp(t,白內(nèi)障) da2(i,11)=1; elseif strcmp(t,白內(nèi)障(雙眼) da2(i,11)=2; elseif

43、 strcmp(t,青光眼) da2(i,11)=3; elseif strcmp(t,視網(wǎng)膜疾病) da2(i,11)=4; else 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,白內(nèi)障) da3(i,11)=1; elseif strcmp(t,白內(nèi)障(雙眼) da3(i,11)=2; elseif strcmp(t,青光眼) da3(i,11)=3; elseif strcmp(t,視網(wǎng)

44、膜疾病) da3(i,11)=4; else da3(i,11)=5; endend%預測已入住者的出院時間%1 白內(nèi)障 2 白內(nèi)障雙眼 3 青光眼 4 視網(wǎng)膜疾病 5 外傷 x=1;day=62;cishu=0;lo(1)=102;f=0;while f102for 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 time=unidr

45、nd(9)+3; da2(i,5)=da2(i,3)+time; elseif da2(i,11)=4 time=unidrnd(11)+4; da2(i,5)=da2(i,3)+time; elseif da2(i,11)=5 time=unidrnd(8)+2; da2(i,5)=da2(i,3)+time; end if da2(i,5)0 if da3(i,11)=5&da3(i,10)=1%不管哪一天首先篩選急診病人 da3(i,2)=day; da3(i,3)=da3(i,2)+1; time=unidrnd(8)+2; da3(i,5)=da3(i,3)+time; da3(i,

46、10)=1; n=n-1; end endendfor i=1:size(da3,1) if n0 if (day-da3(i,1)20&da3(i,10)=1%等待時間超過20天時作急診處理 da3(i,2)=day; if da3(i,11)=1 da3(i,3)=da3(i,2)+1; elseif da3(i,11)=2 da3(i,3)=da3(i,2)+1; da3(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+2; end if da3(i,11)=1 time=unidrnd(3)+1; elseif da3(i,11)=2 time=unid

47、rnd(3)+3; elseif da3(i,11)=3; time=unidrnd(9)+3; elseif da3(i,11)=4 time=unidrnd(11)+4; end %當天是周幾影響其等待時間 if week=1 if da3(i,11)=1 da3(i,3)=da3(i,2)+2; elseif da3(i,11)=2 da3(i,3)=da3(i,2)+7; da3(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+3; end elseif week=2 if da3(i,11)=1 da3(i,3)=da3(i,2)+1; elseif

48、da3(i,11)=2 da3(i,3)=da3(i,2)+6; da3(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+2; end elseif week=3 if da3(i,11)=1 da3(i,3)=da3(i,2)+5; elseif da3(i,11)=2 da3(i,3)=da3(i,2)+5; da3(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+2; end elseif week=4 if da3(i,11)=1 da3(i,3)=da3(i,2)+4; elseif da3(i,11)=2 da3(i,

49、3)=da3(i,2)+4; da3(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+2; end elseif week=5 if da3(i,11)=1 da3(i,3)=da3(i,2)+3; elseif da3(i,11)=2 da3(i,3)=da3(i,2)+3; da3(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+2; end elseif week=6 if da3(i,11)=1 da3(i,3)=da3(i,2)+2; elseif da3(i,11)=2 da3(i,3)=da3(i,2)+2; da3

50、(i,4)=da3(i,3)+2; else da3(i,3)=da3(i,2)+3; end end da3(i,5)=da3(i,3)+time; da3(i,10)=1; n=n-1; end endendweek=mod(day-1,7);if week=0%如果是周日 for i=1:size(da3,1) %首先篩選白內(nèi)障雙眼患者 if n0 if da3(i,11)=2&da3(i,10)=1 da3(i,2)=day; da3(i,10)=1; da3(i,3)=da3(i,2)+1; da3(i,4)=da3(i,2)+3; time=unidrnd(3)+3; da3(i,

51、5)=da3(i,4)+time; n=n-1; end end end for i=1:size(da3,1)%然后篩選白內(nèi)障患者 if n0 if da3(i,11)=1&da3(i,10)=1 da3(i,2)=day; da3(i,3)=da3(i,2)+1; time=unidrnd(3)+1; da3(i,5)=da3(i,3)+time; da3(i,10)=1; n=n-1; end end end for i=1:size(da3,1)%最后篩選其余患者 if n0 if da3(i,10)=0 da3(i,2)=day; da3(i,3)=da3(i,2)+2; if da

52、3(i,11)=3 time=unidrnd(9)+3; else time=unidrnd(11)+4; end da3(i,5)=da3(i,3)+time; da3(i,10)=1; n=n-1; end end endelseif week=6%如果是周六 for i=1:size(da3,1) %首先篩選白內(nèi)障雙眼患者 if n0 if da3(i,11)=2&da3(i,10)=1 da3(i,2)=day; da3(i,10)=1; da3(i,3)=da3(i,2)+2; da3(i,4)=da3(i,2)+4; time=unidrnd(3)+3; da3(i,5)=da3(

53、i,4)+time; n=n-1; end end end for i=1:size(da3,1)%隨后篩選白內(nèi)障患者 if n0 if da3(i,11)=1&da3(i,10)=1 da3(i,2)=day; da3(i,3)=da3(i,2)+2; time=unidrnd(3)+1; da3(i,5)=da3(i,3)+time; da3(i,10)=1; n=n-1; end end end for i=1:size(da3,1)%最后篩選其余患者 if n0 if da3(i,10)=0 da3(i,2)=day; da3(i,3)=da3(i,2)+3; da3(i,10)=1; if da3(i,11)=3 time=unidrnd(9)+3; else time=unidrnd(11)+4; end da3(i,5)=da3(i,3)+time; n=n-1; end end end elseif week=2%如果是周二 for i=1:size(da3,1) %首先篩選白內(nèi)障患者 if n0 if da3(i,11)=1&da3(i,10)=1 da3(i,2)=day; da