A-组合投资收益和风险问题的数学建模-数学建模课程设计报告

1、.XX工程大学  数学建模  课程设计论文组合投资收益和风险问题的数学建模姓 名：何健堃 学 号：3110801216 班 级：数学112 指导教师：周金明 成 绩：完成日期：2021年7月3日摘要在现代商业、金融投资中，投资者总是希望实现收益最大化，然而投资是要承当风险的，收益与风险之间存在难以调和的矛盾，怎样兼顾两者，寻找切实可行的决策思想，是投资的收益和风险决策的一个重要的问题。本论文主要讨论并解决了在组合投资问题中的投资收益与风险的有关问题。分别在不考虑投资工程之间的影响和考虑投资工程之间的影响以及不考虑风险和考虑风险的情况下，建立相应的数学模型

2、，来使得投资所获得的总利润到达最大。问题一是一典型的线性规划问题。根据题目要求，要求第五年末的最大利润，那么建立线性规划模型，在LINDO中编程求得第五年末的最大利润为1418.704万元。第一年投资工程有1，2，3，4，5，6，投资额分别为50000.00，30000.00，40000.00，30000.00，30000.00，20000.00万元；第二年投资工程有1，2，7，投资额分别为10083.00，30000.00，40000.00万元；第三年投资工程有1，2，8，投资额分别为50307.08，30000.00，30000.00万元；第四年投资工程有1，2，3，4，投资额分别为306

3、25.39，30000.00，40000.00，30000.00万元；第五年投资工程有1，2，投资额分别为30689.01，30000.00万元。问题二是在问题一的根底上，增加了约束条件考虑工程间的影响的组合投资问题。建立非线性规划模型，在LINGO中求解得到第五年末的最大利润为231762.8万元。第一年投资工程有1，2，3，4，投资额分别为60000.00，60000.00，35000.00，30000.00万元；第二年投资工程有1，2，5，投资额分别为60000.00，60000.00，30000.00万元；第三年投资工程有1，2，6，投资额分别为60000.00，60000.00，40

4、000.00万元；第四年投资工程有1，2，3，4，投资额分别为60000.00，60000.00，35000.00，30000.00万元；第五年投资工程有1，2，投资额分别为60000.00，60000.00万元。问题三在问题二的根底上考虑投资风险，要求风险最小，收益最大，是一双目标函数问题，求解过程中假设投资者能承受的最大风险度0.01a=，进而将其转化为单目标问题。在LINGO中求解得到第五年末的最大利润为175244.8万元。第一年投资工程有1，2，3，4，6投资额分别为60000.00，60000.00，11670.25，30000.00，7784.336万元；第二年投资工程有1，2，

5、3，5，6投资额分别为57395.06，54729.02，12711.78，2500.000，1342.021万元；第三年投资工程有1，2，3，4，6，投资额分别为60000.00，60000.00，13511.02，30000.00，9012.173万元；第四年投资工程有1，2，3，投资额分别为60000.00，60000.00，17037.64万元；第五年投资工程有1，2，投资额分别为60000.00，60000.00万元。关键词：线性规划 非线性规划 LINGO 风险1、问题重述 由题意知，现有一笔20亿的资金作为未来5年的投资资金，不是将其投资于一个工程上，而是有8个工程可供选择投资，

6、并要求：工程1、工程2每年初投资，当年年末回收本利本金和利润；工程3、工程4每年初投资，要到第二年末才可回收本利；工程5、工程6每年初投资，要到第三年末才可回收本利；工程7只能在第二年年初投资，到第五年末回收本利；工程8只能在第三年年初投资，到第五年末回收本利。1.公司财务分析人员的实验数据，求怎样安排投资，使得第五年末的利润最大，也即是五年累计的利润最大.2.考虑各工程投资的相互影响的同时又考虑了其他的因素：工程5的投资额固定，为500万，可重复投资；对投资工程1，公司管理层争取到一笔资金捐赠，假设在工程1中投资超过20000万，那么同时可获得该笔投资金额的1%的捐赠，用于当年对各工程的投资

