數(shù)學建模競賽承諾書

1、2014高教社杯全國大學生數(shù)學建模競賽承 諾 書我們仔細閱讀了全國大學生數(shù)學建模競賽章程和全國大學生數(shù)學建模競賽參賽規(guī)則（以下簡稱為“競賽章程和參賽規(guī)則”，可從全國大學生數(shù)學建模競賽網(wǎng)站下載）。我們完全明白，在競賽開始后參賽隊員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊外的任何人（包括指導教師）研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反競賽章程和參賽規(guī)則的，如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴格遵守競賽章程和參賽規(guī)則，以保證競賽的公正、公平性。如有違反競賽章程和參賽

2、規(guī)則的行為，我們將受到嚴肅處理。我們授權(quán)全國大學生數(shù)學建模競賽組委會，可將我們的論文以任何形式進行公開展示（包括進行網(wǎng)上公示，在書籍、期刊和其他媒體進行正式或非正式發(fā)表等）。我們參賽選擇的題號是（從a/b/c/d中選擇一項填寫）： d 我們的報名參賽隊號為（8位數(shù)字組成的編號）： 23068006 所屬學校（請?zhí)顚懲暾娜?四川建筑職業(yè)技術(shù)學院 參賽隊員 (打印并簽名) ：1. 王磊 2. 蔡姍姍 3. 蔣國輝 指導教師或指導教師組負責人 (打印并簽名)： 黃磊 （論文紙質(zhì)版與電子版中的以上信息必須一致，只是電子版中無需簽名。以上內(nèi)容請仔細核對，提交后將不再允許做任何修改。如填寫錯誤，論

3、文可能被取消評獎資格。） 日期： 2014 年 9 月 15 日賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）：2014高教社杯全國大學生數(shù)學建模競賽編 號 專 用 頁賽區(qū)評閱編號（由賽區(qū)組委會評閱前進行編號）：賽區(qū)評閱記錄（可供賽區(qū)評閱時使用）：評閱人評分備注全國統(tǒng)一編號（由賽區(qū)組委會送交全國前編號）：全國評閱編號（由全國組委會評閱前進行編號）：儲藥柜的設(shè)計摘要本文針對自動補藥藥柜的設(shè)計進行研究。針對問題一，在只考慮儲藥柜豎向隔板的最小間距種類，在滿足安全送藥的四個條件，即側(cè)間距2mm,無并排，無側(cè)翻，無水平旋轉(zhuǎn)下，建立單目標優(yōu)化模型，并設(shè)計區(qū)間無重疊聚類算法，實現(xiàn)最少間距種類的求解，由程序得到

4、最少四類列寬的分類，分別為19mm，34mm，46mm，58mm針對問題二，我們將總寬度冗余，與列間距類型數(shù)量作為目標，建立雙目標規(guī)劃模型?；诜謱忧蠼舛嗄繕艘?guī)劃模型方法，我們在問題一得到的4個不同類型的基礎(chǔ)上，首先建立冗余權(quán)重模型，首先計算出各中藥盒寬度在原始4種分類基礎(chǔ)上的加權(quán)冗余，并按照其加權(quán)冗余累積貢獻率排序，根據(jù)累積貢獻率，我們討論了90%和95%下，根據(jù)列寬優(yōu)化算法，計算出新的列寬分類，經(jīng)過加權(quán)冗余度和列寬類數(shù)的分析，我們確定在新增3類情況下的最優(yōu)解。列寬分別為19mm，22mm，34mm，37mm，46mm，47mm，58mm并且給出相應的藥盒編號。針對問題三，我們將平面總?cè)哂喽?/p>

5、，與行間距類型最小作為目標，在以藥柜給定規(guī)格為約束條件下，建立雙目標規(guī)劃模型。在問題二的基礎(chǔ)上，我們通過對分布分析法，先按照比列均衡的思想確定藥柜一行放置76個藥槽，在此基礎(chǔ)上為了盡量減少平面冗余，我們按照高相近歸類方法，得到藥柜至少需要26行，并且計算出高大致需要以下9類。34mm，41mm，47mm，54mm，60mm，72mm，85mm，101mm，125mm針對問題四，在藥槽長度1.5米的條件下，我們首先計算出每一種藥盒在藥槽長度方向上能的個數(shù)。因此確定同一種要需要的藥槽數(shù)量。又因為每天僅集中補藥一次，所以設(shè)計的儲藥槽個數(shù)一次性能放藥盒的個數(shù)大于該需求量的最大值才能滿足。 關(guān)鍵詞：雙目

