1、投资的收益与风险问题摘要对市场上的多种风险资产和一种无风险资产（存银行）进行组合投资策略的设计需要考虑两个目标：总体收益尽可能大和总体风险尽可能小，而这两个目标在一定意义上是对立的。本文我们建立了投资收益与风险的双目标优化模型，并通过“最大化策略” ，即控制风险使收益最大，将原模型简化为单目标的线性规划模型一；在保证一定收益水平下，以风险最小为目标，将原模型简化为了极小极大规划模型二；以及引入收益风险偏好系数，将两目标加权，化原模型为单目标非线性模型模型三。然后分别使用 Matlab 的内部函数 linprog，fminmax，fmincon 对不同的风险水平，收益水平，以及偏好系数求解三个模
2、型。关键词：组合投资，两目标优化模型，风险偏好2.问题重述与分析3.市场上有 种资产（如股票、债券、） ( )供投资者选择，某公司有数额为 的一笔相当大的资金可用作一个时期的投资。公司财务分析人员对这 种资产进行了评估，估算出在这一时期内购买 的平均收益率为 ，并预测出购买 的风险损失率为 。考虑到投资越分散，总的风险越小，公司确定，当用这笔资金购买若干种资产时，总体风险可用所投资的 中最大的一个风险来度量。购买 要付交易费，费率为 ，并且当购买额不超过给定值 时，交易费按购买 计算（不买当然无须付费）。另外，假定同期银行存款利率是 , 且既无交易费又无风险。（ ）1、已知 时的相关数据如下：
3、 资产 收益率 (%) 风险率(%)交易费 (%) 阀值 (元)28 2.5 1 10321 1.5 2 19823 5.5 4.5 5225 2.6 6.5 40试给该公司设计一种投资组合方案，即用给定的资金 ，有选择地购买若干种资产或存银行生息，使净收益尽可能大，而总体风险尽可能小。2、试就一般情况对以上问题进行讨论，并利用以下数据进行计算。 资产 收益率 (%) 风险率(%)交易费 (%) 阀值 (元)9.6 42 2.1 18118.5 54 3.2 40749.4 60 6.0 42823.9 42 1.5 5498.1 1.2 7.6 27014 39 3.4 39740.7 68
4、 5.6 17831.2 33.4 3.1 22033.6 53.3 2.7 47536.8 40 2.9 24811.8 31 5.1 1959 5.5 5.7 32035 46 2.7 2679.4 5.3 4.5 32815 23 7.6 131本题需要我们设计一种投资组合方案，使收益尽可能大，而风险尽可能小。并给出对应的盈亏数据，以及一般情况的讨论。这是一个优化问题，要决策的是每种资产的投资额，要达到目标包括两方面的要求:净收益最大和总风险最低，即本题是一个双优化的问题，一般情况下，这两个目标是矛盾的，因为净收益越大则风险也会随着增加，反之也是一样的，所以，我们很难或者不可能提出同时满
5、足这两个目标的决策方案，我们只能做到的是：在收益一定的情况下，使得风险最小的决策，或者在风险一定的情况下，使得净收益最大，或者在收益和风险按确定好的偏好比例的情况下设计出最好的决策方案，这样的话，我们得到的不再是一个方案，而是一个方案的组合，简称组合方案。设购买 Si （i=0,1.n;S 0 表示存入银行， ）的金额为 xi;所支付的交易费为 ci(xi)，则：0()1, 2,()iii iixcxpunc对 Si 投资的净收益为: （i 0，1，n）)()(iiixrxR对 Si 投资的风险为: （i0，1，n） ，q 0=0iqxQ)(对 Si 投资所需资金（即购买金额 xi 与所需的手
6、续费 ci(xi) 之和）是（i 0，1，n）)()(iiicf投资方案用 x=（x 0，x 1，x n）表示，那么，净收益总额为：0()()niRx总风险为：= )(xQ)(mi0inx所需资金为：)()(0inixfF所以，总收益最大，总风险最小的双目标优化模型表示为：0,)(minxMFxRQx但是像这样的双目标模型用一般的方法很难求解出来的，所以经过分析把次模型转化为三种较简单的单目标模型。3.假设与模型假设该公司在这一时期内是一次性投资；除交易费和投资费用外再无其他的费用开支；在这一时期市场发展基本上是稳定的；外界因素对投资的资产无较大影响；无其他的人为干预；社会政策无较大变化；公司
7、的经济发展对投资无较大影响资产投资是在市场中进行的，市场是复杂多变的，是无法用数量或函数进行准确描述的，因此以上的假设是必要的，一般说来物价变化具有一定的周期性，社会政策也并非天天改变，公司自身的发展在稳定的情况下才会用额外的资金进行较大的风险的投资， 市场与社会的系统发展在一个时期内是良性的、稳定的，以上假设也是合理的。3.1 模型 a假设投资的风险水平是 k,即要求总风险 Q（x）限制在 k 内，Q （x） ，则模型可转化为：kmax Rs.t 0,)(,xMFkxQ3.2 模型 b假设投资的收益水平是 h，即净收益总额 )(xR不少于 h： )(xRh，则模型可转化为：)(minxQs.
8、t 0,)(,)(xMFhxR3.3 模型 c假设投资者对风险和收益的相对偏好参数为 （0） ，则模型可转化为：)(1)(minxRxQs.t. 0,MF3.4 模型求解及分析由于交易费 ci(xi)是分段函数，使得上述模型中的目标函数或约束条件相对比较复杂，是一个非线性规划问题，难于求解. 但注意到总投资额 M 相当大，一旦投资资产 Si，其投资额 xi 一般都会超过 ui，于是交易费 ci(xi)可简化为线性函数iixpc)(从而，资金约束简化为 00()()(1)nni iiFfxpxM净收益总额简化为 000()()()()nnniiiiiiRxrcxrpx在实际进行计算时，可设 M=
9、1，此时=（ ） （i 0，1，n）iyip1x可视作投资 Si 的比例. 以下的模型求解都是在上述两个简化条件下进行讨论的.1）模型 a 的求解模型 a 的约束条件 Q(x)k 即k ，00()max()a()i iininQqx所以此约束条件可转化为（i 0，1，n）.kxqi这时模型 a 可化简为如下的线性规划问题： 00max()s.t , =12, ()0niiiniirpxqkn具体到 n=4 的情形，按投资的收益和风险问题中题中给定的数据，模型为： 43210 185.0.9.7.5.max xxs.t kkkk4321 26.,05.,.,2.（i 0，1，4）,6.0450.
10、1432 ixxx利用 matlab7.1 求解模型 a输出结果是0.177638, x0 - 0.158192, x1 - 0.2, x2 - 0.333333, x3 - 0.0909091,x4 - 0.192308这说明投资方案为（0.158192，0.2，0.333333，0.0909091，0.192308）时，可以获得总体风险不超过 0.005 的最大收益是 0.177638M.当 k 取不同的值（00.025） ，风险与收益的关系见图 1.输出结果列表如下：表 1 模型 1 的计算结果风险水平 k 最大收益 0 0.05 1 0 0 0 00.001 0.0755 0.8316
11、 0.04 0.0667 0.0182 0.03850.002 0.1011 0.6633 0.08 0.1333 0.0364 0.07690.003 0.1266 0.4949 0.12 0.2 0.0545 0.11540.004 0.1521 0.3266 0.16 0.2667 0.0727 0.15380.005 0.1776 0.1582 0.2 0.3333 0.0909 0.19230.006 0.2019 0 0.24 0.4 0.1091 0.22120.007 0.2066 0 0.28 0.4667 0.1273 0.10160.008 0.2112 0 0.32 0
12、.5333 0.1271 00.009 0.2155 0 0.36 0.6 0.0233 00.01 0.219 0 0.4 0.5843 0 00.011 0.2223 0 0.44 0.5447 0 00.012 0.2256 0 0.48 0.5051 0 00.013 0.2288 0 0.52 0.4655 0 00.014 0.2321 0 0.56 0.4259 0 00.015 0.2354 0 0.6 0.3863 0 00.016 0.2387 0 0.64 0.3467 0 00.017 0.2419 0 0.68 0.3071 0 00.018 0.2452 0 0.7
13、2 0.2675 0 00.019 0.2485 0 0.76 0.2278 0 00.02 0.2518 0 0.8 0.1882 0 00.021 0.255 0 0.84 0.1486 0 00.022 0.2583 0 0.88 0.109 0 00.023 0.2616 0 0.92 0.0694 0 00.024 0.2649 0 0.96 0.0298 0 00.025 0.2673 0 0.9901 0 0 00 0.005 0.01 0.015 0.02 0.0250.050.10.150.20.250.3偏偏 a偏偏图 1 模型 1 中风险 k 与收益的关系结合图 1，对于
14、风险和收益没有特殊偏好的投资者来说，应该选择图中曲线的拐点（0.006，0.2019） ，这时对 的投资比例见表 1 的黑体所示。从表 1 中的计算结果可以看出，对低风险水平，除了存入银行外，投资首选风险率最低的 S2，然后是 S1 和 S4，总收益较低；对高风险水平，总收益较高，投资方向是选择净收益率（r ipi）较大的 S1 和 S2这些与人们的经验是一致的，这里给出了定量的结果2）模型 b 的求解模型 b 本来是极小极大规划： 0minax()iiqs.t. h x00()niirp0(1)niip但是，可以引进变量 xn+1= ，将它改写为如下的线性规划：0a()iiqx1min()x
15、s.t ,i=0,1,2,n, h, , x01nixq0()iirp0()1niip具体到 n=4 的情形，按投资的收益和风险问题中题中给定的数据，模型为：min x5s.t 54535251 026.,0.,0.,2.0 xxxx ,189.7 43210 h（i 0，1，5）,065.4.2.1. 430 ixxx利用 matlab7.1 求解模型 b，当 h 取不同的值（0.10.25） ，我们计算最小风险和最优决策,收益水平 h 取 ，结果如表 2 所示，风险和收益的关系见图 2.从表 2 看出，对低收益水平，除了存入银行外，投资首选风险率最低的资产 ，然后是 和 ，总收益当然较低。
16、对高收益水平，总风险自然也高，应首选净收益率（ ）最大的 和 。这些与人们的经验是一致的。表 2 模型 2 的计算结果风险水平 k 最大收益 0.002 0.1 0.6702 0.0783 0.1306 0.0356 0.07530.0024 0.11 0.6043 0.094 0.1567 0.0427 0.09040.0027 0.12 0.5383 0.1097 0.1828 0.0499 0.10550.0031 0.13 0.4724 0.1254 0.2089 0.057 0.12050.0035 0.14 0.4064 0.141 0.235 0.0641 0.13560.003
17、9 0.15 0.3405 0.1567 0.2612 0.0712 0.15070.0043 0.16 0.2745 0.1724 0.2873 0.0783 0.16570.0047 0.17 0.2086 0.188 0.3134 0.0855 0.18080.0051 0.18 0.1426 0.2037 0.3395 0.0926 0.19590.0055 0.19 0.0767 0.2194 0.3656 0.0997 0.21090.0059 0.2 0.0107 0.235 0.3917 0.1068 0.2260.0078 0.21 0 0.3114 0.519 0.0569
18、 0.09080.0103 0.22 0 0.412 0.5725 0 00.0134 0.23 0 0.5341 0.4515 0 00.0164 0.24 0 0.6563 0.3305 0 00.0195 0.25 0 0.7784 0.2096 0 00.002 0.004 0.006 0.008 0.01 0.012 0.014 0.016 0.018 0.020.10.150.20.25偏偏 偏偏图 2 模型 2 中风险与收益 h 的关系结合图 2，对于风险和收益没有特殊偏好的投资者来说，应该选择图中曲线的拐点（0.059，0.2） ，这时对 的投资比例见表 2 的黑体所示。3）模
19、型 c 的求解类似模型 b 的求解，我们同样引进变量 xn+1= ，将它改写为如下的线性规划：0ma()iiqxmin xn+1（1） 0()niirps.t ,i=0，1，2，n x0nixq0(1)nii具体到 n=4 的情形，按投资的收益和风险问题题中给定的数据，模型为： )185.0.19.027.05.)(1mi 4325 xxxx s.t 54,535251 6,.,2. （i 0，1，5）,06.14.0 43210 ixxx利用 matlab7.1 求解模型 c，当 取不同的值（0.750.95 ） ，我们计算最小风险和最优决策输出结果列表如下：表 3 模型 3 的计算结果风险
20、水平 最大收益率 0.76 0.0248 0.2673 0 0.9901 0 0 00.77 0.0248 0.2673 0 0.9901 0 0 00.78 0.0092 0.2165 0 0.3693 0.6147 0 00.79 0.0092 0.2165 0 0.3692 0.6148 0 00.8 0.0079 0.2107 0 0.3151 0.5218 0.1431 00.81 0.0092 0.2165 0 0.3687 0.6154 0 00.82 0.0074 0.2084 0 0.2957 0.4928 0.1344 0.05470.83 0.0079 0.2106 0
21、0.3139 0.5235 0.1425 00.84 0.006 0.2017 0 0.2386 0.3973 0.1082 0.2260.85 0.006 0.2016 0 0.2375 0.3967 0.1084 0.22750.86 0.0062 0.2026 0 0.2468 0.4113 0.0753 0.23720.87 0.0066 0.2045 0 0.2653 0.4422 0.0088 0.25520.88 0.0059 0.2017 0 0.2379 0.3964 0.1078 0.22790.89 0.006 0.2017 0 0.239 0.3937 0.1079 0
22、.22940.9 0.006 0.2018 0 0.2392 0.398 0.1029 0.22990.91 0.0059 0.2017 0 0.238 0.3958 0.1075 0.22870.92 0.0059 0.2016 0 0.2376 0.3961 0.108 0.22840.93 0.0059 0.2016 0 0.2375 0.3963 0.108 0.22820.94 0.0061 0.1847 0.135 0.2436 0.3309 0.1106 0.15570.95 0.0059 0.2016 0 0.2376 0.396 0.108 0.22850.96 0.0017
23、 0.0916 0.7256 0.0652 0.1065 0.0302 0.06410.97 0.0054 0.1848 0.114 0.2141 0.3564 0.0974 0.192从图 5 可以看出，模型 3 的风险与收益关系与模型 1 和模型 2 的结果几乎完全一致。0 0.005 0.01 0.015 0.02 0.0250.080.10.120.140.160.180.20.220.240.260.28偏偏 偏偏图 3 模型 3 中风险与收益 的关系0.75 0.8 0.85 0.9 0.95 100.0050.010.0150.020.025偏偏偏偏 偏偏图 4 模型 3 中风险
24、与偏好系数的关系0.75 0.8 0.85 0.9 0.95 10.080.10.120.140.160.180.20.220.240.260.28偏偏偏偏 偏偏图 5 模型 3 中收益与偏好系数的关系四 模型评价与推广本文我们建立了投资收益与风险的双目标优化模型，通过控制风险使收益最大，保证收益使风险最小，以及引入收益风险偏好系数，将两目标模型化为了单目标模型，并使用 matlab7.1 求解，所得结果具有一定的指导意义。但是，本文没有讨论收益和风险的评估方法，在实际应用中还存在资产相关的情形，此时，用最大风险代表组合投资的风险未必合理，因此，对不同风险度量下的最优投资组合进行比较研究是进一
25、步的改进方向。五 总结历经两周的时间终于完成了这次课设，在这次实践课程中，我真的遇到了不少的问题，在同学，老师的帮助以及在图书馆和网站搜集资料，解决了所有遇到的问题。尤其在问题分析的过程中，是难度最大也是问题最多的环节，感觉总是把问题分析的不够全面透彻，经常顾及这个方面而忽视了另一方面，最后我请教了同学，终于完成了问题分析并且建立了模型。在完成这一环节后，接下来的任务都是我独立完成，也遇到了不少的困难，但都是较易解决的。通过这次实践，我确实学到了不少，学会了使用 MATLAB，也知道了分析问题的方法。六 参考文献1MATLAB 程序设计与实例应用。张铮等。北京：中国铁道出版社，2003.102
26、运筹学方法与应用。吴风平。南京：河海大学出版社，2000.123数学模型及方法 。李火林主编。江西高校出版社， 1997.104数学建模教育及竞赛。甘筱青主编。南昌：江西高校出版社。2004.65 萧树铁，面向 21 世纪课程教材：大学数学数学实验，北京：高等教育出版社，1999.7.6赫孝良，戴永红等编著，数学建模竞赛：赛题简析与论文点评，西安：西安交通大学出版社，2002.6.7陈叔平，谭永基，一类投资组合问题的建模与分析，数学的实践与认识， （29）7：45-49，1999.七 附录function result=qiujie()%data为表格数据 data=28 2.5 1 1032
27、1 1.5 2 19823 5.5 4.5 5225 2.6 6.5 40;data1=9.6 42 2.1 18118.5 54 3.2 40749.4 60 6.0 42823.9 42 1.5 5498.1 1.2 7.6 27014 39 3.4 39740.7 68 5.6 17831.2 33.4 3.1 22033.6 53.3 2.7 47536.8 40 2.9 24811.8 31 5.1 1959 5.5 5.7 32035 46 2.7 2679.4 5.3 4.5 32815 23 7.6 131;data=5 0 0 0;data./100;%增加存银行r=data
28、(:,1);q=data(:,2);p=data(:,3);% %模型一求解% result=;% for a=0:0.01:0.5% result=result;moxing1(r,q,p,a);% end% result=round(result.*10000)./10000;% plot(result(:,1),result(:,2)% grid on% xlabel(风险 )% ylabel(收益)% %模型二求解% result=;% for k=0.1:0.01:0.4% result=result;moxing2(r,q,p,k);% end% result=round(resu
29、lt.*10000)./10000;% plot(result(:,1),result(:,2)% grid on% xlabel(风险 )% ylabel(收益)%模型三求解result=;for s=0.76:0.01:0.97result=result;moxing3(r,q,p,s);endresult=round(result.*10000)./10000;figure(1)plot(result(:,2),result(:,3)grid onxlabel(风险 )ylabel(收益 )pausefigure(2)plot(result(:,1),result(:,2)grid on
30、xlabel(偏好系数 )ylabel(风险 )pausefigure(3)plot(result(:,1),result(:,3)grid onxlabel(偏好系数 )ylabel(收益 )function result1=moxing1(r,q,p,a)%线性规划模型%r收益率，为列向量%p交易费率，为列向量%q风险率，为列向量%a风险水平f=(p-r);%转为求极小n=length(q);I=eye(n);for i=2:nA(i-1,:)=q(i)*I(i,:);endb=a*ones(n-1,1);Aeq=(1+p);beq=1;lb=zeros(n,1);ub=;x,fval,e
31、xitflag,output=linprog(f,A,b,Aeq,beq,lb,ub);result1=a,-fval,x;endfunction result2=moxing2(r,q,p,k)%极小极大模型%r收益率，为列向量%p交易费率，为列向量%q风险率，为列向量%k收益水平n=length(q);f=(x)q.*x(1:n);A=(p-r);b=-k;Aeq=1+p;beq=1;lb=zeros(n,1);ub=;x0=rand(n,1);x,fval,maxfval,exitflag=fminimax (f,x0,A,b,Aeq,beq,lb,ub);result2=max(q.*x),k,x;endfunction result3=moxing3(r,q,p,s)%极小极大模型%r收益率，为列向量%p交易费率，为列向量%q风险率，为列向量%s投资偏好系数n=length(q);f=(x)(s*max(q.*x)-(1-s)*sum(r-p).*x);A=;b=;Aeq=1+p;beq=1;lb=zeros(n,1);ub=;x0=rand(n,1);x,fval,exitflag,output = fmincon(f,x0,A,b,Aeq,beq,lb,ub);result3=s,max(q.*x),sum(r-p).*x),x;endend
