数学建模软件lingo示例PPT课件

1、主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-1lindo /lingo软件软件vlingo: linear interactive general optimizer (v10.0)v求解数学规划问题 min z = f (x) s.t x a ( rn )v需要掌握的几个重要方面基本语法结构掌握集合(sets)的应用正确阅读求解报告正确理解求解状态窗口学会设置基本的求解选项(options)主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-2一、一、lingo基本语法基本语法1、定义了目标函数为min=. 2、以一个分号“；”结尾 除sets

2、, endsets, data , enddata, end之外3、可以放在约束条件的右端，同时数字也可放在约束条件的左端。4、假定各变量非负。5、注释：“！”6、为、主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-3建模时需要注意的几个基本问题建模时需要注意的几个基本问题 v尽量使用实数优化，减少整数约束和整数变量v尽量使用光滑优化，减少非光滑约束的个数 v尽量少使用绝对值、符号函数、多个变量求最大/最小值、四舍五入、取整函数等v尽量使用线性模型，减少非线性约束和非线性变量的个数如x/y 5 改为x5yv合理设定变量上下界，尽可能给出变量初始值 v模型中使用的参数数量

3、级要适当如小于103主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-4二、二、lingo集合集合v某公司有6个建筑工地，位置坐标为(ai,bi) (单位：公里),水泥日用量di (单位：吨）例1 选址问题v假设：料场和工地之间有直线道路v（1）现有2料场，位于a(5,1),b(2,7),记(xj,yj),j=1,2, 日储量ej各有20吨。v目标：制定每天的供应计划，即从a, b两料场分别向各工地运送多少吨水泥，使总的吨公里数最小。ia8.750.55.7537.25b0.754.7556.57.75d1.25547611主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模

4、课件module 00-5解：v决策变量：料场j到工地i的运量 cij 12维v线性规划模型vlingo表达式？2622 1/2112161min()() .,1,.,6,1,20,1,.,6,1,2ijjijijiijijijjiijcxaybs tcdicejcijlocation目标：吨公里约束：需求 供应主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-6基本使用1、lingo建模语言建模语言构成：构成：4个段个段v目标与约束段v集合段（sets endsets）v数据段（data enddata）v初始段（init endinit）v（计算段 (calc end

5、calc)）主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-7目标与约束段目标与约束段model:title location problem;sets: demand/1.6/:a,b,d; supply/1.2/:x,y,e; link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25; b=1.25,0.75,4.75,5,6.5,7.75; d=3,5,4,7,6,11; e=20,20; x,y=5,1,2,7; enddatainit:endinitm i n = s u m ( l i n k

6、 ( i , j ) : c ( i , j ) * ( ( x ( j ) - a ( i ) ) 2 + ( y ( j ) -b(i)2)(1/2);for(demand(i):sum(supply(j):c(i,j)=d(i);); for(supply(i):sum(demand(j):c(j,i)=e(i););for(supply: free(x); free(y); );end集合段集合段数据段数据段初始段初始段location供应约束供应约束需求需求需求点的位置需求点的位置供需量供需量供应供应初始点初始点目标目标需求约束需求约束连接连接主讲人：孙云龙主讲人：孙云龙数学建模课件

7、数学建模课件module 00-82、lingo的集合的集合set及其属性及其属性attributev例1v建筑工地位置坐标 (ai,bi) 、水泥日用量di ：对每个建筑工地（6个）都有一个对应的值都是一个由6个元素组成的数组是已知的v料场位置坐标 (xj,yj) 、日储量ej对每个料场（2个）都有一个对应的值都是一个由2个元素组成的数组目前是已知的v料场到建筑工地的供应计划 c i j 对每个料场与建筑工地之间（62）都有一个对应的值是一个62 个元素组成的矩阵是未知数vlindo无数组，每个变量输入麻烦lingo的集合set下标集合100个工地？主讲人：孙云龙主讲人：孙云龙数学建模课件数

8、学建模课件module 00-9集合set及其属性attributev定义数组下标集合demand/1.6/表示6个建筑工地a,b,d称为该集合的属性表示坐标(ai,bi) 、水泥日用量div定义数组下标集合supply/1.2/表示6个建筑工地该集合的属性x,y,e表示坐标(xj,yj) 、日储量ej v定义数组下标集合link(demand,supply)表示62个料场到建筑工地的连接该集合的属性c表示每个料场与建筑工地之间供应计划c i j1到6的整数主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-10v建立下标集合lingo 建模语言集合段 数据段需求点的位置需

9、求点的位置供需量供需量sets: demand/1.6/:a,b,d; supply/1.2/:x,y,e; link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25; b=1.25,0.75,4.75,5,6.5,7.75; d=3,5,4,7,6,11; e=20,20;x,y=5,1,2,7;enddatav赋值需求需求供应供应连接连接主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-113、定义集合定义集合setv直接把元素列举出来v定义格式 集合名 元素列表 属性列表 setname /membe

10、r_list/ : attribute_list; 可选项v元素列表显式列举法列出全部元素, 用逗号或空格分开隐式列举法 1.nv属性列表缺省集合可在程序中作为一循环变量使用，构造更复杂的派生集合v元素列表缺省必须在数据段给出元素列表赋值基本集合基本集合primary set主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-12派生集合派生集合derived setv基于其它集合而派生出来的二维或多维集合v定义格式 集合名 父集合列表 元素列表 属性列表 setname(parent_set_list) /member_list/ : attribute_list;v元素

11、列表缺省所有组合稠密集合稠密集合、或数据段列表赋值v元素列表稀疏集合稀疏集合元素列表法枚举元素过滤法利用过滤条件setname(parent_set_list) |filtrate_condition :attribute_list;主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-13v建立 下标集合例1需求点的位置需求点的位置供需量供需量sets:demand/1.6/:a,b,d;supply/1.2/:x,y,e;link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25; b=1.25,0.75,4

12、.75,5,6.5,7.75; d=3,5,4,7,6,11; e=20,20;x,y=5,1,2,7;enddatav赋值需求需求供应供应连接连接基本基本集合集合派生派生集合集合626152514241323122211211212121654321654321654321,cccccccccccceeyyxxddddddbbbbbbaaaaaa主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-14集合的类型 集合集合 派生集合派生集合 基本集合基本集合 稀疏集合稀疏集合 稠密集合稠密集合 元素列表法元素列表法 元素过滤法元素过滤法 直接列举法直接列举法 隐式列举法隐式

13、列举法setname /member_list/ : attribute_list;setname(parent_set_list) /member_list/ : attribute_list;sets: cities /a1,a2,a3,b1,b2/; roads(cities, cities)/ a1,b1 a1,b2 a2,b1 a3,b2/:d; endsetssets: students /s1.s8/; pairs( students, students) | &2 #gt# &1: benefit, match;endsets主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件

14、module 00-15三、三、lingo运算符和函数运算符和函数v过滤条件使用结果：0 11、运算符及其优先级、运算符及其优先级逻辑运算符#not#否定#eq#相等#ne#不等#and#并且#or#或者#gt#大于#ge#大于等于#lt#小于#le#小于等于算术运算符+ - * / 关系运算符(=)主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-16v运算符的优先级最高最低#not#(负号)* /+(减法)#eq# #ne# #gt# #ge# #lt# #le#and# #or#v先左后右v先括号内，后括号外主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件mod

15、ule 00-17vlingo内部函数”2、lingo基本数学函数基本数学函数abs( x)sin( x)cos( x)tan( x)log( x)exp( x)smax( list )smin( list )sign(x)floor ( x)主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-183、lingo集合循环函数集合循环函数v难点！重点！v循环操作函数集合上的元素下标：集合函数名function(setname (set_index_list)|condition:expression_list);集合名集合索引列表过滤条件表达式formaxminsum pro

16、d主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-19v目标：吨公里for(demand(i):sum(supply(j):c(i,j)=d(i);); for(supply(i):sum(demand(j):c(j,i)=e(i););例12622 1/211min()() ijjijijicxaybv约束：需求 供应2161.,1,.,6,1,2ijijijjis tcdicejmin=sum(link(i,j):c(i,j)*(x(j)-a(i)2+(y(j)- b(i)2)(1/2);主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-20v

17、默认：？4、lingo变量定界函数变量定界函数bnd(l, x,n)l，ufree( x)自由变量bin( x)0-1变量gin( x)整数变量例1v约束：非负for(supply: free(x); free(y); );.0,1,.,6,1,2ijs tcij主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-214、lingo条件判断函数条件判断函数v分段函数if(logical_condition, true_result, false_result)vokv例1主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-22目标与约束段目标与约束段mo

18、del:title location problem;sets: demand/1.6/:a,b,d; supply/1.2/:x,y,e; link(demand,supply):c;endsetsdata:a=1.25,8.75,0.5,5.75,3,7.25; b=1.25,0.75,4.75,5,6.5,7.75; d=3,5,4,7,6,11; e=20,20;x,y=5,1,2,7; enddatainit:endinitm i n = s u m ( l i n k ( i , j ) : c ( i , j ) * ( ( x ( j ) - a ( i ) ) 2 + ( y

19、 ( j ) -b(i)2)(1/2);for(demand(i):sum(supply(j):c(i,j)=d(i);); for(supply(i):sum(demand(j):c(j,i)=e(i););for(supply: free(x); free(y); );end集合段集合段数据段数据段初始段初始段location供应约束供应约束需求需求需求点的位置需求点的位置供需量供需量供应供应初始点初始点目标目标需求约束需求约束连接连接主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-23v2）改建两个新料场，需要确定新料场位置(xj,yj)和运量cij ，在其它条件

20、不变下使总吨公里数最小。决策变量：决策变量：ci j，(xj,yj)16维维非线性规划模型非线性规划模型例1 选址问题2622 1/2112161min()() .,1,.,6,1,20,1,.,6,1,2ijjijijiijijijjiijcxaybs tcdicejcij目标：吨公里约束：需求 供应主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-24lingo集合段集合段数据段数据段初始段初始段目标与目标与约束段约束段 最优：最优：89.8835(吨公里吨公里 ) lp：移到数据段移到数据段location局部局部主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-25边界边界主讲人：孙云龙主讲人：孙云龙数学建模课件数学建模课件module 00-26例2 连续投资问题v设有一笔资金m=10万，未来5年内可以投资4个项目。其中项目1要到第二年末才可回收资金，本利为115%；项目2只能在第三年初投资，到第五年末回收本利125%,但不超过3万；项目3在第二年初投资，第五年末回收本利140%，但不超过4万；项目4每年初投资，年末回收本利106%；试确定5年内如何安排投资？年一二三四五限额回收时间本利项目1第二年末115%项目23万第五年末125%项目34万第五年末140%项目4年末106%主讲人：孙云龙主讲