培训教案——LINGO课件.pptVIP

LINGO实战——整数规划 * LINGO实战——整数规划 * LINGO实战——整数规划 下面是LINGO程序。 MODEL: DATA: M=150; ENDDATA max=3*x1+4*x2+8*x3-100*y1-150*y2-200*y3;!目标函数; 2*x1+4*x2+8*x3=500; 2*x1+3*x2+4*x3=300; x1+2*x2+3*x3=100; 3*x1+5*x2+7*x3=700; x1=M*y1; x2=M*y2; x3=M*y3; @GIN(x1);@GIN(x2);@GIN(x3); !指定产品件数为整数; @BIN(y1);@BIN(y2);@BIN(y3); !指定0-1变量; end 得到的解为x1=100,x2=0,x3=0,y1=1,y2=0,y3=0。最大值为Z=200元。 * LINGO——整数规划 公司在各地有4项业务，选定了4位业务员去处理。由于业务能力、经验和其它情况不同，4业务员去处理4项业务的费用（单位：元）各不相同，见下表： 业务 业务员 1 2 3 4 1 1100 800 1000 700 2 600 500 300 800 3 400 800 1000 900 4 1100 1000 500 700 应当怎样分派任务，才能使总的费用最小？ * LINGO实战——整数规划 * 程序结构——数学算符 二、数学函数 LINGO提供了大量的标准数学函数： @abs(x) 返回x的绝对值 @sin(x) 返回x的正弦值，x采用弧度制 @cos(x) 返回x的余弦值 @tan(x) 返回x的正切值 @exp(x) 返回常数e的x次方 @log(x) 返回x的自然对数 @lgm(x) 返回x的gamma函数的自然对数 @sign(x) 如果x0返回-1；否则，返回1 @floor(x) 返回x的整数部分。当x=0时，返回不超过x的最大整数；当x0时，返回不低于x的最小整数。 @smax(x1,x2,…,xn) 返回x1，x2，…，xn中的最大值 @smin(x1,x2,…,xn) 返回x1，x2，…，xn中的最小值 * 程序结构——应用实例 应用示例： 给定一个直角三角形，求包含该三角形的最小正方形。 解：如图所示。 求最小的正方形就相当于求如下的最优化问题： LINGO代码如下： model: sets: object/1..3/: f; endsets data: a, b = 3, 4; !两个直角边长，修改很方便; enddata f(1) = a * @sin(x); f(2) = b * @cos(x); f(3) = a * @cos(x) + b * @sin(x); min = @smax(f(1),f(2),f(3)); @bnd(0,x,1.57); end A B C D E a b x * 程序结构——概率函数 三、概率函数 1．@pbn(p,n,x)：二项分布的累积分布函数。当n和（或）x不是整数时，用线性插值法进行计算。 2．@pcx(n,x) ：自由度为n的χ2分布的累积分布函数。 3．@peb(a,x) ：当到达负荷为a，服务系统有x个服务器且允许无穷排队时的Erlang繁忙概率。 4．@pel(a,x) ：当到达负荷为a，服务系统有x个服务器且不允许排队时的Erlang繁忙概率。 5．@pfd(n,d,x) ：自由度为n和d的F分布的累积分布函数。 6．@pfs(a,x,c) ：当负荷上限为a，顾客数为c，平行服务器数量为x时，有限源的Poisson服务系统的等待或返修顾客数的期望值。a是顾客数乘以平均服务时间，再除以平均返修时间。当c和（或）x不是整数时，采用线性插值进行计算。 7．@phg(pop,g,n,x) ：超几何（Hypergeometric）分布的累积分布函数。pop表示产品总数，g是正品数。从所有产品中任意取出n（n≤pop）件。pop，g，n和x都可以是非整数，这时采用线性插值进行计算。 8．@ppl(a,x)Poisson ：分布的线性损失函数，即返回max(0,z-x)的期望值，其中随机变量z服从均值为a的Poisson分布。 9．@pps(a,x) ：均值为a的Poisson分布的累积分布函数。当x不是整数时，采用线性插值进行计算。 10．@psl(x) ：单位正态线性损失函数，即返回max(0,z-x)的期望值，其中随机变