信息建模數(shù)學(xué)建模課程論文 奶牛養(yǎng)殖場酮病群體發(fā)病風(fēng)險評估

2014-2015年信息建模數(shù)學(xué)建模課程論文承諾書我們完全明白，參賽隊員不能以任何方式（包括電話、電子郵件、網(wǎng)上咨詢等）與隊外的任何人（包括指導(dǎo)教師）研究、討論與賽題有關(guān)的問題。我們知道，抄襲別人的成果是違反規(guī)則的，如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻的表述方式在正文引用處和參考文獻中明確列出。我們鄭重承諾，嚴(yán)格遵守規(guī)則，以保證公正、公平性。如有違反規(guī)則的行為，我們將受到嚴(yán)肅處理。我們參賽選擇的題號是（從A/B/C/D中選擇一項填寫）：A參賽隊員(打印并簽名)：1.2.3. （論文紙質(zhì)版與電子版中的以上信息必須一致，只是電子版中無需簽名。以上內(nèi)容請仔細(xì)核對，提交后將不再允許做任何修改。如填寫錯誤，論文可能被取消評獎資格。）日期：2014年12月19日A題奶牛養(yǎng)殖場酮病群體發(fā)病風(fēng)險評估摘要酮病是奶牛生產(chǎn)性疾病的主要疾病，主要發(fā)生于圍產(chǎn)期奶牛（又稱過渡期，指奶牛產(chǎn)前2周至產(chǎn)后2周的一段時期）。每年僅酮病造成的經(jīng)濟損失就達數(shù)億美元，因此對酮病發(fā)發(fā)病風(fēng)險預(yù)警顯得至關(guān)重要。本論文的目的就是通過分析一些有可能誘導(dǎo)奶牛患酮病的各個指標(biāo)，建立單個奶牛以及奶牛場的風(fēng)險指數(shù)體系，對酮病發(fā)病的風(fēng)險來進行預(yù)測，來進一步預(yù)防，減少經(jīng)濟損失。第一問中有兩問，即計算各指標(biāo)在酮病發(fā)病率中的權(quán)重，以及建立每頭奶牛酮病發(fā)病風(fēng)險指數(shù)數(shù)學(xué)模型。在計算各指標(biāo)在酮病發(fā)病率的權(quán)重的時候，我們利用的是層次分析法，將發(fā)病風(fēng)險率轉(zhuǎn)化為各個指標(biāo)所占的權(quán)重。這其中最重要的是對比矩陣的構(gòu)建，因為對比矩陣建立的不合理會導(dǎo)致結(jié)果權(quán)重的不真實，甚至不能通過一致性檢驗，所以我們借用發(fā)病風(fēng)險率來構(gòu)建對比矩陣，誘導(dǎo)發(fā)病風(fēng)險率高的我們認(rèn)為該因素所占的權(quán)重比較大。另外在四個主要指標(biāo)進行對比矩陣的構(gòu)建時候，它們之間的重要程度我們是根據(jù)它們的子層各因素誘導(dǎo)發(fā)病率大于%50的個數(shù)來決定的，而同一因素內(nèi)部的重要程度則直接比較風(fēng)險率。另外血液中的葡萄糖（Glu）、β-羥丁酸（BHBA）、游離脂肪酸（NEFA）、門冬氨酸氨基轉(zhuǎn)移酶（AST）數(shù)據(jù)異?？赡軙苯右鹜〉陌l(fā)生，同時也會間接的通過影響能量負(fù)平衡進而來影響誘導(dǎo)發(fā)病風(fēng)險率。通過以上的分析，以及結(jié)合編程計算，我們求得各指標(biāo)在酮病發(fā)病率的權(quán)重：血液生化指標(biāo)、能量負(fù)平衡、日糧營養(yǎng)水平、體況評分（BCS）所占的權(quán)重分別為：0.5601、0.2622、0.1175、0.05529。另外在建立每頭奶牛酮病發(fā)病風(fēng)險指數(shù)數(shù)學(xué)模型時候，我們建立的是風(fēng)險指數(shù)評價函數(shù)，即對每頭牛的各個因素的等級評分加權(quán)求和求得最后的風(fēng)險指數(shù)評分。再對數(shù)據(jù)的利用時，我們需要對每一個數(shù)據(jù)進行處理，將不同變量的值標(biāo)準(zhǔn)化處理，化為統(tǒng)一量綱，即各個權(quán)重（附件中），帶入公式求解。另外日采食量沒有給出標(biāo)準(zhǔn)值，因此我們要利用給出的公式算出每頭奶牛的標(biāo)準(zhǔn)采食量區(qū)間，最后算出來每頭牛的風(fēng)險指數(shù)值，進行排序，使用MATLAB畫出折線圖。根據(jù)第二問對單個奶牛風(fēng)險指數(shù)模型的建立和求解，第三問中我們也可以建立風(fēng)險指數(shù)評價函數(shù)，我們將50頭奶牛的風(fēng)險指數(shù)分為五個等級，每個等級中的奶牛的數(shù)目所占總頭數(shù)的比例為等級權(quán)重，之后我們?yōu)槊總€等級賦予相應(yīng)的經(jīng)驗得分，加權(quán)求和就得到總得分。通過計算結(jié)果數(shù)據(jù)的觀察和分析嗎，我們可以看出該奶牛場患病奶牛的數(shù)目較多，風(fēng)險較高。關(guān)鍵詞酮病層次分析風(fēng)險指數(shù)綜合評價函數(shù)一．問題重述1.問題背景：奶業(yè)發(fā)展水平是衡量一個國家畜牧業(yè)現(xiàn)代化的重要指標(biāo)之一，并已成為我國畜牧業(yè)的支柱產(chǎn)業(yè)。目前,由于我國高產(chǎn)奶牛飼養(yǎng)管理水平錯位，使得奶牛生產(chǎn)性和營養(yǎng)代謝性疾病也日趨突出，每年造成的經(jīng)濟損失高達數(shù)十億美元，并成為嚴(yán)重制約奶業(yè)發(fā)展的瓶頸之一。酮病是奶牛生產(chǎn)性疾病的主要疾病，主要發(fā)生于圍產(chǎn)期奶牛（又稱過渡期，指奶牛產(chǎn)前2周至產(chǎn)后2周的一段時期）。每年僅酮病造成的經(jīng)濟損失就達數(shù)億美元，因此對酮病的發(fā)病風(fēng)險進行預(yù)警顯得至關(guān)重要（詳細(xì)背景見附件1）。為了對酮病的發(fā)病風(fēng)險進行預(yù)警并對奶牛進行健康評估，獸醫(yī)工作者從1000頭奶牛中完全隨機抽取50頭圍產(chǎn)期奶牛，收集以下指標(biāo)：⑴血液生化指標(biāo)對隨機選取的奶牛清晨空腹采集血樣（所選取的奶牛在血樣采集前已確定為近期內(nèi)沒有進行任何藥物治療，并經(jīng)獸醫(yī)人員確認(rèn)無任何異常狀況），進行生化檢測，所檢測的指標(biāo)有：奶牛血漿中葡萄糖（Glu）、β-羥丁酸（BHBA）、游離脂肪酸（NEFA）、門冬氨酸氨基轉(zhuǎn)移酶（AST）等（各血液生化指標(biāo)健康參考范圍及其誘導(dǎo)的酮病發(fā)病風(fēng)險率見附件2）。⑵日糧營養(yǎng)水平于血樣采集當(dāng)日，在奶牛場管理機房收集養(yǎng)殖場不同時期奶牛的日糧營養(yǎng)水平指標(biāo)（各日糧營養(yǎng)水平指標(biāo)正常參考范圍及其誘導(dǎo)的酮病發(fā)病風(fēng)險率見附件3）。⑶體況評分（BCS）體況評分即對奶牛脂肪儲存程度進行評分，由專業(yè)評審員通過手摸和眼觀相互配合，根據(jù)其脂肪儲存變化情況以及眼觀情況，綜合評估被檢奶牛，確定分值。采用5分值，其中1分表示過度消瘦，5分表示過度肥胖（體況評分的健康參考范圍及其誘導(dǎo)的酮病發(fā)病風(fēng)險率見附件2）。⑷能量負(fù)平衡（NEB）的風(fēng)險率（綜合指標(biāo)）能量負(fù)平衡是低能量攝入和高能量輸出的結(jié)果，即奶牛攝入的能量低于維持、消耗和產(chǎn)出的能量，圍產(chǎn)期奶牛能量負(fù)平衡程度的加重會引起臨產(chǎn)時的健康紊亂。評定時，記錄各組奶牛的體況評分、日采食量、年齡、胎次、當(dāng)日泌乳量，乳脂/乳蛋白率，乳蛋白、乳脂、乳糖含量，結(jié)合血液生化指標(biāo)BHBA和NEFA，綜合評定NEB的風(fēng)險率（見附件2表2）。問題提出：請根據(jù)上述資料及附件，通過數(shù)學(xué)建模，解決如下問題：問題1：請根據(jù)上述資料提供的指標(biāo)，計算各指標(biāo)在酮病發(fā)病率中的權(quán)重，并建立每頭奶牛酮病發(fā)病風(fēng)險指數(shù)數(shù)學(xué)模型，以此評估圍產(chǎn)期奶牛酮病的發(fā)病風(fēng)險。問題2：請根據(jù)您所構(gòu)建的酮病發(fā)病風(fēng)險指數(shù)模型，利用附件4中的數(shù)據(jù)計算每頭圍產(chǎn)期奶牛酮病發(fā)病風(fēng)險指數(shù)，并對50頭奶牛的風(fēng)險指數(shù)進行排序（由大到小）。問題3：請您根據(jù)上述單頭奶牛的發(fā)病風(fēng)險指數(shù)數(shù)學(xué)模型，構(gòu)建奶牛養(yǎng)殖場群體酮病發(fā)病風(fēng)險指數(shù)數(shù)學(xué)模型，并進行預(yù)警。二．模型假設(shè)假設(shè)1.假設(shè)各指標(biāo)的權(quán)值和誘導(dǎo)酮病發(fā)病的發(fā)病風(fēng)險率成正比；假設(shè)2.假設(shè)各個因素的評分和其權(quán)重成正比；假設(shè)3.假設(shè)沒有引起風(fēng)險率時的評分值為0（為計算方便，我們記為其權(quán)重為零）；假設(shè)4.忽略圍產(chǎn)期體況評分和泌乳早期體況分?jǐn)?shù)損失單位的影響；三．模型的建立和問題的求解3.1.1符號說明3.1.2問題分析首先我們進行權(quán)值的求解，在權(quán)值的求解過程中，我們主要是利用了各個指標(biāo)誘導(dǎo)的酮病發(fā)病風(fēng)險率，將它轉(zhuǎn)化為各個因素的權(quán)值，而這一過程的主要方法是利用了層次分析法來求解，通過兩兩比較之間的發(fā)病風(fēng)險率來決定兩兩之間的重要程度，進而構(gòu)建對比矩陣來進行求解權(quán)值。另外在建立每頭奶牛的酮病發(fā)病風(fēng)險指數(shù)數(shù)學(xué)模型的時候，我們利用了多指標(biāo)綜合評分方法。多指標(biāo)綜合評分方法是將多個內(nèi)容、量綱、評價方法和評價標(biāo)準(zhǔn)均不統(tǒng)一的指標(biāo)進行標(biāo)準(zhǔn)化處理，使各指標(biāo)的評價結(jié)果或得分值具有可比性，再通過一定的數(shù)學(xué)模型或算法將多個評估指標(biāo)值計算為一個整體性的綜合評估值。每個指標(biāo)的標(biāo)準(zhǔn)分值與其權(quán)重進行加權(quán)平均，就得到風(fēng)險評價的總分值。3.1.3模型建立首先影響發(fā)病風(fēng)險指數(shù)的主因素有四個,如下圖：能量負(fù)平衡（NEB能量負(fù)平衡（NEB）日糧營養(yǎng)水平體況評分（BCS）發(fā)病風(fēng)險指數(shù)血液生化指標(biāo)發(fā)病風(fēng)險指數(shù)血液生化指標(biāo)血液生化指標(biāo)血液生化指標(biāo)血液生化指標(biāo)BHBANEFAGLUAST能量負(fù)平衡乳蛋白率乳蛋白乳脂乳糖預(yù)產(chǎn)前BHBA預(yù)產(chǎn)前NEFA泌乳早期BHBA泌乳早期NEFA日糧營養(yǎng)水平日糧營養(yǎng)水平能量需要蛋白質(zhì)需要粗脂肪粗纖維干物質(zhì)采食量干物質(zhì)采食量干物質(zhì)采食量我們影響酮病的一級因素有四個，我們將其設(shè)置作為準(zhǔn)則層，表示出層次單排序判斷矩陣，并對其一致性檢驗的步驟如下：（=1\*romani）計算一致性指標(biāo) （=2\*romanii）RI：表示平均隨機一致性的指標(biāo)，經(jīng)過大量的事實比較和理論分析，對，Saaty給出了的值，如下表2所示：123456789000.580.901.121.241.321.411.45表表2隨機一致性指標(biāo)RI的數(shù)值RI的值是這樣得到的，用隨機方法構(gòu)造500個樣本矩陣：隨機地從1~9及其倒數(shù)中抽取數(shù)字構(gòu)造正互反矩陣，求得最大特征根的平均值并定義 （ⅲ）計算隨機一致性比率 當(dāng)時，認(rèn)為判斷矩陣的一致性是可以接受的，否則應(yīng)對判斷矩陣作適當(dāng)修正。我們可以構(gòu)建對比矩陣A如下：而針對于血液生化指標(biāo)，它主要包括四部分，即β-羥丁酸（BHBA）、葡萄糖（Glu)、游離脂肪酸（NEFA）、門冬氨酸氨基轉(zhuǎn)移酶（AST）這四部分，我們同樣可以建立對比矩陣B如下：對導(dǎo)致能量負(fù)平衡的，我們選取了八個指標(biāo)來構(gòu)建對比矩陣，對比矩陣D如下：同樣的對日糧營養(yǎng)，我們構(gòu)建了對比矩陣C，如下：以上我們就求出來了這個因素的權(quán)重，而在構(gòu)建風(fēng)險指數(shù)模型的時候我們利用合評價函數(shù)為如下： 因為每個的權(quán)重我們已經(jīng)求得，之后我們評估分級后的分值進行賦值，就求得了風(fēng)險指數(shù)指標(biāo)值。3.1.4模型求解通過編寫程序，運行，得到四個指標(biāo)在酮病發(fā)病率中的權(quán)重CR=0.043327因為CR<0.1,可以通過一致性檢驗所以W1=0.565010.26220.11750.05529即四個指標(biāo)在酮病發(fā)病率中的權(quán)重為下表：指標(biāo)血液生化指標(biāo)能量負(fù)平衡日糧營養(yǎng)水平體況評分（BCS）權(quán)重0.565010.26220.11750.05529CR=0.049685因為CR<0.1,可以通過一致性檢驗所以血液生化指標(biāo)次級指標(biāo)W2=0.522340.24980.176820.05104血液生化指標(biāo)BHBANEFAGluAST權(quán)重0.29510.09980.14110.02901CR=0.060167因為CR<0.1,可以通過一致性檢驗所以次級指標(biāo)W3=0.0305110.0305110.0207850.0305110.120970.120970.322880.32288CR=0.003185因為CR<0.1,可以通過一致性檢驗所以次級指標(biāo)W4=0.648330.229650.12202W3中是其中的八個因素導(dǎo)致負(fù)平衡的權(quán)重，進而繼續(xù)導(dǎo)致奶?？赡芑疾〉臋?quán)重為下表：影響負(fù)平衡的因素預(yù)產(chǎn)前BHBA預(yù)產(chǎn)前NEFA泌乳早期BHBA泌乳早期NEFA乳蛋白率乳蛋白乳脂乳糖直接誘導(dǎo)發(fā)病的權(quán)重0.03170.03170.08460.08460.0080.0080.0080.00563.1.5結(jié)論通過建立模型和計算我們得到權(quán)重為0.56501、0.2622、0.1175、0.05529。3.2.1符號說明3.2.2問題分析問題一中模型已經(jīng)建好，只要是對數(shù)據(jù)進行標(biāo)準(zhǔn)化處理，其中最為重要的就是日采食量，即需要計算每頭奶牛的標(biāo)準(zhǔn)日采食量區(qū)間。3.2.3模型建立模型的建立主要是應(yīng)用上一問中所建立的風(fēng)險指數(shù)評價函數(shù)，如下：，但是我們在利用50頭奶牛的具體數(shù)據(jù)進行求解的時候，日采食量我們不能直接利用，因為沒有標(biāo)準(zhǔn)值，不知道是不是攝入不足，如果是攝入不足的話會誘導(dǎo)患病，因此我們要計算每頭奶牛的標(biāo)準(zhǔn)日采食量區(qū)間，這就要用到附件3所給的計算公式，如下：干物質(zhì)采食量（DMI）⑴泌乳奶牛干物質(zhì)采食量（DMI）預(yù)測方程：⑴4%FCM=（0.4×產(chǎn)奶量（kg）)+(15×乳脂（kg）);BW——為體重（kg）WOL——為泌乳周（泌乳1周以內(nèi)算1周）e=2.71828⑵非泌乳青年牛DMI預(yù)測方程:DMI（kg/d）=（BW0.75×(0.2435×NEM-0.0466NEM2-0.1128)/NEM）⑵BW——為體重（kg）；NEM為飼糧的維持凈能（Mcal/kg），此處NEM取1.24Mcal/kg。能量需要（NEL）⑴維持能量需要：NEL(Mcal/d)=0.08×BW0.75⑶BW——為體重（kg）（注：對第一胎和第二胎奶牛，維持需要應(yīng)在此基礎(chǔ)上分別增加20%和10%。）⑵產(chǎn)奶能量需要：NEL（Mcal/kg）=0.360+0.0969×乳脂%⑷⑶妊娠能量需要量：平均妊娠能量需求從大約從2.5～3.7Mcal/d。⑸（注：泌乳奶牛能量需量=維持能量+產(chǎn)奶能量；非泌乳奶牛能量需要=維持能量+妊娠能量。）蛋白質(zhì)需要粗蛋白量(g)=5.0×BW0.75～5.6×BW0.75⑹BW——為體重（kg）3.2.4模型求解我們通過編程，計算出每頭奶牛的日采食量標(biāo)準(zhǔn)區(qū)間如下：日采食量(kg)正常區(qū)間26.93(21.68,22.88)21.37(22.29,23.49)22.67(23.38,24.58)19.78(23.64,24.84)23.37(23.38,24.58)20.98(22,23.2)21.92(22.58,23.78)20.35(23.1,24.3)21.91(21.92,22.96)24.83(21.92,23.12)21.68(22.29,23.49)20.99(21.98,23.17)25.89(22.79,23.99)24.51(22.03,23.23)19.88(24,25.2)19.58(22.08,23.28)26.47(24.03,25.23)21.58(23.46,24.66)23.66(20.87,20.95)24.93(21.01,21.09)19.98(20.21,20.28)19.44(19.35,19.43)21.94(19.61,19.68)19.82(20.11,20.19)24.96(20.5,20.58)20.69(19.03,19.11)21.95(19.1,19.18)18.45(19,19.07)19.38(19.28,19.35)19.85(19.13,19.21)22.12(19.86,19.93)19.74(20.68,20.76)21.15(19.42,19.5)21.09(20.22,20.3)22.35(19.33,19.41)20.35(20.47,20.55)20.16(19.59,19.67)21.31(19.59,19.68)18.57(21.977,22.05)19.27(20.93,21.01)19.77(21.38,21.46)21.57(21.04,21.11)22.65(21.05,21.13)20.64(21.33,21.4)19.65(21.39,21.46)20.13(22.42,22.51)21.33(22.88,22.96)19.19(22.87,22.95)20.99(22.65,22.73)19.34(22.67,22.75)如果比區(qū)間的左端點值小，那么就可能誘導(dǎo)患病。之后表格中的50頭牛的所有數(shù)據(jù)我們都可以使用了，因為我們把評價標(biāo)準(zhǔn)均不統(tǒng)一的指標(biāo)進行標(biāo)準(zhǔn)化處理，都作為了各自的權(quán)值。接下來我們對評估分級后的分值進行賦值賦值結(jié)果如下：奶牛各個指標(biāo)BCS（x1）血液BHBA(x2)血液NEFAX3血液Glu(x4)預(yù)產(chǎn)前BHBA(x5)預(yù)產(chǎn)前NEFA(x6)泌乳早期BHBA(x7)泌乳早期NEFA(x8)血液AST(x9)日采食量(x10)乳蛋白率(x11)乳蛋白(x12)乳脂(x13)乳糖（x14）分值6141084477592221因此我們得到風(fēng)險指數(shù)評價函數(shù)的最終形式為代入附件中的數(shù)值，我們得到結(jié)果。3.2.5結(jié)論奶牛風(fēng)險指數(shù)排序如下：奶牛序號數(shù)風(fēng)險指標(biāo)值288.17439278.00809307.81754317.81754347.4698337.20934266.15689195.39014174.63614324.27979363.21669243.1597452.8927152.88584492.80679352.71632.4273212.27004181.98144221.8352141.70154441.57474161.5160491.49229461.40524471.40524501.39924481.3892411.3858291.3710951.3698381.243101.1909471.184821.1843121.184341.057561.057581.0575391.0575401.0575411.0575201.05535251.04975430.52229420.34774130.2536110.1268230370其中我們?nèi)?.17439*0.3=2.452317為標(biāo)準(zhǔn)，大于該值就有高風(fēng)險發(fā)病率。所以我們28、27、30、31、34、33、26、19、17、32、36、24、45、15、49、35號奶牛更有可能患亞臨床酮病。做出折線圖如下：3.3.1符號說明3.3.2問題分析根據(jù)第二問對單個奶牛風(fēng)險指數(shù)模型的建立和求解，第三問中我們也可以建立風(fēng)險指數(shù)評價函數(shù)，我們將50頭奶牛的風(fēng)險指數(shù)分為五個等級，每個等級中的奶牛的數(shù)目所占總頭數(shù)的比例為等級權(quán)重，之后我們?yōu)槊總€等級賦予相應(yīng)的經(jīng)驗得分，加權(quán)求和就得到總得分。3.3.3模型建立首先，我們構(gòu)建奶牛場群體的風(fēng)險指數(shù)模型如下：風(fēng)險指數(shù)等級的劃分風(fēng)險指數(shù)等級一二三四五等級區(qū)間[0,1.634878][1.634878,3.269756][3.269756,4.904634][4.904634,6.529512][6.539512,8.17439]奶牛數(shù)2911226權(quán)重0.580.220.040.040.12經(jīng)驗得分123453.3.4模型求解所以最后我們算的奶牛場的風(fēng)險指數(shù)值為：0.58*1+0.22*2+0.04*3+0.04*4+0.12*5=1.93.3.5結(jié)論通過計算結(jié)果數(shù)據(jù)的觀察和分析嗎，我們可以看出該奶牛場患病奶牛的數(shù)目較多，風(fēng)險較高。四模型的評價與改進1.在建立求解各次級權(quán)重的時候，因為題目中沒有給出來BSC各個子層包括哪些內(nèi)容，因此我們沒有作進一步的分析。2.在計算引起各負(fù)平衡的因素的權(quán)重的時候，本來有11項，因為我們考慮到后面50頭牛中沒有其中的3項，因此我們將它忽略掉了。3.我們在設(shè)計各等級評分的時候，只是簡單的設(shè)計成為正比關(guān)系，如果能確定兩者之間，即權(quán)重和評分值具體的函數(shù)關(guān)系，可以建立更為合理的評價函數(shù)。4.我們在利用50頭牛的數(shù)據(jù)時，它的年齡和胎次我們沒有具體的討論，只是認(rèn)為它在BCS之中包括。 五參考文獻[1]王永明張麗娟，動物疾病發(fā)病規(guī)律及防控措施，新疆阿克蘇市沙雅縣蓋孜庫木鄉(xiāng)畜牧獸醫(yī)站，2014。[2]王清艷，動物外來傳染病輸入風(fēng)險評估模型的建立及其應(yīng)用，黑龍江，2008。[3]袁少華張新超，心血管疾病風(fēng)險評估方法及存在問題，中國心血管雜志，2011，第3期。[4]徐磊，農(nóng)業(yè)巨災(zāi)風(fēng)險評估模型的研究，北京，2012。附件（注:此附錄包含論文中計算的所有程序以及處理后的數(shù)據(jù)）%*************************%′求解四個指標(biāo)的權(quán)重%*************************1.A=[1357;1/3135;1/51/313;1/71/51/31];[x,y]=eig(A);n=diag(y);m=n(1);CI=(m-4)/3;CR=CI/0.9w=x(:,1)/sum(x(:,1))2.%*************************%求解血液各次級指標(biāo)的權(quán)重%*************************B=[1337;1/3125;1/31/215;1/71/51/51];[x,y]=eig(B);n=diag(y);m=n(1);CI=(m-4)/3;CR=CI/0.9w=x(:,1)/sum(x(:,1))3.%*************************%求日采食量各次級指標(biāo)的權(quán)重%*************************C=[135;1/312;1/51/21];[x,y]=eig(C);n=diag(y);m=n(1);CI=(m-3)/2;CR=CI/0.58w=x(:,1)/sum(x(:,1))4.%*************************%求引起血液負(fù)平衡的各次級指標(biāo)的權(quán)重%*************************D=[11211/51/51/81/8;11211/51/51/81/8;1/21/211/21/61/61/91/9;11211/51/51/81/8;5565111/31/3;5565111/31/3;88985511;88985511];[x,y]=eig(D);n=diag(y);m=n(1);CI=(m-8)/7;CR=CI/1.41w=x(:,1)/sum(x(:,1))5.%*************************%泌乳期每頭牛的標(biāo)準(zhǔn)采食量區(qū)間%*************************#include<iostream>#include<cmath>usingnamespacestd;intmain(){intBW,WOL;floatDMI,nai,ru,NE,DBmin,DBmax,Zmin,Zmax,NAI;constfloate=2.71828;cin>>BW>>nai>>ru>>WOL;DMI=(0.372*((0.4*nai)+(15*ru))+0.0968*pow(BW,0.75))*(1-pow(e,(-0.192*(WOL+3.62))));NE=0.08*pow(BW,0.75);NAI=0.360+0.0969*ru*0.01;DBmin=5.0*pow(BW,0.75)*0.001;DBmax=5.6*pow(BW,0.75)*0.001;Zmin=DMI+NE+DBmin+NAI;Zmax=DMI+NE+DBmax+NAI;cout<<"Zmin="<<Zmin<<""<<"Zmax="<<Zmax<<endl;return0;}6.%*************************%泌乳期每頭牛的標(biāo)準(zhǔn)采食量區(qū)間%*************************#include<iostream>#include<cmath>usingnamespacestd;intmain(){intBW;floatDMI,NE;cin>>BW;DMI=(pow(BW,0.75)*(0.2435*1.24-0.0466*pow(1.24,2)-0.1128)/1.24);NE=0.08*pow(BW,0.75);cout<<"DMImin="<<DMI+NE+2.5<<""<<"DMImax="<<DMI+NE+3.7<<endl;returnDMI;}處理后50頭奶牛各個指標(biāo)的權(quán)重值表：奶牛體重采樣BCS(x1,0.05529)血液BHBA血液NEFA血液Glu編號(kg)時間(x2,0.2984)（x3,0.1243）（x4,0.1074）160012000.124302623100000366510000.1243046758000056658000.12430661270000763490000865470000960360.05529000.10741060960.05529000.10741162350000126115000013642300001461330.0552900.124301568920.0552900.124301661520.055290001769010.055290.2984001866810.055290001968900.055290.298400.10742069110000.10742162020.05529000.1074226234000.124302363450000246546000.12430.10742567870000.10742662780.055290.298400.10742762890.055290.29840.12430.107428619100.055290.29840.1243029630110.05529000.107430631120.055290.29840.12430.107431663130.055290.29840.12430.107432696140.0552900.12430.107433634150.055290.29840.12430.1074346721600.29840.12430.10743563617000.12430.107436688180.0552900.12430.1074376381900003863920000.124303968321000040641-8000041659-9000042643-90.0552900043646-100.0552900044647-110.0552900.1243045651-12000.1243046649-150.0552900047666-160.0552900048651-160.0552900049652-180.0552900.1243050653-200.05529000預(yù)產(chǎn)前BHBA預(yù)產(chǎn)前NEFA泌乳早期BHBA(x7,0.0846)泌乳早期NEFA(x8,0.0846)血液AST(X9，0.03491)(x5,0.0317)(x6,0.0317)00.031700000.031700000.03170000000000.03170000000000.0317000000000.03170000.03491