基于遗传算法的中国股市波动性研究

提要本文介绍了金融市场波动性过程的长期记忆性特征的分整自回归条件异方差模型——

—FI G A R C H（p，d，q），并介绍了一种建立FI G A R C H模型的新方法——

—遗传算法。对上证综合指数进行实证分析，对其收益率建立自回归模型（A R 模型），由Evi ew s软件可知模型的残差项具有明显的异方差效应，应用遗传算法的基本步骤进行C语言编程，由此求解FI-G A R C H模型的各项参数值，从而建立FI-G A R C H模型，实证结果表明中国股市波动性过程具有明显的长记忆性。

关键词：G A R C H；FI G A R C H；遗传算法

中图分类号：F83文献标识码：A

股票市场是现代市场体系的有机组成部分，也是现代金融市场中最有活力的增长点。我国的股票市场起步较晚，市场制度尚不完善，从而使得股票市场的发展起伏较大。随着市场经济的日趋成熟，股票市场的进一步开放和交易品种的不断丰富，对新形势下的股票市场发展的内在规律和实践进行深入的研究是非常有必要的。

近年来的很多实证研究表明，经济、金融系统中的时间序列大多具有非线形性，也就是说这些时间序列具有长期的记忆性。表现在波动性上的长期记忆型又称持续性。由于波动性不仅是资产风险的决定因素，而且还是衍生证券定价中的一个关键参数，因此很好地理解金融时间序列的波动性具有重要意义。

国外对股票市场价格的波动特征已经进行了大量的实证研究，其中最成功的模拟方差随时间变化的模型是由Engle （1982）首先提出的自回归条件异方差模型，即ARCH模型。ARCH模型实际上是

对时间序列动态模型的推广，它将方差和

条件方差区分开来，并且定义条件方差是

过去误差的函数，为解决异方差问题提供

了新的途径。然而，在实际应用中，当

ARCH模型的阶数过大时，参数的估计则

不再精确；除此之外，为了保证方差为正，

还要求参数值为正，当参数过多时，实际

数据的估计模型往往不能满足这一点，因

此，Bollerslev（1986）在此基础上提出了广

义自回归条件异方差（GARCH）模型。大

量实验表明，金融数据扰动项异方差具有

极大的持续性，这些现象促使Engle和

Bollerslev（1986）提出了IGARCH（q）模

型，给出了单位根的许多特性。Bollerslev

等人近期经过研究又提出了FIGARCH

（p，d，q）模型，这种方法的目的就在于可

以对模型方差灵活定阶，可以更好地解释

金融市场的波动性。目前，ARCH模型和

GARCH模型已经被广泛应用于股票市

场、货币市场、外汇市场、期货市场的研究

中，来描述股票价格、利率、汇率、期货价

格等金融时间序列的波动性特征。

Bollerslev等对美国标准普尔500复合指

数进行FIGARCH建模研究，得到差分阶

数d=0.447，显著不同于0和1，表明美国

股市从波动性上表现出长期的记忆性。我

国学者李汉东、张世英和汤果、何晓群等

人分别从理论方法上和实证分析上对

FIGARCH模型进行了研究，结果表明了

我国股市收益存在长期记忆性。

本文利用遗传算法对上证综合指数

的波动建立FIGARCH（p，d，q）模型，从条

件方差上研究了中国股票市场的长期记

忆性。本文第一部分介绍了FIGARCH模

型；第二部分介绍了遗传算法以及应用遗

传算法进行C语言编程的步骤；第三部

分对上证综合指数进行FIGARCH建模；

第四部分给出了简要的结论。

一、FI G A R C H模型概述和参数估计

（一）模型概述

1、G A R C H模型。ARCH模型介绍了条

件方差的短期记忆性，GARCH模型是对

ARCH模型的扩展，因此GARCH模型具

有ARCH模型的特点，但GARCH模型的

条件方差不仅仅是滞后残差平方项的线

性函数，而且是滞后条件方差的线性函

数，模型如下：

εt=z tσt

E t-1（z t）=0，Var（z t）=1

σ2t=Var（εtψt-1），ψt是到t时刻的信

息集。

σ2t=ω＋α（L）ε2t＋β（L）σ2t式（1）

