【正文】 lon=[lon,long]。 for i=1:length(index) for j=1:11 tt{n2+i,j}=t3{index(i),j}。 t3{nn+i,9}=1。外傷39。 t3{nn+i,10}=1。t3{nn+i,2}=f。 t3{nn+i,11}=11。 t3{nn+i,7}=mod(t3{nn+i,2}1,7)。 end nn=size(t3,1)。 t3{nn+i,8}=1。 for i=1:bs t3{nn+i,1}=39。 t3{nn+i,9}=1。白內(nèi)障39。w=poissrnd(1)。 %產(chǎn)生新的數(shù)據(jù)，補(bǔ)充隊(duì)列 b=poissrnd()。 t3{i,9}=mod(t3{i,4}1,7)。青光眼39。 t3{i,11}=0。視網(wǎng)膜疾病39。 end end elseif mod(f1+2,7)~=1amp。 end t3{i,8}=mod(t3{i,3}1,7)。 if strcmp(t3{i,1},39。n0 index=[index,i]。)|strcmp(t3{i,1},39。 t3{i,10}=mod(t3{i,6}1,7)。) t3{i,6}=f+1+unidrnd(3)+1。 t3{i,3}=f。))amp。 end end elseif mod(f1,7)==0%這天是周日時： for i=1:size(t3,1)%此時確定是安排白內(nèi)障還是（青光眼、視網(wǎng)膜） if (strcmp(t3{i,1},39。 t3{i,8}=mod(t3{i,3}1,7)。 t3{i,11}=0。白內(nèi)障39。 t3{i,9}=mod(t3{i,4}1,7)。 t3{i,3}=f。)amp。 t3{i,10}=mod(t3{i,6}1,7)。 t3{i,4}=f+1。t3{i,11}~=0amp。 end end long=longn。%存儲需要記錄的下標(biāo)號。long=102。endtt=t2。 t3{i,8}=1。 t2{i,9}=mod(t2{i,4}1,7)。 t1{i,10}=mod(t1{i,6}1,7)。 end if t2{i,6}62 t2{i,6}=62。 elseif strcmp(t2{i,1},39。 elseif strcmp(t2{i,1},39。 elseif strcmp(t2{i,1},39。 elseif strcmp(t2{i,1},39。end%計(jì)算機(jī)模擬產(chǎn)生對于t2中的79個人的出院情況%已知到第61天都還沒有出院，所以最早的出院時間是第62天for i=1:size(t2,1) if strcmp(t2{i,1},39。,t2{i,5})。 if ~strcmp(t2{i,5},39。,t2{i,3})。200871239。) t1{i,5}=jsdata(39。200871239。 t1{i,4}=jsdata(39。,t1{i,2})。（2）模擬的程序：clear。else f1=0。使所建模型能準(zhǔn)確的預(yù)測出病人的入院和出院時間。10. 模型的評價(jià)、改進(jìn)及推廣優(yōu)點(diǎn)：（1）根據(jù)我們定下的安排病人入院規(guī)則，建立的模型在一定程度上縮短了病人排隊(duì)的隊(duì)長，因?yàn)樵Ｐ偷娜齻€指標(biāo)，；我們建的模型三個指標(biāo)：，；%，%，%；（2）根據(jù)模型可推算出當(dāng)前病人的出院時間，故我們把表二出院時間的填充了；（3）利用我們建立的模型二，可根據(jù)第二天擬出院的病人確定病人入院的最佳時間；缺點(diǎn)：由于所給數(shù)據(jù)太少以致在統(tǒng)計(jì)數(shù)據(jù)時不是很準(zhǔn)確，又由于計(jì)算機(jī)模擬帶有一定的隨機(jī)性，以致得到模型的三個指標(biāo)不是很讓人滿意。那么 設(shè)病人在在系統(tǒng)內(nèi)的平均逗留時間為，則其中表示患第類病的第個病人的等待入院時間，表示患第類病的第個病人的住院時間，表示一段時間內(nèi)患第類病的人數(shù)）所以，目標(biāo)函數(shù)為因各類病人都分配了相應(yīng)的病床，故每一類病人在可在任意一天住院且是先來先服務(wù)的，這一點(diǎn)和模型一，二不一樣，故約束條件只有：(1) 各類病人的術(shù)前準(zhǔn)備時間(2) 白內(nèi)障的手術(shù)只能安排在周一 三。然后作出隊(duì)列長度的隨時間的變化圖示：如下所示：圖81:隊(duì)長隨時間的變化趨勢此時，我們同時發(fā)現(xiàn):當(dāng)我們忽略外傷病人的等待時間時，隨著時間的推移，病人的等待時間的規(guī)律如下圖所示：圖82:病人的等待時間隨時間的變化趨勢此時，在我們的分配原則下，隊(duì)列的長度不斷減小，并且有趨于穩(wěn)定的趨勢，我們可以認(rèn)為此時的穩(wěn)態(tài)等待人數(shù)是60，而且，病人從就診到第一次手術(shù)的時間也有一個逐步下降的趨勢，所以因此我們認(rèn)為這種分配的方案是比較合理的：據(jù)此我們可以在這種條件下求得的三個指標(biāo)值：現(xiàn)在，我們同樣取得從9月12號后的20天的一個模擬的入院安排情況，然后計(jì)算得我們定義的三個目標(biāo)函數(shù)的最優(yōu)解為：手術(shù)前的平均逗留時間：病人平均術(shù)前準(zhǔn)備時間：；平均每天出院人數(shù)：在這種分配方案下，我們發(fā)現(xiàn)與原模型的手術(shù)前的平均逗留時間：，平均術(shù)前準(zhǔn)備時間：，平均每天出院人數(shù)：，并沒有太大的差別，所以可以同樣按照問題二中的解決辦法，我們會發(fā)現(xiàn)我們的三個評價(jià)指標(biāo)值為：手術(shù)前的平均逗留時間：病人平均術(shù)前準(zhǔn)備時間：；平均每天出院人數(shù)：此時我們就可以發(fā)現(xiàn)，當(dāng)這個排隊(duì)系統(tǒng)在盡量趨于穩(wěn)定狀態(tài)時，它的各個指標(biāo)均比前邊的結(jié)果有了一定的優(yōu)化，從而進(jìn)一步證明我們的排隊(duì)系統(tǒng)比原有的效率更高。（同模型一的目標(biāo)函數(shù)）以三個評價(jià)指標(biāo)最優(yōu)為目標(biāo)函數(shù)：該模型的約束條件除了包含模型一的約束條件外，還有以下幾個：由于白內(nèi)障手術(shù)之后安排在周一和周三，外傷手術(shù)每天都可以安排，所以周六和周日不安排手術(shù)只會影響視網(wǎng)膜和青光眼的手術(shù)安排。66,71,72 ,73 ,789月16日,9月21日,9月22日,9月23日,9月28日10139。我們?nèi)∧M的次數(shù)為10，下面是我們所得的一個近似的結(jié)果：表71：10次模擬后產(chǎn)生的10個模擬的出院時間編號生病類別一二三四五六七八九十1白內(nèi)障(雙眼)636363636363636363632視網(wǎng)膜疾病626262626262626262623青光眼626262626262626262624視網(wǎng)膜疾病626262626262626262625視網(wǎng)膜疾病62626262626262626262… ……………………………………99視網(wǎng)膜疾病76747469747474757469100白內(nèi)障66737272727872717172101視網(wǎng)膜疾病76757469747474757469102視網(wǎng)膜疾病76757474747474757474（注：具體數(shù)據(jù)見附錄一 ）那么我們就可以根據(jù)以上表格中的數(shù)據(jù)確定出病人的大致入院區(qū)間：表72：病人的大致入院區(qū)間編號生病類別最佳入院時間對應(yīng)的日期1白內(nèi)障(雙眼)639132視網(wǎng)膜疾病629123青光眼629124視網(wǎng)膜疾629125視網(wǎng)膜疾病62912…………………9939。然后計(jì)算得我們定義的三個目標(biāo)函數(shù)的最優(yōu)解為：病人手術(shù)前的平均逗留時間：（越小越好）病人平均術(shù)前準(zhǔn)備時間：；（越小越好）平均每天出院人數(shù)：9人（越大越好）在這種分配方案下，我們發(fā)現(xiàn)與原模型的手術(shù)前的平均逗留時間：，平均術(shù)前準(zhǔn)備時間：，平均每天出院人數(shù)：，并沒有太大的差別，我們認(rèn)為這是在隊(duì)列沒有進(jìn)入穩(wěn)定狀態(tài)時的統(tǒng)計(jì)數(shù)據(jù)造成的，我們用以下的方法處理：我們可以認(rèn)為當(dāng)系統(tǒng)服務(wù)了100個顧客后，它已經(jīng)進(jìn)入了穩(wěn)定狀態(tài)，又由于我們要去足夠的數(shù)據(jù)才能具有說服力，因此我們定義一個評價(jià)區(qū)間：從第100個排隊(duì)等待手術(shù)的病人到第30天結(jié)束時最后一個出院的病人。102號以后是[]，表明該行數(shù)據(jù)是我們自己通過possion流對隊(duì)列追加顧客后的效果，是新加入排隊(duì)系統(tǒng)的，方便進(jìn)行連續(xù)模擬。/39。/39。/39。并且可以按照每天到來的各類病人的possion流來補(bǔ)充每天的排隊(duì)隊(duì)列人數(shù)。該模型是為了解決醫(yī)院的病床安排問題，為了使病床的安排更加合理，我們只需使三個評價(jià)指標(biāo)——手術(shù)前的平均逗留時間最小，平均每天出院人數(shù)最大，病人平均術(shù)前的準(zhǔn)備時間最小即可，所以我們建立了如下的目標(biāo)函數(shù)：ⅰ由于白內(nèi)障手術(shù)比較簡單，此類病人的術(shù)前準(zhǔn)備時間只需2天，而且根據(jù)附表一知，各類病人的術(shù)前準(zhǔn)備時間均在17天之內(nèi)（包括1天和7天），所以白內(nèi)障單眼和雙眼病人術(shù)前準(zhǔn)備時間為17天，即ⅱ由于視網(wǎng)膜和青光眼疾病比較復(fù)雜，大致住院以后23天內(nèi)就可以接受手術(shù)，而且根據(jù)附表一知，各類病人的術(shù)前準(zhǔn)備時間均在17天之內(nèi)（包括1天和7天），所以這兩類疾病的術(shù)前準(zhǔn)備時間為27天，即ⅲ外傷疾病有空床時立即安排住院，且住院后第二天便會安排手術(shù)，所以此類病的術(shù)前準(zhǔn)備時間為1天，即ⅳ根據(jù)數(shù)據(jù)分析結(jié)果，為了縮短等待入院病人的隊(duì)列長度，我們制定了如下的病床安排原則：原則1白內(nèi)障單眼患者一般安排在周一、周二、周六、周日入院原則2白內(nèi)障雙眼患者一般安排在周六、周日入院原則3視網(wǎng)膜和青光眼患者安排在周三、周四、周五入院原則4外傷病人當(dāng)天入院，第二天手術(shù)原則5當(dāng)病人等待時間達(dá)到25天時，只要有空病床立即安排入院說明：外傷病人在任意一天優(yōu)先級是最高的，在某一天病人可安排住院此時此類?。ǔ鈧┤说膬?yōu)先級第二高，且若有兩類或兩類以上的病（除外傷）人都可在同一天入院，則這些病人優(yōu)先級第二高且相等，在這天不安排住院的病人的優(yōu)先級最低。首先統(tǒng)計(jì)出外傷病人的住院時間，根據(jù)統(tǒng)計(jì)結(jié)果，我們假定其服從泊松分布，然后通過6SQ統(tǒng)計(jì)軟件進(jìn)行卡方擬合優(yōu)度檢驗(yàn)，.: 外傷病人住院時間服從泊松分布的卡方擬合優(yōu)度檢驗(yàn)結(jié)果假設(shè)檢驗(yàn)零假設(shè)服從泊松分布自由度6卡方統(tǒng)計(jì)量p值顯著性水平結(jié)果接受零假設(shè)5． 問題一的解答本文研究的是某醫(yī)院眼科病床的合理安排的數(shù)學(xué)建模問題，對于病床安排模型的優(yōu)劣，不能憑人們的主觀感受進(jìn)行判斷，而要確定合理的評價(jià)指標(biāo)體系進(jìn)行判斷，為此我們確定了如下的評價(jià)指標(biāo)體系：指標(biāo)1：手術(shù)前的平均逗留時間手術(shù)前的平均逗留時間指門診到第一次手術(shù)的平均時間，其數(shù)學(xué)表達(dá)式為該指標(biāo)值越小，表示病床安排模型越好。（病人到達(dá)時間指病人的門診時間）服從泊松分布由于醫(yī)院就醫(yī)排隊(duì)是典型的排隊(duì)論問題，而一般的排隊(duì)論模型都是泊松輸入，所以我們先假定病人的到達(dá)服從泊松分布，然后根據(jù)附錄一給出的數(shù)據(jù)求出每天到門診看病的各類病人的人數(shù)（統(tǒng)計(jì)結(jié)果見附錄二）,再利用6SQ統(tǒng)計(jì)軟件進(jìn)行卡方擬合優(yōu)度檢驗(yàn)得到如下結(jié)果（）：：各類病人的到達(dá)服從泊松分布的卡方擬合優(yōu)度檢驗(yàn)白內(nèi)障單眼患者到達(dá)時間：白內(nèi)障雙眼患者到達(dá)時間：視網(wǎng)膜患者到達(dá)時間：假設(shè)檢驗(yàn)假設(shè)檢驗(yàn)