简化周永华Lingo入门讲述.ppt

数学建模四步曲 一、分析问题，建立数学模型 （可能是几个小模型的组合） 二、算法（求解方法） 常用算法（思想）：贪婪算法，最短路算法、分支定界 算法等。 特殊算法：单纯形法，共轭梯度法、变尺度法、制约函 数法、广义既约梯度法等 三、编写程序 SAS----统计软件 LINGO----优化软件 MATLAB---数学计算软件（含部分统计 软件、优化软件功能） C语言 四、写作 建模题型--- A题以统计数据分析为主，优化为辅助；-----概率统计 B题以优化为主，可能还包含数据和微分方程。-----运筹学(数学规划、图论、排队论、对策论、存储论）、微分方程 第一节 LINGO的基本用法 例 解数学规划问题： 目标函数 约束条件 LINGO的语法规定： ⑴求目标函数的最大值和最小值分别用max=…和min=… ⑵每个语句必须以分号“；”结束，每行可以有多个语句，语句可以跨行； ⑶变量名必须以字母开头，由字母、数字、下划线组成，长度不超过32个字符，不区分大小写； ⑷可以给语句加上标号，例如[OBJ] [OBJ] max=200*x1+300*x2; ⑸以！开头，以“；”结束的语句是注释语句； ⑹如果对变量的取值范围没有作特殊说明，则默认所有决策变量都非负； ⑺LINGO模型一语句“model：”开头，以“end”结束，对于比较简单的模型，这两个语句可以省略。 选菜单Lingo|Solve（或按Ctrl+S），或用鼠标点击 按钮，对模型求解。 模型2及求解 注意到这个问题中角度和距离的单位是不一致的(角度为弧度，距离为公里)，因此将这4个误差平方和同等对待(相加)不是很合适。并且，4种设备测量的精度(误差限)不同，而上面的方法根本没有考虑测量误差问题。如何利用测量设备的精度信息？这就需要看你对例中给出的设备精度如何理解。 一种可能的理解是：设备的测量误差是均匀分布的。以VOR1为例，目前测得的角度为161.20，测量精度为0.80，所以实际的角度应该位于区间[161.20-0.80，161.20+0.80]内。对其他设备也可以类似理解。由于很小，即测量精度很高，所以在相应区间内正切函数tan的单调性成立。 于是可以得到一组不等式： (i=1, 2, 3) 也就是说，飞机坐标应该位于上述不等式组成的区域内。例如，模型1中得到的目标函数值很小，显然满足测量精度要求，因此坐标(980.6926, 731.5666)肯定位于这个可行区域内。 由于这里假设设备的测量误差是均匀分布的，所以飞机坐标在这个区域内的每个点上的可能性应该也是一样的，我们最好应该给出这个区域的x和y坐标的最大值和最小值。于是我们可以分别以min x、max x、min y、max y、为目标，以上面的区域限制条件为约束，求出x和y坐标的最大值和最小值。 以min x为例，相应的LINGO程序为： MODEL: Title 飞机定位模型2; SETS: VOR/1..3/: x, y, cita, sigma; ENDSETS DATA: x, y, cita, sigma = 746 1393 2.81347 0.0140 629 375 0.78714 0.0105 1571 259 5.39307 0.0227; x4 y4 d4 sigma4 = 155 987 864.3 2.0; ENDDATA ! XX,YY表示飞机坐标; min=xx; @for(VOR: (xx-x)/(yy-y) @tan(cita - sigma) ); @for(VOR: (xx-x)/(yy-y) @tan(cita + sigma) ); d4 - sigma4 @sqrt(@sqr(xx-x4)+@sqr(yy-y4)) ; d4 + sigma4 @sqrt(@sqr(xx-x4)+@sqr(yy-y4)) ; END 用LINGO求解上述模型，LINGO系统返回的信息是这个模型没有可行解。其实这显然是一个不正确的信息，可能只是由于求解空间太大，LINGO没有找到可行解。其实，我们可以想象这个问题的可行解大致就该在模型1中得到的最优解附近，因此可以把这个解作为初始值告诉LINGO。例如，在上面程序中增加以下三行： INIT: xx, yy = 980.6926, 731.