其中：ω为常数，L为推移算子。

α（L）＝α1L＋α2L2＋…＋αq L q

β（L）＝β1L＋β2L2＋…＋βp L p

ω>0；α1，…αq，β1，…βp≥0

GARCH（p，q）不仅反映了短期滞后，

而且考虑到了长期滞后。而且研究表明，

一般只需研究GARCH（1，1）模型。

2、FI G A R C H模型。FIGARCH模型是

对GARCH模型的多项特征的整合和推

广，是其更一般的情形，反映了金融时间

序列的长期记忆性，模型如下：

令v＝ε2

t

－σ2t，代入（1）可得：

［1－α（L）－β（L）］ε2

t

＝ω＋［1－β（L）］v t

式（2）

此时，v

t

是零均值不相关的时间序

列。

当此时的特征多项式1－α（L）－β（L）

＝0有一个单位根时，就得了Engle和

Bollerslev（1986）提出的IGARCH模型。如

同从ARIMA（p，d，q）模型推广到ARFI－

MA（p，d，q）模型来观察经济现象中的长

期记忆性一样，很自然的考虑到从I－

基于遗传算法的中国股市波动性研究

□文/吴蕾

GARCH （p ，q ）到FIGARCH （p ，d ，q ）的推广：

准（L ）（1－L ）d ε2

t =ω＋［1－β（L ）］v t

式

（3）

其中：0

（1－L ）d

＝F （－d ，1，1；L ）＝∞

k =0ΣГ（k-d ）Г

（k+1）-1

Г（－d ）-1L k

=∞

k =0

Σπk L k 式（4）

其中：πk ＝0＜k ≤i

仪（k-1-d ）/k

显然，当d=0时，FIGARCH （p ，d ，q ）模型就是GARCH （p ，q ）模型；当d=1时，FIGARCH （p ，d ，q ）模型就是IGARCH （p ，q ）模型。

（二）参数估计

考虑一般的FIGARCH （p ，d ，q ）：y t =x t ′b+εt εt ＝z t σt

准（L ）（1－L ）d ε2

t =ω＋［1－β（L ）］v t

式

（5）

v t ＝ε2t －σ2t v t φt-1～N （0，σ2t

）由（5

）得：σ2

t

＝ω［1－β（1）］－1

＋｛1-［1－β（L ）］－1

准（L ）（1－L ）d ｝ε

2

t

常用的FIGARCH 模型的参数估计方法是拟极大似然估计法（QMLE ），FI －GARCH （p ，d ，q ）的似然函数如下：

logL （θ；ε1，ε2，…εT ）＝－0.5Tlog （2π）－

0.5T

t =1Σ［log （σ2t ）＋ε2t σ－2t

］其中θ′＝（ω，d ，β1，β2，…，βp ，准1，准2，…，准q ），准k 是准（L ）的系数。

在对FIGARCH 模型进行参数估计时，首先应该确定差分阶数d 的值，常用的方法有四种：尝试法、GPH 方法、周期图法以及重标极差方法（R/S 分析法）。许多经济学家通过选取不同的d 值进行尝试性的研究给出一个近似最优解，这种方法的计算过程是繁琐的，而且是不科学的。GPH 方法是Geweke 与Porder Hudak 在1983年提出的一种最常见的估计阶数d 值的半参数方法，它运用滤波理论，对谱密度对数的函数进行回归而得出d 值。GPH 方法意味着d 可以通过一个简单的回归方程来估计得到，但当样本数足够大时，d 的最小二乘估计才渐进服从正态分

布，而实际中的样本量通常是有限的，尤其是起步较晚的我国股市，很难获得足够大的样本数。除此之外，这种方法所计算出的d 值通常较小，不是明显异于0和

1。周期图法是对加权的周期图量值最小

化而求出d 的估计值的一种方法，它是研究证券市场波动性的有效方法，能过滤大部分序列的相关因素，但不能完全剔除，而且由于没有考虑到宏观政策发布等因素，使得序列不相互独立。经典的R/S 分析是通过计算赫斯特指数来计算d 的一种方法，这种方法计算简单，但序列具有

