尹中江，廖文華，高小麗

(西藏自治區(qū)農牧科學院農業(yè)研究所，西藏 拉薩 850032)

作物種質資源(又稱品種資源、遺傳資源或基因資源)是人類社會生存與發(fā)展的戰(zhàn)略性資源；是作物育種、生物科學研究和農業(yè)生產發(fā)展的物質基礎，是農業(yè)持續(xù)發(fā)展的重要保障，是提高農業(yè)綜合生產能力，是促進生產力發(fā)展的核心，是農業(yè)得以持續(xù)發(fā)展的重要基礎，維系國家食物安全的重要保證和重要資源；同信息資源一樣，必須充分認識作物種質資源的戰(zhàn)略重要性和潛在價值[1]。

西藏農作物各持資源研究工作始于1952年，自1978年開展學科體系以來，主要經(jīng)歷了資源引進、資源收集與考察、農藝性狀鑒定與利用研究、資源貯藏[2]，目前研究工作重點與方向是資源評價利用、建立及適應西藏農業(yè)生產對更多資源需求利用研究。

本研究是基于對所分析青稞資源，研究采用模糊聚類的方法，主要目的是完成區(qū)內地方青稞資源(含野生近緣種)的品質檢測與分析，并建立和完善西藏青稞資源品質性狀信息庫，并做出相關信息的評價，為西藏青稞(及其它作物)育種和生產提供強有力的支撐。

1 材料與方法

1.1 參試青稞資源抗倒性鑒定

項目實施地點資源田間種植部分安排在農業(yè)研究所4號試驗地，占地0.47 hm2。前茬為油菜，土質砂壤，肥力中等，試驗于2017年4月10日播種，內容主要進行田間抗倒性鑒定試驗；安排在農業(yè)研究所品種資源室溫室內進行；實際種植參試800份青稞材料。筆者在數(shù)據(jù)分析實際中，僅選擇了結果中684份材料作為本文分析應用數(shù)據(jù)(表1)。

1.2 模糊聚類分析法

運用模糊數(shù)學方法將模糊性的定義，通過模糊關系運算與合成，將其定量化，確定青稞品種間的親疏關系，通過“物與類聚”，對原有模糊性的概念客觀分型歸類[3]。運用模糊數(shù)學方法之一進行模糊聚類分析。

2 數(shù)據(jù)處理與分析計算

為了對青稞抗倒性進行一個整體評價，主要選擇同抗倒有關株高(cm)、外徑(mm)、壁厚(mm)、莖稈強度4個測量值來衡量青稞莖稈強度與抗倒性。

數(shù)據(jù)處理與分析計算：為了進行評價優(yōu)劣，消除量綱對觀測值的影響，我們對數(shù)據(jù)進行了無量綱化處理計算，為了進行綜合評價，我們采取了模糊聚類分析。通過編制程序，數(shù)據(jù)處理軟件使用：MATLAB 2018 a 和2017 b Windows版上通過運行，計算機硬件使用上：筆記本I7-7700 HQ CPU和臺式I5-2400 CPU 通過運行。

表1 青稞抗倒性數(shù)據(jù)

注：此表限于篇幅只列出部分數(shù)據(jù)。

表2 標準化數(shù)據(jù)矩陣結果

2.1 數(shù)據(jù)標準化

將試驗取得數(shù)據(jù)利用標準差公式中的平移· 極差變換(也稱：極差變換)[4]，對原始數(shù)據(jù)進行數(shù)據(jù)標準化計算處理。

(1)

利用公式(1)標準化后的數(shù)據(jù)范圍介于[0,1]之間，計算得到的數(shù)組矩陣，數(shù)據(jù)如表2所示(限于篇幅，數(shù)據(jù)列出了部分數(shù)據(jù))。

2.2 標定(建立模糊相似矩陣)

標定，也就是計算模糊相似矩陣。其計算公式很多，本文利用歐式距離法(即：一般歐式距離法)[5]對分析數(shù)據(jù)進行標定計算。

rij=1-cd(ui,uj)

(2)

利用公式(2)計算，計算得到的數(shù)組矩陣，數(shù)據(jù)如表3(只保留前5行和后1行部分數(shù)據(jù))。

2.3 模糊聚類分析

