数学建模论文基金使用计划

1、题目基金使用计划摘要本文研究了关于基金使用计划的问题，主要目的在于设计资金的合理安排方法，实 现在一定条件下，使用有限的资金合理投资，达到最大的利润。并且我们建立了相应的 数学模型对该问题进行分析求解。对于第一问，我们在不影响奖学金发放的情况下，对利率较小的银行存款进行排除, 对每年的资金来源进行分析，列出所有可能发生的情况，然后建立一个线性方程组，求 出最大奖学金额度，方程组如下：Sixi ,1 x,2 xi,3 xi,5W (1ri)Xi,iA(i 1)(1ri) Xi,i(12r2)Xi 1,2(1r1) Xi,1(12 r2)Xi 1,2(13r3)Xi 2,3(i 2)(i 3,4)

2、(1 )Xi,1(1 2r2)Xi1,2(13r3)为 2,3(15r5)Xi4,5 A(i 5,6,7,8,9,10)使用Lingo软件对其进行编程求解，最后的得出最大奖学金额数为Z 109.8169（万元）。然而对于第二问，情况与第一问相似，但是又存在不同点，校方允许了购买国库券 这种投资方式。经过分析发现，由于国库券的发行不稳定，会产生三种不同的情况。所 以我们对这三种情况分别进行分析，运用第一题的思路，根据题目要求同样建立了一个 线性方程组（具体方程组见下文）。同样也是使用Lingo软件对其进行编程求解。最后得出在第一种情况下，校方每年能够发放的最大奖学金额数为Z146.8578（万元

3、）；而第二种方案的最大奖学金额数为Z 127.5222（万元）；最后第三种方案的最大奖学金额数为Z 131.7896（万元）第三问比较简单，校方要求第三年的奖金能够多出20% ,但是因为没有规定是只存 款不购国库券还是可存款也可购国库券，所以又要分成两个情况去讨论。我们对第一问和第二问中的方程组加以改进（改进方案见下文），将方程组中第三年支出的奖金额数上调20%,就能够得到满意的答案。最后，还是通过使用Lingo软件对其进行编程求解，在第一种情况下，得出的最大 奖学金额数为Z 107.5524（万元），而第三年百年校庆时的奖学金为 Z 129.0629（万元）； 而在第二种情况下，得出的最大奖

4、学金额数为Z 124.8507（万元），而第三年百年校庆时的奖学金为Z 149.8208（万元）。关键词：线性方程组Lingo求解最优化方案一、问题重述某校基金会有一笔数额为 M元的基金，打算将其存入银行或购买国库券。当前银行 存款及各期国库券的利率见下表。假设国库券每年至少发行一次，发行时间不定。取款 政策参考银行的现行政策。校基金会计划在n年内每年用部分本息奖励优秀师生，要求每年的奖金额大致相同, 且在n年末仍保留原基金数额。校基金会希望获得最佳的基金使用计划，以提高每年的 奖金额。请你帮助校基金会在如下的情况下设计出基金的最合理的使用方案，并且对于 M 5000万元，n 10年给出具体结

5、果：1、只存款不购国库券；2、可存款也可购国库券。3、学校在基金到位后的第3年要举行百年校庆，基金会希望这一年的奖金比其它年度多20%。银行存款税后年利率（%国库券年利率（。%活期0.792半年期1.664一年期1.800二年期1.9442.55三年期2.1602.89五年期2.3043.14二、模型假设1、每年发放的奖金数额都是相同的。2、不考虑通货膨胀或通货紧缩和国家经济政策对存款利率的影响。3、假设基金在年末到位，奖学金在基金到位后发放。4、假设购买国库券后所得的利息不用支付所得税。5、假设资金到位后，立即进入下一轮的存款，中间无时间问隔。6、假设购买国库券时，能够一次满足所有需求。7、

6、假设购买的国库券到期后，通过投资国库券而取得的资金在同年内不能够再购买购 买国库券，必须等到下一年的国库券发行后，才能继续购买。8、假设定期存款如果在没有到期之前取出，就按照活期存款利率计算。三、符号说明xi,j第i年存入银行j年期的存款金额yi,k第i年购买k年期的国库券金额M初始本金5000万A每年发放的固定奖金rjj年期的银行存款税后年利率Si第i年年初投入存款的总金额Wi第i年年末获得的资金总金额Ui第i年年初投入国库券的总金额Pkk年期的国库券年利率四、问题分析4.1 对问题一的分析在第一问中，校方要求只存款而不购买国库券，为了使得每年年终获得的奖学金最 多，那么就必须使得在银行获得

