控制与决策  2019, Vol. 34 Issue (3): 581-590  
0

引用本文 [复制中英文]

常浩, 王春峰, 房振明. 随机利率与随机波动率环境下的DC型养老金计划[J]. 控制与决策, 2019, 34(3): 581-590.
[复制中文]
CHANG Hao, WANG Chun-feng, FANG Zhen-ming. Defined-contribution pension plan with stochastic interest rate and stochastic volatility[J]. Control and Decision, 2019, 34(3): 581-590. DOI: 10.13195/j.kzyjc.2017.1133.
[复制英文]

基金项目

国家自然科学基金面上项目(71671122);中国博士后科学基金项目(2014M560185, 2016T90203);教育部人文社会科学研究基金规划项目(16YJA790004);天津市自然科学基金项目(15JCQNJC04000);天津市高校“中青年骨干创新人才培养计划”项目

作者简介

常浩(1979—), 男, 副教授, 博士后, 从事金融工程与风险管理等研究;
王春峰(1966—), 男, 教授, 博士生导师, 从事金融工程与风险管理等研究。

通讯作者

常浩, E-mail: ch8683897@126.com

文章历史

收稿日期:2017-08-25
修回日期:2018-01-13
随机利率与随机波动率环境下的DC型养老金计划
常浩1,2, 王春峰1, 房振明1    
1. 天津大学 管理与经济学部,天津 300072;
2. 天津工业大学 理学院,天津 300387
摘要:为实现养老金的保值增值, 基金管理人将养老金投资于金融市场.假设金融市场包含一种无风险资产、一种股票和一种零息票债券, 其中, 利率期限结构满足随机仿射利率模型, 而股票价格波动率满足Heston随机波动率模型.基金管理人希望寻找一种最优投资组合以最大化其终端财富的期望效用.假设基金管理人对风险的偏好满足幂效用或指数效用, 运用随机动态规划原理和变量替换方法, 得到幂效用和指数效用下最优投资策略的显式解.最后, 通过数值算例分析主要模型参数对最优投资策略的影响.研究结果表明, 利率风险、股市波动风险以及缴费率都对缴费确定(DC)型养老金的投资决策产生较大的影响.
关键词仿射利率    Heston模型    DC型养老金    幂效用    指数效用    最优投资组合    
Defined-contribution pension plan with stochastic interest rate and stochastic volatility
CHANG Hao1,2, WANG Chun-feng1, FANG Zhen-ming1    
1. College of Management and Economics, Tianjin University, Tianjin~300072, China;
2. School of Science, Tianjin Polytechnic University, Tianjin~300387, China
Abstract: In order to preserve and increase the value of pension fund, a fund manager can invest pension fund into a financial market with a risk-free asset, a stock and a zero-coupon bond. In the financial market, short rate is assumed to follow stochastic affine interest rate model, while stock price is supposed to be driven by Heston's stochastic volatility model. The fund manager expects to obtain an optimal portfolio to maximize the expected utility of terminal utility. Assuming the risk-aversion of the fund manager to satisfy power utility or exponential utility, the explicit solutions to the optimal portfolios with power utility or exponential utility are obtained by using the principle of stochastic dynamic programming and the variable change technique. Finally, a numerical example is given to analyse the impact of main market parameters on the optimal portfolios. Research results show that interest rate risk, volatility risk and contribution rate have a considerable influence on decision-making of fund manager for the defined-contribution pension found.
Keywords: affine interest rate    Heston model    defined-contribution pension fund    power utility    exponential utility    optimal portfolios    
0 引言

养老金是一项重大的公共政策问题.在大多数国家, 养老金占据了政府财政负债的很大比重, 且随着“婴儿潮”一代人的逐渐退休, 这一比例还将继续增大.同时, 由于股票市场的回报低和不断降低的利率以及人们寿命的延长, 许多机构的养老金计划都面临着严重的资金困难.常见的养老金计划有两种:收益确定(DB)型和缴费确定(DC)型. DB型计划中, 养老金支付是确定的, 而且不管养老金如何提供支持, 都要承诺基于确定的收益支付养老金, 这类计划不考虑收益水平或后期发生的通货膨胀. DC型计划中, 缴费的比率是固定的, 养老金支付水平取决于至退休前的期间内累积的缴费值以及投资所得.由此可见, DC型养老金计划比较符合现行的养老金管理机制.

