实验报告模版.docx

《金融风险管理》课程实验报告专业： 班级： 学号: 姓名： 日期： 年 月 日实验项目名称：居中，黑体，四号字）一、实验目的与要求（居中，黑体，小四号字）（简要介绍本项目实验需要解决什么问题，达到什么效果，有哪些规定性要求，比如对 实验内容、样本数据和具体方法等方面的限定；宋体，五号字）二、实验方法与原理（居中，黑体，小四号字）（简要介绍本项目实验所采用的具体方法及其基本原理；宋体，五号字）三、实验过程与步骤（居中，黑体，小四号字）（详细介绍本项目实验的主要过程及具体操作步骤，要求对实验过程的介绍必须完整， 思路清晰，前后步骤不能弄错，关键性步骤必须有截图和文字配合说明；宋体，五号字）四、实验结果与分析（居中，黑体，小四号字）（详细介绍通过前面的操作过程所得到的实验结果，并要求对实验结果进行分析，比如 是否达到了预期目的、可能还存在哪些问题、原因是什么等等，分析时尽可能深入，内容 一定要充实；宋体，五号字）五、心得体会（居中，黑体，小四号字）（这部分内容是对项目实验的进一步延伸，主要围绕实验过程中自己的感受、遇到的问 题及解决、对本实验项目设计的看法和建议等等；宋体，五号字）附：参考模版基于三种方法的VaR估算以中信银行（601998）股票为例一、实验目的与要求【实验目的】本实验项目主要在于通过对具体案例的实验操作，掌握VaR的三种估算 方法，并通过对实验结果的分析，进一步理解这三种方法的基本原理及各自的优劣利弊。

实验要求】投资者于2015年4月21日以收盘价购买一万股中信银行股票，在99%的 置信水平下，预测下一交易日该投资者可能面临的最大损失额要求用历史模拟法、蒙特卡 罗模拟法以及方差协方差法分别实验历史模拟法的样本数据不得少于250个，蒙特卡罗模 拟的次数不得低于10000次二、实验方法与原理风险估值（Value at Risk简称VaR ）是在一定置信水平和一定持有期内，某一进入资产 或组合在正常的市场条件下所面临的最大损失额，是一种用于测量和控制金融风险的量化工 具与传统测量工具相比， VaR 的优点在于其具有简明性和综合性，将市场风险概括为一 个简单的数字，便于高层管理者和监管机构理解VaR 的获取方法主要有三种，即历史模拟法、蒙特卡罗模拟法以及方差协方差法（参 数法），本实验将分别采用这三种方法来估计中信银行股票的VaR数值历史模拟法认为“历史会重演”，从而使用金融资产收益率的历史数据，根据置信度直 接从中选取一个收益率来计算VaR ；蒙特卡罗模拟法是假定金融资产的收益率服从某种分 布，再根据历史数据估计出该分布的参数，然后使用随机数发生器模拟出大量收益率的未来 可能值参数法假定金融资产收益率服从某种分布，再根据修正参数和该资产的期初价值计算 出VaR ；方差一协方差法假定金融资产收益率服从某种分布再根据修正参数和该资产的期 初价值计算出VaR。

三、实验过程与步骤【样本数据】中信银行（601998） 2012年 4月 23日至 2015年 4月 21日的收盘价，共 726 组数据，数据来源为大智慧实验工具】本实验采用的软件及版本为Matlab 7、Stata/ SE相关说明】实验变量选择上对于 PPT 上的内容进行修改了，个人对于统计方面的见 解是，绝对量的波动研究是没有意义的，因为在不同基数下变动相应的数值所需的影响大小 是不同的同时个人觉得股价波动符合这一特征，同时股票收益在于收益率而不是股价绝对 值的变化，因此本实验放弃PPT中以头寸收益来确定VaR值的方法转而使用收益率作为衡 量标准（仅随机游走模型因收益率不满足存在单位根的必要条件，才在该节采用头寸损益作 为变量）一）历史模拟法第一步，建立样本，模拟头寸模型数据选择2012年 4月 23日至2015年4月 21日的 收盘价，共 726组数据模拟头寸在每一个历史交易日中的价值，以当日股价减去前一交易 日股价并除以前一交易日股价作为当日股票的收益率，由此获得模拟头寸的的时间序列第二步，求解VaR使用Matlab的sort命令将收益率序列数据按升序排列，由于所给 定置信水平为99%，找到对应的第726xl%=7.26个数据（本实验采用插值法，对第7、8 个数进行加权平均得到所需数据），为-0.0595，即单日 VaR = 0.0595 x 8.03 x 10000=4721.9 元（99% , 1 天）。

二）蒙特卡罗模拟法1.随机游走模型第一步，单位根检验检验收益率时间序列数据是否存在单位根，本实验采用Dicky-Fuller、 Phillips-Perron 检验，分别对于收益率序列进行单位根鉴定两种检定结果均 为p = 0.0000，即在99%的显著性水平下认为该时间序列不存在单位根因此尝试头寸损 益时间序列是否存在单位根，对于单位根检验本实验采用Dicky-Fuller、Phillips-Perron检验表 单位根鉴定变量变量I (1)过程I (2)过程结论DF检定PP检定DF检定PP检定R值值 z PI (0)过程PZ值 p值I (1)过程其中，R表示收益率，P表示头寸损益第二步，干扰项分布的确定由于设定头寸损益变动服从随机游走过程，因此干扰项与 收益率服从相同的分布函数，两者仅均值不同，因此可以通过估计收益率分布并将其均值变 为零得出干扰项的分布本文证明该收益率数据存在“尖峰厚尾”的情况（如图，证明见节 相关内容），因此考虑使用混合正态分布拟合* t，以减小拟合的误差图 中信银行收益率频率分布实验首先估计收益率的分布本实验使用极大似然估计的方法以Matlab7作为工具，经 过设定初值后重复迭代直至收敛，得出收益率可以由下列混合正态分布近似替代：P 0.8689 x (4.1978,0.45457) + 0.1311x (7.1591,0.52485),从而得出干扰项的混合正态分布服从：j 〜0.8689x (-0.3883,0.45457) + 0.1311x (2.573,0.52485)。

