李冰娜，惠曉峰，李連江

1 東北大學(xué)秦皇島分校 經(jīng)濟(jì)學(xué)院，河北 秦皇島 066004 2 哈爾濱工業(yè)大學(xué) 管理學(xué)院，哈爾濱 150001 3 東北大學(xué)秦皇島分校 控制工程學(xué)院，河北 秦皇島 066004

基于蒙特卡洛RMT去噪法小股票組合風(fēng)險(xiǎn)優(yōu)化研究

李冰娜1，惠曉峰2，李連江3

1 東北大學(xué)秦皇島分校 經(jīng)濟(jì)學(xué)院，河北 秦皇島 066004 2 哈爾濱工業(yè)大學(xué) 管理學(xué)院，哈爾濱 150001 3 東北大學(xué)秦皇島分校 控制工程學(xué)院，河北 秦皇島 066004

在Markowitz證券投資組合理論的框架下，證券收益協(xié)方差矩陣往往受到“維數(shù)災(zāi)禍”的影響而充斥噪聲，這給Markowitz證券投資組合的構(gòu)建及其風(fēng)險(xiǎn)的優(yōu)化帶來(lái)了嚴(yán)重的困擾。在對(duì)證券收益協(xié)方差矩陣去噪進(jìn)而實(shí)現(xiàn)證券投資組合風(fēng)險(xiǎn)優(yōu)化方面，基于隨機(jī)矩陣?yán)碚?RMT)的去噪方法是一種非常具有優(yōu)勢(shì)的有效方法。

從說明Markowitz股票投資組合風(fēng)險(xiǎn)的含義入手，以投資組合風(fēng)險(xiǎn)預(yù)測(cè)的準(zhǔn)確率衡量投資組合風(fēng)險(xiǎn)的優(yōu)劣，繼而對(duì)用于股票收益協(xié)方差矩陣的原有RMT去噪法的原理進(jìn)行闡釋。在此基礎(chǔ)上分析認(rèn)為小股票組合條件下原有RMT去噪法因噪聲特征值邊界界定誤差而會(huì)產(chǎn)生組合風(fēng)險(xiǎn)優(yōu)化作用下降的問題，為解決該問題，采用蒙特卡洛模擬法確定收益協(xié)方差矩陣的最大噪聲特征值，從而建立蒙特卡洛RMT去噪法。對(duì)LCPB去噪法、PG+去噪法和KR去噪法等原有RMT去噪法和蒙特卡洛RMT去噪法的組合風(fēng)險(xiǎn)優(yōu)化作用開展實(shí)證研究。

研究結(jié)果表明，各種RMT方法對(duì)相等加權(quán)協(xié)方差矩陣去噪前后的組合風(fēng)險(xiǎn)結(jié)果與指數(shù)加權(quán)協(xié)方差矩陣的組合風(fēng)險(xiǎn)結(jié)果非常類似；收益協(xié)方差矩陣未去噪時(shí)，隨著股票數(shù)量減小，組合風(fēng)險(xiǎn)的優(yōu)化水平越來(lái)越高，這主要是因?yàn)槭找鎱f(xié)方差矩陣噪聲不斷減小；當(dāng)對(duì)收益協(xié)方差矩陣去噪時(shí)，隨著股票數(shù)量的減小，原有RMT去噪法比蒙特卡洛RMT去噪法具有的去噪優(yōu)勢(shì)因原有方法對(duì)噪聲特征值邊界界定誤差的增大而呈減小趨勢(shì)，因此，組合中股票降至較低數(shù)量時(shí)，原有RMT去噪法的組合風(fēng)險(xiǎn)優(yōu)化作用開始變得低于蒙特卡洛RMT去噪法，這說明蒙特卡洛RMT去噪法是一種解決小股票組合風(fēng)險(xiǎn)優(yōu)化條件下原有RMT去噪法效力下降問題的有效方法。

新的蒙特卡洛RMT去噪法提供了證券投資組合權(quán)重分配的具體方法，提高了組合風(fēng)險(xiǎn)預(yù)測(cè)的準(zhǔn)確率，對(duì)設(shè)計(jì)合理的證券投資策略和優(yōu)化證券投資風(fēng)險(xiǎn)具有參考意義。

Markowitz投資組合；組合風(fēng)險(xiǎn)優(yōu)化；小組合；RMT去噪法；蒙特卡洛模擬

1 引言

Markowitz投資組合理論表明對(duì)證券投資組合進(jìn)行決策往往需要利用證券收益多元波動(dòng)率(協(xié)方差矩陣)[1]。該矩陣一般只能由證券收益率的市場(chǎng)觀測(cè)值來(lái)估計(jì)，但是實(shí)際觀測(cè)值序列長(zhǎng)度的有限(即所謂的“維數(shù)災(zāi)禍”)會(huì)使估計(jì)出的矩陣存在噪聲，尤其是當(dāng)證券數(shù)量較大時(shí)，估計(jì)噪聲甚至?xí)鹬鲗?dǎo)作用[2]。不難想象，在此基礎(chǔ)上測(cè)算的證券投資組合難免會(huì)對(duì)人們的投資決策造成誤導(dǎo)，并且導(dǎo)致組合未來(lái)風(fēng)險(xiǎn)的預(yù)測(cè)也是不準(zhǔn)確的。因此，對(duì)證券收益協(xié)方差矩陣(下文簡(jiǎn)寫為協(xié)方差矩陣)去噪是實(shí)現(xiàn)證券投資組合風(fēng)險(xiǎn)優(yōu)化的前提[3]。

近年來(lái)，源自于經(jīng)濟(jì)物理學(xué)的隨機(jī)矩陣?yán)碚?random matrix theory, RMT)[4]的去噪法憑借其獨(dú)特的優(yōu)勢(shì)在證券協(xié)方差矩陣去噪，進(jìn)而實(shí)現(xiàn)證券投資組合風(fēng)險(xiǎn)優(yōu)化方面發(fā)揮了巨大的作用[5]。然而，相對(duì)于幾百只乃至上千只證券的組合而言，RMT去噪法對(duì)較小數(shù)量證券的組合，特別是百只以下證券組合的協(xié)方差矩陣去噪的作用很小，因此，對(duì)這些規(guī)模較小的組合風(fēng)險(xiǎn)優(yōu)化效果很差。如何改進(jìn)RMT去噪法對(duì)小股票組合協(xié)方差矩陣的去噪效力，以促進(jìn)小股票組合風(fēng)險(xiǎn)的優(yōu)化是值得密切關(guān)注的問題。在美國(guó)，大多數(shù)基金的投資組合囊括了超過100種證券，但個(gè)人投資者的投資組合規(guī)模一般都不大。與美國(guó)相比，中國(guó)證券投資組合的規(guī)模更小，主要原因是中國(guó)資本市場(chǎng)以中小投資者為主，機(jī)構(gòu)投資者所占比例較低，且投資規(guī)模較小[6]?？梢?，現(xiàn)實(shí)經(jīng)濟(jì)生活中小規(guī)模投資組合的使用是非常多的。因此，研究解決小組合條件下RMT去噪法的誤差問題，以期實(shí)現(xiàn)小股票組合風(fēng)險(xiǎn)的優(yōu)化，對(duì)指導(dǎo)證券市場(chǎng)投資風(fēng)險(xiǎn)管理和滿足日益增長(zhǎng)的投資理財(cái)需求有重要的現(xiàn)實(shí)意義。

2 相關(guān)研究評(píng)述

學(xué)者們基于RMT對(duì)金融收益多元波動(dòng)率[7]的噪聲進(jìn)行了大量研究，發(fā)現(xiàn)金融協(xié)方差矩陣的絕大部分特征值處于RMT預(yù)測(cè)的噪聲帶內(nèi)，它們代表協(xié)方差矩陣中的噪聲，而僅有個(gè)別特征值大于最大隨機(jī)特征值，它們包含資產(chǎn)間相關(guān)的信息[8-9]。針對(duì)金融收益多元波動(dòng)率中的噪聲，LCPB去噪法、PG+去噪法和KR去噪法等傳統(tǒng)RMT去噪法被相繼提出，并在證券投資組合風(fēng)險(xiǎn)優(yōu)化方面進(jìn)行了大量的應(yīng)用。也有個(gè)別學(xué)者對(duì)這幾種方法進(jìn)行了適當(dāng)?shù)淖儞Q，不過這方面的探索較少。