6、標規(guī)劃 區(qū)間無重疊聚類 分層法1 問題重述藥柜的結(jié)構(gòu)與書柜相似，若干個橫向隔板和豎向隔板將儲藥柜分割成若干個儲藥槽，橫向隔板決定所放藥品的高度，豎向隔板決定所放藥品的寬度，為了方便使用和保證藥品分揀的準確率，防止發(fā)藥錯誤，一個儲藥槽內(nèi)只能擺放同一種藥品，要求藥盒與兩側(cè)豎向隔板之間、與上下兩層橫向隔板之間應留2mm的間隙，同時還要求藥盒在儲藥槽內(nèi)推送過程中不會出現(xiàn)并排重疊、側(cè)翻或水平旋轉(zhuǎn)。為了更好的在實際中運用，在忽略橫向和豎向隔板厚度的情況下，建立數(shù)學優(yōu)化模型，給出下面幾個問題的解決方案。問題一：因為藥盒尺寸規(guī)格差異較大，根據(jù)提供的數(shù)據(jù)，設(shè)計藥柜的豎向隔板間距類型最小種類數(shù)的數(shù)量和每種類型所對

7、應的藥盒規(guī)格。問題二：寬度冗余是藥盒與兩側(cè)豎向隔板之間的間隙超出2mm的部分，當適當增加豎向隔板間距類型的數(shù)量可以減少寬度冗余，但增加豎向隔板間距類型會增加儲藥柜的加工成本，通過對問題一中的最佳設(shè)計求解方案，設(shè)計出合理的豎向隔板間距類型數(shù)量以及每種類型對應的藥盒編號，使得總寬度冗余盡可能小，同時也希望間距的類型數(shù)量盡可能少。問題三：為了考慮拿藥的方便性和補藥的便利性，儲藥柜的尺寸要具有合理性和可行性，規(guī)定儲藥柜的寬度不超過2.5m，高度不超過2m，儲藥柜允許的最大有效高度為1.5m。藥盒與兩層橫向隔板之間的間隙超過2mm的部分叫做高度冗余，可以得出平面冗余=高度冗余寬度冗余，在問題二中計算結(jié)果

8、的基礎(chǔ)上，確定儲藥柜橫向隔板間距的類型數(shù)量，使得儲藥柜的總平面余量盡可能地小，且橫向隔板間距的類型數(shù)量也盡可能地少。問題四：由附件2可得每一種藥品編號對應的最大日需求量。已知儲藥槽的寬度不超過2.5m，有效高度不超過1.5m，長度為1.5m，每天補藥僅一次，請計算每一種藥品需要的儲藥槽個數(shù)。為了保證藥房儲藥滿足需求，計算稀少需要多少個儲藥柜。2 模型假設(shè)1. 假設(shè)每次從后端放入的藥品都正立平穩(wěn)放入；2. 假設(shè)藥盒水平旋轉(zhuǎn)時中心點在一條直線上；3. 假設(shè)藥盒旋轉(zhuǎn)角度超過90時才為水平旋轉(zhuǎn)；4. 假設(shè)每個藥槽都有藥盒放入；5. 假設(shè)一天中僅有的一次藥品補給是在藥店下班前或者下班后一次性補給完成；3

9、 變量說明：表示豎向隔板的間距；：表示第個型號藥盒的長度；：表示第個型號藥盒的寬度；：表示第個型號藥盒的高度；：表示豎向隔板間距類型數(shù)；：表示定義域為的有效區(qū)間集合；：表示冗余權(quán)重系數(shù)；：表示冗余率；：表示各藥盒尺寸的頻數(shù)；：表示藥盒尺寸出現(xiàn)的總頻次；：表示列寬冗余度；：表示第個藥盒放入的藥槽可能高度；：表示每個藥盒放入對應藥槽時的寬度冗余；:表示每個藥盒放入對應藥槽時的寬度冗余；:表示第個藥盒產(chǎn)生的平面冗余；:表示一排橫向同類型列的藥槽個數(shù)；:表示同類型高的在一列縱向中藥槽個數(shù)；4 模型的建立與求解4.1 問題一的模型建立與求解4.1.1 基于藥盒安全推出下的最小列分類規(guī)劃模型模型分析藥盒為

