Lingo练习题.docVIP

Lingo软件练习题 习题一 分析题目可以知道这是一道简单的运用线性规划来求解最优值的题目。目标函数是工厂的收益最大，约束条件为生产车间的每个加工区域生产时间的限制。题中所提到的晶体管与模块质量控制区域在一个时间内只能加工微型模块或者晶体管，不能同时加工。基于上述分析和假设我们建立如下模型： 符号说明： 1.表示生产的晶体管、微型模块和电路集成器的数量，i=1,2,3 2.表示出售的晶体管、微型模块和电路集成器的数量，i=1,2,3 3.表示一个第i个产品需要占用滴j个加工区域的时间，i=1,2,3, j=1,2,3,4,分别对应晶体管生产线、电路印刷与组装、晶体管与模块质量控制、电路集成器测试与包装。 4..、均为非负整数。 运用lingo编程代码如下： max=2*(x1-3*(x2+x3))+8*(x2-3*x3)+25*x3-(0.7*x1+0.5*x2+2*x3); 0.1*x1=200; 0.5*x1=200; 0.4*x2=200; 0.1*x3=200; 0.5*x3=200; x1-3*(x2+x3)=0; x2-3*x3=0; 0.5*x1+0.4*x2=200; @gin(x1); @gin(x2); @gin(x3); 运行得到最优解如下： Max=568.3 X1=316 X2=105 X3=0 Variable Value Reduced Cost X1 316.0000 -1.300000 X2 105.0000 -1.500000 X3 0.000000 7.000000 送上述结果可以看出，当x3的值增大1的时候，max值就减小7，接下来分析一下7的来源： 一个电路集成器由3个晶体管加上3个微型模块再外加2元成本，售价25元，忽略晶体管和微型模块的成本其利润为23元；而3个的晶体管和3个微型模块单独出售一共为3*2+3*8=30元，忽略晶体管和微型模块成本，30元即为其利润，两者利润相差7元，若都算上晶体管和微型模块的成本，利润差还是7元。这也解释了为什么在这样的约束条件下，最优解中电路集成器产量为0 习题二 查阅资料可以知道：路灯光照强度：I=k*(P*sinα/R^2)。 P处光照强度：+k*3000* sin=5/, sin 要求P处光强最大或最小，将代入可得： k为关照强度系数，为便于计算我们假设k取1,此假设对本题所有方程都成立。 对于式我们用lingo解出其最大最小值，程序如下： 最小值： min=2000*5/@sqrt((5^2+x^2)^3)+3000*6/@sqrt((6^2+(20-x)^2)^3); x=0; x=20; 最大值： max=2000*5/@sqrt((5^2+x^2)^3)+3000*6/@sqrt((6^2+(20-x)^2)^3); x=0; x=20; 运行程序得到如下结果： 最小值解为：x=9.338,=18.24 最大值解为：x=19.977,=84.48 当3kw路灯高度可以从3m变到6m时，我们设3kw路灯高度为h1，此时可以得到方程 分别对x、h1求偏导,令其值为零： 运用mathematica进行求解，联立两式可以得到多组解，其中一组就是我们所要寻找的解，即最暗点的最亮光强，此时x=9.503 h1=7.422, 进一步可以算得Ip=18.56 同理，当2kw路灯高度也可以在3m到9m之间变化时，我们设2kw路灯高度h2，此时可以得到方程： 分别对x、h1、h2求偏导，令其值为零： 联立三式，用mathematica解得多组解，找出我们所需要的解如下： x=9.325 h1=7.548 h2=6.594 进一步可求得Ip=18.99 习题三 本题要求在满足空投条件下的最小费用，那么我们就先要知道对于不同半径r的降落伞的成本和它的最大载重。先计算最大载重： 根据物理知识可以知道，降落伞在降落过程中可以分成两个部分：先做加速运动，后受到的阻力增大，做匀速运动，即在运动一定时间后：G=mg=f阻 题中认为降落伞收到的空气阻力与坠落速度和伞面积的乘积成正比，我们不妨设f阻=kv*2πr2 联立两式可得：mg= kv*2πr2 ,降落伞匀速运动时