短期记忆和非平稳性。为了弥补这种方法的不足，Lo （1991）又提出了修正的R/S 分析，目前已成为实证分析主要采用的方法，但这种方法所计算的d 值通常较小。除此之外，王春峰和张庆翠对中国股市波动性的长期记忆性进行研究时，在OX 统计语言环境下，应用G@RCH2.1软件包，经编程计算，也可以求解d ，并且可以对所求的d 值进行T 统计量检验，这种方法

是比较科学的，而且d 值是显著异于0和1的。我国学者李颖和汤果在理论新探上还提出了BHHH 算法和混合梯度算法，BHHH 算法计算程序比较简单，但迭代次数较多，计算效果较差；混合梯度算法迭代次数少，收敛速度快，耗时较少，可以很大地提高计算效率。本文提出了对金融时间序列建立FIGARCH 模型的一种新方法——

—遗传算法。二、基于遗传算法的C 语言算法设计（一）遗传算法介绍。遗传算法简称GA ，是1962年由美国Michigan 大学的Holland 教授提出的模拟自然界遗传机制和生物进化论而形成的一种并行随机搜索最优化的方法。他将物竞天择的生物进化原则引入优化参数形成的编码群体中，按所选择的适应函数并通过遗传中的复制、交叉和变异对个体进行筛选，使适应性强的个体得到保留，并组成下一代群体，新一代群体既继承了上一代的信息，又优于上一代。这样周而复始，群体中个体适应度不断提高，直到满足所给定的条

件。遗传算法的主要特点是采用群体搜索策略和充分利用群体中个体间的信息交换，具有全局搜索、搜索空间维数较大等

优点，尤其适用于处理传统搜索方法难于解决的复杂问题。其基本操作包括：复制、交叉和变异。

（二）C 语言算法设计

1、理论说明。本文选取模型为FI －ARCH （1，d ，1），如下所示：

σ2t ＝ω＋βσ2t －1＋［1－βL －（1－准L ）（1－L ）d

］

ε2t 其中：ω、β、准都是未知参数。

按公式（4）取（1－L ）d 的20阶泰勒展

开式，将FIGARCH （1，d ，1）模型展成含参数的GARCH 模型：

σ2t ＝ω＋βσ2t －1＋a ［0］ε2t －1＋a ［1］ε2t －2＋…＋a ［20］ε2t －21式（6）

其中：

a ［i ］含有未知参数准，i=0，1，2，…20。由已知的时间序列｛ε2t ｝按公式（6）可计算出｛σ2t ｝的真实值。

由GARCH 模型可知其一步预测为：σ赞2t （1

）＝E ［σ2t －1ε2t ，ε2t －1，…，ε2t －20］＝E ［ω＋βσ2t ＋a ［0］ε2t ＋a ［1］ε2t －1＋

…＋a ［20］ε2t －20ε2t ，ε2t －1，…，ε2t －20］

＝ω＋βσ2t ＋a ［0］ε2t ＋a ［1］ε2t －1＋…＋

a ［20］ε2t －20式（7）

用上述｛ε2t ｝的值根据公式（7）可计算出｛σ2t ｝的预测值｛σ赞2t ｝。

定义：βt ＝σ赞2t －σ2t 2、算法设计步骤步骤1：外层循环

（1

）给出下列参数的取值：种群大小（POP_SIZE ）、交叉率（P_CROSSOVER ）、变异率（P_MUTATION ）、最大进化代数

（GEN

）。（2）选择适应函数为：Q ＝T m

k =1

Σρ2k 式

（8

）其中：ρk 为βk 的自相关函数，T 为ρk

的样本数，m =T 1/2

，ρk =γk /γ0，γk ＝（1/R ）R －k

t =1Σ

（βt+k -（βt -，R 为βk 的样本数。

由于Q ～x （m ），若选取的参数满足Q ＜x 0.05（m ），则说明｛βk ｝是白噪声，即所建立的模型即消除了短期记忆性又消除了长期记忆性。

步骤2：内层循环———遗传算法（复

制、交叉、变异）

（1）初始化。在［0，1］之间随机选取POPSIZE个d值，组成向量d［POPSIZE］，对每个d［i］（1≤i≤POPSIZE），再次使用遗传算法选择参数，具体做法如下：

