投资收益和风险的模型

1、投资收益和风险的模型0 引言在现代商业、金融的投资中，任何理性的投资者总是希望收益能够取得最大化，但是他也面临着不确定性和不确定性所引致的风险。而且，大的收益总是伴随着高的风险。在有很多种资产可供选择，又有很多投资方案的情况下，投资越分散，总的风险就越小。为了同时兼顾收益和风险，追求大的收益和小的风险构成一个两目标决策问题，依据决策者对收益和风险的理解和偏好将其转化为一个单目标最优化问题求解。随着投资者对收益和风险的日益关注,如何选择较好的投资组合方案是提高投资效益的根本保证。传统的投资组合遵循“不要将所有的鸡蛋放在一个蓝子里”的原则, 将投资分散化。 1 问题的提出某公司有数额为M（较大）的

2、资金，可用作一个时期的投资，市场上现有5种资产（Si）(如债券、股票等)可以作为被选的投资项目，投资者对这五种资产进行评估，估算出在这一段时期内购买Si的期望收益率（ri）、交易费率(pi)、风险损失率(qi)以及同期银行存款利率r0（r0=3%）在投资的这一时期内为定值如表1，不受意外因素影响,而净收益和总体风险只受ri,pi,qi影响，不受其他因素干扰 。现要设计出一种投资组合方案, 使净收益尽可能大, 风险尽可能小.表1投资项目Si 存银行S0期望收益率ri(%) 风险损失率qi(%)3 27 22 25 23 210 2.4 1.6 5.2 2.2 1.5交易费率pi(%)0 1 2

3、4.5 6.5 2S1 S2 S3 S4 S5其中i=0,1,2,3,4,5.2 问题假设及符号说明2.1 问题假设(1)总体风险可用投资的这五种中最大的一个风险来度量；(2)在投资中,不考虑通货膨胀因素, 因此所给的Si的期望收益率ri为实际的平均收益率；(3)不考虑系统风险, 即整个资本市场整体性风险, 它依赖于整个经济的运行情况, 投资者无法分散这种风险, 而只考虑非系统风险, 即投资者通过投资种类的选择使风险有所分散；(4)不考虑投资者对于风险的心理承受能力。2.2 符号说明xi：购买第i种资产的资金数额占资金总额的百分比；Mxi：购买第i种资产的资金数额； Mx0：存银行的金额； f

4、(xi)：交易费用； R：净收益；Q：总体风险； i：第i种投资的净收益率。 3 模型的分析与建立令交易费用Mxipi,xi>0(i=0,1, ,5) f(xi)=0,x=0i则净收益为R=M(1+ri)xi-Mi=05总体风险为Q= maxMxiqi 0i5约束条件为f(x)+Mxii=0i=055i=M可以简化约束条件为(1+p)xii=05i=1同时将M=M(1+pi)xi代入，得i=05R=M(1+ri)xi-M(1+pi)xi=M(ri-pi)xii=0i=0i=0555略去M,原问题化为双目标决策问题:max R=xi(ri-pi)i=05minQ= maxxiqi （3.1

5、） 0i55(1+pi)xi=1i=0 s. t .xi0i=0,1, ,5以下设ri-pi>0，否则不对该资产投资。 4 模型的求解4.1 固定R使Q最小的模型固定R使Q最小，将模型（3.1）化为minQ= maxqixi， 0i55(ri-pi)xi=R,(1)i=05s. t . (1+pi)xi=1,(2) （4.1）i=0xi0i=0,1, ,5此模型又可改写为min y(r0-p0)x0+(r1-p1)x1+ +(r5-p5)x5=R(1+p0)x0+(1+p1)x1+ +(1+p5)x5=1s. t . xiqiyxi0,y0i=0,1, ,5令i=(ri-pi)(1+pi

6、)，则i必大于0,否则, i表示第i种投资的净收益率，若10, 则不对Si投资, 因为对该项目投资纯收益率不如存银行, 而风险损失率又大于存银行。将i从小到大排序,设k最大, 则易见对模型（4.1）的可行解必有0.03Rk.当R=0.03时, 所有资金都存银行,Q=0; 当R=k时, 所有资金用于购买Si , Q= qk；当0.03<R<k时,有如下结论7。 1+pk7x1q1= =x5q5。而对于固定收益使风险最小的模型来说，这结论也可换句话说：在前5项投资总额一定的前提下，各项投资的风险损失相等即x1q1=x2q2= =x5q5时，总体风险最小8。证：设y1,y2, ,y5是满