7、。每年对各个工程投资的到期利润率各不一样且受到某些工程同时投资的影响，如同时对工程3和工程4投资时的到期利润率与工程3和工程4独立投资不同。所以应在模型一中加以更多的限制来求解第五年末所得最大利润。3.此题考虑投资风险，并用这笔资金投资假设干种工程时，总体风险可用所投资的工程中最大的一个风险来度量，来求如何投资，使得总体风险最小，第五年末的利润最大。其它约束条件与二问一样。2、条件假设1.假设市场经济在五年内根本稳定。2. 无交易费和投资费用等的费用开支。3. 假设在五年内产生的利润可作为下一年的投资资金。4.对问题一： 1假设不考虑各投资工程间的影响以及存在的风险。2假设在五年内各投资工程的

8、到期利润率不变。5.对问题二： 1假设不考虑投资风险。 2假设仍然考虑一问条件，如工程3、工程4每年初投资，要到第二年末才可回收本利。6. 对问题三： 1某些工程单独投资或共同投资的风险损失率标准差每年根本一样。 2假设仍然考虑一问条件，如工程3、工程4每年初投资，要到第二年末才可回收本利。 3投资者能承受的最大风险度。3、符号说明 下面将会用到以下符号： -第五年末所得的利润。 -第年末回收的本利。-每年年初的可用来投资的资金。-第年初对第个工程的投资金额。-问题一中第个工程的预计到期利润率。-第个工程的投资上限。-问题二中各工程单独投资时，第年第个工程的到期利润率。-问题二中某些两个工程同

9、时投资时，第年第个工程的到期利润率。-问题二中五六八工程同时投资时的到期利润率。-问题二中实际采用的第年第个工程的到期利润率。-问题二中第一工程获得的捐赠资金。 -问题二中第年工程五可以重复投资的次数。 -问题三中的总体投资风险。 -问题三中实际采用的风险损失率。-问题三中第个工程单独投资时的风险损失率。-问题三中某些工程同时投资时的风险损失率。 -问题三中投资者能够承受的最大风险度。4、问题分析 题目背景：随着市场经济的开展，越来越多的人把目光转移到投资上。为了获得最大的利益，组合投资的风险与收益成为投资者必须考虑的问题。投资者在考虑投资工程的预期收益率的同时，也要考虑投资工程存在的风险。如

10、何使利润最大化，风险最小化是投资者首要考虑的问题。下面是对该问题的分析与解答。在分析解答之前给出以下定义：预期收益率：预期收益率也称为期望收益率，是指如果事件不发生的话可以预计到的收益率。此题中到期收益率与之类似，即到期利润率是指对某工程的一次投资中，到期回收利润与本金的比值。风险：风险可以用收益的方差或标准差来进展衡量：方差越大，那么认为风险越大；方差越小，那么认为风险越小。在一定的假设下，有收益的方差或标准差来衡量风险确实是适宜的。风险损失率和风险一样也可以用收益的方差或标准差来进展衡量。1.题一，根据题意，在不考虑各项投资之间的相互影响以及投资存在的风险情况下，求第五年末的最大利润，即五

11、年内利润的和的最大值。那么第五年末的最大利润就等于五年中各项投资工程的金额与该工程的到期利润率的乘积之和。还应考虑的限制条件为：限制条件1：每年可用来投资的最大金额。第一年为20亿元。限制条件2：对每项工程投资的最大金额。2.题二, 根据题意, 在考虑各项投资之间的相互影响且不考虑在投资中的风险情况下，求第五年末的最大利润，即五年内利润的和的最大值。 首先，应根据附表数据求得各年对各工程的单独投资的到期的利润率，以及各年对某些工程同时投资的到期利润率，如同时对工程5、工程6和工程8投资的到期利润率与它们单独投资时的到期利润率不同，应分情况考虑，即它们是否同时投资。然后，对投资工程1应单独考虑，