7、利息最大。因此，活期存款和半年期存款因为利率太低， 而且如果去除，并不会影响奖学金的发放，所以不予考虑。然后我们经过研究，把每年在存款到期日取出的资金分为：发放奖金、转存到一年 期存款、二年期存款、三年期存款、五年期存款这五个款项。通过对这五个款项的分析与求解，求出最大的每年发放的奖学金。4.2 对问题二的分析在第二问中，校方要求可以存款也可以购买国库券，但是国库券的购买与直接存款 有很大的区别，所以要对第一问中的解法做一些调整。因为同年期的国库券的利率要大于同年期的银行存款利率，所以在选择时，我们优 先考虑购买国库券。但是，由于每年发行的国库券的时间和次数不同，所以可能会发生 三种情况：一是

8、，国库券在年初准备存钱的时候发行了，我们就在第一题的基础上，将准备存 入相应年期存款的资金用于购买同样年期的国库券，然后计算出最大的奖金额。二是，国库券没有在存入资金之前发行，为了不让资金发生闲置，我们设立了另一 种解决方案：以二年期的国库券为例：由于在年初年初投放资金时不能购买国库券，我们先将购 买国库券的资金全部用于半年期存款，如果在该半年内发行了国库券，我们就将资金 全部取出购买国库券，在国库券到期的那年将本息全部用于半年期存款，到期后转入活期存款；如果在该半年内没有发行国库券，我们将半年到期的自己全部用于活期存款，用于购买下半年一定会发行的国库券，国库券到期之后再全部转入活期存款。因此

9、, 我们将其运转周期定为三年，在这三年里，不管国库券什么时候发行，该部分资金一 定有两年是用于存国库券，有半年用于存半年期，还有半年是存活期。即采用活期、半 年期、国库券的“组合式”投资。所以同理，三年期国库券和五年期国库券的周期分别 为四年，六年。而第三种则是第二种情况的特殊化，也就是国库券恰好在年中发行了。而在这种情 况下，在国库券发行的前半年，用于购买国库券的资金就可以存入半年期的银行存款； 而在国库券到期后，还会有半年的时间，我们就能够将购买国库券后得到的本息再次存 入半年期的银行存款中去。这样就是半年期和国库券的“组合式”投资，这样的投资方法会比第二种情况获利更多。然后对这三个情况分

10、别求解，得出每年发放的奖学金。4.3 对问题三的分析在第三问中，校方要求第三年的百年校庆中的奖学金能够上涨20%这个问题比较简单。但是因为没有规定是只存款不购国库券还是可存款也可购国库券，所以又要分成 两个情况去讨论。又因为第一种情况与第一问相似，第二种情况与第二问相似，所以只需要在第一问 和第二问的基础上对方程进行一部分的改进，将第三年的奖学金上调 20% ,从而改为 1.2A,再利用Lingo进行求解。五、模型的建立与求解5.1 问题一的模型建立与求解我们对每一年的年末所能够获得的、用于发放奖学金的现金来源进行分析。比如， 在第一年的年末，校方所能够提取的现金，只有可能是在第一年年初投入一

11、年期存款的 基金的本金和所获得的利息；而在第二年年末，校方所能够提取现金的来源就可能有两 个，一是第二年年初投入一年期存款的基金的本金和所获得的利息，二是第一年年初投 入二年期存款的基金本金和所获得的利息；而在第三年的年末，校方所能够提取现金的 方法就有三种。所以根据上述的思想，我们建立了一个线性方程组，用于求解最大的奖学金额数该方程组如下:(i 1)(i 2)(i 3,4)Sixi ,1 x,2 xi,3 xi,5W(1ri)xi,iAW(1ri)xi,i(12r2)xi1,2 A(1 n) xi,1(1 2r2) xi1,2(1 3r3)为 2,3 A(1r1) xi,1(12 r2)xi

