matlab-建模-1998年A题

1风险投资组合的线性规划模型摘要(不要太简单)对市场上的多种风险资产和一种无风险资产（存银行）进行组合投资策略的设计需要考虑两个目标：总体收益尽可能大和总体风险尽可能小，而这两个目标在一定意义上是对立的。本文给出组合投资方案设计的一个线性规划模型。主要思路是通过线性加权综合两个设计目标；假设在投资规模相当大的基础上，将交易费函数近似线性化；通过决策变量的选取化解风险函数的非线性。模型的优点是：计算过程稳定性好，速度快。我们对各种加权因子，求得了最优化决策方案，从而得到问题的有效投资曲线。根据有效投资曲线，投资者可以由自己的主观偏好，直观地选择自己的投资方向。最后通过非线性规划，说明线性规划的结果对于交易费收取的阈值有一定的容忍度。2风险投资组合的线性规划模型一.问题的提出(不要完全重复原题)在风险市场的投资问题中，风险与收益始终是一对矛盾。一般来说想要追求高收益，风险也大;若想风险小，收益也会相应减少。研究表明，大部分的投资者具有以下的行为偏好：对于收益来说，总是越多越好；从风险的角度来说，大部分人都属于风险回避者。我们可以通过选取适当的组合投资方案，在取得良好收益的同时使总体风险减少。设某公司有一笔数额相当大的资金，投资购买若干种风险资产或存银行生息。风险资产收益高但风险大，存银行生息无风险但收益低。公司财务人员对多种资产进行了评估，估算出在这一时期内各种资产的平均收益率和风险损失率，并考虑购买时需付一定的交易费（不买当然无须付费，购买额不超过阈值时，交易费按阈值计算）。现在需要设计一种投资组合方案，以利用好这笔资金使得净收益尽可能大，而总风险尽可能小。二.模型的基本假设及符号说明(一)基本假设H1:只考虑给定时间内的收益和风险，且银行存款利率在给定时间内保持不变;H2:公司用于投资的资金数额相当大，且无贷款或透支;H3:各种资产投资风险相互独立。H4:总体风险可用所投资的资产中最大的一个风险来度量。(二)符号说明Si:第i种资产(i=1，2，...，n，n+1)，其中Sn+1表示存入银行;ri:Si的平均收益率;qi:Si的风险损失率;pi:Si的交易费率;ui:Si购买额阈值;M:资金总额;Xi:投资Si占总额的比重(不含交易费)，以下简称投资;(用比重较间洁！)Yi:投资Si的交易费占总额的比重，以下简称交易费;f1:净收益;f2:总体风险;:权因子;三.模型的建立(一)基本模型我们的目标是对各种资产投资以后，不仅收益尽可能大，同时总体风险还要尽可能小。所以我们的目标函数应为收益和风险两个函数。且由于在一段时间内的各种资产的平均收益3率和风险损失率均已由财务人员分析了出来，因此我们可以建立以下数学模型：（主要模型一定要用数学公式来写！）目标1maxf1=11niiiiYXr目标2minf2=iiniXq1max（风险的定义忠实于原题）s.t.11niiiYX=1其中MuXpXMuXpMuXYiiiiiiiiii000这是一个多目标非线性数学规划模型，且f1不是xi的连续函数，优化求解困难。下面我们将它转化为一个线性规划模型。(二)线性规划模型1目标函数的确定多目标规划有多种方法化为单目标问题解决。我们使用线性加权法。总目标函数minf=f2+(1-)(-f1)反映了风险投资中投资者的主观因素，越小表示投资越冒险。特别地=0表示只顾收益不顾风险，这样的人有可能取得最大收益;=1表示只顾风险而不顾收益，这样的人会将所有资金存入银行。2交易费函数的线性化近似本题难点之一是Yi不是Xi的连续函数。现将Yi近似为Xi的线性函数。Yi=piXi对阈值以下有一定误差(见下图)。但当投资规模充分大时，对优化结果不会有明显影响。一方面，对于i=1，2，...，n，若Si的投资很小，会白白浪费交易费，对优化不利，最优解一般不会出现小Xi;另一方面当投资总额很大时，不足购买费阈值的追加费用对目标函数影响不大。（这部分解释很必要！）YiuXi43风险函数的转化（本文最精彩之处！）令Xn+2=f2，那么必有qiXiXn+2(i=1，2，...n)。由于目标函数优化f，从而最优解必可使iiniXq1max达到Xn+2。这样得到线性规划模型:（主要模型一定要用数学公式来写！）minf=(1-)11niiiiXrpXn+2s.t.2,...2,10,...2,1011211niXniXXqXpiniiniii四.模型的求解(一)求解方法本文采用MATLAB优化工具箱中的线性规划函数lp求解。它优化下列线性规划模型:minCTXs.t.AX或b其使用格式为:X=lp(C，A，b，vlb，vub，X0，N)(现在版本用linprog)其中vlb，vub分别是上下界，X0为初始值，N表示约束条件中前N个约束为等式约束。计算步骤1输入数据，选取权因子;2生成矩阵C，A，b;3根据需要取vlb，vub，X0，N(本问题vlb取零向量，N取1，vub和X0无特殊要求，置为空集)4使用MATLAB函数lp求解;（二）计算结果及分析:(用表格\图形表达很清晰,图表要有编号和标题,文中引用要清楚,主要结果写出来,长表格\图形放在附录中)1.数值计算结果对于附录表1(问题一)和附录表2(问题二)所示数据,使用上述方法分别求解当=0.1，50.2，...，1时的最优决策及风险和收益如附录表3和附录表4(M函数dxz1_1.m和M函数dxz1_2.m)2投资方案分析（1）.从附录表3和附录表4得到问题一的四个典型最优组合，问题二有7个典型最优组合。对于不同风险承受能力，选择该风险水平下的最优投资组合。例如：对问题一，若风险承受水平是0.02，那么取=0.2时的决策方案.（2）净收益和风险都是参数的单调下降函数（见图1）.说明谨慎程度越强，风险越小但受益也越小.具有明确的实际意义。图1收益和风险随参数增大而下降（3）.更详细的计算结果见图2（使用M文件dxz1_1p.m和dxz1_2p.m）。我们用=01内300等分点，求得最优投资组合集及它们形成的有效投资曲线。这条曲线上的任一点都表示该风险水平的最大可能收益和该收益要求的最小风险。实际上在我们计算精度内，问题一只有5个最优方案，问题二只有13个最优方案。其中问题一风险0.0059(即=0.9)的决策(0.2376,0.3960,0.1080,0.2284,0)和问题二风险0.0995(即=0.3)的决策(0,0,0.1658,0,0,0,0.0.1463,0,0.1867,0.2487,0,0,0.2163,0,0,0)具有特别重要的意义，因为它们对应在风险增长较慢情形下最大的收益，可认为是一般意义上的最优解。图2有效投资曲线00.5100.050.10.150.20.250.30.35投投投投投投(n=4)投投投投投投投投00.5100.10.20.30.40.50.60.7投投投投投投(n=15)投投投投投投投投00.010.020.030.050.10.150.20.250.3投投投投投投(n=4)投投投投00.510.050.10.150.20.250.30.350.40.45投投投投投投(n=15)投投投投63适用性分析当0XiMui，线性规划模型可能不是最优解。也就是结果的正确性与M有关，M(iixu)时必最优.但M(iixu)时结果不一定可靠.比如对问题二，当=0.3时，这个临界值为2581；对问题一，=0.9时，这个临界值为500.五.模型的验证与改进(一)非线性规划模型当投资规模较小时，线性规划模型中交易费函数的线性化近似可能造成吸引小投资(低于阈值)，从而使优化结果失去最优性.但由于交易费函数的不连续性，会造成非线性优化算法的不稳定，将交易函数作如下处理YxpXuMupuMxuMXpXuMiiiiiiiiiiiiii100001001..Yi1%Xi可使Yi变成连续函数，这样做的前提假设是H5:极少的投资是不可能出现的(所有投资要超过阈值的1%)，从而得到下列非线性规划模型:7minXn+2-(1-)()rXYiiiin11s.t.2...,2,10...,2,101211niXniXXqYXiniiniii(二)计算结果分析我们使用MATLAB的最优化函数Constr计算(现在版本用linprog)(M文件dxz2_1.m和dxz2_2.m，它的调用函数dxz2f.m)，计算结果发现在线性规划模型M的临界值上方，结果总与线性规划一致.即使对临界值下方线性规划的结果还有相当大的容度(如问题一=0.9,M=100,问题二=0.3,M=2000).对于相当小的M，非线性规划的计算结果出现不稳定，表现出对初值的依赖性.这时，线性规划结果为次优解，而非线性规划结果为局部最优结果.参考文献[1]运筹学，清华大学出版社，北京，1990.[2]赵锡军等，金融投资学，中国人民大学出版社，北京，1996[3]施阳，MATLAB语言工具箱，西北工业大学出版社，西安，19988附录一:数值计算结果表1问题一Siri(%)qi(%)pi(%)ui（元）S1282.51103S2211.52198S3235.54.552S4252.66.540表2问题2Siri(%)qi(%)pi(%)ui(元)S19.6422.1181S218.5543.2407S349.4606.0428S423.9421.5549S58.11.27.6270S614393.4397S740.7685.6178S831.233.43.1220S933.653.32.7475S1036.8402.9248S1111.8315.1195S1295.55.7320S1335462.7267S149.45.34.5328S1515237.6131表3问题1投资组合Si=00.7=0.8=0.9=1S10.99010.36900.23760.0000S20.00000.61500.39600.0000S30.00000.00000.10800S4000.22840.0000存银行0.00000.00000.00001.0000净收益0.26730.21650.20160.0500风险0.02480.00920.00590.00009表4问题2投资组合Si=00.1=0.2=0.3=0.40.5=0.60.7=0.8=0.91S10.00000.00000.00000.00000.00000.08740.0000S20.00000.00000.00000.00000.00000.06800.0000S30.94340.20510.16580.12690.10710.06120.0000S40.00000.00000.00000.00000.15310.08740.0000S50.00000.000000.00000.00000.00000.0000S60.00000.00000.00000.00000.00000.09420.0000S70.00000.18100.14630.11190.09450.05400.0000S80.00000.00000.00000.22790.19250.11000.0000S90.000000.18670.14280.12060.06890.0000S100.00000.30770.24870.19030.16070.09180.0000S110.00000.00000.00000.00000.00000.00000S120.00000.00000.00000.00000.00000.00000.0000S130.00000.26760.21630.16550.13980.07980.0000S140.00000.00