近年来, 许多学者对不同投资环境下的DC型养老金计划进行研究, 取得了很多研究成果.谷爱玲等[1]假设股票瞬时收益率服从Ornstein-Uhlenbeck过程, 运用随机最优控制理论得到了幂效用下最优投资策略的显式解, 并分析了瞬时收益率对最优投资策略的影响.考虑到养老金的投资期限比较长, 一般可以达到20 ~ 40年, 受很多随机因素的影响, 利率一般不是常数或时间的有界确定函数, 而是一种随机过程, 因此, 一些学者对随机利率下的DC型养老金计划进行了研究.如: Boulier等[2]对随机利率环境下带有最低保障水平的DC型计划进行研究, 运用对偶方法得到了幂效用下最优投资策略的显式解, 并分析了投资期限对最优投资策略和财富水平的影响; Deelstra等[3]将Boulier等[2]的模型拓展到仿射利率模型下, 并假设保障水平和缴费率都是随机的, 运用对偶方法进一步研究了DC型养老金计划问题; Gao[4]在Deelstra等[3]模型基础上进一步对仿射利率模型下的DC型计划进行研究, 运用随机最优控制理论得到了对数效用下最优投资策略的显式解; 张初兵等[5]在Gao市场模型[4]基础上进一步对幂效用和指数效用下的DC型养老金计划进行了研究.

另一方面, 在长达20 ~ 40年的投资决策过程中, 单位货币的购买力一般会随着时间的推移而减少, 即长期投资决策过程中的通胀风险不可忽视, 因而一些学者开始对通胀风险环境下的DC型计划进行了研究.如:殷俊等[6]同时考虑利率风险和通胀风险, 研究了幂效用下DC型计划的最优投资策略; 伍慧玲等[7]考虑通胀风险和随机收入的情形, 研究了均值-方差准则下多周期DC型计划的最优投资策略和风险-收益的数量关系, 研究结果表明, 通胀风险和随机收入将对最优投资策略和有效边界产生本质的影响; 常浩等[8]假设计划参与人对风险的偏好满足HARA效用函数, 该效用函数包含幂效用、指数效用和对数效用, 运用随机最优控制理论和Legendre变换-对偶理论研究了通胀风险环境下的DC型养老金计划.这些研究成果为通胀风险环境下的投资决策提供了理论依据.

金融资产价格的波动风险也是一种常见的金融市场风险. Heston模型[9]是刻画金融资产价格波动风险的最常见的模型之一, 该模型结构简单, 已被证实能够很好地刻画金融资产价格的波动率, 是目前连续时间金融中广泛研究的随机波动率模型.目前, 许多学者对Heston随机波动率模型下的各种投资问题进行研究, 取得了很多研究成果.如: Kraft[10]研究了Heston模型下的投资组合选择问题; Li等[11]研究了Heston模型下的时间一致投资-再保险策略问题; Zhao等[12]研究了Heston模型下的超额损失投资-再保险策略问题; 张初兵等[13]研究了Heston模型下的DC型养老金计划问题.这些研究成果分析了Heston模型下各种投资问题中最优投资策略的动态变化过程, 丰富和发展了投资组合选择理论.

上述模型大多是在单一风险环境下取得的, 而实际的金融市场通常是多种风险并存的, 如当利率发生波动时, 金融资产价格也会产生一定的波动.因此, 人们发现同时考虑利率风险和波动风险比较符合实际的投资过程. Liu[14]、Li等[15]以及Noh等[16]根据利率动态对金融资产价格波动的不同影响, 分别构造了3种不同的随机利率与随机波动率模型, 并研究了这些模型下的一些投资问题.在这3种模型中, Liu[14]以及Noh等[16]给出的模型比较一般, 但往往不能得到最优投资组合的显式解; Li等[15]给出的模型虽能得到显式解, 但没有考虑利率变化对股票价格的影响趋势.随后, Chang等[17]在Li等[15]市场模型基础上进一步研究了投资-消费问题; Guan等[18]在Liu[14]的市场模型基础上进一步考虑了最低保障水平, 研究了DC型养老金计划问题.这些成果为随机利率和随机波动率模型的投资决策提供了理论依据, 但这些成果都是在幂效用下取得的.与幂效用相比较, 指数效用属于常系数绝对风险厌恶效用函数, 其对风险的厌恶程度不依赖于当前财富水平.在效用理论中, 幂效用和指数效用经常被用于刻画不同投资者对风险的厌恶程度, 已广泛地应用到投资组合选择理论中.就问题解决的难度而言, 指数效用下最优投资策略的求解更加困难, 原因在于原HJB方程解的形式以及其产生的偏微分方程解的形式都发生了本质的变化, 猜想这些方程解的形式比较困难.