10、了能順利推送過程，且不出現(xiàn)并排重疊、側(cè)翻或水平旋轉(zhuǎn)， 因此每個藥盒存在對應的藥槽列寬區(qū)間，下面根據(jù)三個條件我們分別討論。 設(shè)表示藥槽的寬度，,表示藥盒的長，寬，高。 (a) 順利推出根據(jù)題目，每個藥盒與左右兩側(cè)間距2mm，因此，為了能順利推出，則 (1)(b) 無并排現(xiàn)象為了不發(fā)生并排現(xiàn)象，則藥槽的寬度不能大于兩倍藥盒的寬度，如圖1所示，因此，藥槽寬度和藥盒寬度之間滿足： (2)(c) 無側(cè)翻現(xiàn)象如圖2所示，我們定義當藥盒在藥槽內(nèi)側(cè)倒至藥槽內(nèi)時稱為側(cè)翻?？紤]藥盒在側(cè)翻過程中，橫向最大距離為寬與高的對角線，因此，為避免側(cè)翻，則藥槽的寬度應該大于此對角線長度，即： (3) (d) 無水平旋轉(zhuǎn)現(xiàn)象如

11、圖3所示，我們定義當藥盒在藥槽內(nèi)平面旋轉(zhuǎn)90度時為水平旋轉(zhuǎn)。同樣考慮藥盒在水平旋轉(zhuǎn)過程中，橫向最大距離為寬與長的對角線，因此，為避免水平旋轉(zhuǎn)，則藥槽的寬度應該大于此對角線長度，即： (4）圖1圖2圖3模型建立問題一即在滿足以上四條的情況下，計算最小列寬分類，我們設(shè)為藥槽按照列寬的分類數(shù)，根據(jù)模型分析，建立最小列寬分類規(guī)劃模型: (5) () （6）4.1.2 基于區(qū)間無重疊下最小聚類法的模型求解要求解4.1.1中的模型，我們首先根據(jù)約束條件，根據(jù)附表1中的原始數(shù)據(jù)算出每個藥盒的寬度容許區(qū)間。以每個藥盒尺寸的第一個區(qū)間為例，如下表：（程序見附錄1）表1.藥盒的列寬容許度區(qū)間表寬度有效區(qū)間寬度有效

12、區(qū)間寬度有效區(qū)間寬度有效區(qū)間寬度有效區(qū)間1012,20)2022,40)3032,60)4042,80)5052,82)1113,22)2123,42)3133,43.8)4143,58)5153,72.8)1214,24)2224,38.8)3234,51.2)4244,60)5254,73.5)為了將所有藥盒用最少的藥槽規(guī)格來放置，等價于對個區(qū)間，我們尋找個區(qū)間，使得這個區(qū)間之間交集為空，且這個區(qū)間每一個都至少完全屬于原始區(qū)間中的某一個。區(qū)間無重疊聚類算法第一步，將所有區(qū)間，按照區(qū)間下限從小到大排列，記排序好的區(qū)間為第二步，從第一個區(qū)間開始,比對的上限與的下限，若當時，則為單獨一類，并從分

13、析中暫時剔除，對繼續(xù)聚類當時，交集為，則將代替,，對繼續(xù)聚類。第三步，重復第二步，直至所有原始區(qū)間聚類完畢。聚類過程示意圖如圖4所示：圖4區(qū)間無重疊聚類示意圖按照上述聚類過程，顯然最后得到的個區(qū)間相互間無重疊，且每個區(qū)間一定完全包含于原始的某個區(qū)間中，即所有的藥盒一定能放置于某個聚類后的藥槽中,根據(jù)算法編制程序見附錄2。最后聚類得到的區(qū)間見表2表2：規(guī)格聚類分析表規(guī)格種類序號聚類區(qū)間119,20）234,34.9857113690718)346,46.6690475583121)458,63.6396103067893)考慮到空間節(jié)約，藥槽寬度取每個區(qū)間下限即可，即最少的藥槽列寬規(guī)格為19mm