7、足x1q1=x2q2= =x5q5的一组解，即y1q1=y2q2= =y5q5=Q*。显然此时Q*为总体风险。由于前5项投资总额M是一定的，只要改变其中一项的值，便会导致总体风险增加。（比如说将y1的值增加为y1*会使得y1*q1>Q*，总体风险显然增加；反之，若减小y1的值，必然会导致另外一项或几项的值，总体风险自然增加。）因此,当R(0.03,k)时，可按以下步骤求出最优解：1）将（1）式和（2）5QQ式消去x0；2）将xi=代入解出Q；3）由xi=，1i5，x0=1-(1+pi)xiqiqii=1求出最优解。所以，我们算得如下结果:（1）R=0.03时，x0=1,x1=x2=x3=

8、x4=x5=0,Q=0；（2）R=0.26时，x0=x2=x3=x4=x5=0,x1=,Q=；表1最小风收益R险度Q 投资Si的资金百分比xi （i=0,1,2,3,4,5.） x0 x1 x20.0000 x3 0.0000 x4 0.0000 x5 0.0000 0.03000.04000.05000.06000.07000.08000.09000.10000.11000.12000.13000.14000.15000.16000.17000.18000.19000.20000.2100 0.0000 1.0000 0.0000 0.0002 0.9397 0.0104 0.0156 0.

9、0048 0.0113 0.0166 0.0311 0.0096 0.0226 0.0332 0.0467 0.0144 0.0339 0.0498 0.0005 0.8793 0.0207 0.0007 0.8190 0.0311 0.0010 0.7587 0.0415 0.0622 0.0191 0.0453 0.0664 0.0012 0.6984 0.0519 0.0778 0.0239 0.0566 0.0830 0.0015 0.6380 0.0622 0.0933 0.0287 0.0679 0.0996 0.0017 0.5777 0.0726 0.1089 0.0335 0

10、.0792 0.1162 0.0020 0.5174 0.0830 0.1245 0.0383 0.0905 0.1328 0.0022 0.4571 0.0933 0.1400 0.0431 0.1018 0.1494 0.0025 0.3967 0.1037 0.1556 0.0479 0.1131 0.1660 0.0027 0.3364 0.1141 0.1711 0.0527 0.1245 0.1825 0.0030 0.2761 0.1245 0.1867 0.0574 0.1358 0.1991 0.0032 0.2158 0.1348 0.2023 0.0622 0.1471

11、0.2157 0.0035 0.1554 0.1452 0.2178 0.0670 0.1584 0.2323 0.0037 0.0951 0.1556 0.2334 0.0718 0.1697 0.2489 0.0040 0.0348 0.1660 0.2489 0.0766 0.1810 0.2655 0.0046 0.0000 0.1897 0.2846 0.0876 0.1097 0.3036 0.2589 0.3884 0.1195 0.0000 0.2132 0.0062 0.00000.2200 0.2300 0.2400 0.2500 0.26/1.010.0093 0.000

12、0 0.3858 0.4160 0.1781 0.00000.5471 0.1800 0.2525 0.0000 0.7084 0.00000.2722 0.0000 0.1160 0.0000 0.00000.00000.0000 0.0000 0.0000 0.0000 0.00000.0131 0.0000 0.01700.00000.0209 0.0000 0.8701 0.0000 0.0238 0.00000.99010.0000最小风险度Q随收益R的变化趋势图最小风险度Q收益R4.2 固定Q使R最大的模型max R=(ri-pi)xi，i=05xiqiQ,5s. t .(1+pi

13、)xi=1,i=0xi0,(i=0,1, ,5.)xi=Qqi。证明：反证法。假设i>j，xj>0，而xi<Qqi。选取充分小的正数，使得(xi+)qi<Q，(1+pi)<xj(1+pj)。令xi*=xi+，xj*=xj-(1+pi)(1+pj)，当ki,j时，令xk*=xk，则xk*0，且xk=055*k(1+pk)=ki,j*x*k*k(1+pk)+(xi+)(1+pi)+xj-(1+pi)(1+pj)(1+pj)=1， 5x由此结论, 我们可将i从大到小排序, 使i最大的k 应尽量满足xkqk=Q, 若还有多余资金, 再投资i次大的, 。对于不同的Q ,会有

14、不同的投资方案, 我们可以算出Q的临界值, 从而确定各项目的投资值。因此，设1>2>3>4>5>0 , 则可用下面的方法算出各临界值c1，c2，c3，c4，c5。只有一种投资时, c1(1+p1)=q1,c1=q1(1+p1)=0.023762。当有两种投资时, 将x1=c2q1,x2=c2q2，代入x1(1+p1)+x2(1+p2)=1，得c2=q1q2+p1)q2+(1+p2)q1=0.009449。 同理可得：c3=q1q2q3+p1)q2q3+(1+p2)q1q3+(1+p3)q1q2=0.007941， c4=q1q2q3q4+p1)q2q3q4+(1+

15、p2)q1q3q4+(1+p3)q1q2q4+(1+p4)q1q2q3=0.005736c5=q1q2q3q4q5+p1)q2q3q4q5+(1+p2)q1q3q4q5+(1+p3)q1q2q4q5+(1+p4)q1q2q3q5+(1+p5)q1q2q3q4=0.004131于是得最优解：当Q=0.000000时，x0=1,x1=x2=x3=x4=x5=0。当0<Q0.004131时，x1=Qq1,x2=Qq2,x3=Qq3,x4=Qq4,x5=Qq5,x0=1-(1+pi)xi。i=15当0.004131<Q0.005736时，x1=Qq1,x2=Qq2,x3=Qq3,x4=Qq