最先產(chǎn)生的RMT去噪法是LCPB去噪法，它是由LALOUX et al.[10]建立的，對(duì)該方法的實(shí)證結(jié)果表明去噪使組合風(fēng)險(xiǎn)預(yù)測(cè)誤差大大降低；LEE[11]利用隨機(jī)相關(guān)矩陣譜的中心極限定理識(shí)別出S&P 500的233只代表性股票在1997年1月2日至2001年5月11日的日收益相關(guān)矩陣的噪聲帶后，用LCPB去噪法對(duì)該矩陣進(jìn)行了去噪。RMT去噪法在證券組合選擇方面的作用明顯優(yōu)于其他方法。PAFKA et al.[12]使用模擬方法對(duì)比了樣本相關(guān)矩陣、基于單一指數(shù)法的相關(guān)矩陣和基于RMT的相關(guān)矩陣用于證券組合選擇的表現(xiàn)，結(jié)果顯示，在減弱噪聲影響方面，RMT去噪法的確是一種十分有效的降維方法；TOLA et al.[13]將聚類算法和RMT去噪法用于證券組合優(yōu)化，并從可靠性、風(fēng)險(xiǎn)和有效資產(chǎn)數(shù)量3個(gè)方面對(duì)二者的效果進(jìn)行比較，他們的結(jié)論是在理想化的條件下，即對(duì)未來(lái)收益和波動(dòng)率有完全的預(yù)測(cè)能力和允許賣空時(shí)，聚類法的可靠性優(yōu)于RMT方法，單一連接聚類法的風(fēng)險(xiǎn)水平高于平均連接聚類法和RMT去噪法，而它的有效資產(chǎn)數(shù)量明顯低于后兩種方法。

除LCPB去噪法外，PLEROU et al.[14]也提出了一種基于RMT的去噪法，即PG+去噪法。PG+法與LCPB法的原理非常相似，它們都對(duì)噪聲特征值的替代進(jìn)行簡(jiǎn)單、粗糙的處理。為實(shí)現(xiàn)RMT去噪法算法的改進(jìn)，SHARIFI et al.[15]提出了基于RMT和協(xié)方差矩陣特征向量Krzanowski穩(wěn)定性的KR去噪法。雖然KR去噪法具有非常明顯的理論優(yōu)勢(shì)，但它在實(shí)際應(yīng)用中的組合風(fēng)險(xiǎn)優(yōu)化效果尚需進(jìn)一步驗(yàn)證。對(duì)此，DALY et al.[16]使用LCPB去噪法、PG+去噪法和KR去噪法對(duì)資產(chǎn)數(shù)量為100、250和432的組合協(xié)方差進(jìn)行去噪，發(fā)現(xiàn)KR法能最大程度地降低平均已實(shí)現(xiàn)風(fēng)險(xiǎn)。通過觀測(cè)已實(shí)現(xiàn)風(fēng)險(xiǎn)評(píng)價(jià)組合風(fēng)險(xiǎn)優(yōu)化效果的研究還有很多，如KAROUI[17]應(yīng)用RMT方法研究Markowitz組合的已實(shí)現(xiàn)風(fēng)險(xiǎn)；RUBIO et al.[18]基于RMT研究高維樣本下最小方差組合的已實(shí)現(xiàn)方差；WAKAI et al.[19]分析隨機(jī)矩陣方法構(gòu)建的均值-方差模型的最優(yōu)組合的典型特點(diǎn)。除美國(guó)等發(fā)達(dá)國(guó)家證券市場(chǎng)外，一些學(xué)者基于RMT開展了發(fā)展中經(jīng)濟(jì)體證券市場(chǎng)特征譜的類似研究。TANG et al.[20]發(fā)現(xiàn)上海證券市場(chǎng)523只A股股票和深圳證券市場(chǎng)445只A股股票的收益相關(guān)矩陣存在大量噪聲，基于RMT對(duì)它們?nèi)ピ?，結(jié)果表明去噪后構(gòu)建的投資組合優(yōu)于原始樣本相關(guān)矩陣和Kendallτ相關(guān)矩陣建立的投資組合；OH et al.[21]基于RMT研究韓國(guó)證券市場(chǎng)上證券相關(guān)矩陣的統(tǒng)計(jì)特性及其對(duì)Markowitz組合權(quán)重的影響，發(fā)現(xiàn)由去噪前后相關(guān)矩陣計(jì)算得到的組合風(fēng)險(xiǎn)的熵函數(shù)是一個(gè)冪函數(shù)；韓華等[22]以上海證券市場(chǎng)數(shù)據(jù)為樣本，基于RMT探討金融相關(guān)矩陣特征值的統(tǒng)計(jì)特性，并對(duì)金融相關(guān)矩陣進(jìn)行去噪，以使其更適合于構(gòu)建金融網(wǎng)絡(luò)。隨著RMT在金融領(lǐng)域應(yīng)用的不斷深入，一些研究并非僅單獨(dú)利用RMT方法，而是注重RMT方法與其他理論和方法的有機(jī)整合。ARAI et al.[23]將主成分分析與RMT結(jié)合起來(lái)，以東京股票交易所日收益數(shù)據(jù)為樣本，運(yùn)用所建立的復(fù)雜PCA方法分析股票收益相關(guān)矩陣，結(jié)果表明收益相關(guān)矩陣最大特征值表示一種對(duì)所有股票價(jià)格都會(huì)產(chǎn)生影響的市場(chǎng)力量，可以根據(jù)其他較大特征值對(duì)應(yīng)的特征向量對(duì)股票進(jìn)行分類；YANG et al.[24]應(yīng)用RMT去噪法濾除證券相關(guān)矩陣中的噪聲和共變趨勢(shì)，繼而用加權(quán)多目標(biāo)遺傳算法對(duì)證券進(jìn)行分群；WANG et al.[25]基于RMT研究發(fā)現(xiàn)全球證券指數(shù)間相關(guān)可以由一個(gè)全球因素模型作出解釋，他們進(jìn)一步說明了該模型在降低全球范圍內(nèi)投資組合風(fēng)險(xiǎn)方面的價(jià)值；WANG et al.[26]使用皮爾森相關(guān)系數(shù)法和去趨勢(shì)相關(guān)系數(shù)法構(gòu)建S&P 500指數(shù)的462只成分股票間的相關(guān)矩陣，并采用RMT法分析它們的統(tǒng)計(jì)特性，結(jié)果表明去趨勢(shì)相關(guān)矩陣的一些特點(diǎn)對(duì)最佳組合選擇和風(fēng)險(xiǎn)管理很有幫助。