14、,34mm，46mm，58mm每個規(guī)格下可放置的藥品規(guī)格如下表:表3：藥盒寬度的最佳取值表藥槽規(guī)格d/mm對應的藥盒規(guī)格寬度最佳取值范圍d/mm191017341832463344584556所以綜上所述，豎向隔板間距類型最少的儲藥規(guī)格有4種，每種類型所對應的藥盒規(guī)格如上表。模型二的建立(列寬冗余雙目標規(guī)劃)4.2.1模型分析與建立列寬對于超過規(guī)定尺寸以外的距離，稱為寬度冗余。根據(jù)題目，列寬冗余的計算公式為： （7）整個藥槽的總列寬冗余為問題二即在列分類較少的情況下使得總列寬冗余最小。將列寬冗余的最小作為目標，我們在模型一的基礎(chǔ)上建立列寬冗余雙目標規(guī)劃模型如下： （8）4.2.2 基于分層法求

15、解雙目標規(guī)劃由于模型二是雙目標規(guī)劃，我們根據(jù)分層法進行求解，將冗余度盡可能小作為主要目標，先求解列分類較少這一目標。這一步我們直接應用問題一的結(jié)論，在此基礎(chǔ)上進行冗余度最小的求解。為了減少冗余度，我們首先將問題一結(jié)論中放置方案的每個藥盒的冗余度計算出來，見表4(程序見附件3)表4：問題一結(jié)論基礎(chǔ)上每個藥盒的冗余度表藥品編號藥品規(guī)格d/mm藥槽規(guī)格d/mm冗余度/mm6691019777410197135210197147110197 為下面說明方便，我們定義： 冗余權(quán)重：將同寬的藥盒歸類，考慮每一類在總藥盒數(shù)中的比例，即 (9) 分類加權(quán)冗余度：同寬藥盒的冗余度乘以冗余權(quán)重，即 (10) 分權(quán)

16、貢獻率： 將分類加權(quán)冗余度求和，考慮每一類在總加權(quán)冗余的比列，即 (11)在此基礎(chǔ)上我們計算出分類加權(quán)冗余度表，并按升序排列，加權(quán)冗余度和分權(quán)貢獻率的具體值如下表5（程序見附錄4和附錄5）表5：各種列寬加權(quán)冗余權(quán)度和分權(quán)貢獻率32001700540.1042209480.000144248550.3126628450.0004327744因此，要最盡量減少冗余度，我們首要減少在加權(quán)冗余度較高的那些列。即對此類列增加相應規(guī)格的藥槽。我們?nèi)?.9作為基準，和表5反序，從大到小排列，認為中的各個藥盒，對寬度冗余和豎向隔板間距數(shù)量有較大的影響。因此我們對此考慮優(yōu)化?；趩栴}一基礎(chǔ)上合并優(yōu)化列寬冗余度步

17、驟第一步：從列分類中找到最大的開始，單獨為其設(shè)置一個無冗余的規(guī)格，同時在此規(guī)格里，在其前面的藥盒規(guī)格同樣可以放入這個規(guī)格。第二步：從剩下的列中重復第一步，直至90%之前的列寬都放入新設(shè)置的規(guī)格中。根據(jù)上述步驟，我們很容易找到在90%的基準下，產(chǎn)生7個豎向隔板類型數(shù)時，分別是：19mm，22mm，34mm，37mm，46mm，47mm，58mm。每種藥盒所對應的最佳寬度范圍如下表6：表6：冗余度最佳取值范圍表豎向隔板間距單位（mm）對應的藥盒規(guī)格寬度的最佳取值范圍單位（mm）1910，172218，203421，323733，354636，4447455846，56用excel將7種不同的隔板間

18、距分別求出對應的編號如表7（全部數(shù)據(jù)見附件7）表7：19mm對應的藥品編號（部分）藥品編號藥品編號藥品編號藥品編號669471197129719086147687在新的分列下，加權(quán)冗余度和分權(quán)貢獻率如表8：表8：新的各種列寬加權(quán)冗余權(quán)度和分權(quán)貢獻1042209480.0002550.3126628450.0006靈敏度分析，我們適當修改選取比例由90%變?yōu)?5%，我們?nèi)?.95作為基準，和表5反序，從大到小排列，認為中的各個藥盒，對寬度冗余和豎向隔板間距數(shù)量有較大的影響。因此我們對此考慮優(yōu)化。同理計算出一個新冗余度，如表9表9：冗余度的最佳取值范圍豎向隔板間距單位（m

