波动率于garch模型

1、1.1.波动率波动率是用来描述证券价格、市场指数、利率等在它们均值附近上下波动幅度的术语，是标的资产投资回报率的变化程度的度量。股票的波动率是用于度量股票所提供收益的不确定性。股票通常具有15%-50%之间的波动率。股票价格的波动率可以被定义为按连续复利时股票在1年内所提供收益率的标准差。当t很小时，近似的等于在t时间内股票价格变化百分比的方差。这说明t近似的等于在t时间内股票价格变化百分比的标准差。由标准差来表述股票价格变化不定性的增长速度大约为时间展望期长度的平方根（至少在近似意义下）。1.2.由历史数据来估计波动率为了以实证的方式估计价格的波动率，对股票价格的观察通常是在固定的时间区间内

2、（如每天、每星期或每个月）。定义n+1观测次数；Si 第i个时间区间结束时变量的价格，i=0，1，n；时间区间的长度，以年为单位。令 1.2.1ui的标准差s通常估计为 1.2.2或 1.2.3其中u为的均值。由于的标准差为。因此，变量s是的估计值。所以本身可以被估计，其中 可以证明以上估计式的标准差大约为。在计算中选择一个合适的n值并不很容易。一般来讲，数据越多，估计的精确度也会越高，但确实随时间变化，因此过老的历史数据对于预测将来波动率可能不太相干。一个折中的方法是采用最近90180天内每天的收盘价数据。另外一种约定俗成成俗的方法是将n设定为波动率所用于的天数。因此，如果波动率是用于计算量

3、年期的期权，在计算中我们可以采用最近两年的日收益数据。关于估计波动率表较复杂的方法涉及GARCH模型与EWMA模型，在下文中将进行详细介绍。1.3.隐含波动率首先对于一个无股息股票上看涨期权与看跌期权，它们在时间0时价格的布莱克-斯科尔斯公式为 1.3.11.3.2式中 函数N(x)为标准正态分布变量的累积概率分布函数。式中：c与p分别为欧式看涨期权与看跌期权的价格，S0为股票在时间零的价格，K为执行价格，r为以连续复利的无风险利率，为股票价格的波动率，T为期权的期限。在布莱克-斯科尔斯定价公式中，不能直接观察到的参数只有股票价格的波动率。在前文中已经讨论了如何由股票的历史价格来估计波动率。在

4、实际中，交易员通常使用所谓的隐含波动率（implied volatility）。这一波动率是指由期权的市场价格所隐含的波动率。为了说明隐含波动率的计算思路，假设一个不付股息股票的欧式看涨期权价格为1.875，而S0=21，K=20，r=0.1和T=0.25。隐含波动率是使得式1.3.1所给期权价格c=1.875时对应的值。不幸的是，不能直接通过直接反解式1.3.1来将表示成期权价格与其他变量S0、K、r、T和c的函数，但是可以用迭代的方法求解所隐含的值。例如，开始时令=0.20，对应这一波动率，期权价格c为1.76美元，这一价格太低。由于期权价格为的递增函数，我们需要一个较大的值。再令=0.3

5、0，对应的期权价格c为2.10美元，此值高于市价，这意味着一定介于0.2和0.3之间。接下来，令=0.25，此值对应的期权价格仍太高，所以应在0.20-0.25间。这样继续下去每次迭代都使所在的区间减半，因此我们可以计算出满足任意精确度的近似值。本例中，隐含波动率=0.235，即每年243.5%。隐含波动率可以用来测量市场上对于某一股票波动率的观点。而历史波动率是回望型（backward looking），而隐含波动率则为前瞻型（forward looking）。通常，交易员对于期权所报出的是隐含波动率，而不是期权的价格。这样做会带来许多方便，因为波动率的变化比期权价格变化更加稳定。1.4.估

6、计波动率定义为第n-1天所估计的市场变量在第n天的波动率，第n天波动率的平方为方差率（variance rate），在前面已经对如何从历史数据来估计的标准方法进行了描述。假定市场变量在天末的价格为，变量定义为在第天连续复利收益率（第天末至第天末的收益）： 利用在最近天的观察数据所计算出的每天方差率的无偏估计为 1.4.1其中为的平均值 为了监视每天方差率的变化，式1.4.1中的公式通常会有一些变动：l 被定义为市场变量在第天末与第天末的价格百分比变化 1.4.2l 为假设为零l 为代替以上三个变化对计算结果影响不大，但这些变化会使得方差公式简化成 1.4.3式中由式1.4.2给出。2.模型2.