除將RMT用于資產(chǎn)數(shù)量較大的組合外，也有學(xué)者將其用于小證券組合的研究。CONLON et al.[27]采用LCPB去噪法對(duì)49個(gè)對(duì)沖基金105個(gè)月的月收益協(xié)方差矩陣進(jìn)行去噪，并將去噪后的協(xié)方差矩陣用于構(gòu)建組合基金，結(jié)果表明RMT去噪法對(duì)組合基金風(fēng)險(xiǎn)的改進(jìn)作用較?。籏UMAR et al.[28]利用RMT方法研究20個(gè)金融指數(shù)的相關(guān)和網(wǎng)絡(luò)特性，發(fā)現(xiàn)第二大特征值對(duì)應(yīng)特征向量的成分能夠?qū)⒏髦笖?shù)在正負(fù)方向上分為兩類，這兩類指數(shù)的成分值在金融危機(jī)期間會(huì)向相反方向轉(zhuǎn)換。這些研究表明RMT去噪法用于小證券組合風(fēng)險(xiǎn)優(yōu)化的效果非常不顯著。事實(shí)上，也有一些研究認(rèn)為，當(dāng)將RMT去噪法應(yīng)用于小證券組合時(shí)，甚至無(wú)法產(chǎn)生任何風(fēng)險(xiǎn)優(yōu)化的作用。PAFKA et al.[29]采用RMT方法對(duì)指數(shù)加權(quán)協(xié)方差矩陣去噪，結(jié)果表明當(dāng)組合中資產(chǎn)數(shù)量較小時(shí)，樣本協(xié)方差矩陣的組合風(fēng)險(xiǎn)稍低于采用RMT法去噪的協(xié)方差矩陣；DALY et al.[30]還將LCPB法、PG+法和KR法用于39種外匯資產(chǎn)協(xié)方差矩陣的去噪，結(jié)果表明去噪效果隨資產(chǎn)數(shù)量減少而下降；SANDOVAL et al.[31]發(fā)現(xiàn)用RMT方法去噪并不能產(chǎn)生更好的組合風(fēng)險(xiǎn)預(yù)測(cè)結(jié)果，而將RMT去噪法與消除市場(chǎng)效應(yīng)的回歸模型法結(jié)合起來(lái)卻在62.5%的情況下能帶來(lái)更好的風(fēng)險(xiǎn)預(yù)測(cè)結(jié)果。

與國(guó)外研究相比，中國(guó)學(xué)者基于RMT探討金融市場(chǎng)中投資組合優(yōu)化和風(fēng)險(xiǎn)管理問題的研究很少。羅英等[32]應(yīng)用隨機(jī)矩陣?yán)碚摷胺椒ń馕鲋袊?guó)股票組合協(xié)方差矩陣的信息結(jié)構(gòu)，發(fā)現(xiàn)它存在受市場(chǎng)因素和行業(yè)因素主導(dǎo)的信息結(jié)構(gòu)。韓華等[22]、唐曉清等[33]和惠曉峰等[34]在討論股票協(xié)方差矩陣相關(guān)結(jié)構(gòu)的基礎(chǔ)上，進(jìn)一步開展協(xié)方差矩陣去噪和投資組合構(gòu)建的研究工作。韓華等[22]基于RMT對(duì)金融相關(guān)系數(shù)矩陣去噪，并構(gòu)建了金融網(wǎng)絡(luò)模型；唐曉清等[33]基于RMT改進(jìn)了Markowitz組合投資模型, 結(jié)果表明RMT可以偵測(cè)到“異?！?數(shù)據(jù), 從而改進(jìn)對(duì)模型參數(shù)的估計(jì)；惠曉峰等[34]在Markowitz投資組合理論框架下研究RMT去噪法用于多元GARCH模型的有效性，結(jié)果表明RMT能夠?yàn)槎嘣狦ARCH的降維提供有效的依據(jù)，并且較準(zhǔn)確地確定多元GARCH模型的最佳維度。

從上述研究看，原有RMT去噪法主要包括LCPB法、PG+法和KR法。從這些方法應(yīng)用的對(duì)象看，它們不僅用于對(duì)上百只甚至幾百只資產(chǎn)的大組合進(jìn)行去噪，也用于百只以下資產(chǎn)的小組合；但與大組合的去噪效果相比，它們對(duì)小組合的去噪作用并不明顯，有時(shí)甚至帶來(lái)負(fù)的組合風(fēng)險(xiǎn)優(yōu)化效果。已有研究基于小組合情形下協(xié)方差矩陣噪聲較低的規(guī)律，將RMT去噪法對(duì)小組合作用的不明顯化或“失靈”視為一種十分正常的現(xiàn)象，即這種現(xiàn)象的發(fā)生并不是由RMT去噪法自身的原因?qū)е碌摹Ｈ欢?，?duì)RMT去噪法的原理進(jìn)行細(xì)致的分析發(fā)現(xiàn)，在小組合條件下原有RMT去噪法的算法都有比較嚴(yán)重的缺陷，即最大噪聲特征值的解析解存在較大的計(jì)算誤差，這難免會(huì)對(duì)組合風(fēng)險(xiǎn)優(yōu)化帶來(lái)負(fù)面影響。目前，從RMT去噪法自身存在問題的視角探討它們對(duì)小證券組合風(fēng)險(xiǎn)優(yōu)化作用降低的相關(guān)研究還比較少見。本研究正是從該角度展開研究，闡明Markowitz投資組合風(fēng)險(xiǎn)及其優(yōu)化的有關(guān)基本概念，指出小組合風(fēng)險(xiǎn)優(yōu)化時(shí)原有RMT去噪法的缺陷，即最大噪聲特征值準(zhǔn)確性下降的問題，并分析該問題產(chǎn)生的原因，提出蒙特卡洛RMT去噪法，最后進(jìn)行以驗(yàn)證蒙特卡洛RMT去噪法效力為目的的實(shí)證研究工作。

3 Markowitz投資組合風(fēng)險(xiǎn)及其優(yōu)化

Markowitz均值-方差模型的主要目的是通過對(duì)未來(lái)期的期望收益率和協(xié)方差矩陣的預(yù)測(cè)來(lái)對(duì)未來(lái)期的有效投資組合進(jìn)行預(yù)測(cè)。在實(shí)際應(yīng)用中，人們通常利用樣本協(xié)方差矩陣對(duì)未來(lái)期的真實(shí)協(xié)方差矩陣進(jìn)行預(yù)測(cè)。在此基礎(chǔ)上根據(jù)Markowitz均值-方差模型構(gòu)建的投資組合在本研究中稱為Markowitz投資組合。為了表述方便，本研究將用于預(yù)測(cè)的歷史數(shù)據(jù)所在的時(shí)期稱為預(yù)測(cè)期，將發(fā)生投資的未來(lái)期稱為投資期，Markowitz投資組合在投資期實(shí)現(xiàn)的風(fēng)險(xiǎn)稱為Markowitz投資組合風(fēng)險(xiǎn)。

對(duì)Markowitz投資組合風(fēng)險(xiǎn)預(yù)測(cè)的準(zhǔn)確性越高，則基于該預(yù)測(cè)做出的投資決策就越正確，投資組合給決策者帶來(lái)的利益也就越大。因此，本研究用Markowitz投資組合風(fēng)險(xiǎn)預(yù)測(cè)的準(zhǔn)確率表示Markowitz投資組合風(fēng)險(xiǎn)的優(yōu)化水平。當(dāng)對(duì)Markowitz投資組合風(fēng)險(xiǎn)的預(yù)測(cè)準(zhǔn)確率提高時(shí)，則意味著Markowitz投資組合風(fēng)險(xiǎn)實(shí)現(xiàn)了優(yōu)化。

4 小組合條件下原有RMT去噪法的缺陷

4.1 原有RMT去噪法的原理

到目前為止，還沒有專門針對(duì)小組合而設(shè)計(jì)的RMT去噪法。與大組合相同，能夠用于小組合協(xié)方差的RMT去噪法主要有LCPB法、PG+法和KR法3種，這3種方法對(duì)金融協(xié)方差矩陣的去噪原理類似，都是依據(jù)矩陣的特征分解去噪。

股票樣本協(xié)方差矩陣的計(jì)算方法主要包括相等加權(quán)法和指數(shù)加權(quán)法。假定N個(gè)時(shí)間序列長(zhǎng)度為T、樣本方差為σ2的股票收益序列r1,t，r2,t，…，rN,t的相等加權(quán)樣本協(xié)方差矩陣為Seq，指數(shù)加權(quán)樣本協(xié)方差矩陣為Sex。以Seq和Sex的特征值作為主對(duì)角線元素的對(duì)角矩陣分別記為Deq和Dex，以各特征值對(duì)應(yīng)的特征向量作為相應(yīng)的列形成的矩陣分別為Eeq和Eex，則

(1)

(2)

