陳月霞 陳 龍 查奇芬 景 鵬 謝君平 熊曉夏
(1江蘇大學(xué)汽車與交通工程學(xué)院, 鎮(zhèn)江 212013)(2江蘇大學(xué)財(cái)經(jīng)學(xué)院, 鎮(zhèn)江 212013)
?
基于潛變量SVM的出行方式預(yù)測(cè)模型
陳月霞1陳 龍1查奇芬2景 鵬1謝君平1熊曉夏1
(1江蘇大學(xué)汽車與交通工程學(xué)院, 鎮(zhèn)江 212013)(2江蘇大學(xué)財(cái)經(jīng)學(xué)院, 鎮(zhèn)江 212013)
為提高小樣本下的出行方式選擇模型的預(yù)測(cè)精度,提出了考慮低碳出行心理變量的支持向量機(jī)(SVM)算法.首先基于計(jì)劃行為理論,考慮低碳出行心理因素,建立多原因多指標(biāo)潛變量模型.然后將預(yù)測(cè)后的潛變量帶入SVM分類器,構(gòu)建了帶潛變量的SVM選擇模型.最后,利用交叉驗(yàn)證優(yōu)化所建模型參數(shù),并以長三角地區(qū)城市居民為研究對(duì)象實(shí)證檢驗(yàn)了模型性能.實(shí)證結(jié)果表明,所建帶潛變量的SVM選擇模型具有較好的預(yù)測(cè)效果,比不帶潛變量的SVM選擇模型的精度提高了4.54%,比傳統(tǒng)的帶潛變量的混合選擇模型提高了2.56%,同時(shí)驗(yàn)證了小樣本下模型仍然具有很高的精度.本研究為出行方式選擇模型和低碳出行方式選擇研究提供了一定的理論參考.
混合選擇模型;支持向量機(jī);多原因多指標(biāo);計(jì)劃行為理論;交叉驗(yàn)證算法
出行方式選擇行為模型不斷利用其他學(xué)科領(lǐng)域的理論進(jìn)行完善和更新[1].一些學(xué)者將多原因多指標(biāo)模型(MIMIC)或者結(jié)構(gòu)方程模型(SEM)與傳統(tǒng)離散模型Logit相結(jié)合,建立了潛變量與顯變量共同作用的混合模型[2-3].
雖然預(yù)測(cè)精度得到一定的提高,但是模型估計(jì)過程復(fù)雜而繁瑣.同時(shí),隨著研究的不斷深入,變量增加,模型日趨復(fù)雜,就需要更多的樣本數(shù)據(jù)來滿足傳統(tǒng)統(tǒng)計(jì)模型建立的需求,但實(shí)際上往往不能得到足夠多的數(shù)據(jù),導(dǎo)致模型無法建立.而支持向量機(jī)(SVM)模型是建立在小樣本下的分類算法[4-5],同時(shí)模型的建立比較簡單、智能,只需選擇恰當(dāng)?shù)目刂茀?shù),就能實(shí)現(xiàn)較好的擬合效果,避免了上述傳統(tǒng)統(tǒng)計(jì)方法的缺陷.
目前國外學(xué)者在交通領(lǐng)域主要運(yùn)用SVM模型預(yù)測(cè)短期出行時(shí)間、客流和費(fèi)用[6-7]等因素對(duì)出行選擇的影響,國內(nèi)不少學(xué)者已經(jīng)嘗試將SVM理論運(yùn)用于中國城市居民出行模式識(shí)別和預(yù)測(cè)研究[8-9]中,但單純的SVM模型的精度還有待提高.為了降低建模難度同時(shí)兼顧模型的預(yù)測(cè)精度,本文結(jié)合傳統(tǒng)混合模型和SVM的優(yōu)勢(shì),提出一種新的MIMIC-SVM選擇模型來研究出行方式選擇行為,并將此方法用于低碳出行研究領(lǐng)域.
本文以計(jì)劃行為理論(TPB)為低碳出行心理研究框架,針對(duì)出行方式建立MIMIC模型.首先,將心理潛變量的擬合值作為解釋自變量加入到支持向量機(jī)中,形成帶潛變量的SVM選擇模型.然后,采用實(shí)證數(shù)據(jù)對(duì)傳統(tǒng)離散選擇模型和不帶潛變量的SVM選擇模型的擬合精度進(jìn)行對(duì)比驗(yàn)證.最后,對(duì)小樣本下的模型性能進(jìn)行檢驗(yàn).
1.1 出行方式分類
出行行為研究通常根據(jù)交通工具種類對(duì)出行方式進(jìn)行分類,也有學(xué)者[10]在研究低碳出行時(shí)將出行方式分為低碳出行和非低碳出行兩大類等,但是這些分類方式或過于精細(xì)或過于粗獷,不能滿足本文研究要求.本文根據(jù)各類交通工具的個(gè)人碳排放強(qiáng)度[11]將出行方式重新劃分為3個(gè)等級(jí):① 碳排放強(qiáng)度0~10為低碳出行方式,包含步行、自行車和電動(dòng)車;② 10~100為中碳出行方式,包含摩托車和公交車;③ 100以上為高碳出行方式,包含出租車和私家車.
1.2 數(shù)據(jù)來源
本文采用的數(shù)據(jù)來源于2015年長三角地區(qū)交通低碳出行行為問卷調(diào)查.問卷分為出行者信息調(diào)查、出行信息調(diào)查和出行態(tài)度調(diào)查3個(gè)部分.出行態(tài)度調(diào)查量表主要用于測(cè)量計(jì)劃行為理論框架下通勤者低碳出行的心理變量.參考國內(nèi)外有關(guān)出行方式TPB各變量的問題項(xiàng)[12],采用Likert五級(jí)量表構(gòu)造TPB中4個(gè)基本變量Attitude, SN, PBC,Intension,其中1代表非常不符合,2代表不符合,3代表一般,4代表符合,5代表非常符合.根據(jù)文獻(xiàn)[10,13]關(guān)于中國城市居民低碳出行行為影響因素研究結(jié)果,結(jié)合本文研究內(nèi)容,得出5個(gè)特定變量(見表1),其他常用個(gè)人屬性變量還有kid(有無小孩);bike(有無自行車);elecar(有無電動(dòng)車);car(有無私家車)以及male(性別).
表1 數(shù)據(jù)庫變量
本次問卷共發(fā)出3 200份,實(shí)際回收2 941份,回收率為91.91%,刪除其中有明顯矛盾、填寫不完全的21個(gè)樣本,最后得到2 660個(gè)完整樣本數(shù)據(jù).Stevens[14]建議一個(gè)變量大概需要15個(gè)樣本支持,本文中每個(gè)MIMIC模型有12個(gè)變量,因此樣本量均要求在180以上,而3個(gè)模型樣本總量分別為1 389,726,545,滿足模型建立的樣本量要求.
1.3 數(shù)據(jù)檢驗(yàn)
對(duì)于潛變量數(shù)據(jù),本文從信度和效度2個(gè)角度進(jìn)行檢驗(yàn).采用Cronbachα信度系數(shù)檢驗(yàn)信度,平均方差提取值(AVE)檢驗(yàn)收斂效度,主成分因子分析檢驗(yàn)結(jié)構(gòu)效度.檢驗(yàn)結(jié)果如表2所示.
表2 數(shù)據(jù)檢驗(yàn)結(jié)果
由表2可知,在探索性因子分析中,所有潛變量特征值大于1的主成分因子均唯一,其方差貢獻(xiàn)率都在60%以上,且所有變量問題項(xiàng)的因子載荷都高于0.5.而在一階驗(yàn)證性因子分析中,相應(yīng)的變量與問題項(xiàng)之間均關(guān)系顯著,且標(biāo)準(zhǔn)化路徑系數(shù)都大于0.5,z值也都遠(yuǎn)大于2,證明該量表具有良好的結(jié)構(gòu)效度.各變量的α系數(shù)只有SN為0.64,其他變量都在0.7之上,根據(jù)Nunnally等[15]的評(píng)判標(biāo)準(zhǔn),均在可接受范圍內(nèi),表明各變量的內(nèi)在一致性較高.每個(gè)變量的AVE值均大于0.5,在0.52~0.7 之間,表明該量表收斂效度較好.
本文考慮到模型的復(fù)雜性以及樣本量不足的可能性,將SVM模型運(yùn)用到混合模型中,替代傳統(tǒng)的離散選擇模型.所建模型包含潛變量模型MIMIC和SVM模型,結(jié)構(gòu)如圖1所示.
圖1 帶潛變量的SVM選擇模型
2.1 潛變量模型
定義反映潛變量關(guān)系的結(jié)構(gòu)方程:
η=Λx+ζ
(1)
式中,η為n×1維低碳心理潛變量,本文中n=4,具體為Attitude, SN, PBC,和Intension;x為k×1維外生可觀測(cè)變量向量,k=8,包括male, age, edu, income, kid, bike, elecar和car;Λ為n×k維未知待估計(jì)參數(shù)矩陣;ζ為測(cè)量誤差.
測(cè)量方程為
y=Γη+υ
(2)
式中,y為q×1維潛變量η的可觀測(cè)指標(biāo)向量;Γ為q×n維未知的待估計(jì)參數(shù)矩陣;υ為誤差項(xiàng).對(duì)于誤差項(xiàng)ζ和υ需滿足
E(ζζT)=Ψ,E(υυT)=Θ,E(υζT)=0
(3)
2.2 SVM模型
記n個(gè)樣本點(diǎn)集為{(xi,yi),i=1,2,…,n},其中xi={xi,ηi},xi為可觀測(cè)變量,ηi為潛變量,考慮非線性SVM,引入松弛變量來修正優(yōu)化目標(biāo)和約束項(xiàng),即
(4)
式中,C為懲罰因子.得到最終的最優(yōu)分類函數(shù)如下:
(5)
式中,核函數(shù)選用徑向基核函數(shù)(RBF).
3.1 MIMIC模型
利用統(tǒng)計(jì)軟件Stata14建立3類出行方式的MIMIC模型,擬合優(yōu)度如表3所示.由表可知,3類出行方式選擇模型的卡方自由度比值均在可接受范圍1~3之間,并且近似均方根誤差(RMSEA)、相對(duì)擬合指數(shù)(CFI)、塔克-劉易斯指數(shù)(TLI)和標(biāo)準(zhǔn)化殘差均方根(SRMR)都在模型擬合度可接受范圍內(nèi).
表3 模型擬合指標(biāo)
通勤者社會(huì)經(jīng)濟(jì)屬性與潛變量之間的關(guān)系非常復(fù)雜,在這里僅給出低碳出行方式中兩者的關(guān)系情況,見表4.由表可知,在采用低碳出行方式出行的MIMIC模型中,每個(gè)個(gè)體屬性變量都對(duì)一種或者幾種潛變量產(chǎn)生顯著影響.社會(huì)經(jīng)濟(jì)屬性對(duì)低碳出行“行為態(tài)度”影響最大,受教育程度、月收入、是否有自行車和是否有電動(dòng)車都對(duì)出行者采用低碳方式出行的態(tài)度有顯著性影響,其中月收入為負(fù)影響,即收入越高越不支持低碳出行,這可能是因?yàn)殡S著收入的提高,出行者對(duì)出行方式的要求越高,傾向于舒適、快速的出行方式;相反,受教育程度越高,越認(rèn)同低碳出行,且有低碳出行條件優(yōu)勢(shì)的通勤者比如擁有自行車、電動(dòng)車等也支持低碳出行.此外,社會(huì)經(jīng)濟(jì)屬性對(duì)低碳出行“行為意向”的直接顯著影響最小,只有是否有電動(dòng)車和否是有私家車分別對(duì)其有正、負(fù)顯著影響,即有電動(dòng)車的通勤者更愿意低碳出行,而有私家車的更傾向于開車出行,這可能是由于自身交通工具的限制造成的.社會(huì)經(jīng)濟(jì)屬性對(duì)于其他2個(gè)潛變量的影響程度介于“行為態(tài)度”和“行為意向”之間.
表4 社會(huì)經(jīng)濟(jì)屬性對(duì)潛變量的影響
3.2 SVM模型
考慮到輸入變量之間差異較大,參照文獻(xiàn)[8]將數(shù)據(jù)都規(guī)整到[-3,3]之間.本文的輸出類別可轉(zhuǎn)換成3個(gè)二類劃分問題解決,每個(gè)分類對(duì)應(yīng)的決策函數(shù)為
(6)
若fj(x)=1,則x屬于第j類;若fj(x)=-1,則x不屬于第j類.由此得到多分類問題的總判別函數(shù):
(7)
式中,arg為選取指標(biāo)函數(shù).則本文中的低碳、中碳和高碳出行方式可分別表示如下:
(8)
為了避免盲目湊試,本文利用交叉驗(yàn)證來尋找SVM模型的最佳參數(shù),再利用最佳參數(shù)訓(xùn)練模型進(jìn)行預(yù)測(cè).將分類精度作為評(píng)判模型性能的重要指標(biāo),計(jì)算公式如下:
(9)
式中,Q為模型分類精度;R為預(yù)測(cè)正確樣本個(gè)數(shù);N為預(yù)測(cè)樣本總量.
利用Matlab 2014進(jìn)行編程,將個(gè)人社會(huì)經(jīng)濟(jì)變量和第1階段的MIMIC模型中的心理潛變量作為樣本點(diǎn)帶入建模,樣本總量為2 660,前1 800個(gè)作為訓(xùn)練樣本,其余860個(gè)作為測(cè)試樣本.綜合考慮計(jì)算時(shí)間和分類精度,將SVM參數(shù)設(shè)置如下:交叉驗(yàn)證系數(shù)H=5,懲罰因子-10≤C≤10,核函數(shù)參數(shù)-10≤σ≤10.程序執(zhí)行結(jié)果顯示最佳懲罰系數(shù)C為0.87,最佳核參數(shù)σ為0.21.此模型的訓(xùn)練樣本分類精度為1 577/1 800=87.61%,測(cè)試樣本的分類精度為723/860=84.07%.
為檢驗(yàn)本文所建帶潛變量的SVM選擇模型的可靠性和實(shí)用性,將其與帶潛變量的傳統(tǒng)Logit混合選擇模型和不帶潛變量的SVM選擇模型進(jìn)行比較,預(yù)測(cè)結(jié)果見表5.
表5 模型預(yù)測(cè)性能比較
可看出,帶潛變量的SVM選擇模型的預(yù)測(cè)正確率最高,帶潛變量的傳統(tǒng)Logit混合選擇模型次之,不帶潛變量的SVM選擇模型的預(yù)測(cè)正確率最低.所建模型比不帶潛變量的SVM選擇模型的精度提高了4.54%,比帶潛變量的傳統(tǒng)Logit混合選擇模型提高了2.56%.
表6中數(shù)據(jù)是以月收入敏感性分析為例,分析低碳出行方式比例變化規(guī)律.第1列為原有比例,后面3列分別為帶潛變量SVM選擇模型、帶潛變量Logit混合選擇模型以及不帶潛變量SVM選擇模型進(jìn)行預(yù)測(cè)后的低碳出行在總樣本中的占比,其中比例變化是指收入由低到高時(shí)低碳出行方式占比的減少量.由表可知,隨著月收入提高,低碳出行方式的占比不斷減少,在3個(gè)模型中,帶潛變量SVM選擇模型的減少幅度最大,與原有比例更接近,帶潛變量的Logit混合選擇模型次之,而不帶潛變量SVM選擇模型減少幅度最小.
表6 低碳出行方式占比與月收入間關(guān)系 %
由此可知,在大樣本量下,帶心理潛變量的選擇模型的精確度要高于不帶心理潛變量的選擇模型,說明低碳心理潛變量的補(bǔ)充能進(jìn)一步提高出行選擇模型預(yù)測(cè)能力;而帶心理潛變量的模型中,SVM選擇模型的預(yù)測(cè)效果要優(yōu)于傳統(tǒng)的離散選擇模型,且這種精度的提高有利于后續(xù)變量的敏感性分析研究.
以上對(duì)于模型精度的驗(yàn)證都是建立在樣本量充足的情況下,在樣本量大幅減少時(shí),傳統(tǒng)離散選擇模型由于樣本量要求可能無法完成建模和預(yù)測(cè),而帶潛變量的SVM選擇模型仍然能夠保持較高精確度.本文選擇280個(gè)樣本進(jìn)行驗(yàn)證,其中200個(gè)樣本用于訓(xùn)練,80個(gè)用于預(yù)測(cè),結(jié)果見圖2.
圖2 小樣本下的模型測(cè)試
由圖2可知,在280個(gè)樣本量的情況下,帶潛變量的SVM選擇模型仍然保持較高的訓(xùn)練精度和預(yù)測(cè)精度,分別達(dá)到86.0%和82.5%,與大樣本下的模型性能相差不大.
1) 本文基于計(jì)劃行為理論,考慮低碳出行心理因素,建立了多原因多指標(biāo)潛變量模型,然后將潛變量模型預(yù)測(cè)出的潛變量帶入支持向量機(jī)分類器,構(gòu)建了帶潛變量的SVM選擇模型,并利用交叉驗(yàn)證算法優(yōu)化支持向量機(jī)模型參數(shù).
2) 以長三角地區(qū)2 660個(gè)居民通勤調(diào)查樣本為研究對(duì)象實(shí)證檢驗(yàn)?zāi)P托阅?實(shí)證結(jié)果表明本文所建模型具有較好的預(yù)測(cè)效果,比不帶潛變量的SVM選擇模型的精度提高了4.54%,比傳統(tǒng)的帶潛變量的Logit混合選擇模型提高了2.56%,同時(shí),進(jìn)行了小樣本下的性能測(cè)試,結(jié)果顯示分類正確率仍然較高,所建模型能夠?yàn)槌鲂蟹绞竭x擇模型和低碳出行方式選擇研究的提供一定的理論參考.
3) 本文所建模型在兼顧精度情況下,相對(duì)于傳統(tǒng)預(yù)測(cè)模型,大大降低了樣本量的要求,因此,實(shí)際建模時(shí)難度降低,同時(shí)在進(jìn)行后續(xù)變量敏感性分析時(shí),精度的提高使得分析更加準(zhǔn)確,模型結(jié)果能比較直觀地反映變化情況,有利于低碳出行相關(guān)政策的制定.
References)
[1]Ben-Akiva M, Mcfadden D, Train K, et al. Hybrid choice models: Progress and challenges[J].MarketingLetters, 2002, 13(3): 163-175.
[2]景鵬,雋志才,查奇芬.考慮心理潛變量的出行方式選擇行為模型[J].中國公路學(xué)報(bào),2014,27(11):84-92. Jing Peng, Juan Zhicai, Zha Qifen. Psychological latent variables into travel incorporating choice model[J].ChinaJournalofHighwayandTransport, 2014, 27(11): 84-92. (in Chinese)
[3]Jing Peng, Juan Zhicai, Gao Linjie. Application of the expanded theory of planned behavior in intercity travel behavior[J].DiscreteDynamicsinNatureandSociety, 2014, 2014: 1-10. DOI:10.1155/2014/308674.
[4]劉廣東.基于支持向量機(jī)的地面驅(qū)動(dòng)螺桿泵井工況診斷技術(shù)[J].排灌機(jī)械工程學(xué)報(bào),2014,32(2):125-129. Liu Guangdong. Working conditions diagnosis of surface-driving progressive cavity pump wells based on support vector machine[J].JournalofDrainageandIrrigationMachineryEngineering, 2014, 32(2): 125-129. (in Chinese)
[5]朱周,路小波,衛(wèi)朋,等.基于超像素和支持向量機(jī)的車輛陰影檢測(cè)算法[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,45(3):443-447. DOI:10.3969/j.issn.1001-0505.2015.03.006. Zhu Zhou, Lu Xiaobo, Wei Peng, et al. Vehicle shadow detection algorithm based on superpixel and SVM[J].JournalofSoutheastUniversity(NaturalScienceEdition), 2015, 45(3): 443-447. DOI:10.3969/j.issn.1001-0505.2015.03.006.(in Chinese)
[6]Tsapakis I, Schneider W H. Use of support vector machines to assign short-term counts to seasonal adjustment factor groups[J].TransportationResearchRecord, 2015, 2527: 8-17. DOI: 10.3141/2527-02.
[7]Jiang X S, Zhang L, Chen X Q (Michael). Short-term forecasting of high-speed rail demand: A hybrid approach combining ensemble empirical mode decomposition and gray support vector machine with real-world applications in China[J].TransportationResearchPartC:EmergingTechnologies, 2014, 44: 110-127. DOI:10.1016/j.trc.2014.03.016.
[8]許鐵,高林杰,景鵬,等.基于PSO-SVM的居民出行方式預(yù)測(cè)模型[J].交通運(yùn)輸系統(tǒng)工程與信息,2011,11(5):155-161. DOI:10.3969/j.issn.1009-6744.2011.05.023. Xu Tie, Gao Linjie, Jing Peng, et al. Prediction model of residents’ trip mode based on PSO-SVM[J].JournalofTransportationSystemsEngineeringandInformationTechnology, 2011, 11(5): 155-161. DOI:10.3969/j.issn.1009-6744.2011.05.023.(in Chinese)
[9]楊敏,丁劍,王煒.基于ARIMA-SVM模型的快速公交停站時(shí)間組合預(yù)測(cè)方法[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,46(3):651-656. Yang Min, Ding Jian, Wang Wei. Hybrid dwell time prediction method for bus rapid transit based on ARIMA-SVM model[J].JournalofSoutheastUniversity(NaturalScienceEdition), 2016, 46(3): 651-656. (in Chinese)
[10]金楠.大城市居民低碳出行方式選擇影響因素研究[D].重慶:重慶交通大學(xué)交通運(yùn)輸學(xué)院,2013.
[11]陳月霞,陳龍,查奇芬.鎮(zhèn)江城市交通微觀主體碳排放測(cè)度[J].江蘇大學(xué)學(xué)報(bào):自然科學(xué)版,2015,36(6):645-649. Chen Yuexia, Chen Long, Zha Qifen. Carbon emissions measurement of urban traffic individual in Zhenjiang[J].JournalofJiangsuUniversity:NaturalScienceEditions, 2015, 36(6): 645-649. (in Chinese)
[12]Ajzen I. The theory of planned behaviour: Reactions and reflections[J].Psychology&Health, 2011, 26(9): 1113-11127. DOI:10.1080/08870446.2011.613995.
[13]劉蔚.城市居民低碳出行的影響因素及引導(dǎo)策略研究[D].北京:北京理工大學(xué)管理與經(jīng)濟(jì)學(xué)院,2014.
[14]Stevens J. Applied multivariate statistics for the social sciences[J].JournalofEducationalStatistics, 2015, 57(100): 68-69.
[15]Nunnally J C, Bernstein I H.Psychometrictheory[M]. New York: McGraw-Hill, 1967: 248-292.
Forecasting model of travel mode based on latent variable SVM
Chen Yuexia1Chen Long1Zha Qifen2Jing Peng1Xie Junping1Xiong Xiaoxia1
(1School of Automobile and Traffic Engineering, Jiangsu University, Zhenjiang 212013, China)(2School of Finance & Economics, Jiangsu University, Zhenjiang 212013, China)
In order to improve the prediction accuracy of the travel mode choice model under small samples, a support vector machine (SVM) algorithm considering the low carbon travel psychological variables is proposed. Based on the theory of planned behavior (TPB), considering low carbon travel psychological factors, latent variable models with multiple causes and indicators are established. Substituting the forecasted latent variables into the SVM classifier, a SVM selection model with latent variables is then proposed. The mixed selection parameters are obtained using cross validation optimization, and the model performance is validated based on urban residents’ data in Yangtze River Delta region. Empirical results show that the established SVM selection model with latent variables has a better prediction accuracy, improved by 4.54% compared with the SVM without latent variables, and 2.56% by the traditional model with latent variables. Results prove that the model still has a high precision with small samples. This study provides a theoretical reference for the travel choice model and low carbon travel choice research.
mixed selection model; support vector machine(SVM); multiple indicators and multiple causes; theory of planned behavior; cross validation algorithm
10.3969/j.issn.1001-0505.2016.06.034
2016-05-08. 作者簡介: 陳月霞(1982—),女,博士生;陳龍(聯(lián)系人),男,博士,教授,博士生導(dǎo)師,chenlong@ujs.edu.cn.
國家自然科學(xué)基金資助項(xiàng)目(71373105, 61573171, 51208232)、江蘇省“六大人才高峰”資助項(xiàng)目(2015-JY-025)、江蘇省高??蒲袆?chuàng)新計(jì)劃資助項(xiàng)目(CXZZ12_0663).
陳月霞,陳龍,查奇芬,等.基于潛變量SVM的出行方式預(yù)測(cè)模型[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2016,46(6):1313-1317.
10.3969/j.issn.1001-0505.2016.06.034.
U491.1
A
1001-0505(2016)06-1313-05