数学建模：投资问题.docx

投资的收益与风险问题
摘要
对市场上的多种风险资产和一种无风险资产（存银行）进行组合投资策略的设计需要考虑两个目标：总体收益尽可能大和总体风险尽可能小，而这两个目标在一定意义上是对立的。
本文我们建立了投资收益与风险的双目标优化模型，并通过“最大化策略”，即控制风险使收益最大，将原模型简化为单目标的线性规划模型一；在保证一定收益水平下，以风险最小为目标，将原模型简化为了极小极大规划模型二；以及引入收益——风险偏好系数，将两目标加权，化原模型为单目标非线性模型模型三。然后分别使用Matlab的内部函数linprog，fminmax，fmincon对不同的风险水平，收益水平，以及偏好系数求解三个模型。
关键词：组合投资，两目标优化模型，风险偏好

2. 问题重述与分析
3. 市场上有种资产（如股票、债券、…）()供投资者选择，某公司有数额为的一笔相当大的资金可用作一个时期的投资。公司财务分析人员对这种资产进行了评估，估算出在这一时期内购买的平均收益率为，并预测出购买的风险损失率为。考虑到投资越分散，总的风险越小，公司确定，当用这笔资金购买若干种资产时，总体风险可用所投资的中最大的一个风险来度量。
购买要付交易费，费率为，并且当购买额不超过给定值时，交易费按购买计算（不买当然无须付费）。另外，假定同期银行存款利率是, 且既无交易费又无风险。（）
1、已知时的相关数据如下：
资产
收益率(%)
风险率(%)
交易费(%)
阀值(元)

28
2.5
1
103

21
1.5
2
198

23
5.5
4.5
52

25
2.6
6.5
40
试给该公司设计一种投资组合方案，即用给定的资金，有选择地购买若干种资产或存银行生息，使净收益尽可能大，而总体风险尽可能小。
2、试就一般情况对以上问题进行讨论，并利用以下数据进行计算。
资产
收益率(%)
风险率(%)
交易费(%)
阀值(元)

9.6
42
2.1
181

18.5
54
3.2
407

49.4
60
6.0
428

23.9
42
1.5
549

8.1
1.2
7.6
270

14
39
3.4
397

40.7
68
5.6
178

31.2
33.4
3.1
220

33.6
53.3
2.7
475

36.8
40
2.9
248

11.8
31
5.1
195

9
5.5
5.7
320

35
46
2.7
267

9.4
5.3
4.5
328

15
23
7.6
131

　本题需要我们设计一种投资组合方案，使收益尽可能大，而风险尽可能小。并给出对应的盈亏数据，以及一般情况的讨论。
这是一个优化问题，要决策的是每种资产的投资额，要达到目标包括两方面的要求:净收益最大和总风险最低，即本题是一个双优化的问题，一般情况下，这两个目标是矛盾的，因为净收益越大则风险也会随着增加，反之也是一样的，所以，我们很难或者不可能提出同时满足这两个目标的决策方案，我们只能做到的是：在收益一定的情况下，使得风险最小的决策，或者在风险一定的情况下，使得净收益最大，或者在收益和风险按确定好的偏好比例的情况下设计出最好的决策方案，这样的话，我们得到的不再是一个方案，而是一个方案的组合，简称组合方案。
设购买Si （i=0,1…….n;S0表示存入银行，）的金额为xi;所支付的交易费为ci(xi)，则：

对Si 投资的净收益为:　（i＝0，1，…，n）
对Si投资的风险为: 　（i＝0，1，…，n），q0=0
　对Si 投资所需资金（即购买金额 xi 与所需的手续费 ci(xi) 之和）是
　（i＝0，1，…，n）
　 投资方案用 x=（x0，x1，…，xn）表示，那么，
净收益总额为：

总风险为：
=
所需资金为：

所以，总收益最大，总风险最小的双目标优化模型表示为：

但是像这样的双目标模型用一般的方法很难求解出来的，所以经过分析把次模型转化为三种较简单的单目标模型。
3.假设与模型
假设该公司在这一时期内是一次性投资；除交易费和投资费用外再无其他的费用开支；在这一时期市场发展基本上是稳定的；外界因素对投资的资产无较大影响；无其他的人为干预；社会政策无较大变化；公司的经济发展对投资无较大影响资产投资是在市场中进行的，市场是复杂多变的，是无法用数量或函数进行准确描述的因此以上的假设是必要的，一般说来物价变化具有一定的周期性，社会政策也并非天天改变，公司自身的发展在稳定的情况下才会用额外的资金进行较大的风险的投资， 市场与社会的系统发展在一个时期内是良性的、稳定的，以上假设也是合理的。