①初始化。给出下列参数的取值：种群大小（POP_SIZE1）、交叉率（P_CROSSOVER1）、异率（P_MUTA－TION1）、最大进化代数（GEN1）。

选择适应函数为：L（θ）＝－（T/2）log （2π）-（1/2）

T

t=1

Σ（logσ2t+ε2t/σ2t）式（9）

②复制。随机选取三组POP_SIZE1个［0，1］之间的数ω［POPSIZE1］、β［POP－SIZE1］、准［POPSIZE1］作为初始的种群，把每一组ω［i］、β［i］、准［i］（1≤i≤POP－SIZE1）及相应的d值代入公式（6），将使条件方差为正的ω［i］、β［i］、准［i］保留下来，再根据适应函数（9），按照遗传算法选择的步骤选取一组值（每个值都包含ω、β、准三个参数）作为新一代种群。

③交叉。把新一代种群两两配对，对每一对包含ω、β、准三个参数的向量都随机产生一个［0，1］之间的数，这里设为r，如果r

④变异。对每一个包含ω、β、准三个参数的向量都随机产生一个［0，1］之间的数，如果这个数小于P_MUTATION1，则发生变异，变异的方法为：随机产生一个［0，1］之间的数作为变异点，如果这一点的二进制编码为1就改为0，反之亦然。

⑤将使适应函数依次变大的参数值保留下来。

⑥循环步骤②、③、④、⑤直到达到最大进化代数GEN1，则得到的ω、β、准的值就为d［i］所对应的极大似然估计的参数估计值。

（2）复制。对每一个d［i］（1≤i≤POPSIZE）值及经过上述①-⑥所选取的相应的最优参数ω、β、准的值计算适应函数（8），同样按遗传算法的选择步骤选取一组值作为新一代种群。

（3）交叉。将选择出的新一代种群两两配对，对每一对d［i］（1≤i≤POPSIZE）

随机产生一个［0，1］之间的数，如果小于

P_CROSSOVER，则发生交叉，交叉方法同

上。

（4）变异。对每个d［i］（1≤i≤POP－

SIZE）随机产生一个［0，1］之间的数，如果

小于P_MUTATION，则发生变异，变异方

法同上。

（5）重复上述复制、交叉、变异，直到

所计算的Q值满足Q＜x

0.05

（m）（此时迭代

次数不应超过最大迭代次数GEN），则得

到的d值就为最优的差分阶数。

三、上证综合指数分析

（一）变量说明。T：样本容量，这里选

取2000年1月4日到2006年6月30日

的收盘指数，共1，555个；P

t

：每日收盘指

数，t=1，2，3，…T；R

t

=100（logP t-logP t-1）：每

日收益率，t=1，2，3，…T；ε

t

：对R

t

作确定

性分析后的残差项，t=1，2，3，…T；d：差分

阶数。

（二）数据分析：对上证综合指数每日

收盘指数进行Eviews分析，其走势如图1

所示。（图1）

收益率具有明显的聚类现象，时序图

如图2所示。（图2）

由于股指内部各种股票的非同步交

易会导致股指收益序列的自相关性显著，

为了滤除序列中这种短相关因素而突出

长相关因素，我们建立辅助自回归模型，

分析其残差序列。

对上证收益率建立辅助自回归AR

（2）模型：R

t

=0.0060＋0.0449×R t-1－0.0159×

R t-2＋εt

对模型进行SPSS分析，如表1所示。

（表1）从表1中可看出AR（2）模型中常

数项和R

t-2

的系数的尾概率分别为

0.85784494和0.52724769，说明所建立的

AR模型是不显著的。这是因为AR模型

建立的前提条件是残差项ε

t

必须是白噪

声，这也就说明上证指数收益率的残差项

不是白噪声，事实上，它具有显著的异方

差性，因此，有必要对AR模型的残差项

进行建模分析。

首先，由收益率和AR（2）模型我们可

以计算出残差项ε

t

的值，进而可得出残

差平方项ε2

t

的值，用Eviews分析如图3

所示，具有明显的异方差性。（图3）

由于我国股市具有长期记忆性，我们