12、因为当工程一投资大于20000万时，可以获得1%的捐赠，这是我们必须考虑的。其次，工程5的投资额固定，为500万，可重复投资。所以应考虑工程5在五年内的投资次数。限制条件1：每年可用来投资的最大金额。限制条件2：对每项工程投资的最大金额。限制条件3：同时投资时的到期利润率不同。 3.题三,根据题意考虑各投资工程的风险情况下，求第五年末的最大利润，即五年内利润的和的最大值，并且使得风险最小。此问题为一双目标规划问题。 首先，应该把总的风险表示出来，即总体风险可用所投资的工程中最大的一个风险来度量。风险可以用收益的方差或标准差来衡量。其次，在考虑风险的根底上来求解如何组合投资使得风险最小而利润最大

13、。限制条件1：每年可用来投资的最大金额。限制条件2：对每项工程投资的最大金额。限制条件3：同时投资时的到期利润率不同。限制条件4：同时投资与非同时投资时的风险损失率标准差不同。5、模型建立问题一：不考虑投资风险和各投资工程间的相互影响，建立线形规划模型求解。 1目标函数：第五年末所得利润，即五年内各工程产生的总利润 2约束条件1：每年年初可用来投资的金额应大于等于每年年初投资各工程的总额：第年年初可用来投资的资金等于第年初未投资的资金加上第年末回收的本利。 第一年初可用来投资的资金第一年末回收第一年投资中一二工程的本利 第二年初可用来投资的资金 第二年末回收第一年三四工程和第二年一二工程的本利

14、 第三年初可用来投资的资金 第三年末回收第一年五六工程、第二年三四工程和第三年一二工程的本利 第四年初可用来投资的资金 第四年末回收第二年五六工程、第三年三四工程和第四年一二工程的本利 第五年初可用来投资的资金 第五年末回收第二年七工程、第三年五六八工程、第四年三四工程、第五年一二工程的本利 约束条件2：每项工程在投资期间不能超过工程最大金额的限制 第一二工程每年均可以投资，可以投五次： 第三四工程每年初投资，第二年末才可回收本利，可投资四次：第五六工程每年初投资，第三年末才可回收本利，可投资三次： 第七工程第二年初投资，到第五年末回收本利，可投资一次： 第八工程第三年初投资，到第五年末回收本

15、利，可投资一次：所以可以建立线形规划模型，利用LINGO软件或LINGDO软件求解。问题二：不考虑投资风险时，建立非线性规划模型求解。（1） 根据题目中表二、表三分别求得各投资工程独立投资时的投资额及到期利润率见表一以及一些投资工程同时投资时的投资额及到期利润见表二。工程12345678到期利润率20210.147160.152090.233340.294220.698950.7515318.44862.4939520210.184310.189620.251280.33170.737021.3867-4.27695-1.1981320210.10370.185510.255570.30819

16、-0.161361.49251-4.84479-1.3717420210.144490.210650.519210.259821.039881.299942.829555.1859820210.181710.245810.516380.321011.215511.30994-9.607562.67363表一. 各投资工程独立投资时的到期利润率工程同时投资工程3、4同时投资工程5、6同时投资工程5、6、83456568到期利润率20210.309610.44849-0.016682.1279-0.578280.783271.7063720210.731860.393910.97240.36539

17、1.015870.871912.6029520210.56240.436411.9492-0.65437-0.248452.06485-0.5438520210.263590.391071.006291.694410.880730.387382.3895420210.520860.382091.31418-1.084061.73688-0.74813-1.00638表二. 一些投资工程同时投资时到期利润率 2目标函数：第五年末所得利润，即五年内各工程产生的总利润工程单独投资时，某些工程同时投资时 3约束条件1：某些工程同时投资时的到期利润率。 1当万元时，第年可得捐赠资金，作为第年初的投资资金

