汪月書,李彩娟,許鄭超,凌去非(蘇州大學(xué)醫(yī)學(xué)部基礎(chǔ)醫(yī)學(xué)與生物科學(xué)學(xué)院江蘇 蘇州 215123)
梭鱸不同月齡性狀的主成分分析與判別分析
汪月書,李彩娟,許鄭超,凌去非
(蘇州大學(xué)醫(yī)學(xué)部基礎(chǔ)醫(yī)學(xué)與生物科學(xué)學(xué)院江蘇 蘇州 215123)
為研究梭鱸幼魚期的生長(zhǎng)規(guī)律和生長(zhǎng)特征,對(duì)1—4月齡梭鱸的全長(zhǎng)、體長(zhǎng)、頭長(zhǎng)、體高、體厚、尾柄長(zhǎng)、尾柄高、口裂和體質(zhì)量等9個(gè)性狀進(jìn)行主成分分析與判別分析,并計(jì)算肥滿度。結(jié)果表明:梭鱸各月齡各性狀間多為正相關(guān),體長(zhǎng)與全長(zhǎng)、體質(zhì)量相關(guān)性最顯著,1月齡梭鱸尾柄高與口裂相關(guān)系數(shù)最小。梭鱸各月齡主成分分別為:1、4月齡第一主成分與長(zhǎng)度發(fā)育有關(guān),2、3月齡第一主成分與體質(zhì)量發(fā)育有關(guān);1、2月齡第二主成分與尾柄發(fā)育有關(guān),3、4月齡第二主成分與口裂發(fā)育有關(guān);1、2、3月齡第三主成分別與體高、體厚、尾柄發(fā)育有關(guān)。通過(guò)逐步判別分析,得到各月齡梭鱸判別函數(shù),總判別準(zhǔn)確率為95.6%。4月齡梭鱸肥滿度最大(1.275),3月齡梭鱸肥滿度最小(1.022)。
梭鱸;主成分分析;判別分析;肥滿度
doi:10.3969/j.issn.1004-2091.2016.04.004
梭鱸(Saner lucioperca)屬于鱸形目(Pericformes)、鱸科(Pereidae)、梭吻鱸屬(Lucioperca)[1],原產(chǎn)于歐洲的易北河、咸海、黑海、里海以及波羅的海等水系,在我國(guó)僅分布于新疆伊犁河水系和額爾齊斯河水系。該魚適溫廣(0~33℃)、耐鹽堿,肉質(zhì)細(xì)嫩,肌間刺少,適于加工,營(yíng)養(yǎng)價(jià)值極高,其蛋白質(zhì)含量超過(guò)19.94%,比鱖(17.56%)高,5種致鮮氨基酸含量為9.31%,高于虹鱒(9.11%)[2],有“淡水魚王”之美譽(yù)。梭鱸抗病能力強(qiáng),大大減少了養(yǎng)殖過(guò)程中魚病防治藥物的使用和污染。近年來(lái),梭鱸被譽(yù)為最具養(yǎng)殖潛力的淡水經(jīng)濟(jì)魚類之一,國(guó)際市場(chǎng)對(duì)梭鱸商品魚的需求量增長(zhǎng)迅速[3],且呈供不應(yīng)求態(tài)勢(shì),國(guó)內(nèi)梭鱸商品魚市場(chǎng)需求也極其旺盛,塘口價(jià)已超過(guò)100元/kg,梭鱸已成為國(guó)內(nèi)外市場(chǎng)一致認(rèn)可的名優(yōu)淡水魚類之一。
國(guó)內(nèi)外學(xué)者圍繞梭鱸生物地理學(xué)[4]、群體遺傳學(xué)[5-6]、生理學(xué)[7-8]、人工繁殖[8]、養(yǎng)殖與增殖[9-11]、營(yíng)養(yǎng)與飼料[2,12]等方面進(jìn)行了一定的研究。但是對(duì)于梭鱸生長(zhǎng)尤其是幼魚階段的生長(zhǎng)相關(guān)研究尚未見(jiàn)報(bào)道。
主成分分析和判別分析兩種多元統(tǒng)計(jì)方法在水產(chǎn)基礎(chǔ)生物學(xué)研究方向有較廣泛的應(yīng)用,郭焱等[13]運(yùn)用主成分分析對(duì)伊犁河、塔什庫(kù)爾干河和托什干河3個(gè)群體的斑重唇魚的性狀進(jìn)行了比較分析。王燕飛等[14]對(duì)1—5月齡三疣梭子蟹8個(gè)指標(biāo)進(jìn)行了主成分分析,發(fā)現(xiàn)1—3月齡三疣梭子蟹形態(tài)性狀的增長(zhǎng)優(yōu)先于體質(zhì)量,4—5月齡體質(zhì)量增長(zhǎng)優(yōu)先于形態(tài)性狀的增長(zhǎng)。竇亞琪等[15]運(yùn)用主成分分析、判別分析兩種多元統(tǒng)計(jì)方法,證明了體質(zhì)量、體長(zhǎng)和體高是影響翹嘴鱖早期生長(zhǎng)階段最重要的性狀指標(biāo)。
研究對(duì)1—4月齡梭鱸的全長(zhǎng)等9個(gè)生長(zhǎng)性狀進(jìn)行主成分分析與判別分析,以探究不同月齡梭鱸表型性狀間的相互關(guān)系,同時(shí)結(jié)合肥滿度等生長(zhǎng)指標(biāo),揭示梭鱸幼魚期生長(zhǎng)規(guī)律和生長(zhǎng)特征,并為其科學(xué)養(yǎng)殖管理和生長(zhǎng)選育的指標(biāo)篩選提供理論依據(jù)。
1.1試驗(yàn)材料
2015年4月蘇州沙家浜東湖現(xiàn)代漁業(yè)科技發(fā)展有限公司進(jìn)行了梭鱸人工繁殖,試驗(yàn)用仔魚放入2×667 m2池塘中進(jìn)行幼魚培育。分別于5—8月(間隔30 d)隨機(jī)在試驗(yàn)池中采取梭鱸樣本。
1.2測(cè)量方法
試驗(yàn)魚經(jīng)MS-222(100 g/L)麻醉,測(cè)量全長(zhǎng)、體長(zhǎng)、頭長(zhǎng)、體高、體厚、尾柄長(zhǎng)、尾柄高、口裂、體質(zhì)量等9個(gè)性狀數(shù)據(jù)。使用電子數(shù)顯游標(biāo)卡尺測(cè)量形態(tài)性狀數(shù)據(jù)(精確到0.01 cm),體質(zhì)量使用電子天平測(cè)量(精確到0.01 g)。測(cè)量方法參照Bonar等[16]的淡水魚類形態(tài)特征測(cè)量標(biāo)準(zhǔn)。
1.3數(shù)據(jù)分析方法
1.3.2原始數(shù)據(jù)轉(zhuǎn)換通過(guò)標(biāo)準(zhǔn)化處理以消除原始數(shù)據(jù)中量綱不同的影響,同時(shí)不改變變量間的相關(guān)性。標(biāo)準(zhǔn)化處理公式:
1.3.3主成分選取利用SPSS 19.0進(jìn)行處理數(shù)據(jù)。用Jacobi法求出9項(xiàng)指標(biāo)(體質(zhì)量、體長(zhǎng)、全長(zhǎng)、體高、體厚、頭長(zhǎng)、尾柄長(zhǎng)、尾柄高、口裂)相關(guān)矩陣的特征值、累積貢獻(xiàn)率及特征向量,為保留原多維空間的信息量并簡(jiǎn)化計(jì)算,依據(jù)累積貢獻(xiàn)率>85%的原則選擇主成分。
1.3.4判別分析數(shù)據(jù)利用SPSS 19.0進(jìn)行處理。將各表型性狀數(shù)據(jù)按月齡分為4組,使用逐步判別法建立梭鱸4個(gè)月齡階段的判別函數(shù)。
2.1肥滿度
1—3月齡梭鱸肥滿度呈下降趨勢(shì),3月齡肥滿度最?。?.022),4月齡肥滿度回升,并達(dá)到最大值(1.275)。4月齡梭鱸肥滿度個(gè)體間差異最大,變異系數(shù)達(dá)38.69%(表1)。
表1 不同月齡梭鱸肥滿度
2.2不同月齡梭鱸性狀相關(guān)性比較
利用相關(guān)系數(shù)法,得到各月齡梭鱸的性狀相關(guān)系數(shù)矩陣(表2),各性狀間多為正相關(guān)。1、2、4月齡梭鱸全長(zhǎng)與體長(zhǎng)相關(guān)系數(shù)最大,分別為0.871、0.945、0.992;3月齡梭鱸體長(zhǎng)與體質(zhì)量相關(guān)系數(shù)最大(0.902);1月齡梭鱸尾柄高與口裂相關(guān)系數(shù)最小(0.059)。
2.3不同月齡梭鱸性狀的主成分分析
不同月齡梭鱸9個(gè)性狀參數(shù)的相關(guān)矩陣特征值和累積貢獻(xiàn)率見(jiàn)表3。1月齡梭鱸提取3個(gè)主成分,前3個(gè)主成分累積貢獻(xiàn)率達(dá)89.29%(≥85%)。同理,2月齡提取4個(gè)主成分,累積貢獻(xiàn)率達(dá)89.27%;3月齡提取4個(gè)主成分,累積貢獻(xiàn)率達(dá)86.67%;4月齡提取2個(gè)主成分,累積貢獻(xiàn)率達(dá)91.15%。2、3月齡提取3個(gè)主成分時(shí),累積貢獻(xiàn)率已達(dá)到80%,因此僅列出3個(gè)主成分以便比較。
不同月齡梭鱸入選主成分的特征向量見(jiàn)表4。1月齡梭鱸第一主成分中特征向量絕對(duì)值較大的為全長(zhǎng)X1、體長(zhǎng)X2、尾柄長(zhǎng)X6,故歸納為魚體長(zhǎng)度發(fā)育情況;第二主成分中特征向量絕對(duì)值較大的為尾柄高X7,歸納為尾柄發(fā)育情況;第三主成分中特征向量絕對(duì)值較大的為高X4,歸納為體高發(fā)育情況。不同月齡梭鱸入選主成分結(jié)果見(jiàn)表5,1、4月齡梭鱸第一主成分相同,對(duì)應(yīng)長(zhǎng)度發(fā)育情況,2、3月齡梭鱸第一主成分對(duì)應(yīng)增重情況;1、2月齡第二主成分相同,對(duì)應(yīng)尾柄發(fā)育情況,3、4月齡第二主成分相同,對(duì)應(yīng)口裂發(fā)育情況。
表2 不同月齡梭鱸各性狀的相關(guān)系數(shù)
表3 不同月齡梭鱸各性狀特征值和累積貢獻(xiàn)率
表4 不同月齡梭鱸入選主成分的特征向量
表5 不同月齡梭鱸的主成分
2.4判別分析
運(yùn)用逐步判別分析方法,剔除1項(xiàng)性狀-尾柄高,保留其余8項(xiàng)性狀,建立了不同月齡梭鱸的判別函數(shù),式中X1、X2、X3、X4、X5、X6、X8、X9分別代表全長(zhǎng)、體長(zhǎng)、頭長(zhǎng)、體高、體厚、尾柄長(zhǎng)、口裂、體質(zhì)量。1至4月齡梭鱸的判別式如下:
將記錄的9項(xiàng)性狀數(shù)據(jù)分別代入上述4個(gè)判別式,對(duì)觀測(cè)樣本按判別函數(shù)進(jìn)行預(yù)測(cè)分類,結(jié)果見(jiàn)表6。判別公式總判別準(zhǔn)確率為95.6%,1月齡梭鱸的判別準(zhǔn)確率為100%。
表6 判別函數(shù)對(duì)觀測(cè)樣本的預(yù)測(cè)分類及準(zhǔn)確率
3.1肥滿度
肥滿度是反應(yīng)魚類生理、營(yíng)養(yǎng)狀況的重要指標(biāo),可為漁業(yè)資源評(píng)估及管理提供參考資料[17],捕撈壓力、水溫和餌料數(shù)量等均會(huì)影響肥滿度的大小,對(duì)養(yǎng)殖池中的幼魚進(jìn)行肥滿度監(jiān)測(cè),有利于科學(xué)評(píng)價(jià)幼魚的健康狀況及養(yǎng)殖措施是否得當(dāng)。
該研究中,5、6、8月梭鱸的肥滿度高于7月,且8月梭鱸肥滿度最高,表明5、6、8月梭鱸的生長(zhǎng)優(yōu)于7月,8月份梭鱸的生長(zhǎng)最佳。在幼魚培育早期,梭鱸可攝食池塘中天然餌料生物以滿足生長(zhǎng)需要,但當(dāng)其達(dá)3月齡時(shí),池中天然餌料數(shù)量相對(duì)減少,這對(duì)梭鱸生長(zhǎng)產(chǎn)生一定影響,表現(xiàn)為肥滿度下降。待補(bǔ)充生物餌料后,4月齡梭鱸生長(zhǎng)速度有所加快,肥滿度及其變異系數(shù)均達(dá)到最大值,由此表明及時(shí)配套餌料是培養(yǎng)梭鱸的關(guān)鍵因素之一,同時(shí),對(duì)梭鱸幼魚階段進(jìn)行選育宜從4月齡開始。
3.2主成分分析
主成分分析的基本思想是降維,將多個(gè)彼此相關(guān)的表型性狀指標(biāo)轉(zhuǎn)化為新的、較少的且相互獨(dú)立的綜合指標(biāo)[18],幫助解決多重共線性問(wèn)題,提高信息利用率,便于更加簡(jiǎn)便、直觀地開展分析工作。
主成分分析結(jié)果表明,1、4月齡梭鱸長(zhǎng)度的生長(zhǎng)處于首要位置,2、3月齡梭鱸體質(zhì)量的增長(zhǎng)優(yōu)先于其他性狀的生長(zhǎng)。因全長(zhǎng)、體長(zhǎng)存在較高相關(guān)性,可將全長(zhǎng)作為1、4月齡梭鱸遺傳選育優(yōu)先選擇指標(biāo),將體質(zhì)量作為2、3月齡梭鱸選育優(yōu)先選擇指標(biāo)。唐瞻陽(yáng)等[19]對(duì)2至5月齡尼羅羅非魚進(jìn)行主成分分析后,認(rèn)為應(yīng)始終堅(jiān)持以體質(zhì)量、體長(zhǎng)為尼羅羅非魚選育主選性狀;竇亞琪[15]利用主成分分析的方法證明了在對(duì)廣東清遠(yuǎn)地區(qū)翹嘴鱖的早期選育工作中,應(yīng)重點(diǎn)選擇體質(zhì)量、體長(zhǎng)和體高等3個(gè)性狀。該研究結(jié)果與梭鱸實(shí)際生長(zhǎng)情況相符且與上述學(xué)者的研究結(jié)果類似。
主成分特征向量數(shù)值可反映出各性狀的發(fā)育趨勢(shì),該研究中,1至4月齡梭鱸第一主成分特征向量數(shù)值表明,2、3月齡梭鱸長(zhǎng)度發(fā)育趨勢(shì)較1月齡梭鱸有所減緩,體質(zhì)量增加速度逐漸上升,達(dá)4月齡后,梭鱸長(zhǎng)度發(fā)育速度回升,體質(zhì)量增加速度趨于緩和。
3、4月齡第二主成分對(duì)應(yīng)口裂發(fā)育情況,說(shuō)明7、8月份時(shí)梭鱸口裂發(fā)育速度加快,口裂發(fā)育可提高梭鱸捕食的選擇性,保證其生長(zhǎng)營(yíng)養(yǎng)需求,這與4月齡梭鱸肥滿度最高存在一定相關(guān)性。2、3月齡需要提取四個(gè)主成分后累積貢獻(xiàn)率方達(dá)到85%,推測(cè)可能與6、7月份時(shí)梭鱸各性狀生長(zhǎng)差異較大有關(guān),且3月齡梭鱸恰處于7月份相對(duì)炎熱的環(huán)境中,水溫較高,生長(zhǎng)環(huán)境條件欠佳。
3.3判別分析
判別分析是一種根據(jù)觀測(cè)樣本的若干數(shù)量特征對(duì)樣本進(jìn)行歸屬判別的多元統(tǒng)計(jì)分析方法[20],常用于研究不同種或月齡群體的形態(tài)差異。當(dāng)存在較多的觀測(cè)指標(biāo)時(shí),可利用逐步判別法將判別能力相對(duì)較弱的指標(biāo)剔除,挑選最佳指標(biāo)建立判別函數(shù),以提高判別效果。該研究對(duì)9項(xiàng)性狀進(jìn)行逐步判別后,剔除尾柄高這一性狀,保留其余8項(xiàng)性狀,原因可能是尾柄高在判別分析中的貢獻(xiàn)性較小,與其他性狀間存在較高的相關(guān)性,這與竇亞琪等[15]對(duì)翹嘴鱖進(jìn)行判別分析時(shí)發(fā)現(xiàn)全長(zhǎng)、尾長(zhǎng)、尾柄長(zhǎng)與其他性狀存在著較高的相關(guān)性而這三個(gè)性狀剔除類似。
唐瞻陽(yáng)等[19]對(duì)2至5月齡尼羅羅非魚進(jìn)行判別分析,判別準(zhǔn)確率為99.25%;竇亞琪等[15]對(duì)2至6月齡翹嘴鱖進(jìn)行判別分析,判別準(zhǔn)確率為98.87%。該研究得到的判別公式判別函數(shù)總判別準(zhǔn)確率為95.6%,判別效果較為理想,其中,1月齡梭鱸的判別準(zhǔn)確率為100%,而4月齡梭鱸判別準(zhǔn)確率最低(91.3%),這可能與4月齡時(shí)梭鱸個(gè)體間生長(zhǎng)差異增大有關(guān)。在梭鱸養(yǎng)殖以及育種過(guò)程中,可將該法用于跟蹤探究梭鱸的生長(zhǎng)狀況及生長(zhǎng)差異,以此判定錯(cuò)過(guò)最佳生長(zhǎng)季節(jié)梭鱸的真實(shí)月齡。同時(shí),可以綜合生長(zhǎng)檢測(cè)、肥滿度以及綜合判別分析結(jié)果,評(píng)價(jià)梭鱸魚種培育效果,不斷優(yōu)化管理措施及在幼魚階段的選育策略。
[1]孟慶聞,蘇錦祥,繆學(xué)祖,魚類分類學(xué)[M].北京:中國(guó)農(nóng)業(yè)出版社,1995:643.
[2]陳佳毅,葉元土,郭建林,等.梭鱸、河鱸和加州鱸的肌肉營(yíng)養(yǎng)成分分析[J].飼料研究,2007,(9):52-54.
[3]Fontaine P.Development of European inland fish culture and domestication of new species[J].Cahiers Agricultures. 2009,18(2):144-147.
[4]Haponski A E,Stepine C A.Phylogenetic and biogeographical relationships of the Sander pikeperches(Percidae:Perciformes):patterns across North America and Eurasia[J].Biological J Linnean Society,2013,110:156-179.
[5]黃金善,范兆廷.梭鱸的淋巴細(xì)胞培養(yǎng)及其核型研究[J].東北農(nóng)業(yè)大學(xué)學(xué)報(bào),2012,43(3):104-109.
[6]Poulet N,Balaresque P,Aho T,et al.Genetic structure and dynamics of a small introduced population:the pikeperch,Sander lucioperca,in the Rhone delta[J].Genetica,2009,133:77-86.
[7]周志峰,凌去非.梭鱸夏花魚種耗氧率和窒息點(diǎn)的研究[J].水利漁業(yè),2006,26(4):10-11.
[8]陳羅明,凌去非.三種重金屬離子和兩種消毒劑對(duì)梭鱸胚胎發(fā)育的影響[J].水生態(tài)學(xué)雜志,2008,1(1):109-111.
[9]李永剛,盧自銀,于仕斌,等.梭鱸人工繁殖和苗種培育技術(shù)[J].中國(guó)水產(chǎn),2004(4):82-83.
[10]Szkudlared M,Zakes Z.Effect of stocking density on survival and growth performance of pikeperch,Sander lucioperca (L.),larvae under controlled conditions[J].Aquaculture International,2007,15:67-81.
[11]Policar T,Stejskal V,Kristan J,et al.The effect of fish size and stocking density on the weaning success of pond-cultured pikeperch Sander lucioperca L.juveniles[J].Aquaculture International,2013,21:869-882.
[12]Schulz C,Bohm M,Wirth M,Rennert B.Effect of dietary protein on growth,feed conversion,body composition and survival of pike perch fingerlings(Sander lucioperca)[J].Aquaculture Nutrition,2007,13(5):373-380.
[13]郭焱,孟瑋,劉建,等.新疆不同河流斑重唇魚(Diptychus maculates)形態(tài)特征比較研究 [J].水產(chǎn)學(xué)雜志,2009,22(3):1-5.
[14]王燕飛,王春琳,母昌考,等.三疣梭子蟹不同月齡形態(tài)性狀的主成分分析[J].生態(tài)科學(xué),2014,3(3):527-532
[15]竇亞琪,梁旭方,易提林,等.翹嘴鱖不同月齡性狀的主成分與判別分析 [J].中國(guó)水產(chǎn)科學(xué),2014,6(6):1116-1124.
[16]Bonar S A,Thomas G L,Pauley G B.Evaluation of the separation of triploid and diploid grass carp,Ctenopharyngodon idella(Valenciennes),by external morphlogy[J].Journal of Fish Biology,1988,33(6):895-898
[17]李忠爐,金顯仕,單秀娟,等.小黃魚體長(zhǎng)-體質(zhì)量關(guān)系和肥滿度的年際變化[J].中國(guó)水產(chǎn)科學(xué),2011,18(3):602-610.
[18]廖銳,區(qū)又君,勾效偉,等.黃唇魚、大黃魚、丁氏和棘頭梅童魚的形態(tài)差異和判別分析 [J].大連水產(chǎn)學(xué)院學(xué)報(bào),2009,24(4):305-310.
[19]唐瞻陽(yáng),肖俊,李莉萍,等.尼羅羅非魚(Oreochromis niloticus)不同月齡性狀的主成分與判別分析 [J].海洋與湖沼,2012,43(2):288-293.
[20]陳敏瓊,利用SPSS進(jìn)行判別分析的幾個(gè)問(wèn)題的說(shuō)明[J].現(xiàn)代計(jì)算機(jī),2015(2):34-39.
[21]孫詩(shī)萌,安百勝,蔡明城,等.大瀧六線魚不同群體的形態(tài)差異分析[J].水產(chǎn)學(xué)雜志,2015,28(4):6-11.
[22]何 銅,劉小林,楊長(zhǎng)明,等.凡納濱對(duì)蝦各月齡性狀的主成分與判別分析[J].生態(tài)學(xué)報(bào),2009,29(4):2134-2142.
[23]Sara P S,Bahram F,Ghobad A T.Physiological changes in male and female pikeperch Sander lucioperca(Linnaeus,1758)subjected to different photoperiods and handling stress during the reproductive season[J].Fish Physiology and Biochemistry,2013,33(5):1253-1266.
Principal component and discrimination analysis of Sander lucioperca at differnet ages
Wang Yueshu,Li Caijuan,Xu Zhengchao,Ling Qufei
(School of Preclinical Medicine and Life Science,Soochow University,Suzhou 215123,China)
In order to investigate the growth characteristic of juvenile of pikeperch,Sander lucioperca,data of total length,body length,head length,body height,body width,caudal peduncle length,caudal peduncle height,mouth morphometrics and body weight of pikeperch were collected in 4 months after hatch in this study. The morphological traits were analyzed by principal component and discrimination analysis,and the relative fatness of pikeperch was also evaluated.The results showed that the character parameter of traits of pikeperch at different ages mostly had positive correlation.The relationship between body length,total length and body weight was the most notable while that between caudal peduncle and mouth morphometrics was the lowest at the age of one month.The principal components of pikeperch at different ages were different.The first principal component reflected the growth of total length and body length at one month and four months respectively while that reflects the growth of body weight at two and three months of age.In terms of the second principal component,the caudal peduncle development of one-month-old and two-month-old pikeperch was prominent.While at three and four months old,it was related to the mouth morphometrics.The third principal component reveals the growth of body height,body width and caudal peduncle at one,two and three months old,respectively.After removing caudal peduncle height,the discrimination equations were established with the overall accuracy of 95.6%.The relative fatness at 4 months was the largest(1.275),while 3 months the lowest(1.022).
Sander lucioperca;principal component analysis;discrimination analysis;relative fatness
Q174
A
1004-2091(2016)04-0016-07
江蘇省科技廳項(xiàng)目(BY2015039-10);蘇州大學(xué)大學(xué)生課外學(xué)術(shù)科研基金項(xiàng)目(KY2015110A);常熟市科技局項(xiàng)目(CN201408)
汪月書(1995-),女,大學(xué)生,研究方向:水產(chǎn)食品安全,E-mail:yueshu_wang@163.com
凌去非,博士,副教授,E-mail:Lingqf@suda.edu.cn.
(2015-12-10)