16、4,x5=1-(1+pi)xi+p5),x0=0。i=14当0.005736<Q0.007941时，x1=Qq1,x2=Qq2,x3=Qq3,x4=1-(1+pi)xi+p4),x5=x0=0。i=13当0.007941<Q0.009449时，x1=Qq1,x2=Qq2,x3=1-(1+pi)xi+p3),x4=x5=x0=0。I=12当0.009449<Q0.023762时，x1=Qq1,x2=1-(1+p1)x1(1+p2),x3=x4=x5=x0=0。当Q>0.023762时，x1=+p1),x2=x3=x4=x5=x0=0。当然，我们也可以换个角度来考虑上面这个

17、模型。为了能够给不同风险承受能力的投资者提供某种风险水平下的最优投资组合的决策方案，我们必须确定最优收益值R和最小风险度Q的值之间的对应关系。maxR=(r0-p0)x0+(r1-p1)x1+ +(r5-p5)x5，x0+1.01x1+1.02x2+1.045x3+1.065x4+1.02x5=1,0.024x1Q,0.016x2Q,0.052x3Q, s.t0.022x4Q,0.015x5Q,xi0,(i=0,1,2,3,4,5.)为此编写MATLAB程序（见附录），从风险度Q=0开始，以每次增加0.001的风险度进行搜索5。根据附录中程序一，最优收益值R和最小风险度Q以及投资额分配之间的对

18、应关系计算结果列表如下：风险度最优收益投资Si的资金百分比xi （i=0,1,2,3,4,5.）Q0 0.0010 0.0020 0.0030 0.0040 0.0050 0.0060 0.0070 0.0080 0.0090 0.0100 0.0110 0.0120 0.0130 0.0140 0.0150 0.0160R0.0300 0.0702 0.1103 0.1505 0.1907 0.2044 0.2092 0.2130 0.2167 0.2193 0.2219 0.2245 0.2271 0.2297 0.2322 0.2348 0.2374x01.0000 0.7577 0.5

19、153 0.2730 0.0306 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000x10.0000 0.0417 0.0833 0.1250 0.1667 0.2083 0.2500 0.2917 0.3333 0.3750 0.4167 0.4583 0.5000 0.5417 0.5833 0.6250 0.6667x20.0000 0.0625 0.1250 0.1875 0.2500 0.3125 0.3750 0.4375 0.4927 0.4317 0.3708

20、0.5266 0.2489 0.1879 0.1269 0.0660 0.0051x30.0000 0.0192 0.0385 0.0577 0.0769 0.0962 0.1154 0.1346 0.1538 0.1731 0.1923 0.0000 0.2308 0.2500 0.2692 0.2885 0.3077x40.0000 0.0455 0.0909 0.1364 0.1818 0.0285 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000x50.0000 0.0667 0.1

21、333 0.2000 0.2667 0.3333 0.2396 0.1162 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.00000.0170 0.0180 0.0190 0.0200 0.0210 0.0220 0.0230 0.0240 0.0250 0.09900.2400 0.2426 0.2451 0.2477 0.2503 0.2529 0.2555 0.2574 0.2574 0.25740.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.000

22、0 0.00000.7083 0.7500 0.7917 0.8333 0.8750 0.9167 0.9583 0.9901 0.9901 0.99010.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.00000.2723 0.2321 0.1918 0.1515 0.1112 0.0710 0.0307 0.0000 0.0000 0.00000.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.00000.0000 0.0000 0.

23、0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000最优收益R随风险度Q的变化趋势图最优收益R0.010.020.030.04风险度Q0.070.080.090.1按照收益风险最大原则, 可取模型max，5(1+pi)xi=1,s. ti=0x0,(i=0,1, ,5.)i由于q0=0，因而取x0=1,x1=x2= =x5=0时，max=+。当然，也可取模型minQR，5(1+pi)xi=1,s.ti=0x0,(i=0,1, ,5.)i同上，由于q0=0，因而取x0=1,x1=x2= =x5=0时，minQR=0，从而可知, 全部钱存银行是