根据求出的干扰项分布随机生成十万组数据，在给定置信水平为 99%的前提下，通过分位数计算可以计算出该分位数下干扰项为 * t = -0-5453，从而计算出收益率变动为P = -0.5453 x 10000=5453 元(99%, 1天)t(三) 参数法1. 静态估计法假定O不随时间的变化而变化，因此O =O,即O恒等于一常数由这一假设可得； t t t历史标准差与未来标准差相同，即可根据历史数据求出0.023根据公式VaR p，并假设*服从正态分布的条件下，得出VaR = 3804.61元(99% , 1天)t2. 动态估计法一一GRACH族模型(1)收益率分布特征分析中信银行 2013年 4月 23日至 2015年 4月 21日收益率统计出的频率分布如下图所示， 不难看出其与所对应的正态分布相比具有“尖峰厚尾”的特征，考虑该收益率可能不服从正 态分布，因此对于该数据是否服从正态分布加以检验1 -.05 0 .05 .1R@£ldUdl£图 中信银行收益率频率分布对于收益率是否服从正态分布，本实验采用Q-Q图和Jarque-Bera检验综合判断收益率 是否服从正态分布，结果如下图、表所示。

由 Q-Q 图可以看出该收益率在尾部和正态分布 尾部相差较大，即期分布可能不服从正态分布Jarque-Bera检验结果可以看出其峰度、偏度 以及联合检定的与正态分布是否存在差异的显著性都在 99%的置信水平上拒绝了其服从正 态分布的假设图 中信银行收益率Q-Q图表 Jarque-Bera 检验 统计值 2 p峰度偏度组合对于存在的“尖峰后尾”使得收益率不满足的情况，下文考虑收益率的随机误差项服从 t 分布或 GED 分布2)模型估计 考虑到股票市场所具有的信息对于股价波动的冲击通常具有非对称性，本实验首先尝试 使用EGARCH (1,1)和TGARCH (1,1 )模型采用极大似然估计进行参数估计，其结果如 下GED-TGRACH 结果不收敛舍去由回归结果，TGARCH效应不显著，因此排除模型(4); EGARCH效应在10%的显著性水平下成立，说明不排除存在负向信息和正向信息的影响作用非对称的存在根据 AIC准则确定的模型为t-GARCH （1,1）,而根据BIC准则确定的模型为t-EGARCH （1,1）由于GRACH 类模型成立需要满足残差为白噪声，因此对于两模型进行白噪声检验（见图），发 现残差基本都落在置信区间以内，说明模型的使用是合理的。

对于两模型进行似然比鉴定， 发现由于本实验目的在于确定单日VaR值，因此不再讨论两者的优劣，对两模型均予以接 受，分别计算出其单日VaR值D OO 0 10 0 2D a 3D D 40FrcquencrBarfvfalEilahFtiiitp- i 61 Pn* B« SJiiUOumulsilFje PenolDgrsm ^hite-No-Ee Te^10 00 0.1D OSD 03D 0 4DFrequency出■w*・■«朝 fc B-aa!n]图 t-GARCH （1,1）、t-EGARCH（1,1 ）残差白噪声检验首先，计算出在 t-GARCH（ 1,1 ）、t-EGARCH（ 1,1）下 的分别 2.7449灯0-2、 t2.9967X10-2接着，根据公式VaR = p，在取置信水平为99%的情况下，VaR = aop,分别得出单日 VaR 值为 VaR = 3837.33兀（99% , 1天）、GARCHVaR 二 4189.29元（99% ,1天）其中 t （3） = —3.4819 ）EGARCH 0.98四、实验结果与分析本次实验共采用了三大类五种方法求解VaR值，计算结果汇总如下： 表VaR值实验结果汇总 历史模拟法 蒙特卡罗模拟法VaR参数法（静态估计） 参数法（t-GRACH） 参数法（t-EGARCH）其中VaR值为在99%的置信水平下，下一交易日的最大损失额从VaR值结果来看，预测的结果比较接近，因此VaR结果具有稳定性。

在99%的置信 水平下，参数法（静态估计）所预测的下一交易日VaR值最小，为元（99%，1天），蒙特卡罗 模拟法所预测的VaR值最大为元（99%，1天）根据实验结果可以发现单日VaR值集中在[3500,5500]区间，所以说预测2015年4月22日的VaR值应该在此区间内，由于4月22日 该股上涨，因此已经涵盖了下一交易日投资者的最大损失额，因此在一定程度上说明该结果 是有效的五、心得体会本实验分别采用历史模拟法、蒙特卡罗模拟法和参数法三种方法估算了中信银行股票 的VaR值其中历史模拟法并未采用头寸损益作为变量而是用收益率替代；蒙特卡罗模拟 法从单位根检验中发现头寸收益时间序列含有一阶单位根，认为设定随机游走模型合理，设 定干扰项服从混合正态分布，给出在此设定下的VaR值，考虑到随机游走过程的不合理本 实验也给出以几何布朗运动模型描述得出的VaR值作为参考；参数法中分别采用静态、动 态方法估计 值，动态方法采用的是GARCH族模型通过本项目实验，本人不仅掌握了 VaR的三种方法估算方法，同时，还真切地体会到 VaR 系统的强大正如老师课上所说的， VaR 的最大生命力在于其灵活性，允许操作者可 以根据自。