Seq和Sex的噪聲能夠在它們的特征值中反映出來(lái)。Seq和Sex的一部分特征值能夠反映真實(shí)的相關(guān)信息，但由于用于計(jì)算它們的時(shí)間序列長(zhǎng)度往往是有限(非無(wú)窮)的，因此，Seq和Sex會(huì)受到隨機(jī)因素的干擾而產(chǎn)生一些無(wú)法反映真實(shí)相關(guān)信息的噪聲特征值。由于它們的噪聲特征值具有隨機(jī)性特征，因此，可以依據(jù)相應(yīng)隨機(jī)矩陣的特征值分布來(lái)確定這些噪聲特征值的范圍。

首先說明Seq的噪聲特征值的識(shí)別。Seq對(duì)應(yīng)的隨機(jī)協(xié)方差矩陣為R，即

(3)

其中，A為一個(gè)N·T矩陣，其中元素是均值為0、方差為σ2的獨(dú)立、相同分布的隨機(jī)變量。由于用于計(jì)算R的隨機(jī)序列的個(gè)數(shù)、長(zhǎng)度和方差與用于計(jì)算Seq的時(shí)間序列相同，因此，可以依據(jù)R的特征值分布區(qū)間確定Seq的噪聲特征值的范圍。

(4)

其中，γ為衰減因子，k為其指數(shù)，xi,k、xj,k為均值為0、方差為σ2的正態(tài)分布的隨機(jī)變量，i=1,2,…,N，j=1,2,…,N，k=0,1,…,∞。用于計(jì)算M的隨機(jī)序列的個(gè)數(shù)和方差與計(jì)算Sex的股票收益序列的個(gè)數(shù)和方差對(duì)應(yīng)相等，因此，可以依據(jù)M的特征值分布區(qū)間確定Sex的噪聲特征值的范圍。

接下來(lái)，用[λ-，λ+]表示R和M的特征值區(qū)間，則Seq和Sex的特征值中屬于[λ-，λ+]的部分應(yīng)被識(shí)別為噪聲特征值。然而，由于小于λ-的特征值對(duì)應(yīng)的特征向量并不穩(wěn)定，因此，在應(yīng)用RMT去噪法時(shí)，通常將不大于λ+的特征值都當(dāng)做噪聲特征值而對(duì)其進(jìn)行替換。在替換噪聲特征值時(shí)，LCPB法、PG+法和KR法的做法各不相同。由于λ+界定了噪聲特征值的上限，因此，這里將λ+稱為最大噪聲特征值。由于大于λ+的特征值并不是由噪聲產(chǎn)生的，因此讓它們保持不變。設(shè)對(duì)Deq和Dex中的元素進(jìn)行部分保持和部分替換后形成的矩陣分別為Deq-filtered和Dex-filtered，則對(duì)Seq和Sex去噪后的結(jié)果分別為

(5)

(6)

4.2 最大噪聲特征值的計(jì)算誤差

(7)

其中，Q為股票時(shí)間序列長(zhǎng)度與股票數(shù)量之比。

(8)

(9)

其中，v為函數(shù)F(v)取值為0時(shí)的解。

(10)

(11)

根據(jù)λ+的上述計(jì)算方法，原有RMT去噪法并沒有真正計(jì)算λ+，而是用隨機(jī)協(xié)方差矩陣在極限條件下的最大特征值近似地表示λ+。由于實(shí)際應(yīng)用中N都是有限的，因此，原有RMT去噪法在N→∞時(shí)計(jì)算的λ+必然會(huì)偏離于真正的λ+。不難想象，兩者的偏差會(huì)隨N的減小而增大。當(dāng)N較小時(shí)，對(duì)λ+的計(jì)算準(zhǔn)確性也會(huì)因此降至較低水平。這種在小組合條件下對(duì)λ+的較大計(jì)算誤差會(huì)造成噪聲特征值的錯(cuò)誤識(shí)別，并進(jìn)一步對(duì)噪聲特征值的替換造成不利的影響，導(dǎo)致發(fā)生替換過度或替換不足。無(wú)論是替換過度還是替換不足，都會(huì)令Deq-filtered和Dex-filtered無(wú)法被正確地確定，從而造成RMT方法去噪效力的下降，并進(jìn)一步地對(duì)股票投資組合風(fēng)險(xiǎn)的優(yōu)化帶來(lái)不利的影響。盡管過去多年中RMT去噪法的應(yīng)用較多，不過它們的上述誤差問題卻一直沒有引起學(xué)者們的關(guān)注。

5 蒙特卡洛RMT去噪法的設(shè)計(jì)和實(shí)現(xiàn)

5.1 設(shè)計(jì)思想

為了消除小組合條件下RMT去噪法中λ+準(zhǔn)確性下降的問題對(duì)投資組合風(fēng)險(xiǎn)優(yōu)化的不利影響，本研究用蒙特卡洛模擬法確定λ+，從而設(shè)計(jì)了蒙特卡洛RMT去噪法，簡(jiǎn)稱為MKR去噪法。除了λ+的確定方法外，蒙特卡洛RMT去噪法與原有RMT去噪法的原理相同。蒙特卡洛模擬方法的原理是：根據(jù)隨機(jī)過程，用計(jì)算機(jī)模擬的方法生成時(shí)間序列，并將其作為抽樣結(jié)果，以此計(jì)算參數(shù)估計(jì)值。蒙特卡洛模擬方法對(duì)參數(shù)的估計(jì)精度隨模擬次數(shù)增多而逐漸提高，這保證了蒙特卡洛模擬法一般都能取得較高的估計(jì)精度。由λ+的原始計(jì)算模型可知，λ+的計(jì)算問題具有明顯的概率特征，因此，蒙特卡洛模擬法非常適合于計(jì)算λ+。蒙特卡洛模擬法對(duì)相等加權(quán)協(xié)方差矩陣和指數(shù)加權(quán)協(xié)方差矩陣λ+的計(jì)算方法相同，因此，本研究以相等加權(quán)協(xié)方差矩陣為例介紹用蒙特卡洛模擬法確定λ+的原理。

5.2 蒙特卡洛模型構(gòu)建

(12)

其中，

(13)

(14)

第1步 根據(jù)矩陣A概率模型中各行隨機(jī)變量的分布，在計(jì)算機(jī)上產(chǎn)生隨機(jī)數(shù)。每行產(chǎn)生符合概率模型的T個(gè)隨機(jī)數(shù)，共產(chǎn)生N行，構(gòu)成矩陣A的一次樣本抽樣。

第2步 應(yīng)用A的一次樣本抽樣計(jì)算出矩陣R和它的全部特征值。設(shè)E為單位矩陣，那么|mE-R|=0中m的全部根就是R的全部特征值。設(shè)R的全部特征值為Rλ，Rλ={λ1,λ2,…,λN}，遍歷所有特征值找到最大值λmax1。

第4步λmax的誤差的計(jì)算公式為

ε=λασ

(15)

第5步 取得λ+的最終估計(jì)值。如果ε≤ε0，則第3步得到的λmax就是對(duì)λ+的最終估計(jì)值。否則，進(jìn)行第101模擬，并計(jì)算λmax和其誤差ε。若ε仍大于ε0，則繼續(xù)進(jìn)行第102次模擬。直至獲得ε≤ε0時(shí)的λmax，該λmax就是對(duì)λ+的最終估計(jì)值。

從上面計(jì)算過程能夠看出，蒙特卡洛RMT去噪法是用隨機(jī)模擬法近似計(jì)算λ+的。由于它是在股票數(shù)量N取實(shí)際值的條件下計(jì)算λ+，因此能夠有效解決對(duì)小組合去噪時(shí)原有RMT去噪法因限定N→∞的條件而引起對(duì)λ+計(jì)算準(zhǔn)確性較低的問題，從而在組合資產(chǎn)數(shù)量較少時(shí)取得比原有RMT去噪法更好的去噪和組合風(fēng)險(xiǎn)優(yōu)化效果。