18、，当万元时，第年可得捐赠资金。 2当三四工程同时投资时，此时， 否那么。 3同理：当五六工程同时投资时，此时，否那么。 4当第三年初五六八工程同时投资时，此时，否那么。5工程五为固定投资500万且可以重复投资，。 约束条件2：每年年初可用来投资的金额应大于等于每年年初投资各工程的总额： 第一年初可用来投资的资金第一年末回收第一年投资中一二工程的本利 第二年初可用来投资的资金 第二年末回收第一年三四工程和第二年一二工程的本利 第三年初可用来投资的资金 第三年末回收第一年五六工程、第二年三四工程和第三年一二工程的本利 第四年初可用来投资的资金 第四年末回收第二年五六工程、第三年三四工程和第四年一二

19、工程的本利 第五年初可用来投资的资金 第五年末回收第二年七工程、第三年五六八工程、第四年三四工程、第五年一二工程的本利 约束条件2：每项工程在投资期间不能超过工程最大金额的限制 第一二工程每年均可以投资，可以投五次： 第三四工程每年初投资，第二年末才可回收本利，可投资四次：第五六工程每年初投资，第三年末才可回收本利，可投资三次： 第七工程第二年初投资，到第五年末回收本利，可投资一次： 第八工程第三年初投资，到第五年末回收本利，可投资一次：所以可以建立非线性规划模型，利用LINGO软件或LINDO软件求解。 问题三：考虑投资风险的情况下，求第五年末的最大利润，并使风险最小。为此建立双目标规划，但

20、为计算简便，通过一定的假设条件将其转化为单目标规划问题。（1） 根据题目，我们用预计到期收益率的标准差来衡量风险损失率。根据各投资工程独立投资时的到期利润率表一和一些投资工程同时投资时到期利润率表二，利用EXCEL中的STDEV函数，求得各工程单独投资时的风险损失率表三以及某些工程同时投资的风险损失率表四：工程1工程2工程3工程4工程5工程6工程7工程8标准差0.0294580.0308930.1330060.0246860.4742150.2576979.6889092.508485表三. 各工程单独投资时的风险损失率工程工程3、4影响工程5、6影响工程5、6、8影响工程3工程4工程5工程6

21、工程5工程6工程8标准差0.1718900.0267370.6362341.2598490.8538140.9044511.510221表四.某些工程同时投资时的风险损失率在这里我们定义工程风险度=工程风险损失率*工程投资额/投资资金总额。那么可以得到目标函数：风险最小 利润最大为了简化模型，我们给定投资者能承受的最大风险，那么将双目标函数转化为单目标函数。也就是说：在满足投资者能承受的最大风险的情况下，如何使收益最大。（2） 目标函数：（3） 约束条件： 约束条件1：与问题二中约束条件完全一样，这里不再重复。 约束条件2：总体投资风险不能超过投资者能承受的最大风险 那么工程的风险度应小于投资

22、者能承受的最大风险。即6 模型求解问题一：根据上面的分析和建立的模型知： 模型化解:实际在此题中可以使计算简便。由以上分析并在LINDO中编程源程序见一得到第五年末所得利润的最大值为1418.704万元。各年的投资工程见下表单位：万元：第一年第二年第三年第四年第五年工程一50000.00000010083.00000050307.08203130625.39062530689.015625工程二30000.00000030000.00000030000.00000030000.00000030000.000000工程三40000.0000000040000.000000-工程四30000.00

23、00000030000.000000-工程五30000.000000000-工程六20000.000000000-工程七-40000.000000-工程八-30000.000000-由上表可知：第一年投资工程有1，2，3，4，5，6，投资额分别为50000.00，30000.00，40000.00，30000.00，30000.00，20000.00万元；第二年投资工程有1，2，7，投资额分别为10083.00，30000.00，40000.00万元；第三年投资工程有1，2，8，投资额分别为50307.08，30000.00，30000.00万元；第四年投资工程有1，2，3，4，投资额分别为3

