魏正元,霍艷,李文
(重慶理工大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,重慶 400054)
高頻數(shù)據(jù)下基于VaR模型的我國金融市場(chǎng)研究
魏正元,霍艷,李文
(重慶理工大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,重慶 400054)
基于“國內(nèi)外市場(chǎng)間的異構(gòu)性質(zhì)”和“傳統(tǒng)的風(fēng)險(xiǎn)模型是否能有效地應(yīng)用于今天的高頻金融市場(chǎng)”的考慮,設(shè)計(jì)了歷史模擬法、蒙特卡洛模擬法和極值理論方法,依次對(duì)深市的“寶安地產(chǎn)、長(zhǎng)江證券”與滬市的“滬深300”3支個(gè)股進(jìn)行對(duì)比研究,并進(jìn)行失敗率檢驗(yàn)。結(jié)果發(fā)現(xiàn):國內(nèi)兩大金融市場(chǎng)并無明顯差異;歷史模擬法和蒙特卡洛模擬法并不能對(duì)風(fēng)險(xiǎn)值進(jìn)行有效的估計(jì);極值理論可以對(duì)國內(nèi)高頻金融市場(chǎng)進(jìn)行有效的風(fēng)險(xiǎn)度量。同時(shí)實(shí)驗(yàn)的結(jié)果也反向論證了模型適用的前提條件的重要性。
歷史模擬法;蒙特卡洛模擬法;在值風(fēng)險(xiǎn)(VaR);極值理論
由于金融市場(chǎng)的信息對(duì)證券市場(chǎng)價(jià)格運(yùn)動(dòng)過程的影響是連續(xù)性的,所以數(shù)據(jù)采集的離散程度決定了獲取信息的程度。然而,隨著存儲(chǔ)和計(jì)算技術(shù)的快速發(fā)展,各數(shù)據(jù)庫可以采集到以小時(shí)、分鐘或秒為單位的頻率越來越高的數(shù)據(jù),即高頻數(shù)據(jù)。在研究高頻金融數(shù)據(jù)的具體應(yīng)用時(shí),對(duì)數(shù)據(jù)的正確認(rèn)識(shí)和使用往往以數(shù)據(jù)本身的統(tǒng)計(jì)特征為基本前提。由于金融高頻數(shù)據(jù)與低頻數(shù)據(jù)(即以日、周、月等為采集頻率的數(shù)據(jù))間特性的差異,使先前一些基于低頻數(shù)據(jù)上的“短期價(jià)格波動(dòng)服從高斯隨機(jī)游程的古典經(jīng)濟(jì)假設(shè)”、“不同方法的風(fēng)險(xiǎn)估計(jì)模型”等理論都不再適用。因此,對(duì)高頻金融數(shù)據(jù)的度量與應(yīng)用日益得到國內(nèi)外更多學(xué)者的高度重視,其中代表人物是Andersen和Engle。
本文基于金融市場(chǎng)中微觀結(jié)構(gòu)理論的視角,采集國內(nèi)流動(dòng)性高的5 min高頻交易金融資產(chǎn)價(jià)格數(shù)據(jù),設(shè)計(jì)最優(yōu)的抽樣方案,并基于上述模型做實(shí)證分析。利用參數(shù)方法(歷史模擬法和蒙特卡洛模擬法)和非參方法(極值理論)對(duì)在險(xiǎn)價(jià)值(VaR)進(jìn)行估計(jì),比較滬市的“滬深300”和深市的“長(zhǎng)江證券、寶安地產(chǎn)”的在值風(fēng)險(xiǎn)(VaR)走勢(shì),觀測(cè)我國2個(gè)金融市場(chǎng)的異同性。
1VaR模型的定義
在險(xiǎn)價(jià)值(VaR)是指在市場(chǎng)的正常波動(dòng)下,某一金融資產(chǎn)可能發(fā)生的最大損失。數(shù)學(xué)描述是在給定的置信度水平下,某一金融資產(chǎn)在未來特定持有期內(nèi)的最大可能損失,表達(dá)式為
其中:α是置信水平(95%或99%);L是在持有期Δt內(nèi)的資產(chǎn)損失。由式(1)可知:要確定金融機(jī)構(gòu)的在險(xiǎn)價(jià)值(VaR)模型,首先要確定持有期間的長(zhǎng)短、置信區(qū)間的大小、觀察期間3個(gè)指標(biāo)。本文對(duì)3個(gè)指標(biāo)的選取依據(jù)如下:
1)持有期間的長(zhǎng)短。即確定在某時(shí)間段內(nèi)的持有資產(chǎn)的最大損失值。本文選擇的持有期以每日為周期,原因如下:首先,因?yàn)? min的高頻數(shù)據(jù)具有“包含的市場(chǎng)信息量和市場(chǎng)微觀結(jié)構(gòu)噪聲影響最小”達(dá)到均衡的優(yōu)越性質(zhì);其次,“持有期的選擇依據(jù)所持有資產(chǎn)的特點(diǎn)來確定,對(duì)于流動(dòng)性很強(qiáng)的交易往往以每日為周期計(jì)算在險(xiǎn)價(jià)值(VaR)”。
2)置信區(qū)間的大小。置信區(qū)間的選擇在一定程度上反映的是金融機(jī)構(gòu)對(duì)風(fēng)險(xiǎn)的不同偏好。較大的置信水平意味著對(duì)風(fēng)險(xiǎn)比較厭惡,希望可以達(dá)到把握性較大的預(yù)測(cè)結(jié)果,使模型對(duì)于極端事件的預(yù)測(cè)準(zhǔn)確性較高?;诒疚牡难芯磕康?,選取95%的置信水平進(jìn)行預(yù)測(cè)。
3)觀測(cè)期間。針對(duì)給定持有期限的回報(bào)的波動(dòng)性和關(guān)聯(lián)性考察的整體時(shí)間長(zhǎng)度,是整個(gè)數(shù)據(jù)選取的時(shí)間范圍。例如選擇對(duì)某資產(chǎn)組合在未來6個(gè)月或是1年的觀察期間內(nèi),考察其每周回報(bào)率的波動(dòng)性(風(fēng)險(xiǎn))。這種選擇要在歷史數(shù)據(jù)的可能性和市場(chǎng)發(fā)生結(jié)構(gòu)性變化的危險(xiǎn)之間進(jìn)行權(quán)衡。為克服商業(yè)循環(huán)等周期性變化的影響,歷史數(shù)據(jù)越長(zhǎng)越好,但是時(shí)間越長(zhǎng),收購兼并等市場(chǎng)結(jié)構(gòu)性變化的可能性越大,歷史數(shù)據(jù)因而越難以反映現(xiàn)實(shí)和未來的情況。巴塞爾銀行監(jiān)管委員會(huì)目前要求的觀察期間為1年,因此本文選擇的觀測(cè)期間是2013年1月到12月。
1.2 歷史模擬法
Butler和Schachter[1]提出了基于核估計(jì)的歷史模擬法來估計(jì)VaR值,適用于觀測(cè)周期較長(zhǎng)的情況。Cabedo和Moya[2]提出一般的歷史模擬法,它是通過歷史模擬法基于誤差項(xiàng)對(duì)VaR值進(jìn)行估計(jì)。以上2種方法都給不同時(shí)期的觀測(cè)值賦予了相同的權(quán)重,即默認(rèn)組合的歷史收益率是獨(dú)立同分布的。近年來,Asfaha[3]通過對(duì)數(shù)據(jù)的處理,結(jié)合bootstrap方法對(duì)歷史模擬法進(jìn)行了改進(jìn)。陳玉峰、孫洪祥[4]在應(yīng)用歷史模擬法計(jì)算在值風(fēng)險(xiǎn)(VaR)的基礎(chǔ)上,采取增加權(quán)重和波動(dòng)率等提高在值風(fēng)險(xiǎn)(VaR)的估計(jì)精度。但歷史模擬法的原理之一是:假定以后事件的發(fā)生是之前事件發(fā)生的重復(fù),并沒有考慮到之前事件沒有發(fā)生的情況,也忽略了由于數(shù)據(jù)采集頻率低、數(shù)據(jù)不足而導(dǎo)致的對(duì)在值風(fēng)險(xiǎn)(VaR)的不精確估計(jì)。
歷史模擬法原理:根據(jù)觀測(cè)計(jì)算得到的股市收益率變化值模擬股市的未來損益分布,再通過分位數(shù)給出一定置信度下的在險(xiǎn)價(jià)值(VaR)估計(jì)。
歷史模擬法的優(yōu)點(diǎn):概念直觀,不需要對(duì)資產(chǎn)收益的分布特征做假設(shè),易操作。缺點(diǎn):假定未來的概率分布與金融收益的歷史數(shù)據(jù)的概率分布完全一致,顯然這一假定違背了金融波動(dòng)的易變性;歷史模擬法需要大量的觀測(cè)數(shù)據(jù)(通常至少為1 500個(gè)),否則會(huì)嚴(yán)重地低估風(fēng)險(xiǎn)。
1.3 蒙特卡洛模擬法
蒙特卡洛模擬法是一種隨機(jī)模擬法,它用通過觀測(cè)值的歷史波動(dòng)產(chǎn)生觀測(cè)值的未來波動(dòng)的可能路徑,同時(shí)克服了歷史模擬法正態(tài)性假設(shè)的缺陷。王春峰、萬海輝提出了“馬爾科夫鏈蒙特卡洛模擬法”來估計(jì)在值風(fēng)險(xiǎn)(VaR)。
蒙特卡洛法原理:根據(jù)觀測(cè)數(shù)據(jù)本身的特征來模擬一系列符合要求的偽隨機(jī)數(shù),對(duì)這些偽隨機(jī)數(shù)通過歷史模擬過程得到在險(xiǎn)價(jià)值(VaR)。
蒙特卡洛模擬法的優(yōu)點(diǎn):不需要對(duì)金融收益的歷史數(shù)據(jù)的概率分布做假設(shè),可以較好地處理非線性和厚尾的問題。缺點(diǎn):計(jì)算量大,對(duì)隨機(jī)模型存在依賴,若模型有缺陷,會(huì)導(dǎo)致在險(xiǎn)價(jià)值(VaR)嚴(yán)重失真。
1.4 極值理論
傳統(tǒng)的非參數(shù)風(fēng)險(xiǎn)度量方法往往忽略了分布的尾部信息,極值理論可以表現(xiàn)出良好的尾部時(shí)間預(yù)測(cè)能力。longin[5]通過極值理論對(duì)S&P500的多頭和空頭頭寸進(jìn)行在值風(fēng)險(xiǎn)(VaR)估計(jì),通過與置信水平為5%的正態(tài)分布假設(shè)下的在值風(fēng)險(xiǎn)(VaR)相比,發(fā)現(xiàn)其具有更好的精確度。Ramazan 和Faruk[6]應(yīng)用極值理論方法估計(jì)9個(gè)不同市場(chǎng)的股票在值風(fēng)險(xiǎn)(VaR),發(fā)現(xiàn)在值風(fēng)險(xiǎn)(VaR)的估計(jì)精度隨著置信水平的升高而增加。
無論是歷史模擬法還是蒙特卡洛模擬法,得到的在險(xiǎn)價(jià)值(VaR)都是基準(zhǔn)值,而未能考慮尾部分布,而極值理論將研究重心放到了收益的尾部而非整個(gè)分布。
原理:將收益率分成多個(gè)子區(qū)間,然后選取各個(gè)子區(qū)間的最大和最小值,得到極值序列,依據(jù)最值服從廣義極值分布(GEV)理論進(jìn)行建模。缺點(diǎn):適用于收益率集中的最值,會(huì)損失較多的尾部信息。
主要針對(duì)“不同的風(fēng)險(xiǎn)度量模型在中國市場(chǎng)的表現(xiàn)差異”進(jìn)行實(shí)證分析,分別通過歷史模擬法、蒙特卡洛模擬法、極值方法對(duì)在值風(fēng)險(xiǎn)(VaR)進(jìn)行估計(jì)。通過R軟件繪制出不同條件下2013年的實(shí)時(shí)在值風(fēng)險(xiǎn),使風(fēng)險(xiǎn)值可視化。通過實(shí)證分析,對(duì)比國內(nèi)2個(gè)金融市場(chǎng)的異同。
2.1 歷史模擬法
首先,根據(jù)股票每5 min的交易數(shù)據(jù)計(jì)算第d 天t時(shí)刻的股票收益率Rd,t=ln( Pd,(t+1)Pd,t);將每天的收益率按升序排列,根據(jù)損益分布計(jì)算出置信度分位數(shù)為95%的VaR值。
其次,本文中每支股每天有48個(gè)收益率觀測(cè)值,第24觀測(cè)值的最大損失值即為95%的置信度所對(duì)應(yīng)的分位數(shù)。2013年來自2個(gè)市場(chǎng)的3支股票238天的在值風(fēng)險(xiǎn)(VaR)如圖1所示。
由圖1可以發(fā)現(xiàn)(VaR的負(fù)值情況代表收益):實(shí)際的損失值幾乎全部大于VaR值的個(gè)數(shù),即所估計(jì)的在值風(fēng)險(xiǎn)(VaR)估計(jì)幾乎全部失效。以下采用Kupiec失敗率檢驗(yàn)方法對(duì)在值風(fēng)險(xiǎn)(VaR)的有效性進(jìn)行檢驗(yàn),結(jié)果見表1。
圖1 2013年3支股票238天的在值風(fēng)險(xiǎn)(VaR)
表1 對(duì)模擬結(jié)果的Kupiec失敗率檢驗(yàn)
由表1知:經(jīng)過Kupiec失敗率檢驗(yàn),發(fā)現(xiàn)在95%的置信水平下LR的值大于臨界值。因此,可以初步認(rèn)定對(duì)于以上3支樣本,歷史模擬法并沒有完全通過Kupiec失敗率檢驗(yàn),沒能對(duì)3支樣本的風(fēng)險(xiǎn)值做出有效的度量。
2.2 蒙特卡洛模擬法
Shapiro-Wilk檢驗(yàn)是用順序量W檢驗(yàn)分布的正態(tài)性,檢驗(yàn)樣本容量為3~5 000,經(jīng)研究Shapiro-Wilk效果最高。
首先,通過每天觀測(cè)值的直方圖觀察歷史數(shù)據(jù)的分布特征,判斷是否符合正態(tài)分布。本文選取“滬深300”進(jìn)行Shapiro-Wilk檢驗(yàn)(表2),判斷“滬深300”是否符合正態(tài)分布。結(jié)果顯示W(wǎng)統(tǒng)計(jì)量的值均接近于0,前2個(gè)交易日的p值均大于0.05,第5個(gè)交易日的p小于0.05,即前兩日的數(shù)據(jù)符合正態(tài)分布,第5日的數(shù)據(jù)并不符合正態(tài)分布,即本文滬深300的237個(gè)日內(nèi)交易數(shù)據(jù)并非全部符合正態(tài)分布,但目的之一是“利用高頻數(shù)據(jù)選取非參數(shù)方法對(duì)國內(nèi)兩大金融市場(chǎng)進(jìn)行比較分析”,所以下文僅基于蒙特卡洛方法進(jìn)行在險(xiǎn)值(VaR)計(jì)算。一者對(duì)結(jié)果進(jìn)行Kupiec失敗率檢驗(yàn),從反面論證正態(tài)性條件對(duì)蒙特卡洛方法計(jì)算在險(xiǎn)值的必要性;一者應(yīng)用其結(jié)果與本文其他的方法所計(jì)算的在險(xiǎn)值進(jìn)行量化比較,給出定性的分析,判斷國內(nèi)兩大金融市場(chǎng)的異同性,尋找最合適在險(xiǎn)值模型。
表2 對(duì)‘滬深300’的Shapiro-Wilk檢驗(yàn)
其次,通過R語言函數(shù)rnorm()產(chǎn)生每天相應(yīng)的10 000個(gè)隨機(jī)數(shù)。通過歷史模擬法對(duì)以上隨機(jī)數(shù)進(jìn)行分析,得到2013年來自2個(gè)市場(chǎng)的4支股票238天的在險(xiǎn)值(VaR),如圖2所示。
圖2 蒙特卡洛模擬下的在值風(fēng)險(xiǎn)(VaR)
圖2表明:就方法而言,采用蒙特卡羅模擬法計(jì)算而得的在值風(fēng)險(xiǎn)要小于歷史模擬法計(jì)算得到的在值風(fēng)險(xiǎn);就市場(chǎng)而言,初步識(shí)別國內(nèi)2個(gè)市場(chǎng)并沒有極大的差異性。通過Kupiec失敗率檢驗(yàn)法對(duì)蒙特卡羅模擬法的有效性進(jìn)行檢驗(yàn),結(jié)果見表3。
表3 對(duì)蒙特卡洛模擬結(jié)果的Kupiec失敗率檢驗(yàn)
結(jié)果表明:蒙特卡洛模擬法對(duì)3支股票在險(xiǎn)值的估計(jì)的失敗次數(shù)均大于歷史模擬法所得在險(xiǎn)值的失敗次數(shù),從側(cè)面證實(shí)了采用蒙特卡羅模擬進(jìn)行在險(xiǎn)值(VaR)估計(jì)時(shí)正態(tài)性條件的必要性。
2.3 極值理論法
無論是應(yīng)用歷史模擬法還是蒙特卡洛模擬法,得到的VaR值都是基準(zhǔn)值,沒能考慮尾部分布,而極值理論將研究重心放到了收益的尾部而非整個(gè)分布。
原理:首先假設(shè)Rd,t的值是彼此獨(dú)立且來自于同一分布FXn(x)=1-(1-FR(x))n,因此易知FX∞=0,即當(dāng)n→∞Xn的極限分布是收斂的,設(shè)存在規(guī)模參數(shù)αn>0,位置參數(shù)βn>0使(Xn-βn)αn=X*為標(biāo)準(zhǔn)化且非收斂,依據(jù)Gumbel的定理:標(biāo)準(zhǔn)化后的極值存在的極限分布形式為:
綜合以上分布知:所有不同的收益分布有相同的極限分布(區(qū)別僅在于參數(shù)不同)。
首先,根據(jù)極大似然估計(jì)法估計(jì)廣義極值分布的參數(shù)值(見表4),從而進(jìn)一步估計(jì)極值的漸近分布。
表4 廣義極值分布的參數(shù)估計(jì)值
其次,依據(jù)極大值序列的廣義極值分布計(jì)算得到
最后,計(jì)算出3支股票的在值風(fēng)險(xiǎn)(VaR),如圖3所示。
圖3 極值理論下的在值風(fēng)險(xiǎn)(VaR)
由圖3可見:用極值理論計(jì)算得出的在值風(fēng)險(xiǎn)(VaR)要遠(yuǎn)遠(yuǎn)大于用蒙特卡洛和歷史模擬法計(jì)算得出的在值風(fēng)險(xiǎn)(VaR)。因此,如果在銀行內(nèi)部采用“極值理論”計(jì)算在險(xiǎn)值,則會(huì)提高銀行風(fēng)險(xiǎn)管理的資本,進(jìn)一步從側(cè)面證實(shí)了極值方法在銀行界采用較少的原因。下面通過Kupiec失敗率檢驗(yàn)方法對(duì)長(zhǎng)期頭寸VaR值和短期頭寸VaR值的有效行分別進(jìn)行檢驗(yàn),結(jié)果見表5。
表5 對(duì)極值理論結(jié)果的Kupiec失敗率檢驗(yàn)
檢驗(yàn)結(jié)果顯示:在95%的置信水平下LR統(tǒng)計(jì)量大于卡方統(tǒng)計(jì)量的臨界值。說明:①對(duì)同一支股票而言,不論是長(zhǎng)期頭寸還是短期頭寸,它們的失敗次數(shù)是一致的;②同蒙特卡洛和歷史模擬法相比,極值理論計(jì)算所得的VaR值通過了Kupiec失敗率檢驗(yàn),初步說明與歷史模擬法和蒙特卡洛模擬法相比,應(yīng)用極值理論可以利用高頻數(shù)據(jù)對(duì)我國市場(chǎng)的在值風(fēng)險(xiǎn)(VaR)進(jìn)行有效的估計(jì)。
為了判別國內(nèi)金融市場(chǎng)間是否存在異構(gòu)型,本文選取了國內(nèi)滬深兩大金融市場(chǎng)2013年的“滬深300”、“寶安地產(chǎn)”和“長(zhǎng)江證券”的5 min高頻數(shù)據(jù)作為研究對(duì)象,選取經(jīng)典的歷史模擬法、蒙特卡洛模擬法和極值理論依次對(duì)3支股票的在值風(fēng)險(xiǎn)進(jìn)行預(yù)測(cè),并通過Kupiec失敗率檢驗(yàn)對(duì)國內(nèi)金融市場(chǎng)異構(gòu)性做出判斷。初步得到以下結(jié)果:歷史模擬法并不能很有效地應(yīng)用于高頻金融數(shù)據(jù)的國內(nèi)市場(chǎng);為了考查國內(nèi)市場(chǎng)高頻數(shù)據(jù)的正態(tài)性,采用Shapiro-Wilk檢驗(yàn)方法,依次對(duì)每天的歷史數(shù)據(jù)進(jìn)行分析,發(fā)現(xiàn)全年每天的交易數(shù)據(jù)并不能全部顯示出正態(tài)性,從而初步斷定蒙特卡洛模擬法對(duì)VaR值得估計(jì)會(huì)失效,Kupiec失敗率檢驗(yàn)的結(jié)果驗(yàn)證了這一斷定。為了判定是否數(shù)據(jù)存有后尾特性,本文在蒙特卡洛模擬之后應(yīng)用了極值理論對(duì)3支個(gè)股數(shù)據(jù)進(jìn)行分析。結(jié)果顯示:通過極值理論計(jì)算的VaR值能夠大部分通過失敗率檢驗(yàn),可有效地對(duì)國內(nèi)市場(chǎng)的在值風(fēng)險(xiǎn)進(jìn)行度量。
[1]Butler J S,Schachter B.Estimating value-at-risk with a precision measure by combining kernel estimation with historical simulation[J].Review of Derivatives Research. Springer,1998,1(4):371-390.
[2]David J,Moya.Estimating oil price“Value at Risk”using the historical simulation approach[J].Energy Economics,2003,25(3):239-253.
[3]Asfaha T,Desmond A F,Hailu G.Statistical Evaluation of Value at Risk Models for Estimating Agricultural Risk [J].Journal of Statistical and Econometric Methods,2014,3(1):13-34.
[4]陳玉峰,孫洪祥,溫巧燕.基于歷史模擬法的VaR計(jì)算及其優(yōu)化[J].中國科技論文在線精品論文,2011,4 (13):1231-1237.
[5]Longin F M,F(xiàn)rom value at risk to stress testing:The extreme value approach[J].Journal of Banking and Finance,2000,24(7):1097-1103.
[6]Genay R,Seluk F.Extreme value theory and value-atrisk:relative performance in emerging markets[J].International Journal of Forecasting,2004,20(2):287-303.
(責(zé)任編輯 劉舸)
Study on the High Frequency Data of the Financial Market in China by Using the VaR Method
WEI Zheng-yuan,HUO Yan,LI Wen
(School of Mathematics and Statistics,
C
hongqing University of Technology,Chongqing 400054,China)
With the rapid development of storage and computing technology,high frequency data with real-time transaction acquisition market become true.In order to study the relations between the two major domestic financial market,we choose two stocks named as“Changjian zhengquan”and“Bao’an dichan”from Shanghai market and“Hushen 300”from Shenzhen market.Monte Carlo simulation method,extreme value and historical simulation method will be selected to analysis each of the above three stock,then tested the failure rate.The results show that there is no remarkable difference between Chinese two financial market;the historical simulation method and Monte Carlo simulation method is not effective for risk estimation;extreme value theory can effectively measure the risk of domestic high-frequency Market.
historical simulation method;Monte Carlo simulation method;value at risk(VaR);extreme value theory
O21
A
1674-8425(2014)08-0126-06
10.3969/j.issn.1674-8425(z).2014.08.026
2014-03-28
重慶市自然科學(xué)基金資助項(xiàng)目(cstc2012jjA00018);重慶市教委科學(xué)技術(shù)研究項(xiàng)目(KJ130810);重慶市教委高教研究項(xiàng)目(1203053)
魏正元(1975—),男,博士,副教授,主要從事應(yīng)用概率統(tǒng)計(jì)、金融統(tǒng)計(jì)、金融數(shù)學(xué)研究。
魏正元,霍艷,李文.高頻數(shù)據(jù)下基于VaR模型的我國金融市場(chǎng)研究[J].重慶理工大學(xué)學(xué)報(bào):自然科學(xué)版,2014(8):126-131.
format:WEI Zheng-yuan,HUO Yan,LI Wen.Study on the High Frequency Data of the Financial Market in China by Using the VaR Method[J].Journal of Chongqing University of Technology:Natural Science,2014(8): 126-131.
重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))2014年8期