12、 1,2(13r3)xi 2,3(15r5)xi 4,5 A(i 5,6,7,8,9,10)然后我们将数据代入方程组，利用 Lingo软件对其进行求解（编程代码见附录I ）,求得答案为（输出结果见附录II ）:最大奖学金额数为Z 109.8169（万元） 而存款方案如下表所示：表1 ：问题1第i年存入j年期的金额（万元）第一 年第二 年第三 年第四 年第五 年第六 年第七 年第八 年第九 年第十 年一年 期396.8000000000二年 期200.5000000000三年 期195.6195.600000000五年 期420798.598.598.598.5458200005.2 问题二的模

13、型建立与求解首先对第一种情况进行求解。对于第一种情况，经过分析发现，与第一问十分相似,直接将第一问中的与国库券有着相同年期的银行存款利率改为国库券利率，然后通过运用Lingo软件进行编程求解（编程代码见附录I ,输出结果见附录II ）,得出：最大奖学 金额数为Z 146.8578（万元）。而投资方案如下表所示：表2:问题2.1第i年存入j年期的金额（万元）第一 年第二 年第三 年第四 年第五 年第六 年第七 年第八 年第九 年第十 年一年存款277.0000000000二年 国库260.5000000000三年 国库490.7135.10259.500135.1000五年 国库39710126