19、m）對應的藥盒規(guī)格寬度的最佳取值范圍單位（mm）1910,172218,202721,253426,323733,354636,4447454946,475848,56對比分析減少90%與95%時，冗余度與隔板類型的差別，如表10。表10：隔板數(shù)量與加權(quán)冗余度比較合并累積貢獻度 隔板數(shù)量 加權(quán)冗余度90%7205.395%9367.3原始4722.58從表10中看到，我們以設(shè)計的豎向隔板寬度在90%的基準下，產(chǎn)生7個豎向隔板類型數(shù)時，冗余度減少與隔板類型較優(yōu)，分別是：19mm，22mm，34mm，37mm，46mm，47mm，58mm4.3 問題三的模型建立與求解4.3.1 平面冗余度最小多目

20、標規(guī)劃模型模型的分析問題二我們已經(jīng)求得7個最優(yōu)列寬分類，在此基礎(chǔ)上考慮高度上的冗余。記表示可能的高度，其中根據(jù)題意，首先得到高度冗余計算公式 （12） 表示第個藥盒放入的藥槽可能高度。表示每個藥盒放入對應藥槽時的寬度冗余。因此，每個藥盒產(chǎn)生的平面冗余為 （13）總平面冗余為 （14）同時根據(jù)，根據(jù)題意，藥柜的設(shè)計要滿足藥柜的寬度和高度的要求，同時要將1919種藥品放入柜中，因此滿足一下三個約束條件。橫向約束：設(shè)每類列寬在2.5米規(guī)格的藥柜里， 表示一排橫向同類型列的藥槽個數(shù)，因此有 （15）縱向約束：同時，高度不能超過1.5米，因此有 （16）其中表示同類型高的在一列縱向中藥槽個數(shù)。藥品種類約

21、束： (17)模型的建立根據(jù)上述分析，在冗余度最小的情況下，同時兼顧橫向隔板類型最小，因此建立雙目標規(guī)劃模型如下 目標： (18) (19)模型的求解模型也屬于多目標優(yōu)化，對此我們將減少隔板類型作為主要目標，同時兼顧平面冗余度。由于模型在平面橫向，縱向兩方向都存在未知量，因此我們可以先將列項優(yōu)先考慮，然后再考慮橫向隔板劃分。首先我們，在問題二基礎(chǔ)上統(tǒng)計7類中，每一類所包含的藥品的高度分布表，（程序7）得到表11表11：高度分布表列寬高度282930313222004812342221283700000并且統(tǒng)計出按列每一類的包含的藥的數(shù)量如表12表12列19223437464758包含藥品數(shù)28

22、642766912222054141因為藥柜設(shè)計時為橫列矩形格子，所以，在考慮的個數(shù)時，應當盡可能滿足 ， (20)因此根據(jù)表13可算出，之間的比例關(guān)系大致為表13：的比例關(guān)系表比例關(guān)系5.297.912.42.34.112.6 根據(jù)約束條件，為了節(jié)約空間，盡可能達到2500，并且藥槽列寬大的可以放置部分藥盒列寬較小藥盒，而反過來不行，因此，根據(jù)此原則，上述比例適當修正可以得到表14表14：對應的個數(shù)表個數(shù)1117285726此時，長度上滿足 (21)根據(jù)列的個數(shù)確定，我們估算行的個數(shù)。通過的數(shù)量與的比值我們得到，橫向隔板大致數(shù)量。表15：對應的的數(shù)量表列的類型1922343746475826

23、25.123.824.431.4262723.5從上面數(shù)據(jù)中可以看出，數(shù)據(jù)越小表明此類格子行上比較空閑，例如34列寬的藥槽，在藥柜中按照原始藥盒的放置會有較多剩余。因此可以考慮將22中部分放入34中，同時將19中部分放入22中達到平衡，同理46，47中部分可以放入58中。還可以估計上述的平均值為25.4,因此橫向隔板的數(shù)量為26。此時也滿足條件約束 (22)最后，為了盡量控制平面冗余，在實際設(shè)計時，同行的高度需要盡可能的相等，因此，我們根據(jù)已經(jīng)確定好的一行76個藥槽，對高度進行劃分。我們可以將所有藥盒按照高度排列，同時考慮按列寬分類統(tǒng)計。通過程序（見附錄8）算出76節(jié)點處的高度，為每一行高度的

24、基本值。36 39 43 46 48 51 53 57 62 63 64 67 68 69 71 72 73 74 75 77 78 82 83 87 102 127計算可知該基準值值和大于1500，分析發(fā)現(xiàn)其中部分為最后個別高度所限制，所以可以將該藥品放入下一行中，從而調(diào)節(jié)該行高度，計算出高大致需要以下9類。34mm，41mm，47mm，54mm，60mm，72mm，85mm，101mm，125mm4.4 問題4的模型建立與求解4.4.1 最小槽數(shù)規(guī)劃模型模型分析：設(shè)為不同藥盒的長度，對應的日需求量為，在每天補藥一次的情況下，保證藥房儲藥滿足的需求。首先找出1.5m長的藥槽每一個槽最多放多少

