摘""" 要:為了給高山高原地區(qū)青花菜的良種選育提供理論指導,運用灰色關聯(lián)度分析法對14個青花菜新品種單株產(chǎn)量與表型性狀進行分析。結(jié)果表明,青花菜單株產(chǎn)量與11個表型性狀的關聯(lián)度大小依次為:花球主莖橫徑gt;花球主莖高gt;葉片寬度gt;花球縱徑gt;開展度gt;花球橫徑gt;株高gt;葉片長度gt;分枝數(shù)量gt;外葉數(shù)量gt;生育期。得出對單株產(chǎn)量的影響以花球主莖橫徑和花球主莖高最大,其關聯(lián)度分別為0.795 4和0.759 5。根據(jù)各性狀權重數(shù)據(jù)得到青花菜綜合評價模型:H=0.080 8Z1+0.090 4Z2+0.092 3Z3+0.083 5Z4+0.088 8Z5+0.095 4Z6+0.084 0Z7+0.093 6Z8+0.091 6Z9+0.097 5Z10+0.102 1Z11,篩選出中青518、中青16、中青118和中青193可作為下一年當?shù)卦囼炇痉锻茝V品種。因此,在青花菜雜交育種過程中,以產(chǎn)量為育種目標,可注重花球主莖橫徑和花球主莖高2個性狀的選擇,同時兼顧其他性狀進行綜合篩選。
關鍵詞:青花菜;產(chǎn)量;表型性狀;灰色關聯(lián)度;綜合評價
中圖分類號:S635.3"""""""""" 文獻標志碼:A""""""""""" 文章編號:1673-2871(2024)11-138-07
收稿日期:2024-02-05;修回日期:2024-08-24
基金項目:財政部和農(nóng)業(yè)農(nóng)村部:國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系專項(CARS-23-G07);河北省現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術體系露地蔬菜創(chuàng)新團隊建設項目(HBCT202414)
作者簡介:申領艷,女,副研究員,研究方向為蔬菜栽培與育種。E-mail:jiayou7132212@163.com
共同第一作者:陳子臻,男,在讀本科生,研究方向為十字花科蔬菜性狀與遺傳。E-mail:2789730539@qq.com
通信作者:閆鳳岐,男,研究員,研究方向為蔬菜栽培與育種。E-mail:zjkyfq@126.com
"""""""""""""""陳寶剛,男,高級農(nóng)藝師,研究方向為蔬菜栽培與推廣。E-mail:66658297@qq.com
DOI:10.16861/j.cnki.zggc.2024.0086
Phenotypic traits analysis and comprehensive evaluation of broccoli in northern Hebei plateau
SHEN Lingyan1, CHEN Zizhen2, LI Shufang1, KANG Shaohui1, CAO Xuemei1, PANG Jianwen1, YAO Yuexia3, CHEN Baogang4, YAN Fengqi1
(1. Zhangjiakou Academy of Agricultural Sciences, Zhangjiakou 075000, Hebei, China; 2. College of Horticulture, Hebei Agricultural University, Baoding 071000, Hebei, China; 3. Chicheng County Agricultural and Rural Bureau, Chicheng 075500, Hebei, China; 4. Zhangjiakou Rural Revitalization Promotion Center, Zhangjiakou 075000, Hebei, China)
Abstract: In order to provide theoretical guidance for the breeding of broccoli in the alpine plateau region, in this study, the yield per plant and phenotypic quantitative traits of 14 new varieties of broccoli were analyzed by grey correlation analysis. The results showed the correlation degree between yield per plant and 11 phenotypic quantitative traits of broccoli was as follows: Center column width gt; Center column height gt; Leaf width gt; Flower head vertical diameter gt; Expansion degree gt; Flower head transverse diameter gt; Plant height gt; Leaf length gt; Branch number gt; Outer leaf number gt; Growth period. It was concluded that the width and height of the central column had the greatest influence on the yield per plant, and the correlation degrees were 0.795 4 and 0.759 5, respectively. According to the weight data of each trait, the comprehensive evaluation model of broccoli was obtained : H = 0.080 8Z1 + 0.090 4Z2 + 0.092 3Z3 + 0.083 5Z4 + 0.088 8Z5 + 0.095 4Z6 + 0.084 0Z7 + 0.093 6Z8 + 0.091 6Z9 + 0.097 5Z10 + 0.102 1Z11. Four varieties were selected as Zhongqing 518, Zhongqing16, Zhongqing 118 and Zhongqing 193, which could be used as local test demonstration and promotion varieties. Therefore, in the process of hybrid breeding of broccoli, if the yield was breeding goal, the two traits of center column width and heigh could be paid attention, and comprehensively consider other traits for comprehensive screening.
Key words: Broccoli; Yield; Phenotypic quantitative traits; Grey correlation degree; Comprehensive evaluation
青花菜又名綠菜花、西藍花、綠花椰菜等,為十字花科蕓薹屬甘藍類蔬菜變種之一[1]。青花菜的主要食用部位是花球,營養(yǎng)豐富,富含抗癌活性成分[2-4],是備受國內(nèi)外消費者青睞的高營養(yǎng)及保健蔬菜。中國農(nóng)業(yè)科學院蔬菜花卉研究所是國內(nèi)最早開展青花菜育種的單位之一,早在20世紀80年代初就開始從事青花菜種質(zhì)資源收集、創(chuàng)新和遺傳育種工作[5-6],保存了1000余份種質(zhì)資源,在青花菜育種技術和品種培育方面處于國內(nèi)領先地位。張家口市農(nóng)業(yè)科學院與之合作承擔了國家青花菜重大科研聯(lián)合體一級測試點任務。目前,青花菜在我國的種植面積超過86 000 hm2,總產(chǎn)量高達4 158 000 t,年出口量基本穩(wěn)定在140 000 t左右[7]。然而,我國青花菜育種工作起步晚,國內(nèi)品種在球型、蕾粒、色澤一致性等外觀品相和抗性等方面存在不足;青花菜種子主要從國外進口,不利于中國青花菜產(chǎn)業(yè)的健康發(fā)展。研究青花菜表型性狀是開展青花菜育種的前提,可以加速育種進程并提高育種效率。申領艷等[8]采用通徑分析法對86份青花菜材料進行研究,發(fā)現(xiàn)產(chǎn)量與花球高、外葉長和外葉寬呈極顯著相關;高旭等[9]利用相關性分析發(fā)現(xiàn),青花菜單株產(chǎn)量與球位高、莖徑和花球橫徑呈極顯著相關;崔麗紅[10]利用相關性分析發(fā)現(xiàn),青花菜單株產(chǎn)量與葉面積指數(shù)、花球直徑呈極顯著相關。但有關高山高原地區(qū)特定的環(huán)境條件下利用灰色關聯(lián)度分析法探索青花菜產(chǎn)量與表型性狀相關性的研究較少。筆者采用灰色關聯(lián)度法研究14個青花菜新品種生育期、株高、開展度、外葉數(shù)、葉片長度、葉片寬度、分枝數(shù)、花球縱徑等11個表型數(shù)量性狀與單株產(chǎn)量間的關系,以揭示各表型性狀對單株產(chǎn)量影響的主次關系,進而提高青花菜育種效率,促進青花菜種質(zhì)資源改良和創(chuàng)新,為選育適合市場的新品種提供理論依據(jù)。
1 材料與方法
1.1 材料
參試材料由中國農(nóng)業(yè)科學院蔬菜花卉研究所培育并提供,青花菜品種共13個,以耐寒優(yōu)秀(日本坂田種苗株式會社)為試驗對照,共計14個青花菜品種(表1)。
1.2 方法
田間試驗于2021年6—9月在張家口市農(nóng)業(yè)科學院壩上(張北縣喜順溝)基地進行,地勢平坦、土壤肥力均勻,前茬作物為結(jié)球甘藍。參試材料于2021年4月28日在張家口市農(nóng)業(yè)科學院沙嶺子基地溫室播種育苗,6月2日定植于壩上(張北縣喜順溝)基地試驗田。試驗地采用黑地膜覆蓋起壟,壟寬60 cm,壟高10 cm,壟溝寬50 cm,每壟2行,株距50 cm。采用隨機區(qū)組設計排列,3次重復,小區(qū)10 m2。試驗地四周設保護行,田間管理與當?shù)厍嗷ú舜筇锷a(chǎn)保持一致。
1.3 表型數(shù)量性狀調(diào)查
按照《花椰菜和青花菜種質(zhì)資源描述規(guī)范和數(shù)據(jù)標準》[11]的方法。每個小區(qū)隨機選取5株青花菜,測定生育期、株高、開展度、外葉數(shù)、葉片長度、葉片寬度、分枝數(shù)、花球縱徑、花球橫徑、花球主莖高、花球主莖橫徑、單株產(chǎn)量等12個性狀。
1.4 數(shù)據(jù)分析
根據(jù)灰色系統(tǒng)理論[12],將青花菜新品種單株產(chǎn)量及11個表型數(shù)量性狀看作一個整體,為一個灰色系統(tǒng),設青花菜新品種單株產(chǎn)量(kg)作為參考序列記為Z0;11個主要農(nóng)藝性狀作為比較序列Zi,生育期(Z1)/ d、株高(Z2)/cm、開展度(Z3)/cm、外葉數(shù)量(Z4)/片、葉片長度(Z5)/cm、葉片寬度(Z6)/cm、分枝數(shù)量(Z7)/個、花球縱徑(Z8)/cm、花球橫徑(Z9)/cm、花球主莖高(Z10)/cm、花球主莖橫徑(Z11)/cm。
采用Microsoft Excel 2010計算各性狀指標數(shù)據(jù)取平均值,采用DPS 7.05統(tǒng)計軟件進行灰色關聯(lián)分析。
2 結(jié)果與分析
2.1 數(shù)據(jù)無量綱標準化處理
14個青花菜新品種的單株產(chǎn)量和主要表型數(shù)量性狀調(diào)查結(jié)果求平均值(表2)。由于原始數(shù)據(jù)中各性狀量綱不同,需進行無量綱標準化處理,計算方法見式(1)。數(shù)據(jù)無量綱標準化處理結(jié)果見表3。
[Z(K)=zi(K)-ziSi]。"""""""""""""""""""""""""""""" (1)
式中,[zi(K)]:原始數(shù)據(jù);[zi:]平均值;[Si:]標準差;K:品系(K=1、2、3、……14;[i]=0、1、2、3、……11);[i]:某個表型性狀。
2.2 參考序列與比較序列間的絕對差值序列
求絕對差值序列見式(2)。
ΔZi(K)= | Z0(K)- Zi(K) |。""""""""""""""""""""""" (2)
式中,ΔZi(K):絕對差值序列;Z0(K):參考序列;Zi(K):比較序列。
參考序列與各比較序列的絕對差值序列計算結(jié)果見表4。
2.3 單株產(chǎn)量性狀與表型數(shù)量性狀間的關聯(lián)系數(shù)
由表4可得,最小絕對差值:
△min=minmin|Z0(K)-Zi(K)|=0.007 6。
最大絕對差值:
△max=maxmax|Z0(K)-Zi(K)|=3.198 3。
ρ稱為分辨系數(shù),ρ[∈](0,1),取值0.5[12];ΔZi(K)為參考序列與比較序列間對應的絕對差值,代入公式(3),計算各性狀的關聯(lián)系數(shù)。
求關聯(lián)系數(shù)[[ξi(K)]],計算公式如下。
["""""""""""" ξi(K)=△min+ρ△maxΔZi(K)+ρ△max][""""""""""""""" =0.007 6+0.5×3.198 3ΔZi(K)+0.5× 3.198 3]。"""""" """""""""(3)
將表4式中差值序列[ΔZi(K)]值代入公式(3),求出參試青花菜單株產(chǎn)量與表型數(shù)量性狀的關聯(lián)系數(shù)[ξi(K)]的值,計算結(jié)果見表5。
參試青花菜單株產(chǎn)量與各表型數(shù)量性狀的關聯(lián)系數(shù)不同,與同一表型數(shù)量性狀關聯(lián)系數(shù)也存在很大差異。如耐寒優(yōu)秀花球縱徑與單株產(chǎn)量關聯(lián)系數(shù)最大為1.000 0,表明該品種對花球縱徑性狀進行選擇就易獲得高產(chǎn)的青花菜新品種;中青520的生育期與單株產(chǎn)量關聯(lián)系數(shù)最小,為0.334 9,表明該品種的生育期對單株產(chǎn)量選擇影響較小。由此可見,關聯(lián)系數(shù)越大,說明該性狀對目標性狀貢獻越大;關聯(lián)系數(shù)越小,說明該性狀對目標性狀貢獻越小。
2.4 單株產(chǎn)量與表型數(shù)量性狀的關聯(lián)度分析
單株產(chǎn)量與各表型數(shù)量性狀間的關聯(lián)度為該性狀所有品種關聯(lián)系數(shù)的算術平均值,計算公式如下。
[Ri=1nk=0nξi(K)。""""""""""""""""""""""""""""""""""""""""""""" (4)]
式中,[R]i為單株產(chǎn)量與各表型數(shù)量性狀間的關聯(lián)度,[ξi(]K)為表型數(shù)量性狀與單株產(chǎn)量性狀間的關聯(lián)系數(shù)。計算出單株產(chǎn)量與表型數(shù)量性狀的關聯(lián)度及關聯(lián)度排序,結(jié)果見表6。
灰色權重系數(shù)([Wi]),計算公式:
[Wi]= [RiRi。""""""""""""""""""""""""""""""""""""""""""""""""""""" (5)]
比較青花菜單株產(chǎn)量與各表型數(shù)量性狀的關聯(lián)度并排序(表6),關聯(lián)度從大到小依次為:花球主莖橫徑、花球主莖高、葉片寬度、花球縱徑、開展度、花球橫徑、株高、葉片長度、分枝數(shù)量、外葉數(shù)量、生育期。單株產(chǎn)量與花球主莖橫徑、花球主莖高的關聯(lián)度較高,分別為0.795 4和0.759 5,這說明二者是影響青花菜單株產(chǎn)量形成的重要因素。其余表型性狀對單株產(chǎn)量形成影響較小,其中,生育期與單株產(chǎn)量關聯(lián)度最小,為0.629 8,對單株產(chǎn)量形成的貢獻小,可不作為主要選擇性狀,但在該地域生態(tài)環(huán)境條件下則需作為一個考慮因素。
2.5 青花菜新品種加權關聯(lián)度及其排序
加權關聯(lián)度計算公式:
Gi =[Wi(K)ξi(K)。"""""""""""""""""""""""""""""""""" (6)]
將表5中的表型數(shù)量性狀與單株產(chǎn)量性狀間的關聯(lián)系數(shù)[ξi(K])和表6中的各性狀的權重系數(shù)代入公式(6)中,計算出青花菜新品種加權關聯(lián)度及其排序,結(jié)果見表7。加權關聯(lián)度是對各參試新品種綜合性能的具體反映,表明參試品種間各指標的差異及綜合性能高低。綜合評價13個青花菜新品種,排名前4的品種為中青518、中青16、中青118、中青193,可作為下一年當?shù)卦囼炇痉锻茝V品種。
3 討論與結(jié)論
灰色關聯(lián)度分析法可用于分析品種性狀與目標性狀間的主次關系,計算品種性狀對目標性狀的影響程度,有助于篩選優(yōu)良種質(zhì)資源及品種,提高育種效率,易學易用,操作簡單,已廣泛應用于多種作物,如玉米[13-17]、葵花[18]、大豆[19]、芝麻[20-21]、甘藍[22]等。
筆者利用灰色關聯(lián)度分析法分析14個青花菜品種單株產(chǎn)量與各表型數(shù)量性狀的關聯(lián)度,從大到小依次為花球主莖橫徑、花球主莖高、葉片寬度、花球縱徑、開展度、花球橫徑、株高、葉片長度、分枝數(shù)量、外葉數(shù)量、生育期。若某一性狀對單株產(chǎn)量形成的因素越密切,該性狀與單株產(chǎn)量關聯(lián)度數(shù)值就越大;反之,則該性狀與單株產(chǎn)量關聯(lián)度數(shù)值就越小。在本研究中,花球主莖橫徑對單株產(chǎn)量形成最密切,朱長志等[23]利用相關性及主成分和聚類的方法對76個青花菜品種的16個主要農(nóng)藝性狀進行分析研究,結(jié)果表明,青花菜的莖徑(花球主莖橫徑)與單株產(chǎn)量形成最密切,這與本研究結(jié)果一致。同時筆者發(fā)現(xiàn),生育期對單株產(chǎn)量形成的關聯(lián)度最小,這與申領艷等[8]利用相關性分析、通徑分析和回歸分析的方法對青花菜品種11個主要農(nóng)藝性狀進行研究發(fā)現(xiàn)用生育期來提高單株產(chǎn)量效果不理想的結(jié)果一致。與申領艷等[24]利用聚類分析、相關性和主成分分析的方法,對86個青花菜品種進行表型性狀綜合評價,研究表明生育期對單株產(chǎn)量形成的結(jié)果一致。但高旭等[9]利用相關性及主成分分析以及系統(tǒng)聚類的方法,對26個青花菜品種的13個重要農(nóng)藝性狀進行分析與評價,發(fā)現(xiàn)青花菜的花球橫徑與單株產(chǎn)量形成最密切,其次是莖徑(花球主莖橫徑),與本研究結(jié)果略有差異,這可能是參試材料、土壤條件、栽培技術以及試驗區(qū)生態(tài)環(huán)境不同等原因造成的。雖然,青花菜單株產(chǎn)量隨花球主莖橫徑和花球主莖高增加而提高,但花球主莖高僅能增大到市場收購標準,反之直接影響其經(jīng)濟效益,因此,根據(jù)市場收購標準,盡可能增大中心柱寬度和中心柱高度,確保青花菜單株產(chǎn)量最大化,創(chuàng)造更高的經(jīng)濟價值。在高山高原地區(qū)青花菜良種選育種過程中,以產(chǎn)量為育種目標,注重中心柱寬度和中心柱高度的選擇,同時綜合兼顧其他性狀的選擇,才能選育出適合特定區(qū)域綜合性狀較好的青花菜新品種。
根據(jù)關聯(lián)度計算各性狀的權重系數(shù),使各性狀權重系數(shù)從主觀賦值變?yōu)榭陀^求值[25],以提高青花菜選育的準確性和效率。根據(jù)各性狀權重數(shù)據(jù)可得青花菜綜合評價模型[26]:H=0.080 8Z1+0.090 4Z2+0.092 3 Z3+0.083 5Z4+0.088 8Z5+0.095 4Z6+0.084 0Z7+0.093 6Z8+0.091 6Z9+0.097 5Z10+0.102 1Z11,進一步對參試青花菜新品種進行綜合評價及排序,與對照品種相比,綜合評價前4的品種分別為為中青518、中青16、中青118和中青193,但本研究試驗品種來源于同一個科研單位,造成了青花菜品種的遺傳背景比較窄的局面,品種之間同質(zhì)化程度高,因此,下一步研究應補充更豐富多樣的品種進行比較分析。
筆者利用灰色關聯(lián)度分析14個青花菜品種的11個表型數(shù)量性狀對單株產(chǎn)量形成的影響程度,結(jié)果表明,花球主莖橫徑對單株產(chǎn)量影響最大,生育期對單株產(chǎn)量影響最小。根據(jù)各性狀權重數(shù)據(jù)得出青花菜綜合評價模型,篩選出4個品種,分別為中青518、中青16、中青118和中青193,可作為下一年當?shù)卦囼炇痉锻茝V品種。在青花菜雜交育種過程中,以產(chǎn)量為育種目標,可注重花球主莖橫徑的選擇,同時綜合兼顧其他性狀,這樣才能獲得單株產(chǎn)量及綜合性狀較好的青花菜新品種。
參考文獻
[1]"" 張小麗,劉玉梅,方智遠,等.青花菜及近緣種屬種質(zhì)資源抗根腫病鑒定[J].植物遺傳資源學報,2016,17(6):1106-1115.
[2]"" 姚雪琴,謝祝捷,李光慶,等.青花菜不同器官中4-甲基亞磺酰丁基硫苷及蘿卜硫素含量分析[J].中國農(nóng)業(yè)科學,2011,44(4):851-858.
[3]"" 林俊城,吳秋云,高燦紅,等.青花菜硫、硒代謝競爭及其對保健功能的影響研究進展[J].中國細胞生物學學報,2011,33(4):422-432.
[4]"" FARAG M A,MOTAAL A A.Sulforaphane composition,cytotoxic and antioxidant activity of crucifer vegetables[J].Journal of Advanced Research,2010(1):65-70.
[5]"" 方智遠,孫培田,劉玉梅,等.青花菜自交不親和系選育初報[J].中國蔬菜,1987(1):27-29.
[6]"" 方智遠,孫培田.青花菜雜種優(yōu)勢利用研究初報[J].中國蔬菜,1990(6):2-5.
[7]"" 李占省,戚如詩,劉玉梅,等.我國青花菜生產(chǎn)布局、價格變化及趨勢[J].長江蔬菜,2021(4):1-5.
[8]"" 申領艷,閆鳳岐,栗淑芳,等.青花菜單株產(chǎn)量與主要農(nóng)藝性狀的相關性分析[J].山西農(nóng)業(yè)科學,2021,49(7):860-864.
[9]"" 高旭,檀國印,朱長志,等.青花菜品種資源重要農(nóng)藝性狀的鑒定與評價[J].浙江農(nóng)業(yè)科學,2022,63(5):1098-1105.
[10] 崔麗紅.青花菜花球質(zhì)量與植株主要農(nóng)藝性狀的相關性研究[J].湖南農(nóng)業(yè)科學,2011(21):19-20.
[11] 李錫香,方智遠.花椰菜和青花菜種質(zhì)資源描述規(guī)范和數(shù)據(jù)標準[M].北京:中國農(nóng)業(yè)出版社,2008.
[12] 鄧聚龍.灰色系統(tǒng)基本方法[M].武漢:華中理工大學出版社,1987.
[13] 王美霞,陳保國,張之奇,等.早熟玉米雜交組合主要農(nóng)藝性狀與產(chǎn)量的灰色關聯(lián)度分析[J].種子,2021,40(1):108-111.
[14] 李清超,馬浪浪,文瓊,等.玉米雜交組合主要農(nóng)藝性狀與產(chǎn)量的灰色關聯(lián)度分析[J].中國農(nóng)學通報,2015,31(30):74-78.
[15] 李光發(fā),李忠南,王越人,等.玉米PH6WC組合與產(chǎn)量相關性狀的灰色關聯(lián)分析[J].玉米科學,2013,21(6):45-48.
[16] 王平喜,進茜寧,吳向遠,等.14個玉米品系主要農(nóng)藝性狀與產(chǎn)量的灰色關聯(lián)度分析[J].河南科技學院學報(自然科學版),2021,49(5):12-18.
[17] 陳春艷,楊珊,李清超,等.5個玉米新組合產(chǎn)量與主要農(nóng)藝性狀的灰色關聯(lián)度分析[J].中國農(nóng)學通報,2022,38(15):11-16.
[18] 董紅業(yè),劉勝利,柳延濤,等.26份油葵品種(系)產(chǎn)量與農(nóng)藝性狀分析及綜合評價[J].種子,2023,42(8):63-69.
[19] 張永芳,錢肖娜,王潤梅,等.不同大豆材料的抗旱性鑒定及耐旱品種篩選[J].作物雜志,2019(5):41-45.
[20] 徐東陽,孫妍,李偉峰,等.13個河南白芝麻品種單株產(chǎn)量與相關農(nóng)藝性狀的灰色關聯(lián)分析及綜合評價[J].種子,2022,41(10):73-77.
[21] 呂偉,文飛,韓俊梅,等.芝麻產(chǎn)量與相關農(nóng)藝性狀的灰色關聯(lián)分析[J].種子,2021,40(7):110-114.
[22] 陳寶剛,閆鳳岐,康少輝,等.結(jié)球甘藍單株產(chǎn)量與表型數(shù)量性狀的灰色關聯(lián)度分析[J].中國瓜菜,2024,37(1):39-44.
[23] 朱長志,張志仙,劉君,等.青花菜主要農(nóng)藝性狀相關性、主成分與聚類分析[J].中國農(nóng)學通報,2015,31(4):73-79.
[24] 申領艷,康少輝,閆鳳岐,等.冀西北環(huán)境下青花菜品種表型性狀綜合評價[J].種子,2021,40(10):90-96.
[25] 任喜波,魏毓棠.蘿卜主要性狀與產(chǎn)量性狀間的灰色關聯(lián)度分析[J].沈陽農(nóng)業(yè)大學學報,2007(4):598-601.
[26] 梁萬鵬,何振富,李曉莉,等.玉米新品種主要農(nóng)藝性狀與產(chǎn)量灰色關聯(lián)度分析[J].飼料研究,2022,45(11):80-84.