金融VaR實(shí)驗(yàn)報(bào)告

1、叮叮小文庫(kù)上海金融學(xué)院“實(shí)驗(yàn)超市”實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)項(xiàng)目名稱：金融VaR計(jì)算實(shí)驗(yàn)指導(dǎo)教師：元如林學(xué)生姓名：陳莉?qū)W生所在院系： 保險(xiǎn)學(xué)院學(xué)生專業(yè)：保險(xiǎn)學(xué)實(shí)驗(yàn)時(shí)間：實(shí)驗(yàn)教學(xué)與教育技術(shù)中心制實(shí)驗(yàn)?zāi)康?通過(guò)本實(shí)驗(yàn)，我理解度量金融風(fēng)險(xiǎn)的 VaR模型，了解國(guó)內(nèi)外主要的金融 數(shù)據(jù)庫(kù)，學(xué)習(xí)國(guó)際先進(jìn)的金融計(jì)算軟件的使用方法， 初步掌握金融數(shù)據(jù)采集整理，模型 選 擇，模型參數(shù)確定，VaR計(jì)算，計(jì)算結(jié)果分析的基本方法。實(shí)驗(yàn)過(guò)程（一）數(shù)據(jù)準(zhǔn)備對(duì)2012.01.012014.12.31期間債券代碼為600550的保變電氣股票進(jìn)行測(cè)算。一共在網(wǎng)易（網(wǎng)易首頁(yè)網(wǎng)易財(cái)經(jīng)行情滬深中國(guó)石油 資金流向歷史交易數(shù)據(jù)）下載了 751個(gè)該股票在

