林擎鑫 馬 寧 王一成 成 晨

(西北民族大學(xué)中國(guó)民族信息技術(shù)研究院 甘肅 蘭州 730000)

0 引 言

隨著現(xiàn)代投資組合理論的不斷完善與發(fā)展，如何保證在風(fēng)險(xiǎn)最低化的同時(shí)實(shí)現(xiàn)收益最大化成為了研究的熱點(diǎn)。針對(duì)直接利用Markowitz提出的基于均值-方差決策模型解決上述問題存在的模型中一系列嚴(yán)格假設(shè)導(dǎo)致了無效資產(chǎn)組合以及均值-方差相同的組合，受偏度影響，其風(fēng)險(xiǎn)程度可能存在較大差異等不足之處[1-3]，研究者們對(duì)如何優(yōu)化投資模型開展了深入研究。高岳林等[4]在均值-方差模型的基礎(chǔ)上增加了風(fēng)險(xiǎn)以及條件風(fēng)險(xiǎn)價(jià)值作為約束條件，并建立了相應(yīng)模型，結(jié)合我國(guó)股票市場(chǎng)實(shí)情開展了實(shí)證分析，證明了方法的合理性及優(yōu)勢(shì)性；黃東賓等[5]在M-CVaR模型的基礎(chǔ)上引入偏好強(qiáng)度，建立了基于偏好強(qiáng)度的M-CVaR模型，通過分析驗(yàn)證了所提方法對(duì)于多位數(shù)據(jù)的整合以及資產(chǎn)選擇能力上所具備的優(yōu)勢(shì)。上述方法其主要內(nèi)涵是通過增加約束條件實(shí)現(xiàn)對(duì)于傳統(tǒng)方法的改進(jìn)，這些改進(jìn)雖然起到了積極意義，卻忽略了交易成本帶來的影響，使得模型在一定程度上脫離了我國(guó)證券市場(chǎng)的實(shí)際情況。同時(shí)，交易成本作為市場(chǎng)中重要摩擦因素，對(duì)于投資結(jié)果產(chǎn)生深遠(yuǎn)的影響，因此必須予以考慮。為了達(dá)到考慮成本因素且具備較小計(jì)算代價(jià)的目的，就需要對(duì)組合投資模型中的算法進(jìn)行優(yōu)化。王懷柱[6]結(jié)合改進(jìn)遺傳算法進(jìn)行CVaR模型求解，通過改進(jìn)遺傳算法的變異算子提高了模型收斂速度，從而達(dá)到了節(jié)約計(jì)算時(shí)間的目的；見靜等[7]提出了結(jié)合動(dòng)態(tài)粒子群算法對(duì)于M-VaR進(jìn)行求解的思路，并結(jié)合實(shí)證分析表明該算法相對(duì)于傳統(tǒng)算法在效率以及精度上存在的優(yōu)勢(shì)。一些智能算法替代傳統(tǒng)算法，但容易陷入局部最優(yōu)解，且求解效率較差[8-9]。

綜上所述，為了從根本上解決上述問題，本文提出基于改進(jìn)蟻獅算法優(yōu)化考慮成本因素的組合投資模型。該模型在傳統(tǒng)模型的基礎(chǔ)上引入二次分段凹交易成本，建立了更為具體的組合投資模型；同時(shí)，針對(duì)蟻獅算法提出相應(yīng)的優(yōu)化方案，以提高其探索能力、收斂速度，并避免陷入局部最優(yōu)解。利用上述模型結(jié)合多種證券投資組合進(jìn)行分析，證明了模型在穩(wěn)定求解的同時(shí)，可以保證求解過程的高效、高精度特性，從而為投資者提供更好的決策參考。

1 模型建立

1.1 CVaR模型的建立與計(jì)算

CVaR即條件風(fēng)險(xiǎn)價(jià)值，表示了在一定的置信水平下，投資者可對(duì)于未來某一時(shí)間所承擔(dān)的且超過VaR均值的損失進(jìn)行預(yù)測(cè)[10]。其數(shù)學(xué)表達(dá)式為：

VaRβ+E[L(w,r)-VaRβ|L(w,r)>VaRβ]

(1)

式中：w=(w1,w2，…，wi，…，wn)T，wi表示第i種資產(chǎn)的投資比例；r=(r1,r2，…，ri，…，rn)T，ri表示對(duì)應(yīng)的i種資產(chǎn)的收益率；L(x,r)表示損失函數(shù)；β表示置信度。