6 基于蒙特卡洛RMT去噪法的小組合風(fēng)險(xiǎn)優(yōu)化實(shí)證研究

6.1 研究方案設(shè)計(jì)

為驗(yàn)證蒙特卡洛RMT去噪法對(duì)原有RMT去噪法的改進(jìn)作用僅適用于小組合的特點(diǎn)，本實(shí)證分析考察大小組合在內(nèi)的不同股票數(shù)量下對(duì)組合風(fēng)險(xiǎn)預(yù)測(cè)的準(zhǔn)確率。這樣，當(dāng)股票數(shù)量不斷減小時(shí)，就可以對(duì)原有RMT去噪法和蒙特卡洛RMT去噪法的組合風(fēng)險(xiǎn)優(yōu)化表現(xiàn)進(jìn)行對(duì)比分析和趨勢(shì)分析。具體地，組合內(nèi)的股票數(shù)量N在[5,400]變化。協(xié)方差矩陣的計(jì)算采用日收益數(shù)據(jù)，所利用的原始股價(jià)數(shù)據(jù)來(lái)自wind數(shù)據(jù)庫(kù)。將股價(jià)轉(zhuǎn)換為收益率的方法為ri,t=lnpi,t-lnpi,t-1，pi,t為第i只股票在第t天的收盤價(jià)，ri,t為第i只股票在第t天的收益率，i=1, 2, …,N，t=1, 2, …,T。

圖1 蒙特卡洛模擬流程Figure 1 Monte Carlo Simulation Process

圖2 正態(tài)分布的累積分布函數(shù)Figure 2 Cumulative Distribution Function of Normal Distribution

步驟1 從2011年1月4日至2014年12月31日的所有日期中隨機(jī)選擇1個(gè)測(cè)試期。測(cè)試期后的1天為投資期，測(cè)試期前(包括測(cè)試期)的700天為預(yù)測(cè)期。

步驟2 從備選股票中隨機(jī)選擇N只股票。備選股票是在上海證券交易所上市交易且在該預(yù)測(cè)期和投資期內(nèi)存在交易數(shù)據(jù)的全部股票。

步驟3 用預(yù)測(cè)期股票協(xié)方差矩陣對(duì)投資期協(xié)方差矩陣進(jìn)行預(yù)測(cè)。當(dāng)采用相等加權(quán)法時(shí)，預(yù)測(cè)期中第i只股票與第j只股票收益率的協(xié)方差為

(16)

其中，ri,t和rj,t分別為第i只股票和第j只股票第t天的收益率，t=1, 2, …, 700；〈·〉為700天的均值。當(dāng)采用指數(shù)加權(quán)法時(shí)，預(yù)測(cè)期中第i只股票與第j只股票收益率的協(xié)方差為

(17)

步驟4 對(duì)于S的不同取值，計(jì)算對(duì)最小方差組合風(fēng)險(xiǎn)的預(yù)測(cè)準(zhǔn)確率。根據(jù)S計(jì)算的最小方差組合中第i只股票的投資權(quán)重wi為[35]

(18)

(19)

(20)

(21)

容易看出PR的值越接近于1，表示最小方差組合風(fēng)險(xiǎn)優(yōu)化水平越高。

步驟5 在N不變的條件下重復(fù)20次步驟1～步驟4，然后分別計(jì)算10種協(xié)方差矩陣預(yù)測(cè)方法在20次試驗(yàn)中的平均PR。

步驟6 重復(fù)步驟1～步驟5，考察N不同時(shí)10種協(xié)方差矩陣預(yù)測(cè)方法的平均PR。

通過上述bootstrap方法，在T和γ不變的情況下得到相等加權(quán)協(xié)方差矩陣和指數(shù)加權(quán)協(xié)方差矩陣去噪前后的組合風(fēng)險(xiǎn)預(yù)測(cè)準(zhǔn)確率隨N的變化曲線，見圖3和圖4。

另外，為了給組合風(fēng)險(xiǎn)預(yù)測(cè)準(zhǔn)確率分析提供更多論據(jù)支持，繪制N的不同取值下MKR去噪法基于蒙特卡洛模擬技術(shù)計(jì)算的λ+的模擬解和原有RMT去噪法計(jì)算的λ+的解析解的圖形，見圖5和圖6。

6.2 實(shí)證結(jié)果和分析

圖3 相等加權(quán)協(xié)方差矩陣去噪前后的PRFigure 3 PR Pre- and Post- Filtering Equally Weighted Covariance Matrices

圖4 指數(shù)加權(quán)協(xié)方差矩陣去噪前后的PRFigure 4 PR Pre- and Post- Filtering Exponentially Weighted Covariance Matrices

圖5 相等加權(quán)協(xié)方差矩陣的最大噪聲特征值Figure 5 The Maximum Noisy Eigenvalues of Equally Weighted Covariance Matrices

圖6 指數(shù)加權(quán)協(xié)方差矩陣的最大噪聲特征值Figure 6 The Maximum Noisy Eigenvalues of Exponentially Weighted Covariance Matrices

對(duì)各種去噪法的組合風(fēng)險(xiǎn)優(yōu)化表現(xiàn)進(jìn)行分析。如前所述，隨著股票數(shù)量的減小，相等加權(quán)和指數(shù)加權(quán)協(xié)方差矩陣的噪聲都越來(lái)越小，這必然導(dǎo)致PR受到的噪聲影響越來(lái)越小，于是各種去噪法引起的PR的增加都呈現(xiàn)減小的趨勢(shì)。由圖3和圖4可知，各種去噪法都能通過去噪而提高組合風(fēng)險(xiǎn)預(yù)測(cè)準(zhǔn)確率，但它們對(duì)其改進(jìn)的程度并不相同。

