购买十只股票的最优投资组合

1作业二：购买十只股票的最优投资组合1、理论基础马科维茨讨论了投资者将一笔资金在给定的持有期进行投资的问题，也就是选择一个最优的证券组合。由于每种证券（从而证券组合）未来的收益率是未知的，因此，不可能做出一个保证获得最高收益的决策。尽管可以估计每种证券未来的收益率（期望收益率），仍然不能满足上面的要求，这是因为，基于期望收益率的决策最多只能获得最高平均收益率（组合的期望收益率）。正是因为对收益率的不确定性（风险）在决策中的关注，马柯维茨指出，任何一位投资者在追求“高收益”的同时，还希望“收益尽可能是确定的”。决策目标应该有两个：第一，“尽可能高的收益率”；第二，“尽可能低的不确定性（风险）”。1.1收益与风险的度量有关风险和收益的度量，本文用期望度量收益，用方差(或标准差)度量风险。具体的用历史数据估计期望收益率和方差——样本均值和样本方差。假设收益率的概率分布是恒定的，那么，实际收益率就是来自同一概率分布的抽样样本。因而，可以用样本均值和样本方差对期望收益率和方差进行估计。假设从时刻1t到nt的实际收益率是),...,2,1(niri，这就是由收益率的时间序列所构成的样本，则样本均值和样本方差为：niirnr11niirrn122)(111.2证券之间的关联性--协方差与相关系数用BArr,分别表示证券A和证券B的收益率，则其联合分布通常表示为：ijBjBAiAPrrrrP),(证券A和证券B的协方差由下式计算：minjijBBjAAiBApErrErrrr11))((),cov(2协方差反映两种证券协同变化的数量，数值大小依赖于证券收益率与自身期望收益率的偏离程度。然而，协方差的数值大小并不能完全反映证券间的关联关系。为了对相关程度做出衡量，应将上面的偏离程度进行标准化，标准化后的协方差就是相关系数。数学公式如下：BABAABrr),cov(，其中，BA,分别是证券A和证券B的标准差。1.3证券投资组合的期望方差计算证券投资组合的方差与期望的确定要计算组合P的期望收益率和方差，除了要知道证券A和B的期望收益率和方差外，还必须知道证券A和B之间的相关系数AB（或协方差）。计算公式如下：)()()(BBAAPrEWrEWrEBAABBABBAAP222222显然,选择不同的组合权数AW，就可以得到不同的证券组合，从而得到不同的期望收益率和方差。投资者可以根据自己的偏好，选择符合自己要求的证券组合。1.3马柯维茨均值方差模型要使用数学模型，就需要对客观现实进行假设。虽然社会科学对条件假设不像数学那样严格，必要的说明还是不可少的，马柯维茨模型做了下面的假设。假设1：投资收益率的概率分布是已知的；假设2：投资风险用投资收益率的方差或标准方差标识；假设3：影响投资决策的主要因素为期望收益率和风险；假设4：投资者都遵守占优原则；同一风险水平下，选择收益率较高的证券；同一收益率水平下，选择风险较低的证券。理性的投资者总希望风险尽量小而收益尽量大，或者收益相等风险最小。收益相同风险最小，即在给定的期望收益率下，寻找最小方差的资产组合;风险相同收益最大，即在给定的方差下，寻找最大期望收益的资产组合.以上两种组合中所有有效组合的集合就称作有效前沿。3求最优投资组合就是在一定的预期收益水平下，求方差最小的组合。可表示为如下最优问题：V21min}{s.t.1),~(IWrERWP其中，))~(),...,~(),~((21NrErErER，表示N个风险资产收益率的期望向量，是已知的。I表示每个分量都等于1的N维向量，V表示N个资产的协方差矩阵，W表示权重，前沿证券组合所对应收益率)~(PrE，体现在约束条件中。2.实证研究2.1数据选择本文选取中国A股中十支股票，分别为：中航动控、沧州明珠、海信电器、太原刚玉、五矿稀土、大康牧业、新华龙、上海物贸、大同煤业、南京熊猫。运用matalb计算一天之中收益率的均值、收益率的方差，进一步用matlab计算十只股票的协方差矩阵、相关系数矩阵。2.1.1十只股票的收益率、收益率标准差在EXCEL中计算原始数据的收益率，计算公式为：)ln(1ttpp中。将得到的数据导入matlab中求出每只股票收益率的均值和标准差分别代表该只股票的收益率和收益率标准差。股票收益率(%)收益率标准差(%)中航动控(zh)0.4780.265沧州明珠(cz)0.03720.245海信电器(hx)0.03760.2649太原刚玉(ty)0.03660.3314五矿稀土(wk)0.05440.2742新华龙(xh)0.053160.3029中国建筑(zg)0.010170.26303中信国安(zx)0.051490.27055南京熊猫(nj)0.014640.17907大同煤业(dt)0.049280.491342.1.2十只股票的协方差矩阵在matlab中用原始数据中的十只股票收益率计算十只股票的协方差矩阵，所得到的结果如下表所示：2.1.3十只股票的相关系数矩阵在matlab中运用原始数据的收益率计算十只股票的相关系数矩阵，计算结果如下表所示：2.2最优组合的确定通过以上十只股票的相关数据，分别通过求解最优化问题:V21min}{zhczhxtywkxhzgzxnjdtzh0.0728-0.0004-0.0011-0.00110.0020.00550.0030.00210.0062-0.0008cz-0.00040.05750.01630.00760.0060.00420.00170.0090.00020.0034hx-0.00110.01630.0710.0040.01410.00410.00280.01150.00050.0089ty-0.00110.00760.0040.11220.02290.01040.00130.00920.00320.0038wk0.0020.0060.01410.02290.07660.00530.00720.02230.00290.0161xh0.00550.00420.00410.01040.00530.09260.00510.0091-0.0009-0.0211zg0.0030.00170.00280.00130.00720.00510.06960.00090.00190.0129zx0.00210.0090.01150.00920.02230.00910.00090.07420.00280.0019nj0.00620.00020.00050.00320.0029-0.00090.00190.00280.03210.0013dt-0.00080.00340.00890.00380.0161-0.02110.01290.00190.00130.2425zhczhxtywkxhzgzxnjdtzh1-0.0054-0.0157-0.0120.0270.06760.04220.02880.1275-0.006cz-0.005410.25490.09470.09050.0570.02660.13810.00370.0288hx-0.01570.254910.04520.19070.050.03970.15780.00950.0681ty-0.0120.09470.045210.24740.10240.01470.10120.05420.0229wk0.0270.09050.19070.247410.06320.09840.29580.0590.1181xh0.06760.0570.050.10240.063210.0630.1092-0.0162-0.1407zg0.04220.02660.03970.01470.09840.06310.01320.04050.0996zx0.02880.13810.15780.10120.29580.10920.013210.05660.0145nj0.12750.00370.00950.05420.059-0.01620.04050.056610.0142dt-0.0060.02880.06810.02290.1181-0.14070.09960.01450.014215s.t.1),~(IWrERWP求出均值-方差组合的有效前沿。在给定样本数据相同的情况下，假定不同的期望收益率水平P，利用前面利用matlab已求出的各只股票的收益率(R)、协方差矩阵(V)，使用matlab求解最优化问题，得到的结果如下表所示。根据前面求出的协方差矩阵，在matlab中求解上述二次规划，求得的结果如下表所示，具体的matlab程序见附录。在下表中，最优组合（zh）表示在)()~(zhrrEp；最优组合(E(R))表示)~(prE为十支股票的均值收益率下的最优组合。股票最优组合（zh）最优组合（cz）最优组合（hx）最优组合(ty)最优组合(wk)zh0.98480.01670.01760.01540.0545cz0.05150.13390.13390.1340.1307hx0.04890.08750.08750.08760.086ty0.03520.05250.05250.05250.0518wk0.03440.03970.03970.03970.0395xh0.00560.09650.09640.09660.0929zg-0.04890.13080.13060.1310.1238zx0.03460.07250.07240.07250.071nj-0.17250.33390.33350.33460.3142dt0.02640.0360.0360.0360.0356最优值3.52E-065.60E-075.60E-075.61E-075.39E-07股票最优组合(xh)最优组合(zg)最优组合(zx)最优组合(nj)最优组合(dt)最优组合(E(R))zh0.0518-0.04260.0481-0.03280.04320.1157cz0.13090.1390.13130.13810.13170.1255hx0.08610.08990.08630.08950.08650.0836ty0.05190.05360.05190.05340.0520.0507wk0.03950.040.03950.040.03960.0392xh0.09320.1020.09350.10110.0940.0872zg0.12420.14180.12490.140.12580.1124zx0.07110.07480.07120.07440.07140.0686nj0.31550.36490.31750.35980.320.2822dt0.03570.03660.03570.03650.03580.035最优值5.40E-076.17E-075.42E-076.06E-075.44E-075.27E-076从上表中可以看出当收益率取十只股票的均值时，所得到的组合最优，此时购买的十只股票的组合为：购买中航动控（zh）11.57%，购买沧州明珠（cz）12.55%，购买海信电器（hx）8.36%，购买太原刚玉（ty）5.07%，购买五矿稀土（wk）3.92%，购买新华龙（xh）8.72%，购买中国建筑（zg）11.24%,购买中信国安（zx）6.86%,购买南京熊猫（nj）28.22%，购买大同煤业（dt）3.5%。附录：部分matlab程序1.求协方差、相关系数矩阵mean(zh);--------------------------------------%求中航控股数据的均值cov(zh);-------------------------------------%求中航控股数据的方差w=[zh,cz,hx,ty,wk,xh,zg,zx,nj,dt];-------%将十只股票的收益率数据整合成矩阵cov(w);-------%求十只股票收益率的协方