對(duì)于CVaR而言其計(jì)算的核心在于解決損失函數(shù)帶來的問題，故而依據(jù)損失函數(shù)是否連續(xù)，可以把CVaR的計(jì)算分為連續(xù)性和離散型，具體計(jì)算如下。

(1) 連續(xù)型計(jì)算。記p(r)為概率密度函數(shù)，其表示投資組合的收益率，則：

(2)

(3)

上述方法表征了對(duì)于CVaR的具體計(jì)算。但是，在實(shí)際應(yīng)用過程中還需對(duì)其中的重要參數(shù)進(jìn)行選擇，以求達(dá)到對(duì)于不同風(fēng)險(xiǎn)偏好的投資者在承受低風(fēng)險(xiǎn)的同時(shí)，其收益最大化的目的。在這一過程中，置信水平以及對(duì)應(yīng)產(chǎn)品的投資期限的選擇對(duì)于投資組合問題的解決具備重要意義。前者不僅體現(xiàn)了投資者的風(fēng)險(xiǎn)偏好程度，同時(shí)還表征了計(jì)算結(jié)果的有效性，其選取決定于CVaR的計(jì)算目的；后者表征了計(jì)算期限，主要決定于資產(chǎn)流通性、回報(bào)率的分布特性和約束條件等因素；兩者共同決定了VaR的取值。

1.2 考慮交易成本的M-CVaR模型

交易成本是指實(shí)際證券交易過程中的各種摩擦因素，其中要分為固定成本以及間接成本[11]。對(duì)于中國(guó)證券交易市場(chǎng)而言，固定成本是交易成本最直接的體現(xiàn)。其具體表達(dá)式為：

(4)

通過式(4)可以發(fā)現(xiàn)交易成本是一個(gè)二次分段凹函數(shù)且與資產(chǎn)投資比例呈正相關(guān)。

1.2.1模型假設(shè)

為了使得建立的模型更加符合我國(guó)證券交易市場(chǎng)的實(shí)際情況，本文進(jìn)行如下假設(shè)的條件的建立。

(1) 投資者均為理性投資人假設(shè)。

(2) 價(jià)格對(duì)于價(jià)值準(zhǔn)確反映的假設(shè)，即投資者之間進(jìn)行公平競(jìng)爭(zhēng)，競(jìng)爭(zhēng)過程中不產(chǎn)生間接成本，產(chǎn)品價(jià)格直接體現(xiàn)產(chǎn)品的真實(shí)價(jià)值。

(3) 收益分布假設(shè)。結(jié)合文獻(xiàn)可以發(fā)現(xiàn)，我國(guó)證券交易市場(chǎng)的收益率分布特征符合拉普拉斯分布，即其概率密度與分布性質(zhì)分布表示為：

(5)

(6)

(5) 最小單位以及最大交易量約束；前者以整手為約束，后者按照相關(guān)規(guī)定單個(gè)基金對(duì)于同一家上市公司股票的持有小于其凈資產(chǎn)的10%。

(6) 不允許出現(xiàn)買空賣空的現(xiàn)象。

1.2.2模型建立

設(shè)一個(gè)投資人共投資了n種風(fēng)險(xiǎn)資產(chǎn)以及一種無風(fēng)險(xiǎn)資產(chǎn)，其中單個(gè)風(fēng)險(xiǎn)資產(chǎn)的收益率記為ri(i=1,2,…,n)；第i類資產(chǎn)的購(gòu)買數(shù)Si(i=1,2,…,n)，無風(fēng)險(xiǎn)資產(chǎn)購(gòu)買數(shù)記為S0;無風(fēng)險(xiǎn)資產(chǎn)收益率記為r0；第i類資產(chǎn)的投資比例為；第i類資產(chǎn)的股票單價(jià)記為Pi；投資組合的總收益率記為r；初始財(cái)富記為R0，末期財(cái)富記為Rm，則有:

(7)

(8)

(9)

同時(shí)，第i種資產(chǎn)的權(quán)重以及投資組合的收益率分別表示為:

(10)

(11)

通過上述計(jì)算過程，組合收益的期望與方差可以表示為：

(12)

式中：σij表示不同資產(chǎn)的協(xié)方差。