对扰动项平方项ε2

t

建立FIGARCH（1，d，

1）模型，这里我们对1，555个数中连续的

900个数进行实证研究，查表可知x0.05

（30）=43.8，只要所计算的Q<43.8，就说明

所建立的模型是正确的。

应用遗传算法编程，选取参数值为：

GEN：200POP_SIZE：30

P_CROSSOVER：0.75

P_MUTATION：0.1GEN1：100

POP_SIZE1：30

P_CROSSOVER1：0.75

P_MUTATION1：0.1

对FIGARCH（1，d，1）进行参数估计

和差分阶数的计算得出结果，如表2所

示。（表2）

经过134次迭代后，Q=2.1950<43.8=

x0.05（30）。

2

-

1

-

4

2

-

8

-

1

6

2

1

-

3

-

3

2

1

-

1

1

-

1

2

2

2

-

6

-

2

8

2

3

-

2

-

1

1

2

3

-

9

-

2

4

2

4

-

5

-

7

2

4

-

1

2

-

2

2

5

-

8

-

3

2

6

-

3

-

2

2500

2000

1500

1000

500

收益指数

图1上证收盘指数走势图

2

-

1

-

4

2

-

8

-

1

6

2

1

-

3

-

3

2

1

-

1

1

-

1

2

2

2

-

6

-

2

8

2

3

-

2

-

1

1

2

3

-

9

-

2

4

2

4

-

5

-

7

2

4

-

1

2

-

2

2

5

-

8

-

3

2

6

-

3

-

2

12

8

4

-4

-8

收益指数

图2上证收盘指数收益率时序图

2

-

1

-

4

2

-

8

-

1

6

2

1

-

3

-

3

2

1

-

1

1

-

1

2

2

2

-

6

-

2

8

2

3

-

2

-

1

1

2

3

-

9

-

2

4

2

4

-

5

-

7

2

4

-

1

2

-

2

2

5

-

8

-

3

2

6

-

3

-

2

100

80

60

40

20

残差平方

图3残差平方项时序图

常数R t-1 R t-2B

0.0006

0.0449

-0.0159

SEB

0.0336

0.0253

0.2529

BETA

0.0450

-0.0160

T

0.1791

1.7722

-0.6323

SIG T

0.8578

0.0765

0.5272

表1A R（2）模型的SPSS分析

迭代次数1

2

21

22

43

46

64

134Q

24909.4296

24891.2363

24862.8046

24826.9238

24743.4746

24703.6367

5032.6601

2.1950

d

0.4456

0.9408

0.3814

0.6896

0.7408

0.0833

0.6653

0.6162

ω

0.1301

0.8078

0.9900

0.2512

0.0894

0.3947

0.0249

0.4225

β

0.3707

0.3689

0.8301

0.8101

0.8660

0.8763

0.9922

0.9928

准

0.3849

0.4939

0.8259

0.8421

0.8632

0.6924

0.9553

0.8802

似然函数值

-1902.1263

-2252.7902

-2392.2110

-2456.2311

-2603.5616

-2088.5766

-2127.1269

-2373.9469表2基于遗传算法的C语言程序求解

说明在建立了AR模型和FIGARCH 模型后，所得到的误差项｛β

t

｝是白噪声，不存在短记忆性和长记忆性，所建立的模型是正确的，此时：

差分阶数：d=0.6162

参数：ω=0.4225

β=0.9928准=0.8802

似然函数值为：-2373.9469

由此可得，上证综合指数收益率（2000-01-04～2006-06-30）符合AR（2）-FIGARCH（1，0.6162，1）

R t=0.0060＋0.0449×R t-1－0.0159×R t-2＋εt

εt=z tσt E t-1（z t）＝0E t-1（z2t）＝0

σ2t＝0.4225＋0.9928σ2t－1＋［1－0.9928L－

（1－0.8802L）（1－L）0.6162］ε2

t

＋βt

其中：β

t

是白噪声。

四、结束语

