1、西安铁路职业技术学院第一届数学建模竞赛承 诺 书我们仔细阅读了西安铁路职业技术学院大学生数学建模竞赛的竞赛规则.我们完全明白,在竞赛开始后参赛队员不能以任何方式(包括电话、电子邮件、网上咨询等)与队外的任何人(包括指导教师)研究、讨论与赛题有关的问题。我们知道,抄袭别人的成果是违反竞赛规则的, 如果引用别人的成果或其他公开的资料(包括网上查到的资料) ,必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺,严格遵守竞赛规则,以保证竞赛的公正、公平性。如有违反竞赛规则的行为,我们将受到严肃处理。我们参赛选择的题号是(从 C/D 中选择一项填写): C 我们的参赛报名号为
2、(如果赛区设置报名号的话): 所属学校(请填写完整的全名): 西安铁路职业技术学院(龙首校区) 参赛队员 (打印并签名) :1. 张玉辇 (组长) 2. 张斌 3. 崔超 指导教师或指导教师组负责人 (打印并签名): 日期: 2011 年 6 月 17 日赛区评阅编号(由赛区组委会评阅前进行编号):西安铁路职业技术学院第一届数学建模竞赛编 号 专 用 页赛区评阅编号(由赛区组委会评阅前进行编号):赛区评阅记录(可供赛区评阅时使用):评阅人评分备注学院统一编号(由赛区组委会送交全国前编号):学院评阅编号(由全国组委会评阅前进行编号):1有价证券投资决策优化模型摘要针对有价证券投资的决策问题,我们
3、在充分理解题意的基础上,提出合理的假设并采用了连续线性规划模型对本题进行了讨论和求解,解决了以下问题。对于问题一,该问题的目标是有价证券回收的利息为最高,要做的决策是投资计划,即应购买的各种证券的数量的分配。综合考虑:特定证券购买、资金限制、平均信用等级、平均年限等条件,按照题目所求,依据决策变量、决策目标和约束条件指定 0-1变量建立 0-1 规划模型,并用 LINGO 软件求得应向市政 1(X1=218.18 万元) 、政府1(X3=736.36 万元) 、市政 2(X5=45.45 万元)投资,才能获得最大收益(Z=29.84 万元)。对于问题二,在引入借贷资金的情形下我们进一步假设约束
4、条件,建立了连续性优化模型,并与问题一进行比较得出借贷资金利息比收益金额还要高,所以对于问题二用借贷资金投资是不划算的。对于问题三的税前收益发生的变动,我们又建立了线性规划模型,解得在第一种情形下投资对象不变;在第二种情形下应向市政 1(X1=336 万元)、政府 2(X4=648 万元) 、市政 2(X5=16 万元)投资,所得的收益为(Z=29.42 万元) 。关键字: 证券投资 0-1 规划 决策目标 连续线性规划 加权算术平均值2一 、 问 题 重 述某 银 行 经 理 计 划 用 一 笔 资 金 进 行 有 价 证 券 的 投 资 , 可 供 购 进 的 证 券 以 及 其 信 用
5、等 级 、 到 期年 限 、 收 益 如 下 表 所 示 。 按 照 规 定 , 市 政 证 券 的 收 益 可 以 免 税 , 其 他 证 券 的 收 益 需 按50%的 税 率 纳 税 。此 外 还 有 以 下 限 制 :一 .政 府 及 代 办 机 构 的 证 券 总 共 至 少 需 购 进400万 元 ;二 .所 购 证 券 的 平 均 信 用 等 级 不 超 过1.4( 信 用 等 级 数 字 越 小 , 信 用 程 度 越 高 );三 .所 购 证 券 的 平 均 到 期 年 限 不 超 过5年 ;表 1 证 券 信 息证 券 名 称 证 券 种 类 信 用 等 级 到 期 年 限
6、 到 期 税 前 收 益( %)A 市 政 ( 1) 2 9 4.3B 代 办 机 构 2 15 5.4C 政 府 ( 1) 1 4 5.0D 政 府 ( 2) 1 3 4.4E 市 政 ( 2) 5 2 4.51. 若 该 经 理 有 1000万 元 资 金 , 应 如 何 投 资 ?2. 如 果 能 够 以 2.75%的 利 率 借 到 不 超 过 100万 元 资 金 , 该 经 理 应 如 何 操 作 ?3. 在 1000万 元 资 金 情 况 下 , 若 证 券 A的 税 前 收 益 增 加 为 4.5%, 投 资 应 否 改 变 ? 若 证 券 C的 税 前收 益 减 少 为 4.
7、8%, 投 资 应 否 改 变 ?二、模型假设1、假设证劵税收政策、信用等级、到期年限、到期税前收益在投资年限内稳定不变;2、假设该经理投资不会受到年限过长而导致资金周转困难的影响;3、假设各证劵的风险损失率为零,无相互影响,且投资不需要任何交易费;4、假设所借贷资金的利率不变(利息等于所给的利率乘上借贷资金再乘上借贷时间) 。3三、符号说明表 2四、问题分析根据现有的投资方式,为解决投资方案问题,运用 0-1 规划和连续性规划模型,根据客观条件,来确定各种投资方案。利用线性规划模型,得到最优投资方案来解决已提出的问题,使该经理的投资获得最大的收益。问题 1:在前面假设成立的条件下,由题目给出
8、的限制条件(市政证券的收益可以免税;其他证券的收益按 50%的税率纳税,即政府及代办机构的证券总共至少需购进400 万元,所购证券的平均信用等级不超过 1.4,所购证券的平均到期年限不超过 5 年)假设 Xi (i=1,2,3,4,5)分别为证券 A、证券 B、证券 C、证券 D、证券 E 的投入资金,所获得的纯利润为 Z,假设税前收益为年收益。对于平均信用等级和平均到期年限则需用加权算术平均值的算法求得,即用各个信用等级乘以相应的权值,然后相加,所得之和再除以所有的权之和。然后建立 0-1 规划模型,利用 LINGO 软件求解。问题 2:该题又引入了借贷资金,我们假设借贷资金为 F 万元,借
9、贷时间为 T=1。由问题 2 知借贷利率为 2.75%,则借贷利息为 0.0275F,因此需要用总利润减去借贷利息所得到的值才为该经理的纯利润。再与借贷利息相比,可以得出投入 100 万元所得的收益。问题 3:需要分两种情形,第一种情形,在假设成立的条件下,同问题 1 一样,只是把 A 证券的税前收益增加为 4.5%,同样可以得到线性优化模型。第二种情形,在假设成立的条件下,同问题 1 一样,只是把 C 证券的税前收益减少为 4.8%,同样可以得到线性优化模型。符号 说明 单位Xi 投资证劵的金额(i=1,2,3,4,5)(万元)YJ 是否投资证券(j=1,2,3,4,5) (万元)Z 投资之
10、后所获得的总收益 (万元)F 借贷金额 (万元)T 借贷时间 (年)S 贷款利息 (万元)R 利润变化量 (万元)4五、模型的建立与求解模型建立:对于问题 1:在前面假设成立的条件下,由题目给出的限制条件(市政证券的收益可以免税;其他证券的收益按 50%的税率纳税,即政府及代办机构的证券总共至少需购进 400 万元,所购证券的平均信用等级不超过 1.4,所购证券的平均到期年限不超过 5 年)假设 Xi (i=1,2,3,4,5)分别为证券 A、证券 B、证券 C、证券 D、证券 E 的投入资金,所获得的总利润为 Z,假设税前收益为年收益。对于平均信用等级和平均到期年限则需用加权算术平均值的算法
11、求得,即用各个信用等级乘以相应的权值,然后相加,所得之和再除以所有的权之和。我们建立了如下 0-1 规划模型:即 MaxZ=0.043X1*Y1+(0.054 0.5)X2*Y2+(0.05 0.5)X3*Y3+(0.044 0.5)X4*Y4+0.045X5*Y5X1+X2+X3+X4+X5=400(2X1+2X2+X3+X4+5X5)/( X1+X2+X3+X4+X5)=0,i=1,2,3,4,5Yj=0 或 1,j=1,2,3,4,5 (指定 0-1 变量)模型整理得:Max Z=0.043X1*Y1+0.027X2*Y2+0.025X3*Y3+0.022X4*Y4+0.045X5*Y5
12、X1+X2+X3+X4+X5=4006X1+6X2-4X3-X4+36X5=0,i=1,2,3,4,5Yj=0 或 1,j=1,2,3,4,5 (指定 0-1 变量)模型求解:(LINGO 程序见附录)应用 LINGO 求解可得银行经理的最大收益为 Z=29.84 万元,五种证券类型最优投入的资金方法如下表:5表 3 五种证券类型最优投入的资金方法证券类型 A B C D E Z投资金额 (万元)218.18 0 736.34 0 45.45 29.84数值 Y1 Y2 Y3 Y4 Y5投资对象 1 0 1 0 1由运算结果可以看出对于 B 证券和 D 证券不需要投资。对于 B 证券它的到期年
13、限相对较长而且还需纳 50%的税率,相对于其他证券类型这显然是不理想的;对于 D 证券它的税前收益相对较低也是不理想的;而 A 证券和 C 证券信用等级、到期年限都是较低的税前收益也可以(A 证券还可以免税) ;E 证券信用等级较低,但到期年限较短可以适量投资。所以我们认为这一优化模型是合理的。模型建立:对于问题 2:该题又引入了借贷资金,我们假设借贷资金为 F 万元,借贷时间为 T=1。由问题 2知借贷利率为 2.75%,则借贷利息为 0.0275F,这样需要用总利润减去借贷利息所得值才为该经理的纯利润。我们建立了以下线性优化模型:即:Max Z=0.043X1+(0.054 0.5)X2+
14、(0.05 0.5)X3+(0.044 0.5)X4+0.045X5-0.0275F*TX1+X2+X3+X4+X5=400S.t (2X1+2X2+X3+X4+5X5)/( X1+X2+X3+X4+X5)=0,X2=0,X3=0,X4=0,X5=0F=400S.t 6X1+6X2-4X3-X4+36X5=0,X2=0,X3=0,X4=0,X5=0FR,说明投入 100 万元所得的收益为负值该经理不应贷款投资。模型建立:对于问题 3:第一种情形:在假设成立的条件下,同问题 2 一样,只是把 A 证券的税前收益增加为 4.5%,同证券类型 A B C D E F投资金额 (万元)240.00 0
15、 810.00 0 50.00 100.007样可以得到如下线性优化模型:即:Max Z=0.045X1+0.027X2+0.025X3+0.022X4+0.045X5X1+X2+X3+X4+X5=400S.t (2X1+2X2+X3+X4+5X5)/( X1+X2+X3+X4+X5)=0,X2=0,X3=0,X4=0,X5=0模型整理得: Max Z=0.045X1+0.027X2+0.025X3+0.022X4+0.045X5X1+X2+X3+X4+X5=400S.t 6X1+6X2-4X3-X4+36X5=0,X2=0,X3=0,X4=0,X5=0模型求解:(LINGO 程序见附录)应用
16、 LINGO 求解可得银行经理的最大收益为 Z=30.27 万元,五种证券类型最优投入的资金方法如下表: 表 5 五种证券类型最优投入的资金方法证券类型 A B C D E投资金额 (万元)218.18 0 736.36 0 45.45第二种情形:在假设成立的条件下,同问题 2 一样,只是把 C 证券的税前收益减少为 4.8%,同样可以得到如下线性优化模型:即: Max Z=0.043X1+(0.054 0.5)X2+(0.048 0.5)X3+(0.044 0.5)X4+0.045X5X1+X2+X3+X4+X5=400S.t (2X1+2X2+X3+X4+5X5)/( X1+X2+X3+X
17、4+X5)=0,X2=0,X3=0,X4=0,X5=0模型整理得:Max Z=0.043X1+0.027X2+0.024X3+0.022X4+0.045X5X1+X2+X3+X4+X5=400S.t 6X1+6X2-4X3-X4+36X5=0,X2=0,X3=0,X4=0,X5=0模型求解:(LINGO 程序见附录)应用 LINGO 求解可得银行经理的最大收益为 Z=29.42 万元,五种证券类型最优投入的资金方法如下表: 表 6 五种证券类型最优投入的资金证券类型 A B C D E投资金额 (万元)336.00 0 0 648.00 16.00对于第一种情形由运算结果分析,在证券 A 的税
18、前收益增加为 4.5%的情况下,对于 X1,X3,X5 的投入金额没有发生变化,只是总收益有所增加,所以该投资方案不需要改变。对于第二种运算结果分析,在证券 C 的税前收益减少为 4.8%的情况下,对于证券投资的种类发生了变化,投资对象变为了 X1(市政)、X4(政府 2)、X5(市政 2),所以投资方案应改变。 六、模型的优缺点优点:1.本文采用的方法,针对投资选择可能出现的情况建立了 3 个数学模型,所采用的方案能够对于实际的投资选择给出相应优化,本方案具有较强的通用性。 2.本方案也具有一般性,针对生活中的连续优化问题也得到了推广,如解决农业种植、渔业养殖的优化问题。3.所建线性规划模型
19、简单易懂,方便使用。缺点:1.本文并未考虑到在金融危机发生的条件下能解决问题的方案。2.本文没有对月税前收益予以考虑。9七、模型扩展仿照文中思路,该模型可容易扩展到更多的选择对象,可以应用于资本市场、风险管理,也可以解决如奶制品的生产与销售、农业种植、渔业养殖的优化问题。对于借贷资金也可以扩展到分期借贷和小额借贷的情况。对于在信用等级、到期年限、到期税前收益等因素发生变化时,解决思路相类似,但模型和求解都更复杂,可作进一步研究。八、参考文献1、肖华勇,实用数学建模与软件应用,西安:西北工业大学出版社,2008.112、姜启源等编,数学模型,第 3 版,北京:高等教育出版社,2003.83、ht
20、tp:/ 1:LINGO 程序如下:Max=0.043*x1*y1+0.027*x2*y2+0.025*x3*y3+0.022*x4*y4+0.045*x5*y5;x2+x3+x4=400;x1+x2+x3+x4+x5=0;x2=0;x3=0;x4=0;x5=0;BIN(y1);BIN(y2);BIN(y3);BIN(y4);BIN(y5);运算结果:Rows= 10 Vars= 10 No. integer vars= 5Nonlinear rows= 1 Nonlinear vars= 10 Nonlinear constraints= 0Nonzeros= 35 Constraint n
21、onz= 23 Density=0.318No. : 6, Obj=MAX Single cols= 5Local optimal solution found at step: 6Objective value: 29.83636Branch count: 0Variable Value Reduced CostX1 218.1818 0.0000000Y1 1.000000 2.045455X2 0.0000000 0.0000000Y2 1.000000 0.0000000X3 736.3636 0.0000000Y3 1.000000 0.0000000X4 0.0000000 0.0
22、000000Y4 1.000000 0.0000000X5 45.45455 0.0000000Y5 1.000000 0.0000000Row Slack or Surplus Dual Price1 29.83636 0.6181818E-032 336.3636 0.00000003 0.0000000 0.00000004 -0.4547474E-12 0.00000005 -0.2273737E-12 0.0000000116 218.1818 0.00000007 0.0000000 0.00000008 736.3636 0.00000009 0.0000000 0.000000
23、010 45.45455 0.0000000对于模型 2:LINGO 程序如下:Max=0.043*x1+0.027*x2+0.025*x3+0.022*x4+0.045*x5-t*f*0.0275;x2+x3+x4=400;x1+x2+x3+x4+x5=0;x2=0;x3=0;x4=0;x5=0;f : 1, Obj=MAX, GUBs =400;x1+x2+x3+x4+x5=0;x2=0;x3=0;x4=0;x5=0;运算结果:Rows= 5 Vars= 5 No. integer vars= 0 ( all are linear)Nonzeros= 25 Constraint nonz=
24、 18( 9 are +- 1) DensitZ=0.833Smallest and largest elements in abs value= 0.220000E-01 1000.00No. : 1, Obj=MAX, GUBs =400;x1+x2+x3+x4+x5=0;x2=0;x3=0;x4=0;x5=0;13运算结果:Rows= 5 Vars= 5 No. integer vars= 0 ( all are linear)Nonzeros= 25 Constraint nonz= 18( 9 are +- 1) DensitZ=0.833Smallest and largest e
25、lements in abs value= 0.220000E-01 1000.00No. : 1, Obj=MAX, GUBs = 1Single cols= 0Global optimal solution found at step: 8Objective value: 29.42400Variable Value Reduced CostX1 336.0000 0.0000000X2 0.0000000 0.3064000E-01X3 0.0000000 0.4400005E-03X4 648.0000 0.0000000X5 16.00000 0.0000000Row Slack or Surplus Dual Price1 29.42400 1.0000002 248.0000 0.00000003 0.0000000 0.2942400E-014 0.0000000 0.6360000E-035 0.0000000 0.2440000E-02