由于MKR去噪法是基于模擬技術(shù)計(jì)算λ+的，無(wú)法實(shí)現(xiàn)100%的精確度，所以MKR去噪法計(jì)算的λ+并非一種精確無(wú)誤的解，且λ+的具體計(jì)算誤差并不因N的改變而變化。根據(jù)4.2節(jié)的分析，LCPB法、PG+法和KR法等原有RMT去噪法對(duì)λ+的計(jì)算誤差卻隨N的增大而減小。這導(dǎo)致當(dāng)N足夠大時(shí)，原有RMT去噪法對(duì)λ+計(jì)算的精確度高于MKR去噪法，且兩者的差距呈現(xiàn)隨N增大而擴(kuò)大的趨勢(shì)。如圖5和圖6所示，當(dāng)N取最大值400時(shí)，原有RMT去噪法和MKR去噪法對(duì)λ+計(jì)算結(jié)果之差達(dá)到最大。在N逐漸減小的過程中，原有RMT去噪法和MKR去噪法對(duì)λ+計(jì)算差距出現(xiàn)減小的趨勢(shì)，說明隨著股票數(shù)量的減少，原有RMT去噪法對(duì)λ+計(jì)算的精確度不斷下降，使原有RMT去噪法和MKR去噪法的去噪表現(xiàn)越來(lái)越接近。直至當(dāng)股票數(shù)量減小到一定水平時(shí)，它們的去噪表現(xiàn)變得相同。對(duì)于相等加權(quán)協(xié)方差矩陣和指數(shù)加權(quán)協(xié)方差矩陣而言，原有RMT去噪法和MKR去噪法的去噪表現(xiàn)達(dá)到一致時(shí)的股票數(shù)量分別處于160～200和140～160的區(qū)間內(nèi)。這與圖5和圖6的觀察結(jié)果一致，由圖5和圖6可知，在相等加權(quán)法和指數(shù)加權(quán)法下，當(dāng)N分別近似等于170和140時(shí)，原有RMT去噪法和MKR去噪法對(duì)λ+計(jì)算結(jié)果大致相等。在原有RMT去噪法和MKR去噪法的組合風(fēng)險(xiǎn)優(yōu)化作用表現(xiàn)一致以后，當(dāng)股票數(shù)量進(jìn)一步減小，原有RMT去噪法引起的PR的增加開始變得低于MKR去噪法。此時(shí)原有RMT去噪法對(duì)噪聲特征值的計(jì)算精度開始小于MKR去噪法，并且這種精度上的差距是隨股票數(shù)量減小而不斷增大的。λ+計(jì)算精度的上述變化特點(diǎn)可以由圖5和圖6直觀地反映出來(lái)，由圖5和圖6可知，在原有RMT去噪法和MKR去噪法計(jì)算的λ+達(dá)到相同水平后，二者間出現(xiàn)了新的差距，即MKR去噪法計(jì)算的λ+的模擬解開始變得大于原有RMT去噪法計(jì)算的解析解，并且兩者之差隨N的減小而增大。這樣，在各種去噪法引起的PR的增加都在隨N不斷減小而降低的過程中，由于原有RMT去噪法對(duì)λ+計(jì)算精度越來(lái)越小于MKR法，所以它們引起的PR的增加比MKR法先降到0，即不再具有去噪作用。由圖3和圖4可知，當(dāng)采用相等加權(quán)法計(jì)算協(xié)方差矩陣時(shí)，PG+法、LCPB法和KR法的去噪能力降到0時(shí)的股票數(shù)量分別是91、98和110。對(duì)于指數(shù)加權(quán)協(xié)方差矩陣而言，PG+法、LCPB法和KR法的去噪效果達(dá)到0時(shí)的股票數(shù)量分別為93、107和100。之后，當(dāng)股票數(shù)量進(jìn)一步減小時(shí)，原有RMT去噪法無(wú)法起到去噪作用，反而造成PR的下降。而此時(shí)MKR法的去噪作用仍然存在。不過，類似于原有RMT去噪法，當(dāng)股票數(shù)量減小至某一更低水平時(shí)，MKR法也會(huì)由于無(wú)法滿足低噪聲時(shí)的高精度要求而失去去噪能力。在相等加權(quán)法和指數(shù)加權(quán)法條件下，MKR法的去噪能力達(dá)到0時(shí)的股票數(shù)量分別為10和16。根據(jù)上面的分析，由于MKR法和原有RMT去噪法的差異主要體現(xiàn)在λ+的計(jì)算上，因此，股票數(shù)量較小時(shí)MKR法更好的去噪表現(xiàn)主要是由于λ+更高的計(jì)算精確度。

對(duì)相等加權(quán)協(xié)方差矩陣和指數(shù)加權(quán)協(xié)方差矩陣的PR進(jìn)行比較分析。由圖3和圖4可知，在不同股票數(shù)量下，相等加權(quán)協(xié)方差矩陣去噪前的組合風(fēng)險(xiǎn)預(yù)測(cè)準(zhǔn)確率都高于指數(shù)加權(quán)協(xié)方差矩陣去噪前的組合風(fēng)險(xiǎn)預(yù)測(cè)準(zhǔn)確率。然而，相等加權(quán)協(xié)方差矩陣和指數(shù)加權(quán)協(xié)方差矩陣去噪后的結(jié)果卻和去噪前的結(jié)果不同。當(dāng)股票數(shù)量較大時(shí)，每一種去噪法對(duì)指數(shù)加權(quán)協(xié)方差矩陣去噪后的PR都高于對(duì)相等加權(quán)協(xié)方差矩陣去噪后的PR，而當(dāng)股票數(shù)量降至接近100時(shí)，前者開始變得低于后者。如果將所有股票數(shù)量綜合起來(lái)考慮，每一種去噪法對(duì)指數(shù)加權(quán)協(xié)方差矩陣去噪后的PR都大于對(duì)相等加權(quán)協(xié)方差矩陣去噪后的PR。根本來(lái)講，這是因?yàn)槊糠N去噪法引起的指數(shù)加權(quán)協(xié)方差矩陣的PR的增加都更大的緣故。

以上面分析為基礎(chǔ)，對(duì)相等加權(quán)法和指數(shù)加權(quán)法下PR的最佳結(jié)果進(jìn)行分析。由圖3和圖4可知，如果將兩種協(xié)方差矩陣綜合起來(lái)考慮，會(huì)發(fā)現(xiàn)PR的最高值因股票數(shù)量不同而變化。當(dāng)股票數(shù)量較大時(shí)，LCPB去噪法對(duì)指數(shù)加權(quán)協(xié)方差矩陣去噪后的PR最高。當(dāng)股票數(shù)量較小時(shí)，MKR去噪法對(duì)相等加權(quán)協(xié)方差矩陣去噪后的PR最高。而當(dāng)MKR去噪法對(duì)相等加權(quán)協(xié)方差矩陣沒有去噪能力時(shí)，相等加權(quán)協(xié)方差矩陣未去噪時(shí)的PR結(jié)果是最好的。

7 結(jié)論

本研究指出了小組合風(fēng)險(xiǎn)優(yōu)化條件下對(duì)協(xié)方差矩陣去噪的原有RMT方法因噪聲特征值邊界界定誤差而產(chǎn)生的去噪效力下降的問題，并為解決該問題，采用蒙特卡洛模擬法確定樣本協(xié)方差矩陣的最大噪聲特征值，從而建立了蒙特卡洛RMT去噪法。通過實(shí)證分析，在股票收益序列長(zhǎng)度和衰減因子不變時(shí)，在不同股票數(shù)量下對(duì)LCPB法、PG+法和KR法等原有RMT方法和蒙特卡洛RMT法的組合風(fēng)險(xiǎn)優(yōu)化效果進(jìn)行對(duì)比研究，研究結(jié)果如下。

(1)實(shí)際應(yīng)用中股票數(shù)量N是有限的，因此，原有RMT去噪法在N→∞時(shí)計(jì)算的λ+會(huì)偏離于真正的λ+。由于該誤差隨N的減小而增大，所以當(dāng)對(duì)小組合去噪時(shí)，對(duì)λ+的計(jì)算誤差會(huì)上升至較高水平，這必然造成小組合條件下原有RMT方法去噪效力的下降，并進(jìn)一步對(duì)股票投資組合風(fēng)險(xiǎn)的優(yōu)化帶來(lái)不利的影響。

(2)各種RMT方法對(duì)相等加權(quán)協(xié)方差矩陣去噪前后的組合風(fēng)險(xiǎn)結(jié)果與指數(shù)加權(quán)協(xié)方差矩陣的組合風(fēng)險(xiǎn)結(jié)果非常類似。協(xié)方差矩陣未去噪時(shí)，隨股票數(shù)量減小，組合風(fēng)險(xiǎn)預(yù)測(cè)的準(zhǔn)確率越來(lái)越高，主要是因?yàn)楣善睌?shù)量的減小弱化了“維數(shù)災(zāi)禍”的作用，從而使協(xié)方差矩陣噪聲減小。