股票价格指数的变动反映出了股票市场所在国的政治、经济、社会和其他状况的变化，人们常常称其为“晴雨表”，因此研究股票市场的波动性是很有必要的。FIGARCH模型擅长于反映金融资产的异方差特性以及长记忆型的波动特性，它的主要应用领域是金融资产，包括证券、期权、利率等多方面。从提出至今，它已被许多人成功地应用到证券市场及汇率市场，很好地反映了金融市场的这种波动性。

本文应用遗传算法的思想进行编程，建立FIGARCH模型，模拟了中国股票市场的波动性过程。结果表明，对收益率进行AR建模后，再对FIGARCH模型进行一阶预测的值与真实值的差得到的时间

序列是白噪声，也就是说经过建立自回归

模型和FIGARCH模型后，金融时间序列

已经消除了短记忆性和长期记忆性。上海

股市d=0.6162，显著不同于0和1，说明

过去的冲击对未来股市的影响将会持续

相当长的时间，即中国股票市场的波动性

过程具有长期记忆性。这也就意味着可以

用过去的历史收益和波动情况来预测未

来的收益情况，从而能利用过去的波动性

建立风险控制模型和增加获取投机利润

的机会。除此之外，模型中的β值为：

0.9928，非常接近于1，这意味着条件方差

收敛于无条件方差的速度较慢，因此对条

件方差的冲击经过相当长一段时间才会

消失，也就是说波动性是持久的，进一步

反映了股市的长期记忆性。

遗传算法不同于传统的优化和搜索

方法，它具有智能性和并行性。智能性使

得所选择的子代具有很强的适应性，通过

交叉和变异所得到的后代更适应环境；而

并行性则实现了空间中的多个区域的同

时搜索，保证了大规模计算在短时间内完

成。遗传算法为建立FIGARCH模型提供

了一个平台，我们可以通过优化目标函数

的方法来获得差分阶数d的值，这种做法

是科学的，所计算出的d值是显著异于0

和1的。然而，我们事先无法预测遗传算

法中的最大迭代次数（GEN），只能通过大

量的实验或是预先给出一个比较大的值

来进行定性的检验。遗传算法是通过交叉

和变异来实现结果

优化，当迭代到一定

次数以后，收敛的速

度就会减慢，增加了

运算的时间。除此之

外，对于不同的适应

函数、不同的数据，

应如何确定遗传算

法中的参数，即如何

确定P_CROSSOVER

和P_MUTATION的

值，也有待于进一步

的研究。

（作者单位：中国

人民大学经济学院）

主要参考文献：

[1]R F.A Engle utoregressive conditional

heteroscedasticity with estimates of the

variance of U.K.Inflation[J].Econometrica，

1982.

[2]Bollerslev T.Generalized autoregressive

conditional heteroskedasticity[J].Journal of

Econometrecs，1986.

[3]Engle R F.Bollerslev T，Modeling the

persistence of conditional variances[J].E-

conometric Reviews，1986.

[4]Richard T.Baillie，Tim Bollerslev，Hans

Ole Mikkelsen.Fractionally integrated gen-

eralized autoregressive conditional het-

eroskedasticity[J].Journal of Ecometrics.

1996.Vol.74.

[5]Bollerslev T，Mikkelsen H O.Modeling

and pricing long memory in stock market

volatility[J].Journal of Econometrics，1996.

[6]李汉东，张世英.自回归条件异方差的

持续性研究[J].预测，2000.1.

[7]汤果，何晓群，顾岚.FIGARCH模型对

股市收益长记忆性的实证分析[J].统计研

究，1999.7.

[8]李颖，汤果，陈方正.FIGARCH模型的

参数检验与估计[J].统计与决策，2003.1.

[9]Geweke J.Porter—Hudak.Estimation

and Appli—cation of Long Memory Time

Series Model[J].Journal of Time Series

Analysis，1983.4.

[10]Andrew W.Lo（1991），Long-Term

Memory in Stock Market Prices[J].Econo-

metrica，59（5），September.

[11]王春峰，张庆翠.中国股市波动性过程

中的长期记忆性实证研究[J].系统工程，

2004.1.

[12]王小平，曹立明.遗传算法——

—理论、

应用于软件实现[M].西安：西安交通大学

出版社，2002.

[13]张晓峒.计量经济学软件Eviews使用

指南[M].天津：南开大学出版社，2003.