据笔者对现有文献的查阅情况来看, 尚无其他文献研究随机利率与随机波动率模型下的DC型养老金计划问题, 除了文献[18]以外.但文献[18]考虑了最低保障水平, 模型处理方法比较特殊, 且仅得到了幂效用下的最优投资策略.本文在文献[18]的市场模型基础上, 进一步研究计划参与人退休前养老金投资过程的投资决策问题, 通过假设基金管理人对风险的偏好满足幂效用或指数效用, 运用随机最优控制理论对随机利率和随机波动率环境下的DC型养老金计划问题进行研究, 得到最优投资策略的显式解, 给出算例并分析利率参数和波动率参数对最优投资策略的影响.数值结果表明:利率风险、波动风险和缴费率都会对养老金的投资决策产生较大的影响.与现有文献相比较, 本文有如下几个创新点: 1)~对随机利率和随机波动率模型下的DC型养老金计划问题进行研究; 2)~得到幂效用和指数效用下最优投资策略的显式解; 3)~克服了指数效用下价值函数猜想方面的难点, 为其他更加复杂投资问题的解决提供了理论参考.

1 DC型养老金的连续时间模型

假设[0, T]是固定的有限投资期, 其中T为终端时刻, (Wr(t), Ws(t))是定义在完备概率空间(Ω, F, {Ft}t∈[0, T], P)上相互独立的二维标准布朗运动, {Ft}t∈[0, T]是由(Wr(t), Ws(t))所产生的信息流.假设金融市场是无摩擦的、没有交易成本且是完备的, 存在一种无风险资产、一种风险资产和一种零息债券可连续交易.

第1种资产为无风险资产, 可视为现金或银行账户, 其时刻t的价格记为S0(t), t∈[0, T], 则S0(t)满足微分方程

(1)

其中r(t)表示瞬时利率.本文假设利率r(t)是随机过程, 其期限结构满足随机仿射利率模型[3]

(2)

其中: , 且abk1k2均为正常数, 参数a表示利率的均值回复水平, 参数b表示利率回复到均值的速度, σr表示利率的波动率.

注1  式(2)包含两种常见的利率模型(作为特例):一种是Vasicek利率模型(当k1=0时); 另一种是CIR利率模型(当k2 =0时).当k2 =0时, 式(2)退化为CIR利率模型, 这时需要条件2a > k1来确保r(t) > 0.

第2种资产为风险资产, 可视为股票, 其时刻t的价格记为S1(t), t∈[0, T].考虑到利率波动必然对股票价格产生很大的影响, 本文在原Heston随机波动率模型基础上考虑利率动态对股票价格的影响, 因而可假设S1(t)服从如下随机微分方程[18]:

(3)
(4)

其中: σsλrvαδσl均为正常数且满足条件2αδ > σl2, 该条件可确保对于任意的t∈[0, T], 恒有L(t) > 0.在方程(3)中, σsσr表示由利率波动引起的股票价格波动率, λr σr表示利率风险的市场价格, 表示股市波动风险的市场价格, 表示股票价格的波动率; 在方程(4)中, 参数α表示股票波动率回复到均值的速度, δ表示股票波动率的均值回复水平, 表示股票波动率所引起的波动率.

注2  注意方程(3)中包含两种不同的布朗运动Wr(t)和Ws(t), 它们相互独立, 这表示股票价格的不确定性来源于两部分:一部分是由利率的不确定性所引起的; 另一部分是由股票价格自身的不确定性所引起的.为了能够成功地开展本文工作, 这里假设波动率的不确定性来源于股票的不确定性, 即方程(4)中的布朗运动也用Ws(t)来表示.

第3种资产为到期时刻为T的零息票债券, 其时刻t时的价格记为B(t, T), t∈[0, T], 则B(t, T)服从如下随机微分方程[3]:

(5)

其中

进一步, 解方程(5), 可以得到B(t, T)=eφ0(t)-φ1(t)r, 其中φ0(t)由下式确定:

本文主要考虑退休前养老金投资过程中应采取的最优投资策略问题.为了实现养老金的保值增值, 基金管理者将其财富投资于金融市场.假设投保人的缴费率是常数, 记为c, 时刻t投资在股票和零息票债券的资金数量分别记为πs(t)和πB(t), 则时刻t投资在银行账户的资金数量可记为π0(t), 且π0(t)=X(t)-πs(t)-πB(t), 其中X(t)表示时刻t的财富过程.假设股票市场允许卖空和借贷, 则在投资策略πs(t)和πB(t)下基金管理者的财富过程X(t)可表示为

(6)

其中初始时刻X(0)=x0 > 0.

定义1(可容许策略)  投资策略πs(t)和πB(t)被称为可容许的, 如果πs(t)和πB(t)满足如下3个条件:

1) πs(t)和πB(t)是Ft可测的, 且满足

2)

3) 对于任意的πs(t)和πB(t), 方程(6)的解存在且唯一.

Π ={π (t)|π (t)=(πs(t), πB(t)), 0≤ tT}为养老金计划所有可容许策略形成的集合.基金管理者期望寻找一种最优投资策略π*(t)以最大化其终端财富的期望效用, 即可将目标表示为

(7)

其中U(·)表示效用函数, 是严格凸的连续可微函数且满足Inada条件: U'(+∞)=0, U'(0)=+∞.本文假设基金管理者对风险的厌恶程度满足下列两种效用函数: 1) 幂效用函数, 其U(x)可写为U(x)=xp/p, p < 1且p≠0, p称为风险厌恶因子; 2)指数效用函数, 其U(x)可写为U(x)=-e-qx/q, q>0, q是风险厌恶因子.

注3  与现有文献相比, 本文在市场模型的构建和效用函数的选择方面有如下几个优点:

1) 与文献[18]的工作相比, 本文研究的DC型养老金问题是不同的.文献[18]考虑了最低保障水平, 得到了幂效用下的最优投资策略, 但其结果并不能退化到本文的结果, 且由于模型考虑了最低保障水平, 其模型处理方法比较特殊.

2) 与文献[14-17]的工作相比, 本文研究的投资问题是不同的.文献[14, 16-17]研究的都是投资-消费问题, 文献[15]研究的是投资组合选择问题, 且都假设投资人对风险的偏好满足幂效用函数; 而本文研究的是DC型养老金计划问题, 不仅得到了幂效用下的最优投资策略, 而且得到了指数效用下的最优投资策略, 为基金管理人在风险偏好视角下的最优投资决策提供了理论参考.

3) 与文献[2-5]相比, 本文不仅考虑了利率风险, 而且还考虑了金融资产价格的波动风险; 而文献[2-5]仅考虑了利率风险.

2 最优投资策略

本节运用随机最优控制理论和变量替换方法求解问题(7).

定义价值函数

(8)