25、個藥盒，設(shè)計出個藥槽，使得一次性能放此藥盒的個數(shù)必須要大于總需求量的最大需求量，建立最小槽數(shù)模型找出滿足條件最小的儲藥槽個數(shù)。模型建立如下： 目標： （23） (24)4.4.2 模型的求解：根據(jù)最小槽數(shù)規(guī)劃模型，分別找出每一種藥品在儲藥柜長度方向上最多能放的個數(shù)，用matlab編程得出各個藥品對應的最小儲藥槽個數(shù)，如表16（全部數(shù)據(jù)見附錄8，程序見附錄第四題程序6）表16：不同藥品對應的藥槽個數(shù)藥品編號藥槽個數(shù)藥品編號藥槽個數(shù)1202212591251312091251211789187854.4.3模型的建立藥店藥品的需求量能夠滿足第一天在銷售中不補給，則需在第二天在銷售藥品之前補藥一次。

26、假設(shè)每個藥槽都放有藥盒，則需儲藥柜最少個數(shù)個數(shù)目標：(進一取整)同時應滿足 （25）模型的求解求每種藥盒裝入藥槽中所需要的藥槽個數(shù)，在excel中用roundup命令，得到進1取整后的個數(shù)，再求出藥槽面積乘以該槽的個數(shù)，對所求乘積進行求和，再對得到總和除以藥柜最大面積，從而得到最少的藥柜數(shù)為2個。5 模型的評價與推廣優(yōu)點：在問題一中通過題目已知條件，建立的最小列寬分類規(guī)劃模型方便有效的找出了不同型號對應的列寬容許度的范圍，區(qū)間無重疊聚類算法方便的找出了有共同區(qū)域的區(qū)間，簡單的剔除了獨立區(qū)間。在問題二中建立的列寬冗余度雙目標規(guī)劃模型和分層法，根據(jù)分層法進行求解，將冗余度盡可能小作為主要目標，先求

27、解列分類較少這一目標。這一步我們直接應用問題一的結(jié)論，在此基礎(chǔ)上進行冗余度最小的求解。缺點：最后找出的是一個區(qū)間，達不到是一個數(shù)的目的，我們根據(jù)累計貢獻度的大小，將一些有效的數(shù)據(jù)剔除掉了，從而在所剩下的數(shù)據(jù)中找出最大值來作為所要求得值。參考文獻1張秀蘭，林峰，數(shù)學建模與實驗，北京：化學工業(yè)出版社，20132周品，趙新芬，matlab數(shù)理統(tǒng)計分析，北京：國防工業(yè)出版社，20093任玉杰，數(shù)值分析及matlab實現(xiàn)，北京：高等出版社，20074丁世飛，靲奉祥，趙相偉，現(xiàn)代數(shù)據(jù)分析與信息模式識別，北京：科學出版社，2012附錄附錄1：function a=qujian(b) %求每一個藥品的容許列區(qū)

28、間m,n=size(b);for i=1:m k=b(i,4); %寬 c=b(i,2); %長 g=b(i,3); %高 dkg=sqrt(k2+g2); %寬高對角線 dck=sqrt(c2+k2); %寬長對角線 a(i,1)=b(i,1); %編號 a(i,2)=k+4;%下限 a(i,3)=min(2*k,dkg,dck);%上限end附錄2： function k=julei(b)%k 類 和具體編號m,n=size(b);c=sortrows(b,2);k=;t=1;kt,2=;for i=1:m-1 if c(i,3)=c(i+1,2) kt,1=c(i,2),c(i,3);

