中国资本市场创新与经济发展研讨会会议.docVIP

中国资本市场创新与经济发展研讨会会议 马科维兹资产组合选择模型的计算方法* 张忠桢 （武汉理工大学管理学院，武汉430070） 引言 马科维兹的资产组合选择模型及其分散风险的思想是现代投资理论的奠基石，在经济发达国家和地区的金融业特别是股市应用广泛，它被用于定量地确定有效投资组合，有助于人们形成合理的投资理念，稳定金融市场。同个股走势分析一样，利用马科维兹模型计算投资组合需要历史数据。但后者考虑资产之间的相关性，前者不考虑。不同资产之间的相关关系是一种结构关系。从系统的观点看，系统的结构往往比个体的行为更为稳定。因此，计算出的投资组合应比个股走势分析有更可靠的参考价值。 马科维兹模型又称为均值-方差资产组合选择模型，是一个凸二次规划，其目标函数是投资组合收益率的方差，是投资风险的度量，但约束条件有多种形式。表示资产或资产比例的变量可以有上界限制也可能没有。变量下界可以是零、正数或负数（称为允许卖空）。也可以有“行业约束”之类的一般不等式，或引入融资成本和交易成本。由于篇幅所限，本文仅讨论马科维兹称为标准约束的资产组合选择模型的解法。熟悉了这一方法，不难用其求解一般形式的均值-方差模型。 自从1952年马科维兹提出他的标准约束模型以来，马科维兹本人以及其他许多学者不断研究这一模型的解法，提出了许多算法。例如，在1987年出版的一本书以及2000年的再版本[1]中，马科维兹用大量篇幅介绍了一种计算方法，称为临界线算法。但这种方法非常复杂且难以及时构造大规模投资组合。为了减少计算量，有的学者利用单因素或多因素套利定价模型构造一个稀疏的协方差矩阵进行计算，有的则干脆修改风险的定义，采用线性规划模型。但这些方法得到的只是近似解。 本文介绍如何用解线性不等式组的一种旋转算法求解与马科维兹模型等价的库恩-塔克条件。笔者的这种方法不仅较为简单而且可以快速计算出马科维兹意义下的有效投资组合。我们利用上海证券交易所1998年8月31日至2000年1月10日之间410支股票70期周末收盘价作计算机实验。实验结果表明，这种方法需要的计算机时间少于用单纯形算法解线性规划模型，甚至少与用高斯消元法解同等规模的线性方程组。对这一方法感兴趣的读者可以用笔者所编的程序和实际股市数据进行计算，以验证其有效性。 线性不等式组的旋转算法 任何线性不等式组可以表示为以下形式： ai x = bi, i=1,…,l; ai x≧bi,, i =l+1,…,m. （1） 其中x是n维未知列向量，ai是n维行向量，bi是标量，i=1,…,m。旋转算法用到以下一些基本概念。a1,…，am中一个最大线性无关组称为不等式组(1)的基，基中的量称为基向量，其余向量称为非基向量。基向量相应的（不）等式称为基（不）等式，否则称为非基（不）等式。由基等式和基不等式构成的不等式组称为基本不等式组。基本不等式组相应的方程组（将“≧”改为“=”） 下面推导不等式组解法中的旋转运算公式。给定不等式组（）x1是基本解，即ajx1=bj, j(I0∪I1，并设 (2) 如果对于某个r (I3∪I2和s(I1有wrs≠0, 则从（2）的第r式可解得： (3) 代入（3）其余各式得： (4) 现在，{r}∪I0∪I1\{s}是新基的编号集，相应的基本解记为x2。在(4)式两边右乘x2， 在(2)式两边右乘x1，两式相减，经整理可得 σi’=σi- (wis /wrs)σr， i∈I3∪I2\{r} 其中σi=aix1-bi，σr=arx1-br，σi’=aix2-bi。同样由式（3）可得 σs’= -σr /wrs, 其中σs’=asx2-bs。 以上运算过程称为旋转运算，详细推导过程见[2]。wrs称为枢轴（元素），其所在行称为枢轴行，所在列称为枢轴列.我们说ar入基，as出基，两者位置的交换记为ar?as。旋转运算前后的数据可用以下两张表表示。 表1 初始表 表2 旋转运算的结果 as aj ar aj ar ai wrs wrj wis wij σr σi as ai 1/wrs -wrj /wrs wis /wrs wij’ -σr /wrs σi’ 其中: wij’= wij-(wis /wrs)wrj， iI3∪I2\{r}, j∈I0∪I1\{s}; σi’=σi