3.1模型a
假设投资的风险水平是k,即要求总风险Q（x）限制在k内，Q（x），则模型可转化为：
max
s.t
3.2模型b
假设投资的收益水平是h，即净收益总额不少于 h：≥h，则模型可转化为：

s.t
3.3模型c
假设投资者对风险和收益的相对偏好参数为ρ（≥0），则模型可转化为：

s.t.
3.4 模型求解及分析
由于交易费 ci(xi)是分段函数，使得上述模型中的目标函数或约束条件相对比较复杂，是一个非线性规划问题，难于求解. 但注意到总投资额 M 相当大，一旦投资资产 Si，其投资额 xi 一般都会超过 ui，于是交易费 ci(xi)可简化为线性函数

从而，资金约束简化为

净收益总额简化为

在实际进行计算时，可设 M=1，此时
=（）（i＝0，1，…，n）
可视作投资 Si 的比例.
以下的模型求解都是在上述两个简化条件下进行讨论的.
1）模型 a 的求解
模型 a 的约束条件 Q(x)≤k 即
≤k，
所以此约束条件可转化为
　　（i＝0，1，…，n）.
这时模型 a 可化简为如下的线性规划问题：

具体到 n=4 的情形，按投资的收益和风险问题中题中给定的数据，模型为：

s.t
（i＝0，1，…，4）
利用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 取不同的值（0~0.025），风险与收益的关系见图1.
输出结果列表如下：
表1 模型1的计算结果
风险水平k
最大收益

0
0.05
1
0
0
0
0
0.001
0.0755
0.8316
0.04
0.0667
0.0182
0.0385
0.002
0.1011
0.6633
0.08
0.1333
0.0364
0.0769
0.003
0.1266
0.4949
0.12
0.2
0.0545
0.1154
0.004
0.1521
0.3266
0.16
0.2667
0.0727
0.1538
0.005
0.1776
0.1582
0.2
0.3333
0.0909
0.1923
0.006
0.2019
0
0.24
0.4
0.1091
0.2212
0.007
0.2066
0
0.28
0.4667
0.1273
0.1016
0.008
0.2112
0
0.32
0.5333
0.1271
0
0.009
0.2155
0
0.36
0.6
0.0233
0
0.01
0.219
0
0.4
0.5843
0
0
0.011
0.2223
0
0.44
0.5447
0
0

0.012
0.2256
0
0.48
0.5051
0
0
0.013
0.2288
0
0.52
0.4655
0
0
0.014
0.2321
0
0.56
0.4259
0
0
0.015
0.2354
0
0.6
0.3863
0
0
0.016
0.2387
0
0.64
0.3467
0
0
0.017
0.2419
0
0.68
0.3071
0
0
0.018
0.2452
0
0.72
0.2675
0
0
0.019
0.2485
0
0.76
0.2278
0

0
0.02
0.2518
0
0.8
0.1882
0
0
0.021
0.255
0
0.84
0.1486
0
0
0.022
0.2583
0
0.88
0.109
0
0
0.023
0.2616
0
0.92
0.0694
0
0
0.024
0.2649
0
0.96
0.0298
0
0
0.025
0.2673
0
0.9901
0
0
0

图1 模型1中风险k与收益的关系
结合图1，对于风险和收益没有特殊偏好的投资者来说，应该选择图中曲线的拐点（0.006，0.2019），这时对的投资比例见表1的黑体所示。
从表1中的计算结果可以看出，对低风险水平，除了存入银行外，投资首选风险率最低的 S2，然后是 S1 和 S4，总收益较低；对高风险水平，总收益较高，投资方向是选择净收益率（ri–pi）较大的 S1 和 S2．这些与人们的经验是一致的，这里给出了定量的结果．
2）模型 b 的求解
模型 b 本来是极小极大规划：

s.t. ≥h　　　　x≥0
但是，可以引进变量 xn+1=，将它改写为如下的线性规划：

s.t ,i=0,1,2,…,n,　≥h,　,　x≥0
具体到 n=4 的情形，按投资的收益和风险问题中题中给定的数据，模型为：
min x5
s.t

