1、 上海金融學(xué)院“實(shí)驗(yàn)超市”實(shí)驗(yàn)報(bào)告實(shí)驗(yàn)項(xiàng)目名稱： 金融VaR計(jì)算 實(shí)驗(yàn)指導(dǎo)教師： 元如林 學(xué)生姓名： 陳莉 學(xué)生所在院系： 保險(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í)間的開盤價(jià)，留下250個(gè)數(shù)據(jù)（2012年12月9日至2014年12月31日的數(shù)據(jù)）作為檢驗(yàn)數(shù)據(jù)及建立模型。收盤價(jià)與收益率的圖形如圖1和圖2。圖1圖2（二）計(jì)算實(shí)驗(yàn)和實(shí)驗(yàn)結(jié)果1、直接法： 對(duì)2012年12月9日至2014年12月31日的數(shù)據(jù)進(jìn)行測(cè)算。均值：-0.000555943標(biāo)準(zhǔn)差：0.025888746直接法測(cè)算結(jié)果如圖3：圖32、移動(dòng)平均法：對(duì)2012年12月9日至2014年12月31日的數(shù)據(jù)進(jìn)行測(cè)算。（1）用office進(jìn)行測(cè)算，測(cè)算結(jié)果如圖4：圖4（2）用Mathlab進(jìn)行測(cè)算：對(duì)2012年1月1日至2014年12月31日的數(shù)據(jù)進(jìn)行測(cè)算。使用代碼如下：data=xlsr

3、ead(D:chen.xls);n=size(data,1);d=data(1:n);m=100;for i=1:n-1 x(i) = (d(i+1)-d(i)/d(i);end y1=0;for i=1:my1=y1+x(i);endmu(1)=y1/m;for i = 2:n-m-1 mu(i) = mu(i-1)-(x(i-1)/m)+(x(m+i-1)/m);endfor i=1:n-m-1 xigma1=0; for j=1:m xigma1=xigma1+(x(i+j-1)-mu(i)*(x(i+j-1)-mu(i); end xigma1=xigma1/(m-1); xgm(i)

4、=sqrt(xigma1); var(i)=mu(i)-1.96*xgm(i);end mt=1:n-m-1;xx=x(m+1:n-1);plot(t,xx,k-,t,mu,r-,t,mu+var,b-)l 置信度為99%，m=160時(shí)，測(cè)算結(jié)果如圖5：圖5l l置信度為97.5%，m=100時(shí)，測(cè)算結(jié)果如圖6：圖6l 置信度為95%，m=160時(shí)，測(cè)算結(jié)果如圖7：圖73、蒙特卡洛模擬法：對(duì)2012年12月9日至2014年12月31日的數(shù)據(jù)進(jìn)行測(cè)算。 測(cè)算代碼如下：data=xlsread(d:chen.xls);n=size(data,1);d=data(1:n);r=price2ret(d

5、);arf=0.025;kn=10000x=r(1:n-251);spec=garchset(R,1,M,1,P,1,Q,1,Display,off);coeff=garchfit(spec,x)y=garchsim(coeff,250,kn,60);yy=y;yyyy=sort(yy);kk=arf*kn;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:250 if rr(i)v1(i) flag=flag+1; b

6、v(flag)=n-251+i; bv(flag)=i; endendflagbv（1）置信度為99%，模擬次數(shù)kn=10000, 用ARMAX(1,1,0)和GARCH(2,2)模型，正態(tài)分布。結(jié)果如下：kn = 10000coeff = Comment: Mean: ARMAX(1,1,0); Variance: GARCH(1,1) Distribution: Gaussian R: 1 M: 1 C: -1.4754e-004 AR: -0.3097 MA: 0.4051 VarianceModel: GARCH P: 1 Q: 1 K: 3.9032e-004 GARCH: 0.25

7、07 ARCH: 0.1236 Display: offflag = 0bv = 0 23 88 243 247 圖形如圖8：圖8（2）（1）置信度為97.5%，模擬次數(shù)kn=10000, 用ARMAX(1,1,0)和GARCH(2,2)模型，正態(tài)分布。結(jié)果如下：kn = 10000coeff = Comment: Mean: ARMAX(1,1,0); Variance: GARCH(1,1) Distribution: Gaussian R: 1 M: 1 C: -1.4754e-004 AR: -0.3097 MA: 0.4051 VarianceModel: GARCH P: 1 Q:

8、 1 K: 3.9032e-004 GARCH: 0.2507 ARCH: 0.1236 Display: offflag = 2bv = 7 88 88 243 247 圖形如圖9：圖9（3）置信度為95%，模擬次數(shù)kn=10000, 用ARMAX(1,1,0)和GARCH(2,2)模型，正態(tài)分布。結(jié)果如下：kn = 10000coeff = Comment: Mean: ARMAX(1,1,0); Variance: GARCH(1,1) Distribution: Gaussian R: 1 M: 1 C: -1.4754e-004 AR: -0.3097 MA: 0.4051 Vari

9、anceModel: GARCH P: 1 Q: 1 K: 3.9032e-004 GARCH: 0.2507 ARCH: 0.1236 Display: offflag = 5bv = 7 23 88 243 247 圖形如圖10：圖10（三）結(jié)果的比較分析下表是巴塞爾委員會(huì)和國(guó)際清算銀行（BCBS）規(guī)定的懲罰區(qū)。如表2：區(qū) 域超限次數(shù)擴(kuò)大因子提高比例綠燈區(qū)0 - 40黃燈區(qū)50.460.570.6580.7590.85紅燈區(qū)10及以上1表2各種模型方法的超限次數(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ò)課程開始的舉例論證，了解了運(yùn)用Var模型進(jìn)行風(fēng)險(xiǎn)測(cè)量的重要性。在實(shí)驗(yàn)中接觸到了resset，新華08等多種數(shù)據(jù)庫(kù)，并學(xué)會(huì)運(yùn)