2、相應(yīng)時(shí)間的開(kāi)盤價(jià)，留下 250個(gè)數(shù)據(jù)（2012年12月9 日至2014年12月31日的數(shù)據(jù)）作為檢驗(yàn)數(shù)據(jù)及建立模型。收盤價(jià)與收益率的圖形如圖1和圖厶收盤價(jià)（J計(jì)算實(shí)驗(yàn)和實(shí)驗(yàn)纟占杲1、直接法：對(duì)2012年12月9日至2014年12月31日的數(shù)據(jù)進(jìn)行測(cè)算。均值：-0.000555943 標(biāo)準(zhǔn)差：0.025888746星信水平1* CLot(X分位點(diǎn)Zot99%0.012.3397.5%0.0251,9695%0.051,65直接法測(cè)算結(jié)果如圖3:5駕0 5 00ZLW& rM黒-51亠4-ZCHl甘6325收益率99% VaR97.5%VaR9S%VaR2、移動(dòng)平均法：對(duì) 2012年12月9日至2

3、014年1231日的數(shù)據(jù)進(jìn)行測(cè)算。(1)用office進(jìn)行測(cè)算，測(cè)算結(jié)果如圖 4:U.15 at0.10.053r*jo CM 015(2)用Mathlab進(jìn)行測(cè)算：對(duì)2012年1月1日至2014年12月31日的數(shù)據(jù)進(jìn) 行測(cè)算。使用代碼如下：data=xlsread(D:che n .xls);n=size(data,1);d=data(1: n);m=100;for i=1: n-1x(i) = (d(i+1)-d(i)/d(i);endy仁0;for i=1:my1=y1+x(i);endmu(1)=y1/m;for i = 2:n-m-1mu(i) = mu(i-1)-(x(i-1)/m

4、)+(x(m+i-1)/m);endfor i=1: n-m-1 xigma1=0;for j=1:mxigmal 二 xigma1+(x(i+j-1)-mu(i)*(x(i+j-1)-mu(i); endxigma1=xigma1/(m-1);xgm(i)=sqrt(xigma1);var(i)=mu(i)-1.96*xgm(i);17endmt=1: n-m-1; xx=x(m+1: n-1);plot(t,xx,k-,t,mu,r-,t,mu+var,b-)置信度為99%, m=160時(shí)，測(cè)算結(jié)果如圖5:0.150.10.050-0.05600置信度為97.5%, m=100時(shí)，測(cè)算結(jié)果

5、如圖6:0.15-0.10.10.050-0.050 100 200300400500600700圖6置信度為95%, m=160時(shí)，測(cè)算結(jié)果如圖7:0.150.10.05-0.05-0.1 J0 100 200f【3004005006003、蒙特卡洛模擬法：對(duì) 2012年12月9日至2014年12月31日的數(shù)據(jù)進(jìn)行測(cè)算。 測(cè)算代碼如下：data=xlsread(d:che n. xls);n=size(data,1);d=data(1: n);r=price2ret(d);arf=0.025;kn=10000x=r(1: n-251);spec=garchset(R,1,M,1,P,1,Q,

6、1,Display,off);coeff=garchfit(spec,x)y=garchsim(coeff,250,k n ,60);yy=y;yyyy=sort(yy);kk=arf*k n;var1=yyyy(kk:kk,:);v1=var1;rr=r( n-250: n-1);u(1:250)=0;t=1:250;plot(t,rr,b-,t,v1,r-,t,u,g-)flag=0;bv(1)=0;for i=1:250if rr(i)v1(i)flag=flag+1;bv(flag)=n-251+i;bv(flag)=i;endendflagbv(1)置信度為 99%，模擬次數(shù) kn=

7、10000,用 ARMAX(1,1,0)和 GARCH(2,2)模型，正態(tài)分布。結(jié)果如下：kn =10000coeff =Comment: Mean: ARMAX(1,1,0); Varianee: GARCH(1,1)Distributi on: Gaussia nR: 1M: 1C: -1.4754e-004AR: -0.3097MA: 0.4051Varia n ceModel: GARCHP: 1Q: 1K: 3.9032e-004GARCH: 0.2507ARCH: 0.1236Display: off*flag = bv =02388243247250圖形如圖8:bHfe-T.1i

8、i111，.Ii.1J1111 1 1h 1L1； rif0.080.060.040.02 0-0.02-0.04-0.06-0.08050100150200圖8(2) ( 1)置信度為 97.5%,模擬次數(shù) kn=10000,用 ARMAX(1,1,0)和 GARCH(2,2)模 型，正態(tài)分布。結(jié)果如下：kn =10000coeff =Comme nt: Mea n: ARMAX(1,1,0); Varia nee: GARCH(1,1)Distributi on: Gaussia nR: 1M: 1C: -1.4754e-004AR: -0.3097MA: 0.4051Varia n ce

9、Model: GARCHP: 1Q: 1K: 3.9032e-004GARCH: 0.2507ARCH: 0.1236Display: offflag = bv =8888243247圖形如圖 9(3)置信度為 95%，模擬次數(shù) kn=10000,用 ARMAX(1,1,0)和 GARCH(2,2)模 型，正態(tài)分布。結(jié)果如下：kn =10000coeff =Comme nt: Mea n: ARMAX(1,1,0); Varia nee: GARCH(1,1)Distributi on: Gaussia nR: 1M: 1C: -1.4754e-004AR: -0.3097MA: 0.4051

10、Varia n ceModel: GARCHP: 1Q: 1K: 3.9032e-004GARCH: 0.2507ARCH: 0.1236Display: offflag = bv =72388243247圖形如圖10:0.080.060.040.02-0.02-0.04-0.06 1050100150200250圖10（三）結(jié)果的比較分析 下表是巴塞爾委員會(huì)和國(guó)際清算銀行（BCBS ）規(guī)定的懲罰區(qū)如表2:區(qū)域超限次數(shù)擴(kuò)大因子提高比例綠燈區(qū)0 - 40黃燈區(qū)50.460.570.6580.7590.85紅燈區(qū)10及以上1各種模型方法的超限次數(shù)比較，如表3:保變電氣置信水平95%600550模型參數(shù)m超限次數(shù)區(qū)域參數(shù)法直接法10紅燈移動(dòng)平均法1608黃燈蒙特卡羅法5黃燈表3:回顧測(cè)試結(jié)果的分區(qū)由表可知，使用蒙特卡羅法計(jì)算金融VaR更為精確，使用性更強(qiáng)三、實(shí)驗(yàn)總結(jié)通過(guò)課程開(kāi)始的舉例論證，了解了運(yùn)用Var模型進(jìn)行風(fēng)險(xiǎn)測(cè)量的重要性。在 實(shí)驗(yàn)中接觸到了 ress