金融工程大作业(股票资产组合在险价值单位根检验格兰杰因果脉冲响应等).docVIP

金融工程大作业

三只股票、板块以及大盘数据选取：

X1=宜华地产；年收益率u1=0.310；收益率日波动sgm1=0.0221；

X2=中青宝；年收益率u2=4.688；收益率日波动sgm2=0.0479；

X3=金鹰股份；年收益率u3=0.164；收益率日波动sgm3=0.0266；

HS300=沪深300；年收益率u4=0.003；收益率日波动sgm4=0.0142;

FDCKF=房地产开发；年收益率u5=0.010231；收益率日波动sgm5=0.01585。

—2013.12.26，共250个交易日。

1、根据一个实际的股票交易账户，计算其30天后的，总资产的期望值、标准差、VaR〔30,99%〕

n=10000;

u1=0.310/250;

u2=4.688/250;

u3=0.164/250;

sgm1=0.0221;

sgm2=0.0479;

sgm3=0.0266;

s10=160640;

s20=82260;

s30=28300;

t=30;

d={10.070.32,0.071-0.03,0.32-0.031};c=chol(d);

y=rndn(n,3);x=(c*y);

z1=s10*exp((u1-0.5*sgm1*sgm1)*t+sgm1*sqrt(t)*x[.,1]);

z2=s20*exp((u2-0.5*sgm2*sgm2)*t+sgm2*sqrt(t)*x[.,2]);

z3=s30*exp((u3-0.5*sgm3*sgm3)*t+sgm3*sqrt(t)*x[.,3]);

z=z1+z2+z3;

e=sumc(z)/n;sg=stdc(z);

printesg;

340039.3744789.828

计算VaR〔30，99%〕：

y=sortc(z,1);

var=e-y[100];

printvary[100]e;

85689.749254349.62340039.37

2、根据上面的计算结果，给出一个资产组合的优化方案，并说明理由。

n=10000;

u1=0.310/250;

u2=4.688/250;

u3=0.164/250;

sgm1=0.0221;

sgm2=0.0479;

sgm3=0.0266;

s10=187100;

s20=85000;

s30=0;

t=30;

d={10.070.32,0.071-0.03,0.32-0.031};c=chol(d);

y=rndn(n,3);x=(c*y);

z1=s10*exp((u1-0.5*sgm1*sgm1)*t+sgm1*sqrt(t)*x[.,1]);

z2=s20*exp((u2-0.5*sgm2*sgm2)*t+sgm2*sqrt(t)*x[.,2]);

z3=s30*exp((u3-0.5*sgm3*sgm3)*t+sgm3*sqrt(t)*x[.,3]);

z=z1+z2+z3;

e=sumc(z)/n;sg=stdc(z);

printesg;

343484.7044800.512

由于u1u3，sgm1sgm3，因此可以通过降低金鹰股份的比重提高宜华地产的比重，使得方差变动很小的情况下较大的提高收益。

3、假设该投资者承受不了上面方案的风险，打算把的VaR〔30,99%〕值降低一半，需要建立几手〔可以小数〕股指期货的空头。

n=10000;

r=0.03/250;

u1=0.003/250;

u2=0.310/250;

u3=4.688/250;

u4=0.164/250;

sgm1=0.0142;

sgm2=0.0221;

sgm3=0.0479;

sgm4=0.0266;

s10=600000;

s20=160640;

s30=82260;

s40=28300;

k=615000;

tt=60;t=30;

d={10.50-0.030.50,0.5010.070.32,-0.030.071-0.03,0.500.32-0.031};c=chol(d);

y=rndn(n,4);x=(c*y);

z1=s10*exp((u1-0.5*sgm1*sgm1)*t+sgm1*sqrt(t)*x[.,1]);

z0=k*exp(-r*(tt-t))-z1;

z2=s20*exp((u2-0.5*sgm2*sgm2)*t+sgm2*sqrt(t)*x[.,2]);

z3=s30*exp((u3-0.5*sgm3*sgm3)*t+sgm3*sqrt(t)*x[.,3]);

z4=s40*exp((u4-0.5*sgm4*s