根据随机最优控制原理, 如果V(t, r, l, x)∈C1, 2, 2, 2([0, TR+× R+× R+), 则V(t, r, l, x)可视为下列Hamilton-Jacobi-Bellman(HJB)方程的一个连续解:

(9)

其中: V(T, r, l, x)=U(x), 且VtVrVlVrrVllVxVxxVrxVlx分别表示V(t, r, l, x)对变量trlx的一阶和二阶偏导数.

根据极值的一阶条件可得最优策略πs*(t)和πB*(t)分别为

(10)
(11)

将式(10)和(11)代入(9), 得

(12)

其中边界条件为V(T, r, l, x)=U(x).

下面分别在幂效用和指数效用下求解方程(12), 从而得到两种效用下问题(7)的最优投资策略.

2.1 幂效用下的最优投资策略

幂效用下, 假设方程(12)的解具有如下结构:

(13)

计算V(t, r, l, x)关于各个变量的一阶和二阶偏导如下:

将上述各阶偏导数代入方程(12), 可得

消去对x的依赖, 得

(14)
(15)

为了方便阅读, 将上述偏微分方程的求解过程表达成下面3个引理.

引理1  假设方程(14)解的结构为f(t, r, l)=eA1(t)l+A2(t)r+A3(t), 其中边界条件为A1(T)=0和A2(T)=A3(T)=0, 则在

条件下, A1(t)、A2(t)和A3(t)分别为

(16)
(17)
(18)

其中

证明  将f(t, r, l)=eA1(t)l+A2(t)r+A3(t)代入方程(14), 整理后可得

消去对lr的依赖, 得

(19)
(20)
(21)

解方程(19)和(20)可得式(16)和(17).进一步, 通过式(21)× k1-式(20)× k2可得式(18).

引理2  假设方程(15)解的结构为g(t, r)= , 则满足方程

(22)

证明  对于任意函数g(t, r), 定义如下变分算子:

则方程(15)可写成如下形式:

(23)

另一方面, 有

(24)
(25)

将式(24)和(25)代入(23), 得

即得

引理3  假设方程(22)解的结构为eB2(t)r+B1(t), 其中边界条件为B2(T)=0和B1(T)=0, 则B2(t)和B1(t)分别为

(26)
(27)

其中

证明  将eB2(t)r+B1(t)代入方程(22), 整理后得

消除对lr的依赖, 得

(28)
(29)

解方程(28)可得式(26).进一步, 通过式(29)× k1-(28)× k2, 可得式(27).

应用式(13)和其各阶偏导数, 可以得到

综上所述, 考虑到式(10)和(11), 将幂效用函数下问题(7)的最优投资策略总结如下.

定理1  当效用函数为U(x)=xp/p, p < 1且p≠0时, 则在

p≠0的条件下, 问题(7)的最优投资策略为

(30)
(31)

其中: A1(t)和A2(t)由引理1确定; g(t, r)和gr分别由下列两式确定:

注4  由式(30)和(31)可以得出: 1)投资于股票的资金数量πs*(t)依赖于参数abk1k2λrvασlpTc, 而与σsδ无关; 2)投资于零息票债券的资金数量πB*(t)依赖于参数abk1k2σSλrvασlpTc, 而与δ无关; 3)经比较发现, B2(t)=-φ1(t), B1(t)=φ0(t), 即B(t, T), 因此, 幂效用下πs*(t)和πB*(t)依赖于财富过程X(t)和零息票债券价格过程B(t, T).

注5  考虑到文献[18]研究的是带有最低收益保障的DC型养老金计划问题, 也得到了幂效用偏好下的最优投资策略.与文献[18]的结论相比, 本文结论有很大的不同: 1)本文的πs*(t)依赖于利率动态, 从而更进一步与零息票债券价格密切相关, 而文献[18]的πs*(t)不依赖于利率动态; 2)本文的πB*(t)直接依赖于财富过程和零息票债券价格过程B(t, T), 其结构和形式都发生了根本的变化; 3)本文的价值函数的结构和形式发生了根本的变化.

如果不考虑缴费率, 即c=0, 则本文将退化为随机利率与随机波动率模型下的动态资产分配问题.此时有g(t, r)=0和gr =0, 从而定理1的结论退化为

由此可见, 如果不考虑缴费率c, 则幂效用下的最优投资策略仅依赖于财富过程X(t), 这与定理1的结论有着本质的区别.

2.2 指数效用下的最优投资策略

指数效用下, 本文猜想方程(12)的解的形式如下:

(32)

其中边界条件为: J(T, r)=1, K(T, r)=0和ϕ(T, r, l)=0.

计算式(32)的各阶偏导数如下:

将上述偏导数代入方程(12), 分离变量后可得如下3个偏微分方程:

(33)
(34)
(35)

假设方程(33)的解具有结构J(t, r)=eD1(t)+D2(t)r, 其边界条件为D1(T)=0和D2(T)=0.将J(t, r)代入方程(33), 分离变量后可得

(36)
(37)

解方程(36)和(37)可得

(38)
(39)

其中

观察方程(34)和(15), 可以发现二者同解, 因此有K(t, r)=g(t, r).

观察方程(35), 可将其解假设为

(40)

将式(40)代入方程(35), 分离变量后可得

(41)
(42)
(43)

求解方程(41) ~ (43)可得

(44)
(45)
(46)

应用式(32)和其各阶偏导数, 可得

综上所述, 可将指数效用下问题(7)的最优投资策略总结如下.

定理2  假设效用函数为U(x)=, q > 0, 则问题(7)的最优投资策略可表示为

(47)
(48)

其中: D2(t)、D4(t)和D5(t)分别由式(38)、(45)和(44)所确定; K(t, r)和Kr分别由下列两式确定:

注6  由式(47)和(48)可以看出: 1)投资于股票的资金数量πs*(t)依赖于参数abk1k2λrvασlqT, 而与cσsδ无关; 2)投资于零息票债券的资金数量πB*(t)依赖于参数abk1k2σSλrvασlqTc, 而与δ无关; 3)经观察发现, D2(t)=φ1(t), D1(t)=-φ0(t), 即J(t, r)=B-1(t, T), 因此, 指数效用下πs*(t)仅依赖于零息票债券价格过程B(t, T), 而πB*(t)不仅依赖于B(t, T), 而且依赖于财富过程X(t).

注7  由定理1和定理2可以看出, 指数效用与幂效用下最优投资策略之间的区别: 1)幂效用下的πs*(t)依赖于当前财富水平X(t), 而指数效用下的πs*(t)不依赖于X(t), 这与两种效用函数所表示的对风险的厌恶程度是一致的, 幂效用下投资者对风险的绝对厌恶程度随着瞬时财富水平的增加而递减, 而指数效用下投资者对风险的绝对厌恶程度是一个常数; 2)与幂效用相比, 指数效用下的πB*(t)其结构和投资数量都发生了很大的变化; 3)与幂效用下最优投资策略的求解相比, 指数效用下价值函数的结构发生了很大的变化, 这给准确猜想价值函数的结构带来了很大困难.

如果不考虑缴费率c, 则有K(t, r)=0和Kr =0.因而, 定理2的结论将退化为

(49)
(50)

与定理2的结论相比, πs*(t)是相同的, 而πB*(t)有着本质的区别.

综合定理1和定理2可以发现, 无论管理者对风险的偏好是满足幂效用还是指数效用, 缴费率c都对DC型养老金计划的最优资产配置有着本质的影响.

3 算例分析

本节将给出数值算例以分析一些主要模型参数对最优投资策略的影响.假定参数取值为: a=0.018 712, b=0.233 9, k1 =0.007 293 16, k2 =0, r(0)=0.05, σs =0.84, λr =0.46, v=0.68, α =0.2, δ =0.6, σl =0.6, l0 =0.2, t=0, T=40, X(0)=100, p=-0.5, q=0.02, c=0.15.

上述取值中, abk1k2的取值来源于文献[3]中对股市利率期限结构的实证结果.幂效用下, 本文运用定理1的结论分析参数bk1αT的变化对最优策略的影响, 具体见图 1 ~ 图 4.指数效用下, 运用定理2的结论分析参数σlσsvc对最优策略的影响, 具体见图 5 ~ 图 8.

图 1 幂效用下b对最优投资策略的影响
图 2 幂效用下k1对最优投资策略的影响
图 3 幂效用下α对最优投资策略的影响
图 4 幂效用下T对最优投资策略的影响
图 5 指数效用下σl对最优投资策略的影响
图 6 指数效用下σs对最优投资策略的影响
图 7 指数效用下v对最优投资策略的影响
图 8 指数效用下c对最优投资策略的影响

图 1 ~ 图 4描绘了幂效用下最优投资策略随参数变化的动态变化过程, 从中可以得到如下一些重要结论:

1) πs*(t)和π0*(t)随参数b递增, 而πB*(t)随参数b递减.从参数b所代表的经济含义来看, 当参数b的值递增时, 利率r(t)的期望值随之递减, 这导致利率风险的市场价格随之递减, 即利率风险的风险补偿也随之递减, 这时, 零息票债券的投资收益随之递减, 这意味着基金管理人应减少投资在零息票债券的资金额, 而增加投资在无风险资产上的资金额.同时, 由于利率的期望值递减导致股票的波动率降低, 基金管理人应逐步增加投资在股票上的资金额.

2) πs*(t)和π0*(t)随参数k1递减, 而πB*(t)随参数k1递增.由方程(2)可以看出, 参数k1决定着利率动态的波动率.当参数k1的值递增时, 利率的波动率随之递增, 同时导致股票的波动率也随之递增, 这意味着基金管理人应减少投资在股票和零息票债券上的资金数量.为了对冲利率风险, 基金管理人应逐步增加投资在零息票债券上的资金数量.

3) πs*(t)随参数α递增, 而πB*(t)和π0*(t)随参数α递减.由方程(4)可以看出, 当参数α的值递增时, 波动率L(t)的期望值随之递减, 这导致股票波动率随之降低, 此时, 基金管理人更愿意将更多资金投资于股票以获取更多的收益, 同时, 基金管理人应减少投资于零息票债券和无风险资产的资金额.