7、1.Arch模型有一种这样的模型为： （2-1）其中为第i天以前观察值所对应的权重，取正值。当选择这些变量的时候，如果，则，也就是我们将较少的权重给予较旧的数据。权重之和必须为一，即对于式（2-1）可以做一推广。假定存在某一长期平均方差，并且应当给予该方差一定权重，这将导致以下形式的模型 （2-2）其中为长期方差率，为所对应的权重，因为权重之和仍为1，我们有此模型就是最先由Engle提出来的ARCH模型。在这一模型中，方差的估计值是基于长期平均方差以及m个观察值，观察数据越陈旧所对应的权重就越小。令，我们可以将式（2-2）写成 （2-3）2.2.指数加权移动平均（EWMA）模型 指数加权移动平

8、均模型是式（2-1）的一个特殊形式，其中权重随着时间以指数速读递减，具体地讲，其中是介于0与1之间的某一常数。 在以上特殊假设下，更新波动率公式被简化为 （2-4）一个变量第n天的波动率估计值（在第n-1天估算）由第n-1天波动率估计值（在第n-2天估算）和变量在最近一天变化百分比决定。 为了说明式（2-4）的权重以指数速读下降，我们将式（2-4）所算出的代入公式中即代入项，进一步得出重复计算，得出当m很大时，项数很小可以忽略，所以当时，式（2-4）与（2-1） 等价，对应于的权重以速度随时间向前推移而递减，每一项的权重是前一项权重与的乘积。EWMA方法的诱人之处是其仅需要相对较少的数据。对于

9、任一时刻，我们只需要记忆对当前波动率的估计以及市场变量的最新观察值。当我们得到市场变量最新观察值后，就可以计算当天价格变化的百分比，然后利用式（2-4）就可以更新方差估计。旧的方差估计与旧的市场变量可以被舍弃。2.3.GARCH(1，1)模型我们现在讨论Bollerslev于1986年提出的GARCH(1，1)模型，GARCH(1，1)模型与EWMA模型的不同就好比式（2-1）与（2-2）的不同。在GARCH(1，1)中，是由长期平均方差以及和计算得出，GARCH(1，1)的表达式为式中为对应于的权重，为对应于的权重，为对应于的权重。因为权重之和仍为1，我们有EWMA模型是GARCH(1，1)

10、模型对应于，及的特例。GARCH(1，1)模型的(1，1)表示是由最近的的观察值以及最新的方差率估计而得出。在更广义的GARCH(p，q)模型中，是最近的p个观察值及q个最新方差率估计而得出的，GARCH(1，1)是最流行的GARCH模型。令，我们可以将GARCH(1，1)模型写成在估计模型的参数时，通常会采用这种形式，一旦、和被估算，我们可以由来计算，长期方差。为了保证GARCH(1，1)模型的稳定，我们需要，否则对应于长期方差的权重会是负值。2.4.GARCH(p，q)模型由GARCH(1，1)模型我们可以推广到一般的GARCH(p，q)模型，即GARCH(p，q)模型被广泛应用于金融资产

11、收益和风险的预测，相比于ARCH模型，GARCH模型更能反映实际数据中的长期记忆性质。由于GARCH(p，q)模型是ARCH模型的扩展,因此GARCH(p，q)同样具有ARCH(q)模型的特点。GARCH模型适合在计算量不大时,方便地描述了高阶的ARCH过程,因而具有更大的适用性。3.实证部分3.1.沪股通指数收益率与上证综合指数收益率的统计性分析（1）沪股通指数收益率与上证综合指数收益率的比较图1 从图1可观察到收益率波动的“集群”现象：波动在一段时间内较小（例如从第70个观察值到第80个观察值），在另一段时间内非常大（例如从第10个观察值到第40个观察值）。本文只收集了102个数据，若数据