24、0625.39，30000.00，40000.00，30000.00万元；第五年投资工程有1，2，投资额分别为30689.01，30000.00万元。问题二：根据上面的分析和建立的模型知：由以上分析并在LINGO中编程源程序见二得到第五年末所得利润的最大值为231762.8万元。各年的投资工程见下表单位：万元：第一年第二年第三年第四年第五年工程一60000.0060000.0060000.0060000.0060000.00工程二60000.0060000.0060000.0060000.0060000.00工程三35000.000035000.00-工程四30000.000030000.00

25、-工程五030000.000-工程六0040000.00-工程七-0-工程八-0-由上表可知：第一年投资工程有1，2，3，4，投资额分别为60000.00，60000.00，35000.00，30000.00万元；第二年投资工程有1，2，5，投资额分别为60000.00，60000.00， 30000.00万元；第三年投资工程有1，2，6，投资额分别为60000.00，60000.00，40000.00万元；第四年投资工程有1，2，3，4，投资额分别为60000.00，60000.00，35000.00，30000.00万元；第五年投资工程有1，2，投资额分别为60000.00，60000.0

26、0万元。问题三：根据上面的分析和建立的模型知：由以上分析并在LINGO中编程源程序见三得到第五年末所得利润的最大值为175244.8万元。各年的投资工程见下表单位：万元：第一年第二年第三年第四年第五年工程一60000.0057395.0660000.0060000.00 60000.00工程二60000.0054729.0260000.0060000.0060000.00工程三11670.2512711.7813511.0217037.64-工程四30000.00030000.000-工程五02500.0000-工程六7784.3361342.0219012.173-工程七-0-工程八-0-

27、由上表可知：第一年投资工程有1，2，3，4，6投资额分别为60000.00，60000.00，11670.25，30000.00，7784.336万元；第二年投资工程有1，2，3，5，6投资额分别为57395.06，54729.02，12711.78，2500.000，1342.021万元；第三年投资工程有1，2，3，4，6，投资额分别为60000.00，60000.00，13511.02，30000.00，9012.173万元；第四年投资工程有1，2，3，投资额分别为60000.00，60000.00，17037.64万元；第五年投资工程有1，2，投资额分别为60000.00，60000.0

28、0万元。7、模型检验与评价本文均采用规划模型实现了组合投资的收益和风险问题的优化投资方案。问题一：问题一是一典型的线性规划问题，在不考虑投资工程间的影响和投资风险时，不难得到目标函数和约束条件，在LINDO中很快得到全局最优解。问题二：问题二在问题一的根底上，把投资工程间的影响考虑了进来，增加了约束条件。由于各工程各年的收益率不同，增加了求解模型的难度。且将模型变得非线性化。在本模型中用LINGO的if选择语句实现其在组合投资中对个工程是否同时投资的选择，在一定程度上实现了非线性模型的建立,并借助第一问的约束条件，在LINGO中编程，求得局部最优解。问题三：问题三在以上两问的根底上，又把投资风

29、险考虑进来，要求在风险最小的情况下建立模型求得投资方案以获得最大利润，为双目标规划模型，但在求解上有一定的难度，于是在此转化为单目标规划模型。即假设投资者能承受的最大风险度，这样将双目标问题转化为与二问类似的问题，并且投资越分散说明投资风险越小。进而得到各年的投资方案，从投资方案中看出，该投资方案相比照一二问的投资工程分散，这说明该模型的假设还是比拟合理的。当改变假设的最大风险度时，第五年末的最大利润也发生相应的变化。时，最大利润；时，最大利润；时，最大利润；由上面的数据分析可知，投资者能承受的最大风险度越大，所得的利润就越大。但这还是有一定的局限性，对此还可以考虑以下方案：1在投资者要求的最