2.3.1 計算模糊等價關系矩陣 按照計算所得的模糊相似矩陣，依據(jù)傳遞閉包法計算出模糊等價矩陣。

在模糊相似矩陣的基礎上，經(jīng)過7次計算，可以計算出模糊等價關系矩陣，計算得到的數(shù)組矩陣，數(shù)據(jù)如表4(只保留部分數(shù)據(jù))。

表3 模糊相似矩陣

表4 模糊等價關系矩陣

2.3.2 根據(jù)不同的λ截距水平進行分類 按照模糊等價矩陣和計算的λ截距水平對所得青稞品種進行聚類分析。

模糊等價矩陣基礎上，可以計算出λ截距矩陣：數(shù)據(jù)如表5。

表5 λ截距矩陣及λ截距值

續(xù)表5 Continued table 5

序號123456450.95660.95630.95590.95560.95520.9546460.95450.95430.95420.95300.95200.9517470.95160.95080.95030.94920.94910.9481480.94690.94660.94620.94580.94510.9446490.94440.94330.94320.94200.94140.9398500.93940.93840.93820.93770.93740.9368510.93570.93510.93440.93410.93380.9337520.93260.92720.92490.92320.92240.9088530.89920.87190.86530.81510.81080.8096

注：318個λ截距值。①說明參試青稞材料，據(jù)λ截距值可以分為318個類型。②但實際應用分類時，也要結合實際情況可再分不同需求類型。

2.3.3 畫模糊聚類圖 有了λ截距數(shù)據(jù)值，就可以將青稞材料，依據(jù)表6的λ截距水平，將材料劃歸于不同的類型。結合農業(yè)生產實際意義，聚類則更有意義。圖1中左邊豎排紅色部分為λ值，頂部粉紅色為聚類后的材料序號，右邊豎排藍色部為分類后的序號。

圖1 青稞模糊聚類圖

圖2 青稞糊聚類圖(圖1λ值和排序后數(shù)組局部放大效果圖)

圖3 青稞糊聚類圖(圖1分類值和排序后數(shù)組局部放大效果圖)

序號123456789101112131415113732925525956215445555916026116036346532658664672618563629642394478274328329355381400346148628830034631634114218516414916520615518942292401672022042332362822872852182442892932695357401296330344393402348351368360371375384422641040441243534539645147629139741642921126527172782993653673664054474554775015105594673105248487511539453521527583625374379406430466378424930535341546550832738535228335936429532239242110497230248308389420499533545506496565593640637116674134254344795075436796416635095173095022231226828033739541444143831536144249553145850057213528306377399445456489557606514243302284326184141972092461621884592906165235405535375746176131561563565951853253660063657062358060524925430116320354335376423498253267380382437463471522473175054745965495665946243503564094084754192325731863164364457663933439845249445748154857740744619488139179551614150177208214261225226314383597206273914113534538550607628655542581622587462214844931191344685192102342762982592624335295542259962063314317247049156221322022165421925764623272286317590632652666670325347403443338369427244645354495125155605646692122793184925922556512566867668235842848064957923958236356815121525226227235258238520530609612247270312332386584370

續(xù)表6 Continued table 6

序號12345678910111213141527645365071697944476592102811101251292880971111281411148913115816315913617519019329859910810712339264266439490578241362426661305585032074316481335711044825164241855619513316505891351694855465985885475566603264436132759312014615317610011213814431934067729460833275638541306646747224303159586249342526349094919510911813218645201242277353243492505851701802971998411712713715716818336195174196178187121263203217231281390205147154371616746756806831139873656761018710688963817118222414866263116448145292284060433337239339525273313444450454342569601115126206571944021619823725132310330762625638730461057706741105665684160228673222140567166191774404174694246035413361833686716783435265651601224358623371561738327481156131191012200441241526811924543613173863068714164578192260303321388432472483331215045411812246825357552451586647

注：參試的青稞，模糊聚類排序后的數(shù)組矩陣的排序值。

3 結果與分析

3.1 青稞抗倒性鑒定聚類結果分析

從收獲的青稞材料整體分析，從模糊聚類結果上講，依據(jù)λ值可以分為318個不同的類型，按照平均數(shù)簡單的統(tǒng)計結果，在684份資源里則可以說每2.15個品種可歸為1個相同的類型，即其品種的特征與特性是一樣的，可以看作是同一個品種資源類型。