4) πs*(t)和π0*(t)随参数T递增, 而πB*(t)随参数T递减.模型中, 参数T代表投资期限, T越小表示投保人越临近退休.因此, 投保人在刚开始缴费时(T的值比较大), 基金管理人应将更多资金投资于股票和无风险资产; 而随着年龄的增加, 投保人越临近退休(T的值比较小), 基金管理人投资于股票和无风险资产的资金将逐步减少, 同时, 应增加投资于零息票债券的资金数量.

图 5 ~ 图 8给出了指数效用下最优投资策略随主要市场参数变化的动态变化过程, 从中可以得出如下重要结论:

1) πs*(t)和πB*(t)随参数σl递增, 而π0*(t)随参数σl递减.从本文的结论看, 当参数σl的波动率增大时, 基金管理人应增加投资于股票和零息票债券上的资金数量, 而减少投资于无风险资产上的资金数量.

2) πB*(t)随参数σs递增, π0*(t)随参数σs递减, 而πs*(t)保持不变.由方程(3)可以看出, 当参数σs的值增大时, 由利率波动引起的利率风险将增大, 为了对冲利率风险, 基金管理人应增加投资于零息票债券的资金数量, 而减少投资于无风险的资金数量, 投资于股票的数量保持不变.

3) πs*(t)和πB*(t)随参数v递增, 而π0*(t)随参数v递减.由方程(3)可以看出, 当参数v的值递增时, 股市波动风险的市场价格随之递增, 这导致股票收益率也随之递增.为了获取更多收益, 基金管理人更愿意将更多资金投资于股票和零息票债券, 同时减少投资于无风险资产的资金数量.

4) πB*(t)随参数c递减, 而π0*(t)随参数c递增, πs*(t)则保持不变.从本文的结论看, 当缴费率c较低时, 基金管理人应将大量的资金投资到股票和零息票债券上以获取更多的收益.随着缴费率的递增, 养老金变得比较富裕, 这时基金管理人变得比较保守, 投资于无风险资产的资产数量将递增, 而购买零息票债券的资金数量将递减, 投资于股票的资金数量保持不变.

4 结论和展望

本文假设利率和股票价格波动率都是随机过程, 在文献[18]的市场模型基础上, 进一步对DC型养老金计划的最优投资策略进行研究, 得到了幂效用和指数效用下最优投资策略的显式解.通过数值算例分析了利率参数和波动率参数对最优投资策略的影响, 为基金管理人在风险偏好视角下的最优投资决策提供了理论依据.研究结果表明:利率风险、波动风险和缴费率都对DC型养老金计划的投资决策产生了较大的影响, 且不同风险偏好的投资人在确定最优投资策略时对模型参数、财富过程、利率风险和波动率风险的敏感性也不同.

养老金计划问题正在变得越来越复杂.首先, 随机收入、保费返还条款、死亡风险和长寿风险也是影响养老金投资的最重要因素, 这些因素的引入会给养老金投资策略的解决带来很大困难; 其次, 养老金投资过程的风险很难准确度量, 其投资过程中应采取的风险准则很难确定, 有待大量的理论和实证研究; 最后, 除均值-方差准则外, 其他风险准则下的投资组合选择理论尚未取得突破性进展, 如VaR风险准则或均值-CVaR风险准则等, 这也给多种风险环境下基于不同风险测度的DC型养老金计划的研究带来了很大困难.这些问题将是本文未来的研究内容.

