《数学建模与数据学实验》课件第9章.ppt

图9.21二次规划（例9.8）的输入模型输入中的第1行（目标函数）只用于给出模型中相应变量的出现顺序：X、Y、RT、ONE、UL，用加号连接。输入中的第2、3行约束是在实际约束前增加的有关变量的一阶最优条件，即Lagrange函数：

3x2+y2－xy－RT(1.2x+0.9y－1.1)+ONE(x+y－1)+UL(y－0.7)

(9.43)

分别对x、y求偏导数，令其大于0可得第2、3行约束（一阶最优条件）。

“end”后面的语句“QCP4”表示原来的二次规划问题真正的约束是从输入的第4行开始的。求解结果如图9.22所示。图9.22二次规划（例9.8）的输出结果这样经过5次迭代，就得到了最优解x=0.666667，y=0.333333，最优值为1.222222。同整数规划时的情况类似，二次规划也没有敏感性分析结果，因此LINDO对应的敏感性分析输出结果也没有意义，不能被利用。【例9.9】求解如下二次规划问题：我们需要用对偶变量（Lagrange乘子）写出其Lagrange函数，设两个“≤”约束的Lagrange乘子分别为LAG1、LAG2，则其Lagrange函数为图9.23二次规划（例9.9）的输入模型最后得到的结果见图9.24，最优整数解x=(0，100)，最小值为-17700。

可见，LINDO用于求解二次规划（QP）问题比较复杂，不容易非常直观地理解。特别要注意的是，原约束是“≥”约束，则对应的Lagrange乘子应该取0或负值（对最小化问题）。

对于二次规划问题，建议直接用LINGO软件进行求解，因为LINGO软件中的输入模型更接近于二次规划的数学表达式，不易出错，而计算效果同样很好。注意：要想学好和灵活应用LINDO软件，首先要多练习使用LINDO来解决问题，熟能生巧。LINDO中的显示报告完全是英文，大家要熟悉其含义，不要太拘泥于书本或他人的方法，要学会举一反三，综合使用，才能用得巧而精。图9.24二次规划（例9.9）的输出结果“DUALPRICES（对偶价格）”表示当对应约束有微小变动时，目标函数的变化率。输出结果中对应于每一个约束有一个对偶价格，若其值为p，表示对应约束中不等式右端项若增加1个，目标函数将增加p个（对max型问题）。显然对于在最优解处约束正好取等号（“紧约束”，即起作用约束）时，对偶价格值才可能不是0。本例中第3）行对偶价格值为2.5，表示当紧约束2A+2B+1.5C=110变为2A+2B+1.5C=111时，目标函数值=1175+2.5=1177.5。对第4）行也可作类似解释。对于非紧约束（本例中的第2）行），DUALPRICES的值为0，表示对应约束中不等式右端项的微小变动不影响目标函数。有时，通过分析DUALPRICES，可以对产生不可行问题的原因有所了解。

图9.10后半部分的输出结果是敏感性分析结果（若在原来求解时没有要求LINDO作敏感性分析，可直接用菜单命令“Reports/Range”）。敏感性分析的作用是给出“RANGESINWHICHTHEBASISISUNCHANGED”，即研究当目标函数的系数和约束右端项在什么范围变化时（此时假定其他系数保持不变），最优基（矩阵）保持不变（报告中的INFINITY表示正无穷），这包括两方面的敏感性分析内容：（1）目标函数中系数变化的范围（OBJCOEFFICIENTRANGES）。

如本例中，目标函数中A变量当前的系数（CURRENTCOEF）为35、允许增加（ALLOWABLEINCREASE）为5、允许减少（ALLOWABLEDECREASE）为15，说明当这个系数在［35-15，35+5］=［20，40］范围变化时，最优基保持不变。对B、C变量，可以作类似解释。由于此时约束没有变化(只是目标函数中某个系数发生了变化)，所以最优基保持不变的意思也就是最优解不变（由于目标函数中的系数发生了变化，所以最优值也会发生变化）。（2）约束右端项变化的范围（RIGHTHANDSIDERANGES）。

本例中第2）行约束中当前右端项（CURRENTRHS）为180、允许增加（ALLOWABLEINCREASE）为INFI