24、最优解。对于此问题, 其他投资的收益与风险损失率都不影响该最优解, 故这种模型不够好。由偏好系数法, 我们选取偏好系数(01),建立模型max(1-)R-y，5xi(1+pi)=1,i=0txiqiy, s.x0,(i=0,1,2,3,4,5.)i具体数据可应用参数规划法进行计算。权重r最小风 险度Q投资Si的资金百分比xi （i=0,1,2,3,4,5）x00.0000 0.0000 0.0000 0.0000x10.9901 0.3309 0.2149 0.1719x20.0000 0.4963 0.3223 0.2579x30.0000 0.1527 0.0992 0.0794x40.0

25、000 0.0000 0.0000 0.1876x50.0000 0.0000 0.3438 0.27510,0.7200 0.02380.7210,0.7920 0.0079 0.7930,0.9070 0.0052 0.9090,0.9750 0.00410.9760,10.0000 1.0000 0.0000 0.0000 0.0000 0.0000 0.0000 风险度Q随权重r的变化趋势图附录一模型一Lingo 语句min y(r0-p0)x0+(r1-p1)x1+ +(r5-p5)x5=R(1+p0)x0+(1+p1)x1+ +(1+p5)x5=1s. t . xiqiyxi0,y

26、0,(i=0,1, ,5.)min=y;0.03*x0+(0.27-0.01)*x1+(0.22-0.02)*x2+(0.25-0.045)*x3+(0.23-0.065)*x4+(0.21-0.02)*x5=0.03;x0+1.01*x1+1.02*x2+1.045*x3+1.065*x4+1.02*x5=1;0.024*x1<=y;0.016*x2<=y;0.052*x3<=y;0.022*x4<=y;0.015*x5<=y;模型一Matlab 程序>> R=0.03>> while R<0.26/1.01;C= 0 0 0 0

27、0 0 1;A= 0 0.024 0 0 0 0 -1;0 0 0.016 0 0 0 -1;0 0 0 0.052 0 0 -1;0 0 0 0 0.022 0 -1;0 0 0 0 0 0.015 -1;B= 0;0;0;0;0;Aeq= 0.03 0.26 0.2 0.205 0.165 0.19 0;1 1.01 1.02 1.045 1.065 1.02 0; Beq= R;1;Vlb= 0;0;0;0;0;0;0;% or Vlb= zeros(7,1);Vub= ;x,fval= linprog(C,A,B,Aeq,Beq,Vlb,Vub);RQ=fvalx=x'plot

28、(R, Q, 'm.')axis(0 0.3 0 0.03)xlabel('收益R')ylabel('最小风险度Q')title('最小风险度Q随收益R的变化趋势图')hold onR=R+0.01;grid onendR=0.26/1.01;C= 0 0 0 0 0 0 1;A= 0 0.024 0 0 0 0 -1;0 0 0.016 0 0 0 -1;0 0 0 0.052 0 0 -1;0 0 0 0 0.022 0 -1;0 0 0 0 0 0.015 -1;B= 0;0;0;0;0;Aeq= 0.03 0.26 0.2

29、 0.205 0.165 0.19 0;1 1.01 1.02 1.045 1.065 1.02 0; Beq= R;1;Vlb= 0;0;0;0;0;0;0;% or Vlb= zeros(7,1);Vub= ;x,fval= linprog(C,A,B,Aeq,Beq,Vlb,Vub)程序二模型二Matlab 程序>> Q=0>> while (1.1-Q)>1 % or Q<0.1;C= -0.03 -0.26 -0.20 -0.205 -0.165 -0.19;A= 0 0.024 0 0 0 0;0 0 0.016 0 0 0;0 0 0 0.05

30、2 0 0;0 0 0 0 0.022 0;0 0 0 0 0 0.015;B= Q;Q;Q;Q;Q;Aeq= 1 1.01 1.02 1.045 1.065 1.02;Beq= 1;Vlb= 0;0;0;0;0;0;% or Vlb= zeros(5,1);Vub= ;x,fval= linprog(C,A,B,Aeq,Beq,Vlb,Vub);QR=-fvalx=x'plot(Q,R,'m.')axis(0 0.1 0 0.5)xlabel('风险度Q')ylabel('最优收益R')title('最优收益R随风险度Q的变化趋

31、势图')hold onQ=Q+0.001;grid onenda=0;while(1.1-a)>1c=-0.05 -0.27 -0.19 -0.185 -0.185;Aeq=1 1.01 1.02 1.045 1.065;beq=1;A=0 0.025 0 0 0;0 0 0.015 0 0;0 0 0 0.055 0;0 0 0 0 0.026;b=a;a;a;a;vlb=0,0,0,0,0;vub=;x,val=linprog(c,A,b,Aeq,beq,vlb,vub);ax=x'Q=-valplot(a,Q,'.')axis(0 0.1 0 0.5)hold ona=a+0.001;endxlabel('a'),ylabel('Q')模型三Lingo 语句max(1-)R-y，5xi(1+pi)=1,i=0txiqiy, s.x0,(i=0,1,2,3,4,5.)imax=(1-0.2)*(