查 靜 虞文美

(1.安徽財貿(mào)職業(yè)學(xué)院財會金融學(xué)院， 安徽合肥 230000；2.安徽財經(jīng)大學(xué)金融學(xué)院， 安徽蚌埠 233000)

一、引言

全球各類經(jīng)濟體系內(nèi)部相互之間的金融關(guān)系日漸增多，金融風(fēng)險和危機也表現(xiàn)出系統(tǒng)性關(guān)聯(lián)。金融市場的開放性與融合性，在增加資本運轉(zhuǎn)、加強資源配合效率的同時，也給金融市場帶來了機遇和挑戰(zhàn)[1]。在金融一體化階段，若市場經(jīng)濟發(fā)生金融危機時，其蔓延速度極快，系統(tǒng)性的金融風(fēng)險很有可能導(dǎo)致價值損失，致使人們對金融系統(tǒng)失去信心，發(fā)生支付危機、貨幣貶值危機，繼而使全球性經(jīng)濟下滑，所呈現(xiàn)的嚴重后果不堪設(shè)想[2]?，F(xiàn)階段，我國潛在的金融風(fēng)險實質(zhì)上是體制風(fēng)險，涉及跨經(jīng)濟領(lǐng)域、跨市場的金融活動。簡單地依靠機構(gòu)監(jiān)管已經(jīng)很難使金融風(fēng)險得到妥善解決，如若不做出根本性改變，建立一個穩(wěn)定長效的監(jiān)管機制，定會增加突發(fā)性金融風(fēng)險的價值成本和社會負面效應(yīng)[3]。

審慎評估的出發(fā)點是要改善傳統(tǒng)貨幣舉措和微觀監(jiān)管在抵抗風(fēng)險方面的不足，把金融業(yè)當(dāng)成一個有機整體，評估金融業(yè)存在的未知風(fēng)險，評估的目標(biāo)要注重金融體系在目前及未來時段經(jīng)濟運轉(zhuǎn)的深層面問題。在實行經(jīng)濟可持續(xù)發(fā)展的同時，預(yù)防和監(jiān)管跨行業(yè)、跨期限的整體金融風(fēng)險，維系金融平穩(wěn)發(fā)展是審慎評估的最終目標(biāo)。

可預(yù)測的風(fēng)險亟須密切監(jiān)視，不可預(yù)測的風(fēng)險同樣不能疏忽遺漏[4]。由于金融市場對國民經(jīng)濟占據(jù)舉足輕重的地位，因此怎樣有效地權(quán)衡、控制風(fēng)險已成為學(xué)術(shù)界共同研討的核心內(nèi)容。

本文提出一種基于納入概率分布的金融風(fēng)險審慎評估模型設(shè)計。首先建立納入概率分布的風(fēng)險評估模型，其次使用基于小波分析與非參數(shù)估計的參數(shù)厘定，進一步確定金融風(fēng)險的損失數(shù)目，使用貝葉斯后驗概率分布得到馬爾科夫鏈，對未知變量模擬，用穩(wěn)定分布抽樣點推算蒙特卡洛積分，增強金融風(fēng)險審慎評估的有效性。

二、基于納入概率分布的風(fēng)險評估模型

在本文中，提出兩種納入概率分布族[5]。使用Q1(α，β，δ，μ)代表第一類納入分布概率族，其具體表達式可描述為：

(1)

其中，α，β，δ，μ表示參數(shù)，且滿足α，β，δR+，μR，上標(biāo)識′(α)′在α>0的情況下可用下列算子進行表達：

(2)

在式(1)中，全部參數(shù)均擁有直接的現(xiàn)實含義，μ表示位置參數(shù)，δ為尺度參數(shù)，β參數(shù)大小可衡量尾部對稱性。在β等于1的情況下，表明該納入分布為均衡狀態(tài)，在β處于β<(>)1的狀態(tài)下，表明納入分布的右尾厚于左尾。參數(shù)α能夠立體地展現(xiàn)出納入分布的尾部厚度，若α的值較小，則表示納入分布的尾部更厚[6]。

第一類納入概率分布族QⅠ(α，β，δ，μ)具備明顯的分布密度，可將其表示為：

(3)

式中，x(-∞,μ)U(μ+∞)，由于第一類納入概率分布族擁有明顯的分布密度方程式，因此能夠直接使用極大似然估算方法推導(dǎo)出參數(shù)α，β，δ，μ。

第二類納入概率分布族使用QⅡ(y;α-，α+，β-，β+，μ)進行表達，其核心作用是擬合具備特殊不平穩(wěn)尾部特征的數(shù)據(jù)分布狀態(tài)，其函數(shù)表達式可描述為：

qⅡ(y;α-，α+，β-，β+，μ)

(4)

式中，模型參數(shù)α-，α+，β-，β+R+，且μR，同時最為關(guān)鍵的是參數(shù)α-≤1、α+≤1的概率分布狀態(tài)。和第一類分布相同，參數(shù)α-、α+也是用來衡量左右尾部的厚度情況。