按照當λ=1時，參試青稞品種各自屬于獨立一類。即：{1,373,29,2,552,595,621,544,555,591,602,611,603,634,653,658,664,672,……,575,52,4,5,15,86,647}；

當λ=0.9957時，分為兩類，即：{1,29,2,552,595,621,544,555,591,602,611,603,634,653,658,664,672,……, 575,52,4,5,15,86,647}和{373}。

按照表3的λ值，依次類推，將參試青稞資源材料歸類分為318個類型。

圖1～3的模糊聚類結果圖顯示，青稞材料的有明顯的差異性。從品種序號排序順序上看，前面的青稞品種如：1,373,29,2,552,595,621,544,555,和排在后面青稞資源如：575,52,4,5,15,86,647有著明顯的差異性。

通過模糊聚類分析后，對應的青稞資源名，經(jīng)過重新排序后，可以摸清其資源親緣關系的遠近程度，對于明確資源的類型是極其重要的。

從表7可以得到一個較為明確的結果類型。例如：前10個資源即：ZDM5766、ZDM5685、ZDM4381、ZDM4611、ZDM5497、ZDM4340、ZDM4684、ZDM44 78、ZDM4387、ZDM4281，從類型上更相似和接近為同親緣類型；后10個資源即：ZDM4821、ZDM4423、ZDM5479、ZDM4513、ZDM5559、ZDM5575、ZDM55 00、ZDM4819、ZDM5781、ZDM4426從類型上更相似和接近為同親緣類型；但前后兩者從親緣和類型上則更加遠離。從資源的利用上，該特性尤其重要，需要人們有針對性的、根據(jù)要求去選擇親緣關系：近、遠及差異大小，才能有明確的目標。

3.2 結合表型數(shù)據(jù)的抗倒性評價

3.2.1 莖粗 模糊聚類的結果，表型數(shù)據(jù)較好是：模糊序號82，種植序號：244，種植名稱：ZDM4423，同聚類結果較為一致。

3.2.2 莖稈強度 模糊聚類的結果，表型數(shù)據(jù)較好是：模糊序號4，種植序號：537，種植名稱：ZDM5575，同聚類結果較為一致。

表7 模糊聚類后的結果排序表

3.2.3 莖稈壁厚 模糊聚類的結果，表型數(shù)據(jù)較好是：模糊序號647，種植序號：313，種植名稱：ZDM4426，同聚類結果較為一致。

3.2.4 株高 結合模糊聚類排序情況：可知道，排序號為682號，模糊聚類后的序號是：15，實際田間編號是：525，種植品種名稱：ZDM4819。其ZDM4819最高121 cm。ZDM5479的株高105.25 cm，其親緣程度同分類結果吻合,可歸屬于高稈類型。參試表型株高最高的種質編號為ZDM4819。青稞株高對于西藏青稞生產具有很好的實際意義，尤其是對于農區(qū)生產和農區(qū)發(fā)展農區(qū)畜牧業(yè)。

3.3 模糊聚類分析的結果綜述

從模糊聚類分析的結果上看，綜合抗倒性優(yōu)質的青稞資源，主要集中模糊排序的最后一些在資源材料中，其整體體現(xiàn)較為抗倒性較為理想和優(yōu)良。

4 小結與討論

(1)通過Matlab程序]編制[6-7]，能夠很容易克服大量的數(shù)組矩陣數(shù)學運算任務，節(jié)省大量的計算時間，避免人為計算出錯。

(2)品種資源數(shù)量大，模糊聚類時，計算費時耗費。雖然分類容易，但出圖的效果比品種資源數(shù)量小的情況差[8]。

(3)由于資源數(shù)量大，評價分析工作和評價方法還需要進一步探討和完善。

(4)對于青稞的抗倒性分析，還應該進行更多的指標進行綜合評價分析；本文所涉及的數(shù)據(jù)僅是綜合指標中的一個部分，對于青稞抗倒和青稞綜合資源整體分析，還需求進一步開展多指標和多數(shù)據(jù)評價分析工作，才能更加有利于向從事青稞育種工作提供更多，更有價值的信息和實際利用的供體材料資源。