29、kt,2=kt,2,c(i,1); t=t+1; kt,2=; else c(i+1,3)=min(c(i,3),c(i+1,3); kt,1=c(i+1,2),c(i,3); kt,2=kt,2,c(i,1); endend 附錄3：求冗余function a=rongyu(b,c)m,n=size(b)t=1;for i=1:m xia=min(bi,1); bh=bi,2; m1,n1=size(bh); for j=1:n1 a(t,1)=bh(j); a(t,2)=c(c(:,1)=bh(j),4); a(t,3)=xia; a(t,4)=xia-a(t,2)-2; t=t+1; e

30、ndend附錄4：求冗余權(quán)重系數(shù)function dd,d=rl(b)m,n=size(b);d=tabulate(b(:,2);d=d(10:end,:); %累計貢獻率數(shù)據(jù)dd=sortrows(d,3);附錄5：程序3：求分類加權(quán)冗余度function b = jryl( a )%untitled5 summary of this function goes here% detailed explanation goes herem,n=size(a); for i=1:m if a(i,1)=10 b(i,2)=(19-a(i,1)-2)*a(i,3); b(i,1)=a(i,1);

31、elseif a(i,1)=18 b(i,2)=(34-a(i,1)-2)*a(i,3); b(i,1)=a(i,1); elseif a(i,1)=33 b(i,2)=(46-a(i,1)-2)*a(i,3); b(i,1)=a(i,1); elseif a(i,1)=45 b(i,2)=(58-a(i,1)-2)*a(i,3); b(i,1)=a(i,1); endendend附件6：function n=cs(x)%untitled2 summary of this function goes here% detailed explanation goes heren=x(:,3)./(

32、1500./x(:,2) end附件7：function c,jug,r=lkjulei(a,b) %a表示原始數(shù)據(jù)，b表示列寬類別,分類，構(gòu)造按列的高度統(tǒng)計表程序m,n=size(a);m1,n1=size(b);for i=1:m p=0; t=1; while p=0 if a(i,4)=b(t)-2 b=b(t); r(i,1)=a(i,1); r(i,2)=a(i,4); r(i,3)=a(i,3); r(i,4)=b; p=1; else t=t+1; end endendfor i=1:n1jugi,1=r(r(:,4)=b(1,i),3);endfor i=1:7 a=jugi

33、,1; b=tabulate(a); b=b(b(:,2)=0,:); m,n=size(b); for j=1:m c(i,b(j)=b(j,2); endendc=c(:,28:end);c=c;28:125;c=b;0,c;附件8：function t=gfl(c)%高分類程序m,n=size(c);t=0;i=2;while i99 a=a+c(1:7,i); bb=sum(a.*c(1:7,1); if sum(a.*c(1:7,1)=2500 i=i+1; else t=t+1; g(t)=c(8,i-1); p=0; endendend附錄7：藥槽規(guī)格編號藥槽規(guī)格編號藥槽規(guī)格編號

34、藥槽規(guī)格編號藥槽規(guī)格編號藥槽規(guī)格編號藥槽規(guī)格編號藥槽規(guī)格19419927224222707221430341053458034191819928224822712221432341063458234192519934225022713221433341103458334193419936225222714221437341143458634196119950225722733221444341193458734196219953226822735221450341243459034196719954227022737221451341273459234198019962227122742221

35、452341343459434198419971227522745221453341373459734198719975227822753221456341443460434199119990227922758221457341463460534199319100422103227612214583416434610341997191012221082276522146034166346163419991910162211322766221461341743461834191071910182211522769221463341753462134191111910222211622770221

36、472341783462534191121910232211822771221473341793462634191171910302212122779221476341923462734191201910322212522786221477341933463234191221910492212622787221483341943463534191231910512213022788221491341953464434191281910532213122790221507342133465134191411910692213222791221509342183465234191511910702

37、213322792221510342193465334191681910712213522797221511342223465434191771910722213622798221512342273465634191841910762214022815221524342343465934191851910792214522816221525342373466134191961910802214722823221530342393466334191971910812216322826221538342463466434191991910822216522831221546342473466734

38、192301910832216922832221550342483467034192411910852218622838221554342493467534192521910922219822846221556342503467634192531910972220322848221571342513467734192541911002220522849221573342593468134192551911112220922852221574342673468234192691911322221022861221584342713468934192741911332221422869221590

39、342813470534192781911352221522877221600342823470634192871911522221622878221602342843470834192981911532221722886221608342863471034193031911692222122910221610342893471534193081911712222522911221614342903471634193091911722222622912221617342913472034193101911732222822913221620342933472134193121911762223

40、32294222162434295347223419317191177222362294322162634297347263419318191179222402296722162834301347273419332191185222442296922163034305347283419333191188222602297722163134311347293419334191195222622298722163234314347303419335191200222632299822163634315347323419348191204222642299922163734322347383419354191209222