使用上述納入概率分布模型，可簡單地估算出金融資產(chǎn)獲益與相應(yīng)的風(fēng)險系數(shù)[7]。假設(shè)用R(t)代表某項資產(chǎn)組合在持有期T中的收益率，VaRθ為其資產(chǎn)組合在置信水準(zhǔn)1-θ(0<θ<1)下的風(fēng)險系數(shù)，那么可得到：

VaRθ=-q(θ)

(5)

式中，q(y)是R(t)的樣本納入概率函數(shù)。

三、基于小波分析與非參數(shù)估計的參數(shù)厘定

在構(gòu)建納入概率分布風(fēng)險評估模型后，為明確金融風(fēng)險造成的價值損失，繼而保證金融產(chǎn)業(yè)的平穩(wěn)運行，在此提出基于小波分析與非參數(shù)估計的參數(shù)厘定方法，具體實現(xiàn)內(nèi)容如下：

小波分析是時域分析的一種，其在時域與頻域均具備優(yōu)秀的局部化特征，擁有多分辨率解析優(yōu)勢。小波轉(zhuǎn)換在低頻區(qū)域擁有較高的頻度分辨率及較低的時段分辨率，在高頻區(qū)域擁有較高的時段分辨率及較低的頻度分辨率，是信號分析顯微鏡[8]。使用小波分析原則，可以將信號分解，將其劃分為多層傳輸至頻度信道內(nèi)。由于分解之后信號的頻度因子比初始信號更加單一，并且小波分解將信號進行分解，分解后信號的穩(wěn)定性要高于初始信號，證明小波分析具有優(yōu)越性。

多分辨分析是一種將信號空間分解的方法，在此前提下，獲得Mallat算法。利用該算法能夠?qū)⑿盘枌嵭袑蛹壏纸?，令各個層分解的結(jié)果為上一次分解獲得的低頻信號，然后再將其分解成低頻與高頻，具體方法為：

(6)

(7)

其中，j表示分解尺度，k、m表示位移系數(shù)，Aj,m表示低頻區(qū)域尺度系數(shù)，Dj+1,k表示高頻區(qū)域小波系數(shù)，h0、h1依次為低通與高通濾波器。

使用分解后的小波系數(shù)能夠重建初始序列，希波系數(shù)的重建表達式可描述為：

(8)

從第一層進行分解，其結(jié)果用于高頻區(qū)域D1與低頻區(qū)域A1，其次針對低頻區(qū)域進行更深入分解，其結(jié)果擁有高頻區(qū)域D2與低頻區(qū)域A2。以此類推，持續(xù)將信號采取分解，通過四次分解后，初始信號A被分解成：

A=A4+D4+D3+D2+D1

(9)

式中，D1、D2、D3、D4依次表示第一層至第四層通過分解所獲得的高頻信號，A4表示第四層分解獲取的低頻信號，分解層數(shù)的個數(shù)可按照實際情況決定。

若X1，X2，…，Xn為一元持續(xù)全局樣本，在隨機點x位置的核密度函數(shù)f(x)可表示為：

(10)

(11)

(12)

(13)

多數(shù)情況下，σ是使用樣本標(biāo)準(zhǔn)差S替換，同時，由于真實數(shù)據(jù)通常與正態(tài)分布呈偏離狀態(tài)，因此系數(shù)1.06降為0.8會更加合理。

四、基于蒙特卡洛算法的模型求解

計算金融風(fēng)險價值損失，可以直接觀測金融風(fēng)險程度。將風(fēng)險模型進行求解可進一步提升金融風(fēng)險審慎評估，將資金利用率發(fā)揮到最優(yōu)。

(14)

(15)

式中，St為t時段的資產(chǎn)數(shù)目，μ為平均值，σ是資產(chǎn)的浮動率，ε是隨機變化參數(shù)。處于正態(tài)分布的設(shè)想前提下，σ可描述為標(biāo)準(zhǔn)差，ε是隨機變化參數(shù)。在蒙特卡洛模擬算法內(nèi)，能運用σ代表收益率的浮動性。

使用極大似然估計，經(jīng)過迭代所獲得的估算值要符合漸進性與相同性，這一結(jié)果要在大樣本情況才會成立。

蒙特卡洛算法的出現(xiàn)為處理貝葉斯后驗積分難以計算的問題，提供新的運算思路。MCMC方法把隨機過程內(nèi)的馬爾科夫鏈加入蒙特卡洛算法內(nèi)，將抽樣概率分布隨蒙特卡洛算法的計算過程而不斷發(fā)生變化，把動態(tài)模擬轉(zhuǎn)化為真實運算[12]。MCMC方法能利用下面推導(dǎo)步驟進行演算。

將隨機變化參數(shù)θ處于t時段描述為θt，θt也是狀態(tài)空間Θ內(nèi)的一個值，若θ在Θ中不同數(shù)值之間的移動概率只依靠于θ的目前狀態(tài)，和之前時段的空間狀態(tài)無關(guān)，若處于式(16)狀態(tài)下，那么就將隨機序列{θt，t≥0}當(dāng)成馬爾科夫鏈。