綜合上述內(nèi)容并結(jié)合式(2)和式(3)以及拉普拉斯假設(shè)，引入ε作為投資者的風(fēng)險(xiǎn)承受極限，同時(shí)考慮本文研究目的在于求解投資者低風(fēng)險(xiǎn)下收益最大化模型，故而基于交易成本的M-CVaR模型為：

(13)

(14)

0≤Si≤μi,且Si?Z

式中：μi表示第i類資產(chǎn)的最大交易約束。

Crystal將抒情詩(shī)定義為“以一種短小、如歌的形式，抒發(fā)作者內(nèi)心主觀思想和情感的詩(shī)歌”。[5]從抒情詩(shī)歌的定義來看，抒情詩(shī)歌有明顯的文體特點(diǎn)：體裁短小精悍、韻律節(jié)奏感強(qiáng)、表達(dá)詩(shī)人的內(nèi)心情感和主觀感受。除此之外，抒情詩(shī)歌還有一些較明顯的特點(diǎn)，如修辭的使用。

上述模型是一個(gè)非線性模型，求解時(shí)存在精度和效率的雙重困擾，故而本文將引入改進(jìn)蟻獅算法對(duì)其進(jìn)行求解，以求高效、高精度實(shí)現(xiàn)對(duì)于模型的求解。

2 改進(jìn)蟻獅算法

蟻獅算法相對(duì)于蟻群算法、遺傳算法等智能算法由于其考慮了螞蟻游走與蟻獅捕食螞蟻的雙重行為，因此對(duì)于復(fù)雜問題具備更強(qiáng)的適應(yīng)性[12-13]；同時(shí)，由于該算法具備少量調(diào)節(jié)參數(shù)，因此使得其實(shí)現(xiàn)更為容易。但是，蟻獅算法也存在固有缺陷，主要體現(xiàn)為容易陷入局部最優(yōu)解、搜索能力不足、求解時(shí)具備較慢的收斂速度[14-15]。對(duì)上述缺點(diǎn)，本文提出改進(jìn)蟻獅算法，其主要思路是在傳統(tǒng)蟻獅算法的基礎(chǔ)上引入相應(yīng)策略，從而提高算法的整體性能。

2.1 傳統(tǒng)蟻獅算法的原理

傳統(tǒng)的蟻獅算法其基本原理是模擬蟻獅捕食幼蟻的自然行為。該算法按如下過程開展[16-17]：

(1) 利用式(15)進(jìn)行初始解的隨機(jī)生成，同時(shí)，對(duì)所有蟻獅的適應(yīng)度進(jìn)行計(jì)算和排序，其中適應(yīng)度最高的視為最優(yōu)解即為精英蟻獅，記為RE。

Xn.m=La+rand(Lb-La)

(15)

式中：Xn.m表示螞蟻位置且初始時(shí)刻螞蟻位置與蟻獅位置相同；La、Lb分別表示螞蟻游走可行域的上下限。

(2) 根據(jù)式(16)計(jì)算螞蟻游走的位置。

Xi=[0,cumsum(2r(1)-1),cumsum(2r(2)-1),…，

cumsum(2r(n)-1)]

(16)

式中：Xi表示螞蟻位置；cumsum表示螞蟻游走位置的累積和；n表示最大迭代步數(shù)；r(n)表示隨機(jī)函數(shù)，按式(17)計(jì)算。

(17)

為了確保每一只螞蟻的隨機(jī)游走過程都在可行域的范圍內(nèi)進(jìn)行，還需要對(duì)螞蟻的位置進(jìn)行離差，其具體計(jì)算如下：

(18)

(19)

(20)

式中：w為常數(shù)，按式(21)取值；t與T分別表示當(dāng)前迭代次數(shù)與最大迭代次數(shù)。

(21)

(3) 對(duì)圍繞蟻獅RA與精英蟻獅RE游走的螞蟻產(chǎn)生的隨機(jī)位置進(jìn)行更新；

(22)

(4) 如果在隨機(jī)游走的螞蟻中，出現(xiàn)某只螞蟻的適應(yīng)度高于蟻獅，則視作該螞蟻被蟻獅捕獲，此時(shí)，蟻獅會(huì)根據(jù)該螞蟻的位置更新其自身位置，重構(gòu)陷阱，其數(shù)學(xué)表達(dá)式為：

(23)

(5) 檢驗(yàn)迭代次數(shù)是否達(dá)到最大，若是則輸出最優(yōu)解，反之重復(fù)上述步驟直至最優(yōu)解產(chǎn)生。

