数学建模投资问题.docx
《数学建模投资问题.docx》由会员分享，可在线阅读，更多相关《数学建模投资问题.docx（22页珍藏版）》请在冰点文库上搜索。
数学建模投资问题
投资的收益与风险问题
摘要
对市场上的多种风险资产和一种无风险资产（存银行）进行组合投资策略的设计需要考虑两个目标：
总体收益尽可能大和总体风险尽可能小，而这两个目标在一定意义上是对立的。
本文我们建立了投资收益与风险的双目标优化模型，并通过“最大化策略”，即控制风险使收益最大，将原模型简化为单目标的线性规划模型一；在保证一定收益水平下，以风险最小为目标，将原模型简化为了极小极大规划模型二；以及引入收益——风险偏好系数，将两目标加权，化原模型为单目标非线性模型模型三。
然后分别使用Matlab的内部函数linprog，fminmax，fmincon对不同的风险水平，收益水平，以及偏好系数求解三个模型。
关键词：
组合投资，两目标优化模型，风险偏好
2•问题重述与分析
3.市场上有”种资产（如股票、债券、，）.：
0丨.小供投资者选择，某公司有数额为匸的
一笔相当大的资金可用作一个时期的投资。
公司财务分析人员对这种资产进行了评估，估算出在
这一时期内购买•「的平均收益率为c，并预测出购买t的风险损失率为％。
考虑到投资越分散，总的风险越小，公司确定，当用这笔资金购买若干种资产时，总体风险可用所投资的：
中最大的一个风
险来度量。
购买」要付交易费，费率为;■.，并且当购买额不超过给定值•；..时，交易费按购买■；.计算（不买当然无须付费）。
另外，假定同期银行存款利率是：
，且既无交易费又无风险。
（•
1、已知"；时的相关数据如下:
资产Si
收益率7\（%）
风险率缶
（%）
交易费刊（％）
阀值％（元）
S1
28
2.5
1
103
21
1.5
2
198
S3
23
5.5
4.5
52
S,
25
2.6
6.5
40
试给该公司设计一种投资组合方案，即用给定的资金有选择地购买若干种资产或存银行生息,
使净收益尽可能大，而总体风险尽可能小。
2、试就一般情况对以上问题进行讨论，并利用以下数据进行计算。
资产$
收益率T\（%）
风险率彊
（%）
交易费p（（%）
阀值％（元）
Si
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
St
40.7
68
5.6
178
s8
31.2
33.4
3.1
220
Sg
33.6
53.3
2.7
475
Sio
36.8
40
2.9
248
5ii
11.8
31
5.1
195
9
5.5
5.7
320
Sl3
35
46
2.7
267
S14
9.4
5.3
4.5
328
15
23
7.6
131
本题需要我们设计一种投资组合方案，使收益尽可能大，而风险尽可能小。
并给出对应的盈亏数
据，以及一般情况的讨论。
这是一个优化问题，要决策的是每种资产的投资额，要达到目标包括两方面的要求：
净收益最大和总
风险最低，即本题是一个双优化的问题，一般情况下，这两个目标是矛盾的，因为净收益越大则风险也会随着增加，反之也是一样的，所以，我们很难或者不可能提出同时满足这两个目标的决策方案，我们只能做到的是：
在收益一定的情况下，使得风险最小的决策，或者在风险一定的情况下，使得净收益最大，或者在收益和风险按确定好的偏好比例的情况下设计出最好的决策方案，这样的
话，我们得到的不再是一个方案，而是一个方案的组合，简称组合方案。
设购买Si（i=0,1…….n;So表示存入银行，）的金额为Xi；所支付的交易费为Ci（Xi），则：
0x=o
Cj（Xj）二〒PiUi0..Xi:
:
:
Uii-1,2,…，n,c°（x°）=0
I
PiXiX-a
对Si投资的净收益为：
R匕）=*Xi-q（xj（i=0,1,,,n）
对Si投资的风险为：
QdXjnqjXi（i=0,1,,,n）,q°=0
对S投资所需资金（即购买金额Xi与所需的手续费Ci（Xi）之和）是
fi（Xi）=XiCi（Xi）（i=0,1,,,n）
投资方案用x=（X0,X1,,,Xn）表示，那么，
净收益总额为：
n
R（x）八R（x）
i=0
总风险为：
Q（X）=rmiinQi（xJ
所需资金为：
n
F（x）八fi（xi）
i=0
所以，总收益最大，总风险最小的双目标优化模型表示为：
Yq（x））
min」F（x）=M,x兰0>
xR（x）丿J
但是像这样的双目标模型用一般的方法很难求解出来的，所以经过分析把次模型转化为三种较简单的单目标模型。
3.假设与模型
假设该公司在这一时期内是一次性投资；除交易费和投资费用外再无其他的费用开支；在这一时期市场发展基本上是稳定的；外界因素对投资的资产无较大影响；无其他的人为干预；社会政策无较大变化；公司的经济发展对投资无较大影响资产投资是在市场中进行的，市场是复杂多变的，是无法用数量或函数进行准确描述的，因此以上的假设是必要的，一般说来物价变化具有一定的周期性,
社会政策也并非天天改变，公司自身的发展在稳定的情况下才会用额外的资金进行较大的风险的投资，市场与社会的系统发展在一个时期内是良性的、稳定的，以上假设也是合理的。
3.1模型a
假设投资的风险水平是k,即要求总风险Q（x）限制在k内，Q（x）乞k，则模型可转化为:
maxRx
s.tQx乞k,F（x）二M,x_0
3.2模型b
假设投资的收益水平是h,即净收益总额R（x）不少于h：
R（x）>h,则模型可转化为:
minQ（x）
s.tR（x）_h,F（x）=M,x_0
3.3模型c
假设投资者对风险和收益的相对偏好参数为p（>0）,则模型可转化为:
ming（x）-（1-?
）R（x）
s.t.F（x）=M,x_0
3.4模型求解及分析
由于交易费Ci（xi）是分段函数，使得上述模型中的目标函数或约束条件相对比较复杂，是一个非线性
规划问题，难于求解.但注意到总投资额M相当大，一旦投资资产Si,其投资额Xi—般都会超过Ui,于是交易费Cig）可简化为线性函数
Ci（Xi）二PiXi
从而，资金约束简化为
nn
F（x）=迟fi（x）=迟（1+pJXj=M
i=0i=0
净收益总额简化为
nnn
R（x）=»R（Xi）=»[nx-G（Xi）]==（n-Pi）x
i=0i=0i=0
在实际进行计算时，可设M=1，此时
yi=（1Pi）Xi（i=0,1,,,n）
可视作投资Si的比例•
以下的模型求解都是在上述两个简化条件下进行讨论的
1）模型a的求解
模型a的约束条件Q（x）Q（x）=maxQi（Xj）=max（qixi）0丈翅0弍疋
所以此约束条件可转化为
qiK乞k（i=0,1,,,n）.
这时模型a可化简为如下的线性规划问题：
n
max'（斤_Pi）Xi
i国
s.t.qXj乞k,i=1,2,…，n
n
'、（1Pi）A",X一0
i=Q
具体到n=4的情形，按投资的收益和风险问题中题中给定的数据，模型为：
max0.05x00.27X!
0.19x20.185x30.185x4
s.t0.025x^k,0.015x2込k,0.055x3三k,0.026x4込k
X。
1.01x11.02x21.045x31.065x4=1,Xi_0（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
最大收益
易
S\
员
鬲
54
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
风险a
图1模型1中风险k与收益的关系
结合图1,对于风险和收益没有特殊偏好的投资者来说，应该选择图中曲线的拐点（0.006,0.2019）,
这时对：
的投资比例见表1的黑体所示。
从表1中的计算结果可以看出，对低风险水平，除了存入银行外，投资首选风险率最低的然后
是S1和S4，总收益较低；对高风险水平，总收益较高，投资方向是选择净收益率（ri-pj较大的
S1和&•这些与人们的经验是一致的，这里给出了定量的结果.
2）模型b的求解
模型b本来是极小极大规划：
minmaX（qiXi）
nn
s.t.'（r—pJXj》h'（1pix）=1x》0
i=0i=0
但是，可以引进变量Xn+1=max（qixi），将它改写为如下的线性规划:
mi门（人1）
s.tqiXi乞Xn1,i=0,1,2,,,n,\（A-pJx》h,（1pjx—l,x>0
i-0i」
具体到n=4的情形，按投资的收益和风险问题中题中给定的数据，模型为：
minX5
s.t0.025xi_X5,0.015x2-Xs,0.055x3_X5Q.O26X4_X5
0.05x00.27x10.19x20.185x30.185x4_h,
x01.01x11.02x21.045x31.065x4=1必-0,（i=0,1,,,5）
利用matlab7.1求解模型b，当h取不同的值（0.1~0.25）,我们计算最小风险和最优决策，收益水
平h取-U.2-；,结果如表2所示，风险和收益的关系见图2.
从表2看出，对低收益水平，除了存入银行外，投资首选风险率最低的资产然后是二和二，总
收益当然较低。
对高收益水平，总风险自然也高，应首选净收益率（-）最大的：
■［和二。
这些
与人们的经验是一致的。
表2模型2的计算结果
风险水平k
最大收益
So
Si
S2
£
£
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的关系
0.2）,
结合图2，对于风险和收益没有特殊偏好的投资者来说，应该选择图中曲线的拐点（0.059,
这时对，的投资比例见表2的黑体所示。
3）模型c的求解
类似模型b的求解，我们同样引进变量
将它改写为如下的线性规划:
n
minpx+i-（1-p）Z（斤一pJXi
i=0
n
s.tqiXi/ni,i=0,1,2,,,nx（1Pix）=1x>0
i=0
具体到n=4的情形，按投资的收益和风险问题题中给定的数据，模型为：
min'x5-（1-;?
）（0.05x00.27x10.19x20.185x30.185x4）
s.t0.025为兰x5,0.015x2兰x5,0.055x3兰x5,0.026x4兰x5
x01.01x11.02x21.045x31.065x4=1,x^0（i=0,1,,,5）
利用matlab7.1求解模型c,当p取不同的值（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.
七附录
functionresult=qiujie（）
%data为表格数据
data=[282.51103
211.52198
235.54.552
252.66.540];data1=[9.6422.1181
18.5543.2407
49.4606.0428
23.9421.5549
8.11.27.6270
14393.4397
40.7685.6178
31.233.43.1220
33.653.32.7475
36.8402.9248
11.8315.1195
95.55.7320
35462.7267
9.45.34.5328
15237.6131];
data=[[5000];data]./100;%增加存银行r=data（:
1）;
q=data（:
2）;
p=data（:
3）;%%模型一求解%result=[];
%fora=0:
0.01:
0.5
%result=[result;moxing1（r,q,p,a）];%end
%result=round（result.*10000）./10000;%plot（result（:
1）,result（:
2））
%gridon
%xlabel（'风险）
%ylabel（'收益'）
%%模型二求解%result=[];
%fork=0.1:
0.01:
0.4
%result=[result;moxing2（r,q,p,k）];%end
%result=round（result.*10000）./1000
升级会员 