30、低利润条件下，来求风险的最小值。即给定一个最小利润，在以此为约束的条件下来求风险的最小值，这样也能在一定程度上将双目标函数优化为单目标函数。2采用加权法，对风险和投资利润加于权重因子构造出目标函数，再以一定的约束条件以及权重因子的约束为条件，用LINGO软件编程求解相应的投资方案。8、参考文献1韩中庚。数学建模方法及其应用。：高等教育，2006.2谢金星，薛毅。优化建模与LINDO/LINGO软件。：清华大学，2006.9、附录附源程序一：LINDO中求解max 0.001x11+0.001x21+0.001x31+0.001x41+0.001x51+0.0011x12+0.0011x22+0

31、.0011x32+0.0011x42+0.0011x52+0.0025x13+0.0025x23+0.0025x33+0.0025x43+0.0027x14+0.0027x24+0.0027x34+0.0027x44+0.0045x15+0.0045x25+0.0045x35+0.005x16+0.005x26+0.005x36+0.008x27+0.0055x38STx11<60000x21<60000x31<60000x41<60000x51<60000x12<30000x22<30000x32<30000x42<30000x52<

32、;30000x23+x13<40000x33+x23<40000x43+x33<40000x24+x14<30000x34+x24<30000x44+x34<30000x25+x15+x35<30000x26+x16+x36<20000x27<40000x38<30000x11+x12+x13+x14+x15+x16<200000x21+x22+x23+x24+x25+x26+x27-0.001x11-0.0011x12+x13+x14+x15+x16<200000x31+x32+x33+x34+x35+x36+x38-0

33、.001x21-0.0011x22-0.001x11-0.0011x12-0.0025x13-0.0027x14+x23+x24+x15+x16+x25+x26+x27<200000x41+x42+x43+x44-0.001x31-0.0011x32-0.001x21-0.0011x22-0.001x11-0.0011x12-0.0025x23-0.0027x24-0.0025x13-0.0027x14-0.0045x15-0.005x16+x33+x34+x25+x26+x35+x36+x27+x38<200000x51+x52-0.001x41-0.0011x42-0.001x

34、31-0.0011x32-0.001x21-0.0011x22-0.001x11-0.0011x12-0.0025x33-0.0027x34-0.0025x23-0.0027x24-0.0025x13-0.0027x14-0.0045x25-0.005x26-0.0045x15-0.005x16+x43+x44+x35+x36+x27+x38<200000END附源程序二：LINGO中求解max=0.14716*x11+0.15209*x12+b13*x13+b14*x14+b15*x15+b16*x16+ 0.18431*x21+0.18962*x22+b23*x23+b24*x24+