2.2 基于不同策略下蟻獅算法的優(yōu)化

本文主要通過以下兩種策略進(jìn)行蟻獅算法的優(yōu)化：

(1) 邊界策略。通過對(duì)傳統(tǒng)蟻獅方法的計(jì)算過程進(jìn)行分析，我們可以發(fā)現(xiàn)隨著迭代次數(shù)的增加，螞蟻圍繞蟻獅走的可行域在逐漸減小，而對(duì)于同一輪迭代中呈線性分段變化趨勢(shì)的螞蟻，其游走邊界完全相同。這種方法一定程度上降低了算法的多樣性與隨機(jī)性，不易于全局最優(yōu)解的獲取。本文對(duì)于式(20)進(jìn)行如下改進(jìn)：

(24)

這種方法使得I值隨著迭代次數(shù)的增加由線性遞增轉(zhuǎn)變?yōu)榉蔷€性隨機(jī)自適應(yīng)遞增，由式(24)可知，邊界將按非線性隨機(jī)自適應(yīng)遞減，這就使得相對(duì)于傳統(tǒng)方法，螞蟻圍繞蟻獅游走更具備多樣性與隨機(jī)性，從而提高了算法的全局搜索能力。

(2) 選擇策略。傳統(tǒng)的蟻獅方法其選擇策略為輪盤賭策略，但這種策略難以保證所有螞蟻都圍繞著較優(yōu)蟻獅游走。因此，本文將預(yù)備參加輪賭選擇的蟻獅按適應(yīng)度值從大到小排序，同時(shí)設(shè)定只有當(dāng)適應(yīng)度值F滿足式(25)時(shí)，蟻獅參與輪盤賭。

F

(25)

式中：P為預(yù)備蟻獅適應(yīng)度的平均值。

這種方法，提高了輪賭的門檻，在一定程度上保留蟻獅的同時(shí)，加大了螞蟻圍繞較優(yōu)螞蟻?zhàn)叩母怕?，從而提高了算法的性能。在?yōu)化考慮交易成本的M-CVaR模型時(shí)，適應(yīng)度函數(shù)與目標(biāo)函數(shù)的關(guān)系可以用式(26)表示。

(26)

3 實(shí)證分析

通過對(duì)我國(guó)證券市場(chǎng)的實(shí)情進(jìn)行分析可以發(fā)現(xiàn)多數(shù)投資者會(huì)選擇支付成本更大的積極策略進(jìn)行相關(guān)的投資以求在控制風(fēng)險(xiǎn)的同時(shí)，獲得更大的利益。故而，本節(jié)將建立考慮交易成本的M-CVaR模型并結(jié)合改進(jìn)蟻獅算法求解最優(yōu)資產(chǎn)組合，同時(shí)，在這一過程中驗(yàn)證考慮成本因素的積極意義以及算法的有效性。

3.1 數(shù)據(jù)預(yù)處理

為了更有效地證明求解模型的有效性，以近期熱門的券商、科技、新能源和農(nóng)業(yè)等作為股票選取的題材，具體如下：京東方A(000725)、西部證券(002673)、牧原股份(002714)、比亞迪(002594)、格力電器(000651)、美的集團(tuán)(000333)、中興通訊(000063)、五糧液(000858)和一種無風(fēng)險(xiǎn)資產(chǎn)，本文選擇國(guó)債為例；同時(shí)，所有股票數(shù)據(jù)來源于“萬得咨詢”。

假設(shè)初始投資為1 000萬元，數(shù)據(jù)取2019年8月16日至2019年10月18日共計(jì)10個(gè)交易周來進(jìn)行數(shù)據(jù)監(jiān)測(cè)，得T=10。通過實(shí)際數(shù)據(jù)調(diào)研，不存在停盤以及復(fù)盤現(xiàn)象。

3.2 參數(shù)設(shè)定

按照M-CVaR模型求解的相關(guān)需求，對(duì)以下參數(shù)進(jìn)行詳細(xì)說明：

(1) 風(fēng)險(xiǎn)值。設(shè)置最大承受風(fēng)險(xiǎn)系數(shù)δ=0.05。

(2) 置信水平。其選擇受多因素約束，為了保證計(jì)算的有效性要求置信水平較低；為了更準(zhǔn)確反映投資者的風(fēng)險(xiǎn)偏好，則需要選擇相對(duì)較高的置信水平。為了更有效地對(duì)問題進(jìn)行說明，同時(shí)為了展開相應(yīng)的比較分析，本文取置信水平β為0.85，0.9，0.95。