(3)當(dāng)股票數(shù)量足夠大時(shí)，由于原有RMT去噪法對(duì)噪聲特征值邊界界定的精確度高于蒙特卡洛RMT去噪法，因此，原有RMT去噪法引起的組合風(fēng)險(xiǎn)優(yōu)化水平的增加更大。隨著股票數(shù)量的減小，相對(duì)于蒙特卡洛RMT去噪法，原有RMT去噪法的去噪優(yōu)勢(shì)因?qū)υ肼曁卣髦颠吔缃缍ㄕ`差的增大而呈減小趨勢(shì)。直至當(dāng)原有RMT去噪法的誤差開始超過蒙特卡洛RMT去噪法的誤差時(shí)，原有RMT去噪法的組合風(fēng)險(xiǎn)優(yōu)化作用也開始變得低于蒙特卡洛RMT法。這說明蒙特卡洛RMT去噪法是一種解決小股票組合條件下原有RMT方法去噪效力下降的有效方法。

股票收益樣本協(xié)方差矩陣含有的噪聲會(huì)造成Markowitz股票投資組合風(fēng)險(xiǎn)優(yōu)化水平的下降，而原有RMT去噪法在小組合條件下的效力卻很小，甚至為0或負(fù)值。針對(duì)原有RMT去噪法的上述缺陷，本研究設(shè)計(jì)了蒙特卡洛RMT去噪法。本研究能夠?yàn)樽C券市場(chǎng)投資中股票投資組合的權(quán)重分配提供具體的方法，對(duì)證券投資決策及有關(guān)政策制定均有一定的借鑒意義。由于本研究的實(shí)證研究?jī)H以中國(guó)股票數(shù)據(jù)為樣本，沒有考慮其他國(guó)家，如一些發(fā)達(dá)經(jīng)濟(jì)體的證券市場(chǎng)情況，所以可能會(huì)影響實(shí)證結(jié)果的普適性，這是后續(xù)研究中需要進(jìn)一步解決的問題。

[1]李洋,余麗霞.基于馬科維茨理論的最優(yōu)證券組合分析.財(cái)會(huì)月刊,2013(22):53-55.

LI Yang,YU Lixia.Analysis of optimal stock portfolio based on Markowitz theory.FinanceandAccountingMonthly,2013(22):53-55.(in Chinese)

[2]ZHANG M.Finite-sample linear filter optimization in wireless communications and financial systems.IEEETransactionsonSignalProcessing,2013,61(20):5014-5025.

[3]黃飛雪.基于單鏈接聚類過濾法的均值方差模型.預(yù)測(cè),2011,30(1):66-70.

HUANG Feixue.Mean-variance based on filtered single linkage cluster.Forecasting,2011,30(1):66-70.(in Chinese)

[4]曾杏元.有關(guān)隨機(jī)矩陣領(lǐng)域最新研究動(dòng)態(tài)與進(jìn)展的綜述報(bào)告.數(shù)學(xué)理論與應(yīng)用,2011,31(3):7-19.

ZENG Xingyuan.A survey on the latest dynamic and progress of random matrix field.MathematicalTheoryandApplications,2011,31(3):7-19.(in Chinese)

[5]CHU Chen,FANG Zhaoben.Optimal portfolio project with modified covariance matrix and its stability.JournalofUniversityofScienceandTechnologyofChina,2011,41(12):1035-1041.

[6]毛磊,王宗軍,王玲玲.機(jī)構(gòu)投資者持股偏好、篩選策略與企業(yè)社會(huì)績(jī)效.管理科學(xué),2012,25(3):21-33.

MAO Lei,WANG Zongjun,WANG Lingling.Institutional investors preferences,screening and corporate social performance.JournalofManagementScience,2012,25(3):21-33.(in Chinese)

[7]王明進(jìn).多元波動(dòng)率模型的一些新進(jìn)展.數(shù)理統(tǒng)計(jì)與管理,2010,29(2):232-247.

WANG Mingjin.Recent developments in multivariate volatility modeling.JournalofAppliedStatisticsandManagement,2010,29(2):232-247.(in Chinese)

[8]NGUYEN Q.One-factor model for the cross-correlation matrix in the Vietnamese stock market.PhysicaA:StatisticalMechanicsandItsApplications,2013,392(13):2915-2923.

[9]LIM K,KIM M J,KIM S,et al.Statistical properties of the stock and credit market:RMT and network topology.PhysicaA:StatisticalMechanicsandItsApplications,2014,407:66-75.

[10] LALOUX L,CIZEAU P,POTTERS M,et al.Random matrix theory and financial correlations.InternationalJournalofTheoreticalandAppliedFinance,2000,3(3):391-397.

[11] LEE Y H.Noisedetectionfromfinancialcorrelationmatrices.(2001-05-19)[2015-05-14].http:∥web.mit.edu/bazant/www/teach/18.325/projects.lee.ps.

[12] PAFKA S,KONDOR I.Estimated correlation matrices and portfolio optimization.PhysicaA:StatisticalMechanicsandItsApplications,2004,343:623-634.

[13] TOLA V,LILLO F.Cluster analysis for portfolio optimization.JournalofEconomicDynamics&Control,2008,32(1):235-258.

[14] PLEROU V,GOPIKRISHNAN P,ROSENOW B,et al.Random matrix approach to cross correlations in financial data.PhysicalReviewE,2002,65(6):066126-1-066126-18.

[15] SHARIFI S,CRANE M,SHAMAIE A,et al.Random matrix theory for portfolio optimization:a stability approach.PhysicaA:StatisticalMechanicsandItsApplications,2004,335(3/4):629-643.

[16] DALY J,CRANE M,RUSKIN H J.Random matrix theory filters in portfolio optimisation:a stability and risk assessment.PhysicaA:StatisticalMechanicsandItsApplications,2008,387(16/17):4248-4260.

[17] KAROUI N E.On the realized risk of high-dimensional Markowitz portfolios.SIAMJournalonFinancialMathematics,2013,4(1):737-783.

[18] RUBIO F,MESTRE X,PALOMAR D P.Performance analysis and optimal selection of large minimum variance portfolios under estimation risk.IEEEJournalofSelectedTopicsinSignalProcessing,2012,6(4):337-350.

[19] WAKAI R,SHINZATO T,SHIMAZAKI Y.Random matrix approach for portfolio optimization problem.JournalofJapanIndustrialManagementAssociation,2014,65(1):17-28.

[20] TANG L C,CHEN Y N.Noise dressing of ECM and investment portfolio selection based on RMT in China financial market∥ProceedingsofInternationalConferenceonManagement&ServiceScience.Piscataway:IEEE,2010:1-4.

[21] OH G,EOM C,WANG F,et al.Statistical properties of cross-correlation in the Korean stock market.TheEuropeanPhysicalJournalB,2011,79(1):55-60.

[22] 韓華,吳翎燕,宋寧寧.基于隨機(jī)矩陣的金融網(wǎng)絡(luò)模型.物理學(xué)報(bào),2014,63(13):138901-1-138901-10.

HAN Hua,WU Lingyan,SONG Ningning.Financial networks model based on random matrix.ActaPhysicalSinica,2014,63(13): 138901-1-138901-10.(in Chinese)

[23] ARAI Y,YOSHIKAWA T,IYETOMI H.Complex principal component analysis of dynamic correlations in financial markets∥Proceedingsofthe5thKESInternationalConferenceonIntelligentDecisionTechnologies.IOS Press,2013:111-119.

[24] YANG J,WANG Z,LIU X,et al.Group identification in stock correlation networks with an application to portfolio selection.JournalofComputationalInformationSystems,2014,10(7):2937-2949.

[25] WANG D,PODOBNIK B,HORVATIC D,et al.Quantifying and modeling long-range cross correlations in multiple time series with applications to world stock indices.PhysicalReviewE,2011,83(4):046121-1-046121-9.

[26] WANG G J,XIE C,CHEN S,et al.Random matrix theory analysis of cross-correlations in the US stock market:evidence from Pearson′s correlation coefficient and detrended cross-correlation coefficient.PhysicaA:StatisticalMechanicsandItsApplications,2013,392(17): 3715-3730.

[27] CONLON T,RUSKIN H J,CRANE M.Random matrix theory and fund of funds portfolio optimisation.PhysicaA:StatisticalMechanicsandItsApplications,2007,382(2):565-576.

[28] KUMAR S,DEO N.Correlation and network analysis of global financial indices.PhysicalReviewE,2012,86(2):026101-1-026101-8.

[29] PAFKA S, POTTERS M,KONDOR I. Exponential weighting and random-matrix-theory-based denoising of financial covariance matrices for portfolio optimization.(2004-02-24)[2015-05-14].http:∥eprintweb.org/S/article/cond-mat/0402573/cited.

[30] DALY J,CRANE M,RUSKIN H J.Random matrix theory filters and currency portfolio opimisation.JournalofPhysics:ConferenceSeries,2010,221:1-14.

[31] SANDOVAL L,Jr,BORTOLUZZO A B,VENEZUELA M K.Not all that glitters is RMT in the forecasting of risk of portfolios in the Brazilian stock market.PhysicaA:StatisticalMechanicsandItsApplications,2014,410:94-109.

[32] 羅英,蔡玉梅,崔小梅,等.資產(chǎn)組合協(xié)方差矩陣的信息結(jié)構(gòu).預(yù)測(cè),2013,32(4):26-30.

LUO Ying,CAI Yumei,CUI Xiaomei,et al.The information structure of the covariances between financial returns.Forecasting,2013,32(4):26-30.(in Chinese)

[33] 唐曉清,白延琴,劉念祖,等.基于隨機(jī)矩陣?yán)碚摰腗arkowitz組合投資模型.上海大學(xué)學(xué)報(bào):自然科學(xué)版,2013,19(3):293-297.

TANG Xiaoqing,BAI Yanqin,LIU Nianzu,et al.Markowitz portfolio model based on random matrix theory.JournalofShanghaiUniversity:NaturalScience,2013,19(3):293-297.(in Chinese)

[34] 惠曉峰,李冰娜.基于隨機(jī)矩陣?yán)碚摏Q定多元GARCH模型最佳維度研究.運(yùn)籌與管理,2011,20(4):141-148.

HUI Xiaofeng,LI Bingna.Study on determining the optimal dimension of multi-GARCH Models based on random matrix theory.OperationsResearchandManagementScience,2011,20(4):141-148.(in Chinese)

[35] CLARKE R,DESILVA H,THORLEY S.Minimum-variance portfolio composition.TheJournalofPortfolioManagement,2011,37(2):31-45.

[36] 瞿慧,劉燁.滬深300指數(shù)收益率及已實(shí)現(xiàn)波動(dòng)聯(lián)合建模研究.管理科學(xué),2012,25(6):101-110.

QU Hui,LIU Ye.A joint model for CSI300 index return and realized volatility.JournalofManagementScience,2012,25(6):101-110.(in Chinese)

ResearchonRiskOptimizationofSmallStockPortfolioBasedontheFilteringMethodofRMTUsingMonteCarloSimulation

LI Bingna1，HUI Xiaofeng2，LI Lianjiang3

1 School of Economics, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China2 School of Management, Harbin Institute of Technology, Harbin 150001, China3 School of Control Engineering, Northeastern University at Qinhuangdao, Qinhuangdao 066004, China

In the framework of classical Markowitz portfolio optimization problem, covariance matrices of stock returns obtained from historical data contain a high amount of noise, caused “curse of dimensions”, which will forbid the perfect construction of Markowitz portfolio and the optimization of its risk seriously. According to numerous studies, the filtering approaches based on random matrix theory are effective in reducing the estimation noise of financial covariance matrices and then achieving portfolio risk optimization.

Starting from the meaning of Markowitz stock portfolio risk, this work focuses on the accuracy of portfolio risk prediction as the measure for assessing optimal performance of portfolio risk and interprets the principle of the existing RMT-based filtering methods for stock return covariances. On this basis, we point out the calculation error of the maximum noisy eigenvalue in the existing RMT methods, which will turn out to reduce the ectiveness of these methods in the context of risk optimization of small portfolios. For any practical use of Markowitz portfolio theory, it would therefore be necessary to reduce this calculation error of the existing methods of RMT in the case of a small number of stocks. We introduce a new filtering method based on random matrix theory, which determines the theoretical edge of the random part of the eigenvalue distribution by a Monte Carlo simulation. Finally, we illustrate, on empirical data, the effect of the newly introduced method and the existing random-matrix-theory-based methods including LCPB, PG+ and KR on portfolio risk optimization.

Our empirical research has returned the following results. Firstly, whether covariance matrices of stock returns are filtered or unfiltered, it can be seen that, in general, the experimental results for equally weighted covariance matrices are very similar to those for exponentially weighted covariance matrices. Secondly, for unfiltered forecasts, the performance of portfolio risk is better when a smaller number of stocks is used for the estimation. The reason for this could be that the fewer stocks are used, the smaller amount of estimation noise can be in return covariance matrices. Thirdly, if enough stocks are taken into account, the existing RMT-based filters are superior to the newly proposed filter in terms of the risk optimization performance. However, the fewer the stocks are, the more inaccurately the maximum noisy eigenvalue in the existing filtering methods is calculated. As a consequence, the ectiveness of the existing filtering methods in the optimization of portfolio risk deteriorates as the portfolio size decreases. When the portfolio size is reduced to less than 200 stocks, the newly proposed filter begins to outperform the existing RMT filters. The results show that our novel filter is of great help in improving the ectiveness of the existing RMT-based filters in the context of risk optimization of small portfolios.

Combining the filtering procedure based on random matrix theory with the Monte Carlo simulation technique, the newly introduced filtering method can be very powerful in constructing portfolios with better weights allocation and improving the accuracy of portfolio risk prediction. In addition, our method also helps substantially in reasonable design of stock investment strategies and the optimization of stock investment risk.

Markowitz investment portfolio；portfolio risk optimization；small portfolio；RMT based filtering method；Monte Carlo simulation

Date:June 24th, 2015

DateMarch 4th, 2016

FundedProject:Supported by the National Natural Science Foundation of China(71401028), the Fundamental Research Funds for the Central Universities(N130323008) and the Funded by Northeastern University at Qinhuangdao(XNB201418)

Biography:LI Bingna, doctor in management, is a lecturer in the School of Economics at Northeastern University at Qinhuangdao. Her research interests include financial investment portfolio and risk management. Her representative paper titled “Optimization of stock portfolio risks based on correlation filtering using minimum perturbations of eigenvectors” was published inTheJournalofSystemsEngineering-theory&Practice(Issue 10, 2013). E-mail：libingna519@163.com

HUI Xiaofeng, doctor in management, is a professor in the School of Management at Harbin Institute of Technology. His research interests include exchange rate forecast, financial security and financial market efficiency. He is the principal researcher for the project titled “Disclosure of listed firms and capital costs” supported by the National Natural Science Foundation of China(70573030). E-mail：xfhui@hit.edu.cn

LI Lianjiang, doctor in engineering, is a lecturer in the School of Control Engineering at Northeastern University at Qinhuangdao. His research interests include multidisciplinary combined simulation. His representative paper titled “Research on improvement of ladar imaging quality based on fusion technology” was published in theLasersinEngineering(Issue 1-2, 2012). E-mail：lilianjiangpro@163.com

F830.91

A

10.3969/j.issn.1672-0334.2016.02.011

1672-0334(2016)02-0134-12

2015-06-24修返日期2016-03-04

國(guó)家自然科學(xué)基金(71401028)；中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金(N130323008)；東北大學(xué)秦皇島分校校內(nèi)基金(XNB201418)

李冰娜，管理學(xué)博士，東北大學(xué)秦皇島分校經(jīng)濟(jì)學(xué)院講師，研究方向?yàn)榻鹑谕顿Y組合和風(fēng)險(xiǎn)管理等，代表性學(xué)術(shù)成果為“基于最小擾動(dòng)相關(guān)去噪法股票投資組合風(fēng)險(xiǎn)優(yōu)化”，發(fā)表在2013年第10期《系統(tǒng)工程理論與實(shí)踐》，E-mail：libingna519@163.com

惠曉峰，管理學(xué)博士，哈爾濱工業(yè)大學(xué)管理學(xué)院教授，研究方向?yàn)閰R率預(yù)測(cè)、金融安全和金融市場(chǎng)效率等，曾主持國(guó)家自然科學(xué)基金項(xiàng)目“上市公司信息披露與資本成本研究”(70573030)，E-mail：xfhui@hit.edu.cn

李連江，工學(xué)博士，東北大學(xué)秦皇島分校控制工程學(xué)院講師，研究方向?yàn)槎鄬W(xué)科聯(lián)合仿真等，代表性學(xué)術(shù)成果為“Research on improvement of ladar imaging quality based on fusion technology”，發(fā)表在2012年第1-2期《Lasers in Engineering》，E-mail：lilianjiangpro@163.com