致谢: 衷心感谢审稿专家提出的富有启发性的修改意见!
参考文献
[1]
谷爱玲, 李仲飞, 曾燕. Ornstein-Uhlenbeck模型下DC养老金计划的最优投资策略[J]. 应用数学学报, 2013, 36(3): 715-726.
(Gu A L, Li Z F, Zeng Y. Optimal investment strategy under Ornstein-Uhlenbeck model for a DC pension plan[J]. Acta Mathematicae Applicatae Sinica, 2013, 36(3): 715-726.)
[2]
Boulier J F, Huang S J, Taillard G. Optimal management under stochastic interest rates: The case of a protected defined contribution pension fund[J]. Insurance: Mathematics and Economics, 2001, 28(1): 173-189.
[3]
Deelstra G, Grasselli M, Koehl P F. Optimal investment strategies in the presence of a minimum guarantee[J]. Insurance: Mathematics and Economics, 2003, 33(1): 189-207.
[4]
Gao J W. Stochastic optimal control of DC pension funds[J]. Insurance: Mathematics and Economics, 2008, 42(3): 1159-1164.
[5]
张初兵, 荣喜民. 仿射利率模型下确定缴费型养老金的最优投资[J]. 系统工程理论与实践, 2012, 32(5): 1048-1056.
(Zhang C B, Rong X M. Optimal investment for DC pension under the affine interest rate model[J]. Systems Engineering - Theory & Practice, 2012, 32(5): 1048-1056. DOI:10.3969/j.issn.1000-6788.2012.05.017)
[6]
殷俊, 李媛媛. 基于随机利率和通货膨胀的缴费确定型养老金计划最优资产配置策略[J]. 当代经济科学, 2013, 35(2): 11-20.
(Yin J, Li Y Y. Stochastic interest rate and inflation based optimal asset allocation tactics for defined-contribution pension plans[J]. Modern Economic Science, 2013, 35(2): 11-20. DOI:10.3969/j.issn.1002-2848.2013.02.002)
[7]
伍慧玲, 董洪斌. 带有通胀风险和随机收入的确定缴费养老计划[J]. 系统工程理论与实践, 2016, 36(3): 545-558.
(Wu H L, Dong H B. Multi-period mean-variance defined contribution pension management with inflation and stochastic income[J]. Systems Engineering - Theory & Practice, 2016, 36(3): 545-558.)
[8]
常浩, 王春峰, 房振明. 通胀风险下基于HARA效用的DC型养老金计划[J]. 运筹学学报, 2016, 20(4): 39-51.
(Chang H, Wang C F, Fang Z M. Defined contribution pension fund scheme with HARA preference under inflation risk[J]. Operations Research Transactions, 2016, 20(4): 39-51.)
[9]
Heston S L. A close-form solution for options with stochastic volatility with applications to bond and currency options[J]. Review of Financial Studies, 1993, 6(2): 327-343. DOI:10.1093/rfs/6.2.327
[10]
Kraft H. Optimal portfolios and Heston's stochastic volatility model: An explicit solution for power utility[J]. Quantitative Finance, 2005, 5(3): 303-313. DOI:10.1080/14697680500149503
[11]
Li Z F, Zeng Y, Lai Y Z. Optimal time-consistent investment and reinsurance strategies for insurers under Heston's SV model[J]. Insurance: Mathematics and Economics, 2012, 50(1): 191-203.
[12]
Zhao H, Rong X M, Zhao Y G. Optimal excess-of-loss reinsurance and investment problem for an insurer with jump-diffusion risk process under the Heston model[J]. Insurance: Mathematics and Economics, 2013, 53(3): 504-514. DOI:10.1016/j.insmatheco.2013.08.004
[13]
张初兵, 荣喜民, 侯如靖, 等. Heston模型下确定缴费型养老金的投资组合优化[J]. 系统工程, 2012, 30(12): 39-44.
(Zhang C B, Rong X M, Hou R J, et al. Investment portfolio optimization for defined-contribution pension under a Heston model[J]. Systems Engineering, 2012, 30(12): 39-44.)
[14]
Liu J. Portfolio selection in stochastic environments[J]. The Review of Financial Studies, 2007, 20(1): 1-39. DOI:10.1093/rfs/hhl001
[15]
Li J Z, Wu R. Optimal investment problem with stochastic interest rate and stochastic volatility: Maximizing a power utility[J]. Applied Stochastic Models in Business and Industry, 2009, 25(3): 407-420. DOI:10.1002/asmb.v25:3
[16]
Noh E J, Kim J H. An optimal portfolio model with stochastic volatility and stochastic interest rate[J]. J of Mathematical Analysis and Applications, 2011, 375(2): 510-522. DOI:10.1016/j.jmaa.2010.09.055
[17]
Chang H, Rong X M. An investment and consumption problem with CIR interest rate and stochastic volatility[J]. Abstract and Applied Analysis, 2013(1): 1-12.
[18]
Guan G H, Liang Z X. Optimal management of DC pension plan in a stochastic interest rates and stochastic volatility framework[J]. Insurance: Mathematics and Economics, 2014, 57(1): 58-66.