（i＝0，1，…，5）
利用 matlab7.1 求解模型 b，当 h 取不同的值（0.1~0.25），我们计算最小风险和最优决策,收益水平h取，结果如表2所示，风险和收益的关系见图2.
从表2看出，对低收益水平，除了存入银行外，投资首选风险率最低的资产，然后是和，总收益当然较低。对高收益水平，总风险自然也高，应首选净收益率（）最大的和。这些与人们的经验是一致的。
表2 模型2的计算结果
风险水平k
最大收益

0.002
0.1
0.6702
0.0783
0.1306
0.0356
0.0753
0.0024
0.11
0.6043
0.094
0.1567
0.0427
0.0904
0.0027
0.12
0.5383
0.1097
0.1828
0.0499
0.1055
0.0031
0.13
0.4724
0.1254
0.2089
0.057
0.1205
0.0035
0.14
0.4064
0.141
0.235
0.0641
0.1356
0.0039
0.15
0.3405
0.1567
0.2612
0.0712
0.1507
0.0043
0.16
0.2745
0.1724
0.2873
0.0783
0.1657
0.0047
0.17
0.2086
0.188
0.3134
0.0855
0.1808
0.0051
0.18
0.1426
0.2037
0.3395
0.0926
0.1959
0.0055
0.19
0.0767
0.2194
0.3656
0.0997
0.2109
0.0059
0.2
0.0107
0.235
0.3917
0.1068
0.226
0.0078
0.21
0
0.3114
0.519
0.0569
0.0908
0.0103
0.22
0
0.412
0.5725
0
0
0.0134
0.23
0
0.5341
0.4515
0
0
0.0164
0.24
0
0.6563
0.3305
0
0
0.0195
0.25
0
0.7784
0.2096
0
0

图2 模型2中风险与收益h的关系
结合图2，对于风险和收益没有特殊偏好的投资者来说，应该选择图中曲线的拐点（0.059，0.2），这时对的投资比例见表2的黑体所示。

3）模型 c 的求解
类似模型 b 的求解，我们同样引进变量 xn+1=，将它改写为如下的线性规划：
min ρxn+1–（1–ρ）
s.t ,i=0，1，2，…，n　　　　x≥0
具体到 n=4 的情形，按投资的收益和风险问题题中给定的数据，模型为：

s.t
（i＝0，1，…，5）
利用 matlab7.1 求解模型 c，当 ρ 取不同的值（0.75~0.95），我们计算最小风险和最优决策
输出结果列表如下：
表3 模型3的计算结果
风险水平
最大收益率

0.76
0.0248
0.2673
0
0.9901
0
0
0
0.77
0.0248
0.2673
0
0.9901
0
0
0
0.78
0.0092
0.2165
0
0.3693
0.6147
0
0
0.79
0.0092
0.2165
0
0.3692
0.6148
0
0
0.8
0.0079
0.2107
0
0.3151
0.5218
0.1431
0
0.81
0.0092
0.2165
0
0.3687
0.6154
0
0
0.82
0.0074
0.2084
0
0.2957
0.4928
0.1344
0.0547
0.83
0.0079
0.2106
0
0.3139
0.5235
0.1425
0
0.84
0.006
0.2017
0
0.2386
0.3973
0.1082
0.226
0.85
0.006
0.2016
0
0.2375
0.3967
0.1084
0.2275
0.86
0.0062
0.2026
0
0.2468
0.4113
0.0753
0.2372
0.87
0.0066
0.2045
0
0.2653
0.4422
0.0088
0.2552
0.88
0.0059
0.2017
0
0.2379
0.3964
0.1078
0.2279
0.89
0.006
0.2017
0
0.239
0.3937
0.1079
0.2294
0.9
0.006

0.2018
0
0.2392
0.398
0.1029
0.2299
0.91
0.0059
0.2017
0
0.238
0.3958
0.1075
0.2287
0.92
0.0059
0.2016
0
0.2376
0.3961
0.108
0.2284
0.93
0.0059
0.2016
0
0.2375
0.3963
0.108
0.2282
0.94
0.0061
0.1847
0.135
0.2436
0.3309
0.1106
0.1557
0.95
0.0059
0.2016
0
0.2376
0.396
0.108
0.2285
0.96
0.0017
0.0916
0.7256
0.0652
0.1065
0.0302
0.0641
0.97
0.0054
0.1848
0.114
0.2141
0.3564
0.0974
0.192

