哈工大数学建模2013考查题

1、某银行经理计划用一笔资金进行证券投资业务，可供购进的证券及其相应信息如下表所示，且有如下规定和限制：（1）市政证券的收益可以免税，其它证券的收益需要按50%的税率纳税；（2）政府及代办机构的证券总共至少购进400万元；（3）所购证券的平均信用等级不超过1.4（信用等级越小，信用程度越高）；（4）所购证券的平均到期年限不超过5年；证券名称证券种类信用等级到期年限到期税前收益率(%)A市政294.3B代办机构2155.4C政府145.0D政府134.4E市政524.5请回答下列问题：（1）若该经理有1000万资金，应如何投资？（2）如果能够以2.75%的利率借到不超过100万元，该经理应该如何操作？（3）在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C的税前收益减少为4.8%，投资应否改变？注：为简化问题起见，题中的税前收益率和利率都与年限无关，即都为固定值。问题分析这个问题主要是为了实现收益最大化，通过合理投资证券实现，受到信用等级、到期年限、到期税前收益率以及每种证劵的纳税税率的约束。决策变量：设每种证劵分别投资x1、x2、x3、x4、x5（万元），平均信用等级为u，平均到期年限为v。目标函数：设投资总金额为Q，投资的利润为y（万元），则有y=0.043x1+0.027x2+0.025x3+0.022x4+0.045x5约束条件：平均信用等级u=(2×x1+2×x2+x3+x4+5×x5)/Q≤1.4平均到期年限v=(9×x1+15×x2+4×x3+3×x4+2×x5)/Q≤5非负约束xi=0,j=1,2,3,4,5附加约束x1+x2+x3≥400模型分析与假设xi对目标函数的“贡献”与xi取值成正比xi对约束条件的“贡献”与xi取值成正比xi对目标函数的“贡献”与xj取值无关xi对约束条件的“贡献”与xj取值无关xi取值连续变量xi符合比例性、可加性、连续性。模型求解目标函数y=0.043x1+0.027x2+0.025x3+0.022x4+0.045x5约束条件x1+x2+x3+x4+x5=1000=Qx2+x3+x4≥4002×x1+2×x2+x3+x4+5×x5≤1.4Q=14009×x1+15×x2+4×x3+3×x4+2×x5≤5Q=50000≤x1≤10000≤x2≤10000≤x3≤10000≤x4≤10000≤x5≤1000模型求解用LINGO软件求解，程序如下max=0.043*x1+0.027*x2+0.025*x3+0.022*x4+0.045*x5;x1+x2+x3+x4+x5=1000;x2+x3+x4=400;0.6*x1+0.6*x2-0.4*x3-0.4*x4+3.6*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;运行结果Globaloptimalsolutionfound.Objectivevalue:29.83636Totalsolveriterations:4VariableValueReducedCostX1218.18180.000000X20.0000000.3018182E-01X3736.36360.000000X40.0000000.6363636E-03X545.454550.000000RowSlackorSurplusDualPrice129.836361.00000020.0000000.2983636E-013336.36360.00000040.0000000.6181818E-0250.0000000.2363636E-02从程序运行的结果中可以得到最优解为：x1=218.1818，x2=0，x3=736.3636，x4=0，x5=45.45455,y=29.83636。即对A中证劵投资218.1818万元，对C中证劵投资736.3636万元，对E中证劵投资45.45455万元，对B和D中证劵不进行投资，此种投资方案可使投资利润最大为29.83636万元。（2）从问题一程序运行的结果中可以看到，程序的第二条所对应的影子价格为2.983636％，大于2.75％，所以应该借足100万元进行投资。用LINGO软件求解，程序如下max=0.043*x1+0.027*x2+0.025*x3+0.022*x4+0.045*x5;x1+x2+x3+x4+x5=1100;x2+x3+x4=400;0.6*x1+0.6*x2-0.4*x3-0.4*x4+3.6*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;运行结果：Globaloptimalsolutionfound.Objectivevalue:32.82000Totalsolveriterations:4VariableValueReducedCostCostX1240.00000.0000000.000000X20.0000000.3018182E-010.3018182E-01X3810.00000.0000000.000000X40.0000000.6363636E-030.6363636E-03X550.000000.0000000.000000RowSlackorSurplusDualPrice132.820001.0000001.00000020.0000000.2983636E-010.2983636E-013410.00000.0000000.00000040.0000000.6181818E-020.6181818E-0250.0000000.2363636E-020.2363636E-02从程序运行的结果中可以得到最优解为：x1=240，x2=0，x3=810，x4=0，x5=50,y=32.82。也就是对A中证劵投资240万元，对C中证劵投资810万元，对E中证劵投资50万元，对B和D中证劵不进行投资，此种投资方案可使投利润最大为32.82万元。（3）若证券A的税前收益增加为4.5%时，程序如下：max=0.045*x1+0.027*x2+0.025*x3+0.022*x4+0.045*x5;x1+x2+x3+x4+x5=1000;x2+x3+x4=400;0.6*x1+0.6*x2-0.4*x3-0.4*x4+3.6*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;运行结果如下：Globaloptimalsolutionfound.Objectivevalue:30.27273Totalsolveriterations:4VariableValueReducedCostX1218.18180.000000X20.0000000.3436364E-01X3736.36360.000000X40.0000000.2727273E-03X545.454550.000000RowSlackorSurplusDualPrice130.272731.00000020.0000000.3027273E-013336.36360.00000040.0000000.6363636E-0250.0000000.2727273E-02投资结果不改变。若证券C的税前收益减少为4.8%，则程序如下：max=0.045*x1+0.027*x2+0.024*x3+0.022*x4+0.045*x5;x1+x2+x3+x4+x5=1000;x2+x3+x4=400;0.6*x1+0.6*x2-0.4*x3-0.4*x4+3.6*x5=0;4*x1+10*x2-x3-2*x4-3*x5=0;运行结果如下：Globaloptimalsolutionfound.Objectivevalue:30.09600Totalsolveriterations:5VariableValueReducedCostX1336.00000.000000X20.0000000.3456000E-01X30.0000000.7600000E-03X4648.00000.000000X516.000000.000000RowSlackorSurplusDualPrice130.096001.00000020.0000000.3009600E-013248.00000.00000040.0000000.6440000E-0250.0000000.2760000E-02投资应改变。3、某公司有三个加工厂，某月公司接到4份分别来自于A，B，C，D客户的订单依次为15,17,22,12件，三个加工厂开工费用（固定值，与生产量无关）和生产量及各加工厂的至各个客户的每件运输费用如下面两个表格所示，请制定生产和运输费用，使得总费用最小。单件运输费用（元）开工费用（元）生产量第一加工厂9139第二加工厂7035第三家工厂2431问题分析：这个问题是为了实现生产与运输费用总和最少，决策需考虑每个加工厂的开工费用、每个工厂生产量、每个工厂到每个客户单件运输费用和每个客户的订单量这四个因素。基本模型决策变量：设A、B、C、D四个客户分别给三个工厂的分配变量分别为x11、x12、x13；x21、x22、x23；x31、x32、x33；x41、x42、x43。目标函数：设开工费用为y1（元），运输费为y2（元），生产与运输总费用为设生产与运输总费用为y（元），y=y1+y2。约束条件：1.工厂生产量：x11+x21+x311+x41≤39，x12+x22+x32+x42≤35，x13+x23+x33+x43≤31。2.总订货量为15+17+22+12=66件，又由生产量可知至少有两家工厂开工。3.总费用y2=6×x11+4×x12+8×x13+2×x21+9×x22+8×x23+6×x31+5×x32+x33+7×x41+3×x42+5×x43。4.订单量x11+x12+x13=15，x21+x22+x23=17，x31+x32+x33=22，x41+x42+x43=12。5.结果为正整数。模型求解方案：一．只开两个工厂1、只开工厂一、工厂二则有：x13=x23=x33=x43=0x11+x12=15x21+x22=17x31+x32=22客户A客户B客户C客户D第一加工厂6267第二加工厂4953第三家工厂8815x41+x42=12x11+x21+x31+x41≤39x12+x22+x32+x42≤35y1=91+70=161y2=6×x11+4×x12+2×x21+9×x22+6×x31+5×x32+7×x41+3×x42将上述的y2变为如下式子y2=6×x11+4×(15-x11+2×x21+9×(17-x21)+6×x31+5×(22-x31)+7×x41+3×(12-x41)=2×x11-7×x21+x31+4×x41+359即变为如下条件y=y1+y2=2×x11-7×x21+x31+4×x41+520x11+x21+x31+x41≤39x11+x21+x31+x41≥310≤x11≤150≤x21≤170≤x31≤220≤x41≤12模型求解：程序如下：min=2*x11-7*x21+x31+4*x41+520;x11+x21+x31+x41=39;x11+x21+x31+x41=31;x11=15;x21=17;x31=22;x41=12;运行结果如下：Globaloptimalsolutionfound.Objectivevalue:415.0000Totalsolveriterations:1VariableValueReducedCostX110.0000001.000000X2117.000000.000000X3114.000000.000000X410.0000003.000000RowSlackorSurplusDualPrice1415.0000-1.00000028.0000000.00000030.000000-1.000000415.000000.00000050.0000008.00000068.00