(3) 交易成本。根據(jù)假設(shè)，交易成本僅考慮固定交易成本，按2.5計(jì)算。

(4) 無風(fēng)險(xiǎn)資產(chǎn)。以短期貼現(xiàn)的國(guó)家債券為例，其2019年利率水平為0.013%，故而r0=0.013%(以天為計(jì)算單位)。

(5) 收益率。前人的分析中僅考慮算術(shù)平均收益率，這種方式忽略了資本的時(shí)間價(jià)值同時(shí)會(huì)出現(xiàn)上傾現(xiàn)象，因此本文提出使用幾何平均收益率克服上述問題。通過計(jì)算，資產(chǎn)的期望收益率如表1所示。

表1 資產(chǎn)期望收益率(%)

3.3 數(shù)據(jù)結(jié)果分析

3.3.1模型計(jì)算

通過計(jì)算，考慮成本因素的計(jì)算結(jié)果如表2所示。

表2 考慮成本因素的最優(yōu)投資組合

可以看出，改進(jìn)蟻群算法對(duì)于求解考慮成本因素的M-CVaR模型有效性。同時(shí)，由于置信水平的選擇體現(xiàn)了投資者的風(fēng)險(xiǎn)偏好程度，投資者可以根據(jù)自身情況選擇合適的置信度組合進(jìn)行投資。如果投資者厭惡風(fēng)險(xiǎn)，則選擇置信度高的組合進(jìn)行投資，如果投資者對(duì)風(fēng)險(xiǎn)厭惡程度低，可選擇置信度低的組合。結(jié)合我國(guó)實(shí)際證券市場(chǎng)，在后面的分析中置信度均取0.95。因此，在置信水平為0.95時(shí)，不考慮成本因素的投資最優(yōu)組合如表3所示。

表3 不考慮成本因素的最優(yōu)投資組合

3.3.2數(shù)據(jù)分析

在置信度為0.95時(shí)，考慮成本因素的實(shí)際收益為451 802元；不考慮成本因素的實(shí)際收益為393 145元；因此可以說明在風(fēng)險(xiǎn)系數(shù)δ=0.05時(shí)，考慮成本因素計(jì)算出來的最優(yōu)投資組合收益有95%的可能要高于理論綜合效益；而不考慮成本因素的實(shí)際收益值則小于理論綜合收益。綜上所述考慮交易成本的模型對(duì)于最優(yōu)組合的求取更具備意義。

3.4 方法對(duì)比

傳統(tǒng)智能算法如遺傳算法、蟻群算法等在求解模型時(shí)也被廣泛應(yīng)用，基于遺傳算法與蟻群算法對(duì)于模型求解的結(jié)果如表4所示。

表4 遺傳算法與蟻群算法的計(jì)算結(jié)果

由不同方法在置信度為0.95時(shí)的計(jì)算結(jié)果可以發(fā)現(xiàn)，改進(jìn)蟻獅子算法的求解結(jié)果優(yōu)于遺傳算法與蟻獅算法，且其求解效率也高于遺傳算法與蟻獅子算法，因此對(duì)于考慮成本因素的M-CVaR模型該算法具備更好的適用性。同時(shí)，該算法不需要對(duì)金融變量的分布以及市場(chǎng)參數(shù)做其余假設(shè)，避免了求解時(shí)的情景限制，因此其適用于對(duì)于多種風(fēng)險(xiǎn)度量下帶現(xiàn)實(shí)交易條件的投資組合問題，相對(duì)于其余算法更具備普遍性。

4 結(jié) 語

本文為了解決有限資產(chǎn)的條件下，對(duì)具有不同風(fēng)險(xiǎn)收益特征的證券進(jìn)行配置的問題，首先建立了考慮交易成本的M-CVaR投資組合模型；其次，引入了改進(jìn)蟻獅算法對(duì)于該模型進(jìn)行求解；最后，通過對(duì)多種證券展開實(shí)證分析，證明了改進(jìn)蟻獅算法求解的可行性以及考慮考慮交易成本對(duì)于求解最優(yōu)投資組合求解的積極性與必要性。目前的研究都假定投資者為理性經(jīng)濟(jì)人，如何應(yīng)對(duì)投資人為非理性情況的研究將是未來繼續(xù)研究和努力的方向。