从图5可以看出，模型3的风险与收益关系与模型1和模型2的结果几乎完全一致。

图3 模型3中风险与收益的关系

图4模型3中风险与偏好系数的关系

图5 模型3中收益与偏好系数的关系
四 模型评价与推广
本文我们建立了投资收益与风险的双目标优化模型，通过控制风险使收益最大，保证收益使风险最小，以及引入收益——风险偏好系数，将两目标模型化为了单目标模型，并使用matlab7.1求解，所得结果具有一定的指导意义。
但是，本文没有讨论收益和风险的评估方法，在实际应用中还存在资产相关的情形，此时，用最大风险代表组合投资的风险未必合理，因此，对不同风险度量下的最优投资组合进行比较研究是进一步的改进方向。
五 总结
历经两周的时间终于完成了这次课设，在这次实践课程中，我真的遇到了不少的问题，在同学，老师的帮助以及在图书馆和网站搜集资料，解决了所有遇到的问题。尤其在问题分析的过程中，是难度最大也是问题最多的环节，感觉总是把问题分析的不够全面透彻，经常顾及这个方面而忽视了另一方面，最后我请教了同学，终于完成了问题分析并且建立了模型。在完成这一环节后，接下来的任务都是我独立完成，也遇到了不少的困难，但都是较易解决的。通过这次实践，我确实学到了不少，学会了使用MATLAB，也知道了分析问题的方法。
六 参考文献
[1]MATLAB程序设计与实例应用。张铮等。北京：中国铁道出版社，2003.10
[2]运筹学—方法与应用。吴风平。南京：河海大学出版社，2000.12
[3]《数学模型及方法》。李火林主编。江西高校出版社，1997.10
[4]数学建模教育及竞赛。甘筱青主编。南昌：江西高校出版社。2004.6
[5] 萧树铁，面向21世纪课程教材大学数学数学实验，北京：高等教育出版社，1999.7.
[6]赫孝良，戴永红等编著，数学建模竞赛：赛题简析与论文点评，西安：西安交通大学出版社，2002.6.
[7]陈叔平，谭永基，一类投资组合问题的建模与分析，数学的实践与认识，（29）7：45-49，1999.
七 附录
function result=qiujie()
%data为表格数据
data=[28 2.5 1 103
21 1.5 2 198
23 5.5 4.5 52
25 2.6 6.5 40];
data1=[9.6 42 2.1 181
18.5 54 3.2 407
49.4 60 6.0 428
23.9 42 1.5 549
8.1 1.2 7.6 270
14 39 3.4 397
40.7 68 5.6 178
31.2 33.4 3.1 220
33.6 53.3 2.7 475
36.8 40 2.9 248
11.8 31 5.1 195
9 5.5 5.7 320
35 46 2.7 267
9.4 5.3 4.5 328
15 23 7.6 131];
data=[[5 0 0 0];data]./100;%增加存银行
r=data(:,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(result.*10000)./10000;
% plot(result(:,1),result(:,2))
% grid on
% xlabel('风险 ')
% ylabel('收益')

%模型三求解
result=[];
for s=0.76:0.01:0.97
result=[result;moxing3(r,q,p,s)];
end
result=round(result.*10000)./10000;
figure(1)
plot(result(:,2),result(:,3))
grid on
xlabel('风险 ')
ylabel('收益')
pause
figure(2)
plot(result(:,1),result(:,2))
grid on
xlabel('偏好系数 ')
ylabel('风险')
pause
figure(3)
plot(result(:,1),result(:,3))
grid on
xlabel('偏好系数 ')
ylabel('收益')

function result1=moxing1(r,q,p,a)
%线性规划模型
%r收益率，为列向量
%p交易费率，为列向量
%q风险率，为列向量
%a风险水平
f=(p-r)';%转为求极小
n=length(q);
I=eye(n);
for i=2:n
A(i-1,:)=q(i)*I(i,:);
end
b=a*ones(n-1,1);
Aeq=(1+p');
beq=1;
lb=zeros(n,1);
ub=[];
[x,fval,exitflag,output]=linprog(f,A,b,Aeq,beq,lb,ub);
result1=[a,-fval,x'];
end

function 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'];
end

function 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'];
end
end