14、.9126.9044480000然后对第二种情况进行求解。因为国库券的周期与存款不同，所以需要重新进行计 算，根据不同的存款利率和国库券利率列出如下的“组合式”购买利率公式：PiP2P31r1 -) 1.0639421r1 -) 1.1000821r1) 1.1712521(1 2ri)(1 r。-)(121(1 3r2)(1r。-)(121(15r5)(1r。)(12然后根据第一问中的思想建立一个新的线性方程组用于求解最大的奖学金额数。该方程组如下：Sixi ,1 x,2 xi,3 xi,5U k yi,1 yi,2yi ,3W(1h)xi,1A(i1)W (1 h) xi,1(12r2)

15、xi 1,2 A(i 2)W(11)xi,1(12r2)xi 1,2(1 3r3)为 2,3(1 R) yi2,1 A(i3)W(1r1)xi,1(12r2)xi1,2(13r3)为2,3(1P1)yi 2,1(1 P2)yi3,2 A (i 4)W(1r1)xi,1(12r2)xi1,2(1%)X2,3(15r5)为 4,5(1P1) Yi2,1(1P2) Yi 3,2 A(i 5)W(1r1)xi,1(12r2)xi1,2(13r3)x2,3(15r5)xi 4,5(1P1)yi2,1 (1P2)yi3,2(1P3)yi5,3 A(i 6,7,8,9,10)然后我们将数据代入方程组，利用

16、Lingo软件对其进行求解（编程代码见附录I ）,求得答案为（输出结果见附录II ）:最大奖学金额数为Z 127.5222（万元）而投资方案如下表所示:表3: 2.2第i年存入j年期的金额(万元)第一 年第二 年第三 年第四 年第五 年第六 年第七 年第八 年第九 年第十 年一年存款232.2000125.200000二年 存款227.5000000000三年 仔演222.0000000000五年 存款0000000000二年 国库0000000000三年 国库4209000000000五年 国库108.8108.8108.8108.8437700000最后对第三种情况进行分析求解，因为是第二

17、种情况的特殊化，所以只需要在第二问中进行一些改动。由于在第三种情况中，没有活期银行存款的存在，而是变成了两个 半年期，所以“组合式”购买利率发生了改变，更改后的公式为：R (1 2ri)(1 ri) 1.06849P2 (1 3r2 )(1 ri) 1.10478访(1 5r5)(1 A) 1.17625然后还是根据第二种情况中的线性方程组求解最大的奖学金额数。最后求出的最大奖学金额数为(编程代码见附录I ,输出结果见附录II ): Z 131.7896(万元)表4: 2.3第i年存入j年期的金额(万元)第一 年第二 年第三 年第四 年第五 年第六 年第七 年第八 年第九 年第十 年一年存款3

18、56.7000000000二年 存款234.7000000000三年 仔演0000000000五年0000000000存款二年 国库228.2000000000三年 国库4068119.200000000五年 国库112.0112.0112.0112.04362000005.3 问题三的模型建立与求解对于第三问，由于第三年为百年校庆，所以校方要求当年的奖学金要上调 20%。但 是因为没有规定是只存款不购国库券还是可存款也可购国库券，所以又要分成两个情况去讨论。在第一种情况中，我们只存款而不购买国库券，我们将第三年的奖学金发放额改为 原来的120%,对第一问方程组中的i 3,4的方程进行修改：W

19、i (1 ri) Xi,i (1 2r2) x 1,2 (1 3r3)为 2,3 1.2A(i 3)W (1 1) Xi,1 (1 2r2) x 1,2 (1 3r3) Xi 2,3 A(i 4)然后再通过利用Lingo软件对其进行求解(编程代码见附录I),求得答案为(输出 结果见附录II ):最大奖学金额数为Z 107.5524(万元)。而投资方案如下表所示：表5: 3.1第i年存入j年期的金额(万元)第第第第四 年第五 年第六 年第七 年第八 年第九 年第十 年一年存款388.6000122.600000二年 存款196.4000000000三年 仔演211.8191.600000000五

20、年 仔演420396.4496.4496.4496.4445790000而在第二种情况中，我们可以存款也可以购买国库券，我们对第二问中第二种情况的求解方程组中的i 3的方程进行修改：W (1r1)Xi,1(12也)x 1,2 (13r3)Xi2,3(1P1)yi2,1 1.2A(i 3)将A改为了 1.2A,再通过利用Lingo软件对其进行求解(编程代码见附录I),求得 答案为(输出结果见附录II ):最大奖学金额数为Z 124.8507(万元)。而投资方案如下表所示：表6: 3.2第i年存入j年期的金额（万元）第一 年第二 年第三 年第四 年第五 年第六 年第七 年第八 年第九 年第十 年一

21、年存款227.3000122.600000二年 存款222.7000000000三年 仔演240.0000000000五年 仔演0000000000年库 二国0000000000三年 国库4202000000000年库 五国106.6106.6106.6106.6437500000六、模型的评价与改进6.1 模型的优点1、运用线性规划的思想，使用Lingo进行编程求解，方便运算，简单实用。2、模型细致，分析具体，将各种情况都考虑了进去。3、运用图表，结果更加直观，可操作性高。6.2 模型的缺点1、假设过于理想化，与现实生活有一定的差异。2、重复计算较多，不借助软件无法求解。七、模型的推广这个模

22、型的假设过于理想化，是建立在不受社会经济情况的影响的情况下。比如， 通货膨胀、通货紧缩、国家颁布的经济政策等，都没有考虑到。这个模型可以适用于短期的一些较为复杂的资金分配问题中，或者是一些经济较为稳定的，银行利率不会产生大幅度波动的地区所产生的一些经济问题中，从而提高资金 的利用率，能够在固定的时间中获得最大的利益。八、参考文献1马君儿，李东明，资金的最优分配，工科数学，18卷5期：66-70, 20022姜启源，数学模型M,北京：高等教育出版社，19893沧浦，最优控制的理论与方法M,北京：国防工业出版社，1989附录 I问题 1 的 Lingo 软件编程代码：MODEL:SETS:YEAR

23、1/1.10/:S;! 题目的 1-10 年;YEAR2/1.10/:W;NIANQI/1.6/:r;! 题目所给的年利率;LINK(YEAR1,NIANQI):x;! 第 i 年投给 j 年期的金额;ENDSETSMAX=A; ! 每年发给优秀师生的奖金;for(YEAR1(i):S(i)=x(i,3)+x(i,4)+x(i,5)+x(i,6);for(year2(i)|i#ge#5#and#i#le#10:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1+3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)-A);for(year2(i)|

24、i#ge#3#and#i#le#4:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1+3*r(5)*x(i-2,5)-A);for(year2(i)|i#ge#2#and#i#le#2:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)-A);for(year2(i)|i#ge#1#and#i#le#1:W(i)=(1+r(3)x)*x(i,3)-A);S(2)=W(1);S(3)=W(2);S(4)=W(3);S(5)=W(4);S(6)=W(5);S(7)=W(6);S(8)=W(7);S(9)=W(8);S(10)=W(9);S

25、(1)=M;W(10)=5000;DATA:M=5000;R=0.00792 0.01664 0.01800 0.01944 0.02160 0.02304;ENDDATAEND问题 2.1 的 Lingo 软件编程代码：MODEL:SETS:YEAR1/1.10/:S;! 题目的 1-10 年;YEAR2/1.10/:W;NIANQI/1.6/:r;! 题目所给的年利率;LINK(YEAR1,NIANQI):x;! 第 i 年投给 j 年期的金额;ENDSETSMAX=A; ! 每年发给优秀师生的奖金;for(YEAR1(i):S(i)=x(i,3)+x(i,4)+x(i,5)+x(i,6)

26、;for(year2(i)|i#ge#5#and#i#le#10:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1+3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)-A);for(year2(i)|i#ge#3#and#i#le#4:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1+3*r(5)*x(i-2,5)-A);for(year2(i)|i#ge#2#and#i#le#2:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)-A)for(year2(i)|i#ge#1#and

27、#i#le#1:W(i)=(1+r(3)*x(i,3)-A);S(2)=W(1);S(3)=W(2);S(4)=W(3);S(5)=W(4);S(6)=W(5);S(7)=W(6);S(8)=W(7);S(9)= W(8);S(10)=W(9);S(1)=M;W(10)=5000;DATA:M=5000;R=0.00792 0.01664 0.01800 0.02550 0.02890 0.03140;ENDDATAEND问题 2.2 的 Lingo 软件编程代码：MODEL:SETS:YEAR1/1.10/:S,U;YEAR2/1.10/:W;NIANQI1/1.6/:r;NIANQI2/1

28、.3/:P;LINK(YEAR1,NIANQI1):x;LINKS(YEAR1,NIANQI2):Y;endsetsMAX=A; !每年发给优秀师生的奖金;for(YEAR1(i):S(i)=x(i,3)+x(i,4)+x(i,5)+x(i,6);for(YEAR1(i):U(i)=Y(i,1)+Y(i,2)+Y(i,3);for(YEAR2(i)|i#ge#6#and#i#le#10:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+( 1+3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(

29、i-3,2)+ (1+P(3)*Y(i-5,3)-A);for(YEAR2(i)|i#ge#5#and#i#le#5:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)-A );for(YEAR2(i)|i#ge#4#and#i#le#4:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)-A

30、);for(YEAR2(i)|i#ge#3#and#i#le#3:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1+3*r(5)*x(i-2,5)+(1+P(1)*Y(i-2,1)-A);for(YEAR2(i)|i#ge#2#and#i#le#2:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)-A);for(YEAR2(i)|i#ge#1#and#i#le#1:W(i)=(1+r(3)*x(i,3)-A);S(2)+U(2)=W(1);S(3)+U(3)=W(2);S(4)+U(4)=W(3);S(5)+U(5)=W(4);S(

31、6)+U(6)=W(5);S(7)+U(7)=W(6);S(8)+U(8)=W(7);S(9)+U(9)=W(8);S(10)+U(10)=W(9);S(1)+U(1)=M;W(10)=5000;DATA:M=5000;r=0.00792 0.01664 0.01800 0.01944 0.02160 0.02304;P=0.06394 0.10008 0.17125;ENDDATAEND问题 2.3 的 Lingo 软件编程代码：MODEL:SETS:YEAR1/1.10/:S,U;YEAR2/1.10/:W;NIANQI1/1.6/:r;NIANQI2/1.3/:P;LINK(YEAR1,

32、NIANQI1):x;LINKS(YEAR1,NIANQI2):Y;endsetsMAX=A; !每年发给优秀师生的奖金;for(YEAR1(i):S(i)=x(i,3)+x(i,4)+x(i,5)+x(i,6);for(YEAR1(i):U(i)=Y(i,1)+Y(i,2)+Y(i,3);for(YEAR2(i)|i#ge#6#and#i#le#10:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+( 1+3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)+ (1+P(3)*Y

33、(i-5,3)-A);for(YEAR2(i)|i#ge#5#and#i#le#5:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)-A );for(YEAR2(i)|i#ge#4#and#i#le#4:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)-A);for(YEAR2(i)|i#

34、ge#3#and#i#le#3:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)+(1+P(1)*Y(i-2,1)-A);for(YEAR2(i)|i#ge#2#and#i#le#2:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)-A) ;for(YEAR2(i)|i#ge#1#and#i#le#1:W(i)=(1+r(3)*x(i,3)-A);for(YEAR2(i)|i#ge#1#and#i#le#9:W(i)=S(i+1)+U(i+1);S(1)+U(1)=M;W(10)=5000;DATA

35、:M=5000;r=0.00792 0.01664 0.01800 0.01944 0.02160 0.02304;P=0.06849 0.10478 0.17625;ENDDATAEND问题 3.1 的 Lingo 软件编程代码：MODEL:SETS:YEAR1/1.10/:S;! 题目的 1-10 年;YEAR2/1.10/:W;NIANQI/1.6/:r;! 题目所给的年利率;LINK(YEAR1,NIANQI):x;! 第 i 年投给 j 年期的金额;ENDSETSMAX=A; ! 每年发给优秀师生的奖金;for(YEAR1(i):S(i)=x(i,3)+x(i,4)+x(i,5)+x

36、(i,6);for(year2(i)|i#ge#5#and#i#le#10:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+( 1+3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)-A);for(year2(i)|i#ge#4#and#i#le#4:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)-A);for(year2(i)|i#ge#3#and#i#le#3:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5

37、)-1.2*A);for(year2(i)|i#ge#2#and#i#le#2:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)-A);for(year2(i)|i#ge#1#and#i#le#1:W(i)=(1+r(3)*x(i,3)-A);for(YEAR2(i)|i#ge#1#and#i#le#9:W(i)=S(i+1);S(1)=M;W(10)=5000;DATA:M=5000;R=0.00792 0.01664 0.01800 0.01944 0.02160 0.02304;ENDDATAEND问题 3.2 的 Lingo 软件编程代码：MODEL:SET

38、S:YEAR1/1.10/:S,U;YEAR2/1.10/:W;NIANQI1/1.6/:r;NIANQI2/1.3/:P;LINK(YEAR1,NIANQI1):x;LINKS(YEAR1,NIANQI2):Y;endsetsMAX=A; !每年发给优秀师生的奖金;for(YEAR1(i):S(i)=x(i,3)+x(i,4)+x(i,5)+x(i,6);for(YEAR1(i):U(i)=Y(i,1)+Y(i,2)+Y(i,3);for(YEAR2(i)|i#ge#6#and#i#le#10:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+( 1+3*r(5)

39、*x(i-2,5)+(1+5*r(6)*x(i-4,6)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)+ (1+P(3)*Y(i-5,3)-A);for(YEAR2(i)|i#ge#5#and#i#le#5:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1 +3*r(5)*x(i-2,5)+(1+5*r(6)*x(i-4,6)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)-A );for(YEAR2(i)|i#ge#4#and#i#le#4:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)

40、+(1+3*r(5)*x(i-2,5)+(1+P(1)*Y(i-2,1)+(1+P(2)*Y(i-3,2)-A);for(YEAR2(i)|i#ge#3#and#i#le#3:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)+(1+3*r(5)*x(i-2,5)+(1+P(1)*Y(i-2,1)-1.2A);for(YEAR2(i)|i#ge#2#and#i#le#2:W(i)=(1+r(3)*x(i,3)+(1+2*r(4)*x(i-1,4)-A);for(YEAR2(i)|i#ge#1#and#i#le#1:W(i)=(1+r(3)*x(i,3)-A);S(2)

41、+U(2)=W(1);S(3)+U(3)=W(2);S(4)+U(4)=W(3);S(5)+U(5)=W(4);S(6)+U(6)=W(5);S(7)+U(7)=W(6);S(8)+U(8)=W(7);S(9)+U(9)=W(8);S(10)+U(10)=W(9);S(1)+U(1)=M;W(10)=5000;DATA:M=5000;r=0.00792 0.01664 0.01800 0.01944 0.02160 0.02304;P=0.06394 0.10008 0.17125;ENDDATAEND附录 II问题 1 的 Lingo 编程输出结果：Global optimal soluti

42、on found.Objective value:109.8169X( 1, 1)X( 1, 2)X( 1, 3)X( 1, 4)X( 1, 5)X( 1, 6)X( 2, 1)X( 2, 2)X( 2, 3)X( 2, 4)X( 2, 5)X( 2, 6)X( 3, 1)X( 3, 2)X( 3, 3)X( 3, 4)X( 3, 5)X( 3, 6)X( 4, 1)X( 4, 2)X( 4, 3)X( 4, 4)0.000000 0.000000 396.7621 200.4946 195.6140 4207.129 0.000000 0.000000 0.000000 0.000000 1

43、95.6140 98.47287 0.000000 0.000000 0.000000 0.000000 0.000000 98.47287 0.000000 0.000000 0.000000 0.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.2709247E-030.7466756E-030.0000000.0000000.0000000.0000000.7316685E-030.4673388E-030.8708696E-030.0000000.0000000.0000000.0000000

44、.8496704E-030.1367398E-03 0.000000 0.000000 0.0000000.1087030E-02 0.8346467E-03 0.3884287E-03 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000000.2429382E-03 0.6695441E-03 0.0000000.9874166E-010.0000000.0000000.6560873E-03 0.4190628E-03 0.7809089E-03 0.9675710E-01 0.000

45、000 0.000000 0.0000000.7618995E-03 0.9440178E-01 0.9440178E-01 0.000000 0.0000000.9747403E-03 0.9273259E-01 0.9273259E-01 0.9273259E-01X( 4, 5)X( 4, 6)X( 5, 1)X( 5, 2)X( 5, 3)X( 5, 4)X( 5, 5)X( 5, 6)X( 6, 1)X( 6, 2)X( 6, 3)X( 6, 4)X( 6, 5)X( 6, 6)X( 7, 1)X( 7, 2)X( 7, 3)X( 7, 4)X( 7, 5)X( 7, 6)X( 8,

46、 1)X( 8, 2)X( 8, 3)X( 8, 4)X( 8, 5)X( 8, 6)X( 9, 1)X( 9, 2)X( 9, 3)X( 9, 4)X( 9, 5)X( 9, 6)X( 10, 1)X( 10, 2)X( 10, 3)X( 10, 4)X( 10, 5)X( 10, 6) 0.000000 98.47287 0.000000 0.000000 0.000000 0.000000 0.000000 98.47287 0.000000 0.000000 0.000000 0.000000 0.000000 4581.974 0.000000 0.000000 0.000000 0

47、.000000 0.000000 0.000000 0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.0000000.000000 0.000000 0.000000 0.000000 0.000000 0.0000002.1 的 Lingo 编程输出结果：Global optimal solution found.Objective value:146.8578X( 1, 1)0.0000000.000000X( 1, 2)X( 1, 3)X(

48、 1, 4)X( 1, 5)X( 1, 6)X( 2, 1)X( 2, 2)X( 2, 3)X( 2, 4)X( 2, 5)X( 2, 6)X( 3, 1)X( 3, 2)X( 3, 3)X( 3, 4)X( 3, 5)X( 3, 6)X( 4, 1)X( 4, 2)X( 4, 3)X( 4, 4)X( 4, 5)X( 4, 6)X( 5, 1)X( 5, 2)X( 5, 3)X( 5, 4)X( 5, 5)X( 5, 6)X( 6, 1)X( 6, 2)X( 6, 3)X( 6, 4)X( 6, 5)X( 6, 6)X( 7, 1)X( 7, 2)X( 7, 3)X( 7, 4)X( 7,

49、5)X( 7, 6)X( 8, 1)X( 8, 2)X( 8, 3)0.000000 277.0126 260.5020 490.7406 3971.745 0.000000 0.000000 0.000000 0.000000 135.1411 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 126.9298 0.000000 0.000000 0.000000 0.000000 259.5002 126.9298 0.000000 0.000000 0.000000 0.000000 0.000000 0.000000 0.000

50、000 0.000000 0.000000 0.000000 0.000000 4448.451 0.000000 0.000000 0.000000 0.000000 135.1411 0.000000 0.000000 0.000000 0.0000000.0000000.0000000.0000000.0000000.0000000.0000000.0000000.1592579E-020.1761287E-020.0000000.2985110E-030.0000000.0000000.1705985E-020.1657247E-030.1420565E-020.0000000.000

51、0000.0000000.0000000.1373897E-020.0000000.0000000.0000000.0000000.2796277E-020.1620767E-020.1193949E-020.2746949E-030.0000000.0000000.2626484E-030.0000000.0000000.0000000.0000000.0000000.1118469E-020.1264284E-020.0000000.9831608E-010.0000000.0000000.1474490E-020.3927650E-030.1227800E-020.9547899E-010.0000000.0000000.2416936E-030.1187465E-020.9234234E-010.9234234E-010.0000000.0000000.21