35、b25*x25+b26*x26+(-4.27695)*x27+ 0.1037*x31+0.18551*x32+b33*x33+b34*x34+b35*x35+b36*x36+b38*x38+ 0.14449*x41+0.21065*x42+b43*x43+b44*x44+ 0.18171*x51+0.24581*x52;x15=500*n1;x25=500*n2;x35=500*n3;b13=if(x13*x14*EQ*0,0.23334,0.30961);b14=if(x13*x14*EQ*0,0.29422,0.44849);b15=if(x15*x16*EQ*0,0.69895,-0.0

36、1668);b16=if(x15*x16*EQ*0,0.75153,2.1279);b23=if(x23*x24*EQ*0,0.25128,0.73186);b24=if(x23*x24*EQ*0,0.3317,0.39391);b25=if(x25*x26*EQ*0,0.73702,0.9724);b26=if(x25*x26*EQ*0,1.3867,0.36539);b33=if(x33*x34*EQ*0,0.25557,0.5624);b34=if(x33*x34*EQ*0,0.30819,0.43641);b35=if(x35*x36*EQ*0,-0.16136,if(x38*EQ*0

37、,1.9492,-0.24845);b36=if(x35*x36*EQ*0,1.49251,if(x38*EQ*0,-0.65437,2.06485);b38=if(x35*x36*x38*NE*0,-0.54385,-1.37174);b43=if(x43*x44*EQ*0,0.51921,0.26359);b44=if(x43*x44*EQ*0,0.25982,0.39107);x11<60000;x21<60000;x31<60000;x41<60000;x51<60000;x12<60000;x22<60000;x32<60000;x42

38、<60000;x52<60000;x23+x13<35000;x33+x23<35000;x43+x33<35000;x24+x14<30000;x34+x24<30000;x44+x34<30000;x25+x15+x35<30000;x26+x16+x36<40000;x27<30000;x38<30000;g11=if(x11*LT*20000,0,x11*0.01);g21=if(x21*LT*20000,0,x21*0.01);g31=if(x31*LT*20000,0,x31*0.01);g41=if(x41*

39、LT*20000,0,x41*0.01);g51=if(x51*LT*20000,0,x51*0.01);x11+x12+x13+x14+x15+x16-g11<200000;x21+x22+x23+x24+x25+x26+x27-0.14716*x11-0.15209*x12+x13+x14+x15+x16-g21<200000;x31+x32+x33+x34+x35+x36+x38-0.18431*x21-0.18962*x22-0.14716*x11-0.15209*x12-b13*x13-b14*x14+x23+x24+x15+x16+x25+x26+x27-g31<

40、200000;x41+x42+x43+x44-0.1037*x31-0.18551*x32-0.18431*x21-0.18962*x22-0.14716*x11-0.15209*x12-b23*x23-b24*x24-b13*x13-b14*x14-b15*x15-b16*x16+x33+x34+x25+x26+x35+x36+x27+x38-g41<200000;x51+x52-0.14449*x41-0.21065*x42-0.1037*x31-0.18551*x32-0.18431*x21-0.18962*x22-0.14716*x11-0.15209*x12-x33*b33-x

41、34*b34-b23*x23-b24*x24-b13*x13-b14*x14-x25*b25-x26*b26-x15*b15-x16*b16+x43+x44+x35+x36+x27+x38-g51<200000;gin(n1);gin(n2);gin(n3);free(b13);free(b23);free(b33);free(b43);free(b14);free(b24);free(b34);free(b44);free(b15);free(b25);free(b35);free(b16);free(b26);free(b36);free(b38);free(b27);end附源程序

42、三：LINGO中求解max=0.14716*x11+0.15209*x12+b13*x13+b14*x14+b15*x15+b16*x16+ 0.18431*x21+0.18962*x22+b23*x23+b24*x24+b25*x25+b26*x26+(-4.27695)*x27+ 0.1037*x31+0.18551*x32+b33*x33+b34*x34+b35*x35+b36*x36+b38*x38+ 0.14449*x41+0.21065*x42+b43*x43+b44*x44+ 0.18171*x51+0.24581*x52;y0=200000;x15=500*n1;x25=500

43、*n2;x35=500*n3;b13=if(x13*x14*EQ*0,0.23334,0.30961);b14=if(x13*x14*EQ*0,0.29422,0.44849);b15=if(x15*x16*EQ*0,0.69895,-0.01668);b16=if(x15*x16*EQ*0,0.75153,2.1279);b23=if(x23*x24*EQ*0,0.25128,0.73186);b24=if(x23*x24*EQ*0,0.3317,0.39391);b25=if(x25*x26*EQ*0,0.73702,0.9724);b26=if(x25*x26*EQ*0,1.3867,0

44、.36539);b33=if(x33*x34*EQ*0,0.25557,0.5624);b34=if(x33*x34*EQ*0,0.30819,0.43641);b35=if(x35*x36*EQ*0,-0.16136,if(x38*EQ*0,1.9492,-0.24845);b36=if(x35*x36*EQ*0,1.49251,if(x38*EQ*0,-0.65437,2.06485);b38=if(x35*x36*x38*NE*0,-0.54385,-1.37174);b43=if(x43*x44*EQ*0,0.51921,0.26359);b44=if(x43*x44*EQ*0,0.2

45、5982,0.39107);o13=if(x13*x14*EQ*0,0.133006,0.171890);o23=if(x23*x24*EQ*0,0.133006,0.171890);o33=if(x33*x34*EQ*0,0.133006,0.171890);o43=if(x43*x44*EQ*0,0.133006,0.171890);o14=if(x13*x14*EQ*0,0.024686,0.026737);o24=if(x23*x24*EQ*0,0.024686,0.026737);o34=if(x33*x34*EQ*0,0.024686,0.026737);o44=if(x43*x4

46、4*EQ*0,0.024686,0.026737);o15=if(x15*x16*EQ*0,0.474215,0.636234);o25=if(x25*x26*EQ*0,0.474215,0.636234);o35=if(x35*x36*EQ*0,0.474215,if(x38*EQ*0,0.636234,0.853814);o16=if(x15*x16*EQ*0,0.257697,1.259849);o26=if(x25*x26*EQ*0,0.257697,1.259849);o36=if(x35*x36*EQ*0,0.257697,if(x38*EQ*0,1.259849,0.904451

47、);o38=if(x35*x36*x38*NE*0,1.510221,2.508485);o11=0.029458;o21=0.029458;o31=0.029458;o41=0.029458;o51=0.029458;o12=0.030893;o22=0.030893;o32=0.030893;o42=0.030893;o52=0.030893;o27=9.688909;a=0.01;o11*x11/(y0+g11)<a;o21*x21/(y0+0.14716*x11+0.15209*x12-x13-x14-x15-x16+g21)<a;o31*x31/(y0+0.18431*x

48、21+0.18962*x22+0.14716*x11+0.15209*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o41*x41/(y0+0.1037*x31+0.18551*x32+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36-x27-x38+g41)<a;o51*x51/(y0+0.14449*x41+0.21065*x4

49、2+0.1037*x31+0.18551*x32+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+x33*b33+x34*b34+b23*x23+b24*x24+b13*x13+b14*x14+x25*b25+x26*b26+x15*b15+x16*b16-x43-x44-x35-x36-x27-x38+g51)<a;o12*x12/(y0+g11)<a;o22*x22/(y0+0.14716*x11+0.15209*x12-x13-x14-x15-x16+g21)<a;o32*x32/(y0+0.18431*x21+0.189

50、62*x22+0.14716*x11+0.15209*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o42*x42/(y0+0.1037*x31+0.18551*x32+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36- x27-x38+g41)<a;o52*x52/(y0+0.14449*x41+0.21065*x42+0.103

51、7*x31+0.18551*x32+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+x33*b33+x34*b34+b23*x23+b24*x24+b13*x13+b14*x14+x25*b25+x26*b26+x15*b15+x16*b16-x43-x44-x35-x36-x27-x38+g51)<a;o13*x13/(y0+g11)<a;o23*x23/(y0+0.14716*x11+0.15209*x12-x13-x14-x15-x16+g21)<a;o33*x33/(y0+0.18431*x21+0.18962*x22+

52、0.14716*x11+0.15209*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o43*x43/(y0+0.1037*x31+0.18551*x32+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36- x27-x38+g41)<a;o14*x14/(y0+g11)<a;o24*x24/(y0+0.14716*x11+0.

53、15209*x12-x13-x14-x15-x16+g21)<a;o34*x34/(y0+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+b13*x13+b14*x14-x23-x24-x15-x16-x25-x26-x27+g31)<a;o44*x44/(y0+0.1037*x31+0.18551*x32+0.18431*x21+0.18962*x22+0.14716*x11+0.15209*x12+b23*x23+b24*x24+b13*x13+b14*x14+b15*x15+b16*x16-x33-x34-x25-x26-x35-x36- x27-x38+g41)<a;o15*x15/(y0+g11)<a;o25*x25/(y0+0.14716*x11+0.15209*x12-x13-x14-x15-x16+g21)<a;o35*x35/(y0+0.