12、更多，则现象更显著。从下图可以看出，上证综合指数收益率曲线有与之相似的变化趋势。图2 为了能更清楚地比较二者的变化情况，现将二者在同一曲线图中画出来，结果如下：图3 从图3可以看出：在整体上，沪股通指数收益率的变化趋势与上证综合指数收益率的变化趋势并无显著差异，但是当波动幅度较大时（例如第10个观察值到第45个观察值），沪股通指数收益率的振荡幅度明显小于上证综合指数收益率；当当波动幅度较小时（例如第50个观察值到第100个观察值），沪股通指数收益率的振荡幅度明显大于上证综合指数收益率。（2）“尖峰厚尾”现象沪股通收益率时间序列柱形统计图图4图5 由图4可知，上证沪港通指数收益率时间序列均值（M

13、ean）为0.，标准差（Std.Dev.）为0.，偏度（Skewness）为-0.，接近于-1，说明该收益率时间序列分布有非常长的左拖尾性。峰度（Kurtosis）为6.，高于正态分布的峰度值3，说明收益率时间序列具有尖峰和厚尾的特征。Jarque-Bera统计量为64.18696，P值为0.，故拒绝该收益率时间序列服从正态分布的假设。 再与上证综合指数收益率进行比较（如图5），从图5中可以看出，上证沪港通指数收益率时间序列的偏度（Skewness）为-1.，峰度（Kurtosis）为6.，说明其“尖峰厚尾”的现象更加严重。并且，其Jarque-Bera统计量为75.40078，P值为0.，因

14、而拒绝该收益率时间序列服从正态分布的假设。3.2.基于GARCH(p，q)模型对沪股通收益率的波动性进行估计 以上证沪股通指数在2014年11月17日至2015年4月17日期间的日收益率为研究对象，基于GARCH(p，q)模型对沪股通收益率的波动性进行预测。通过多次对样本进行拟合，发现GARCH(1，4)模型对沪股通收益率的波动性进行估计时效果最好。（1） 考察收益率时间序列的平稳性 通过Eviews对上证沪股通指数收益率时间序列r进行ADF检验，结果如下：图6 由ADF结果（图6）可知，t统计量的值为-10.39138，对应的P值与0没有显著性差别，说明该时间序列r平稳。（2） 序列自相关性

15、和偏相关性检验 通过Eviews对上证沪股通指数收益率时间序列r进行序列自相关性和偏相关性检验，结果如下：图7从图7中可以看出，收益率时间序列的自相关和偏相关系数均落入两倍的估计标准差内，且Q-统计量对应的P值均大于置信度0.05，故该序列在5%的显著性水平上不存在显著的相关性。（3）检验ARCH效应由于序列r不存在显著的相关性，因此将均值方程设定为白噪声，设立回归模型为：再将r去均值化得到残差w,即本文利用残差的平方相关图来检验ARCH效应，令 再通过Eviews对时间序列rt进行序列自相关性和偏相关性检验,结果如下：图8 由图8可知，Q-统计量对应的P值小于置信度0.05，故该序列在5%的

16、显著性水平上存在显著的自相关性，所以有ARCH效应。(4) GARCH(p,q)模型比较表1：GARCH(p,q)模型估计参数对应的P值GARCH(1，1）GARCH(1，2）GARCH(1，3）GARCH(1，4）GARCH(1，5）GARCH(1，6）C0.51110.08620.291600.34720.1969RESID(-1)20.45340.29110.05970.05090.00430.1672GARCH(-1)0.39950.0592000.00020.2379GARCH(-2)0.2212000.00020.4706GARCH(-3)0.0006000.7771GARCH(-

17、4)00.00010.8925GARCH(-5)0.09070.8175GARCH(-6)0.6364AIC值-5.-5.-5.10792-5.14277-5.-5.SC值-5.-5.-4.-4.-4.-4.表1表2：GARCH(p,q)模型估计参数对应的P值GARCH(1，1）GARCH(2，1）GARCH(3，1）GARCH(2，2）C0.51110.10410.6770.1536RESID(-1)20.45340.93870.60370.9405RESID(-2)20.39950.01580.13440.0168GARCH(-1)0.56710.11320.604GARCH(-2)0.1510.966