Pr(θt+1|θ1…θt)=Pr(θt+1|θ1)

(16)

假設(shè){θ(t)}t≥0是狀態(tài)空間Θ內(nèi)的齊次馬爾科夫鏈，則能夠得到移動納入概率函數(shù)：

(17)

式中，p(θ，θ*)表示移動核。

針對某個分布π(θ)，經(jīng)過較長時間后，馬爾科夫鏈會漸漸忽視其原始狀態(tài)θ0，表明π(θ)既不依靠原始狀態(tài)，也不依靠時段t，則θt的邊際分布為π(θ)。

因此，不管原始數(shù)值是什么，θt會慢慢約束至同一個分布內(nèi)，也就是穩(wěn)定分布。而在約束沒有出現(xiàn)的前m次迭代中，π(θ)還不是任意狀態(tài)的邊際分布，應(yīng)將前面m次迭代數(shù)值進行剔除，僅使用后面n-m個迭代值采取估算，可將其描述為：

(18)

對于模擬推導(dǎo)角度而言，建立的移動核將已知的納入概率分布π(θ)變?yōu)榉€(wěn)定分布。若t維隨機矢量θ=(θ1，…，θt)擁有聯(lián)合分布π(θ1，…，θt)，那么函數(shù)[f(θ)]的后驗數(shù)學(xué)期望可表示為：

(19)

式(19)內(nèi)的積分運算過程較為繁雜，因此可利用MCMC方法進行簡化計算，具體步驟為：

(1)在θ內(nèi)設(shè)計出恰當(dāng)?shù)鸟R爾科夫鏈，移動核是P( )，讓其對照穩(wěn)定分布π(θ)；

(2)通過θ內(nèi)某一點θ0，使用步驟(1)內(nèi)的馬爾科夫鏈生成點序列θ1，…，θn；

(3)剔除估算E[f(θ)]之前m個預(yù)判迭代數(shù)值，那么隨機函數(shù)f(θ)的期望估計可描述成：

(20)

使用MCMC方法對未知變化參數(shù)采取抽樣模擬，通過一定的抽樣次數(shù)后，馬爾科夫鏈逐漸趨向于平穩(wěn)狀態(tài)，此時的穩(wěn)定狀態(tài)分布為求解出的后驗分布，再使用穩(wěn)定狀態(tài)分布內(nèi)的抽樣點推算蒙特卡洛積分，完成精確的金融風(fēng)險審慎評估模型設(shè)計，對風(fēng)險擁有正確估判。

五、仿真實驗

使用基于MATLAB平臺，對本文方法、EVT-POT方法和AR(1)-GARCH(1，1)-SKST方法(簡稱SKST方法)采取仿真實驗，以此驗證本文方法的優(yōu)越性。

圖1 風(fēng)險控制對比

在圖1中展現(xiàn)了近十幾年的深成指數(shù)收益率浮動情況，在3 562個測量指數(shù)比較中，EVT-POT、SKST方法的突破次數(shù)依次是58次和62次，突破率是1.3685%和1.4896%。而使用本文方法的突破數(shù)量為29次，其突破率為0.5267%，明顯優(yōu)于其他兩種方法，可有效對金融風(fēng)險采取精準(zhǔn)控制，將風(fēng)險系數(shù)降到最低。

除突破率外，對三種方法進行損失函數(shù)對比，驗證其金融風(fēng)險概率，具體如表1所示。

表1 損失函數(shù)分析

從表1中可以看出，由于EVT-POT、SKST方法均為現(xiàn)階段較好的參數(shù)方法，所以都能通過無條件覆蓋測試、獨立性測試和條件覆蓋測試，但本文方法的損失函數(shù)值最低，證明本文方法為最優(yōu)損失函數(shù)，而其他兩種方法具備一定的局限性，不能兼顧突破次數(shù)和突破過程中的損失規(guī)模，也由此證明本文方法的實用性。

六、結(jié)論

金融領(lǐng)域近年來在全球發(fā)展迅猛，但其中也包含大量不確定性及風(fēng)險，為提升金融風(fēng)險預(yù)判精度，本文提出一種基于納入概率分布的金融風(fēng)險審慎評估模型設(shè)計。審慎評估在金融風(fēng)險機制中發(fā)揮重要作用，可以極大改善經(jīng)濟風(fēng)險帶來的損失。通過建立納入概率分布風(fēng)險評估模型，確定風(fēng)險相關(guān)系數(shù)；利用小波分析與非參數(shù)估計相結(jié)合的參數(shù)厘定方法，保證金融產(chǎn)業(yè)穩(wěn)步運轉(zhuǎn)；采取蒙特卡洛算法對模型進行求解，實現(xiàn)金融風(fēng)險的有效監(jiān)管。通過仿真實驗，證明本文方法擁有較強的實用性，魯棒性較高。