數(shù)學(xué)建模軟件lingo示例PPT課件

1、主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-1LINDO /LINGO軟件軟件vLINGO: Linear INteractive General Optimizer (V10.0)v求解數(shù)學(xué)規(guī)劃問題 Min Z = f (x) s.t x A ( Rn )v需要掌握的幾個重要方面基本語法結(jié)構(gòu)掌握集合(SETS)的應(yīng)用正確閱讀求解報告正確理解求解狀態(tài)窗口學(xué)會設(shè)置基本的求解選項(OPTIONS)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-2一、一、Lingo基本語法基本語法1、定義了目標(biāo)函數(shù)為MIN=. 2、以一個分號“；”結(jié)尾 除SETS

2、, ENDSETS, DATA , ENDDATA, END之外3、可以放在約束條件的右端，同時數(shù)字也可放在約束條件的左端。4、假定各變量非負(fù)。5、注釋：“！”6、為、主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-3建模時需要注意的幾個基本問題建模時需要注意的幾個基本問題 v盡量使用實數(shù)優(yōu)化，減少整數(shù)約束和整數(shù)變量v盡量使用光滑優(yōu)化，減少非光滑約束的個數(shù) v盡量少使用絕對值、符號函數(shù)、多個變量求最大/最小值、四舍五入、取整函數(shù)等v盡量使用線性模型，減少非線性約束和非線性變量的個數(shù)如x/y 5 改為x5yv合理設(shè)定變量上下界，盡可能給出變量初始值 v模型中使用的參數(shù)數(shù)量

3、級要適當(dāng)如小于103主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-4二、二、Lingo集合集合v某公司有6個建筑工地，位置坐標(biāo)為(ai,bi) (單位：公里),水泥日用量di (單位：噸）例1 選址問題v假設(shè)：料場和工地之間有直線道路v（1）現(xiàn)有2料場，位于A(5,1),B(2,7),記(xj,yj),j=1,2, 日儲量ej各有20噸。v目標(biāo)：制定每天的供應(yīng)計劃，即從A, B兩料場分別向各工地運(yùn)送多少噸水泥，使總的噸公里數(shù)最小。ia8.750.55.7537.25b0.754.7556.57.75d1.25547611主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模

4、課件Module 00-5解：v決策變量：料場j到工地i的運(yùn)量 cij 12維v線性規(guī)劃模型vlingo表達(dá)式？2622 1/2112161min()() .,1,.,6,1,20,1,.,6,1,2ijjijijiijijijjiijcxaybs tcdicejcijlocation目標(biāo)：噸公里約束：需求 供應(yīng)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-6基本使用1、Lingo建模語言建模語言構(gòu)成：構(gòu)成：4個段個段v目標(biāo)與約束段v集合段（SETS ENDSETS）v數(shù)據(jù)段（DATA ENDDATA）v初始段（INIT ENDINIT）v（計算段 (CALC END

5、CALC)）主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-7目標(biāo)與約束段目標(biāo)與約束段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集合段集合段數(shù)據(jù)段數(shù)據(jù)段初始段初始段location供應(yīng)約束供應(yīng)約束需求需求需求點的位置需求點的位置供需量供需量供應(yīng)供應(yīng)初始點初始點目標(biāo)目標(biāo)需求約束需求約束連接連接主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件

7、數(shù)學(xué)建模課件Module 00-82、Lingo的集合的集合Set及其屬性及其屬性Attributev例1v建筑工地位置坐標(biāo) (ai,bi) 、水泥日用量di ：對每個建筑工地（6個）都有一個對應(yīng)的值都是一個由6個元素組成的數(shù)組是已知的v料場位置坐標(biāo) (xj,yj) 、日儲量ej對每個料場（2個）都有一個對應(yīng)的值都是一個由2個元素組成的數(shù)組目前是已知的v料場到建筑工地的供應(yīng)計劃 c i j 對每個料場與建筑工地之間（62）都有一個對應(yīng)的值是一個62 個元素組成的矩陣是未知數(shù)vLINDO無數(shù)組，每個變量輸入麻煩Lingo的集合Set下標(biāo)集合100個工地？主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)

8、學(xué)建模課件Module 00-9集合Set及其屬性Attributev定義數(shù)組下標(biāo)集合demand/1.6/表示6個建筑工地a,b,d稱為該集合的屬性表示坐標(biāo)(ai,bi) 、水泥日用量div定義數(shù)組下標(biāo)集合supply/1.2/表示6個建筑工地該集合的屬性x,y,e表示坐標(biāo)(xj,yj) 、日儲量ej v定義數(shù)組下標(biāo)集合link(demand,supply)表示62個料場到建筑工地的連接該集合的屬性c表示每個料場與建筑工地之間供應(yīng)計劃c i j1到6的整數(shù)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-10v建立下標(biāo)集合Lingo 建模語言集合段 數(shù)據(jù)段需求點的位置需

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賦值需求需求供應(yīng)供應(yīng)連接連接主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-113、定義集合定義集合Setv直接把元素列舉出來v定義格式 集合名 元素列表 屬性列表 setname /membe

10、r_list/ : attribute_list; 可選項v元素列表顯式列舉法列出全部元素, 用逗號或空格分開隱式列舉法 1.nv屬性列表缺省集合可在程序中作為一循環(huán)變量使用，構(gòu)造更復(fù)雜的派生集合v元素列表缺省必須在數(shù)據(jù)段給出元素列表賦值基本集合基本集合primary set主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-12派生集合派生集合derived setv基于其它集合而派生出來的二維或多維集合v定義格式 集合名 父集合列表 元素列表 屬性列表 setname(parent_set_list) /member_list/ : attribute_list;v元素

11、列表缺省所有組合稠密集合稠密集合、或數(shù)據(jù)段列表賦值v元素列表稀疏集合稀疏集合元素列表法枚舉元素過濾法利用過濾條件setname(parent_set_list) |filtrate_condition :attribute_list;主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-13v建立 下標(biāo)集合例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賦值需求需求供應(yīng)供應(yīng)連接連接基本基本集合集合派生派生集合集合626152514241323122211211212121654321654321654321,cccccccccccceeyyxxddddddbbbbbbaaaaaa主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件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主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模

14、課件數(shù)學(xué)建模課件Module 00-15三、三、Lingo運(yùn)算符和函數(shù)運(yùn)算符和函數(shù)v過濾條件使用結(jié)果：0 11、運(yùn)算符及其優(yōu)先級、運(yùn)算符及其優(yōu)先級邏輯運(yùn)算符#NOT#否定#EQ#相等#NE#不等#AND#并且#OR#或者#GT#大于#GE#大于等于#LT#小于#LE#小于等于算術(shù)運(yùn)算符+ - * / 關(guān)系運(yùn)算符(=)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-16v運(yùn)算符的優(yōu)先級最高最低#NOT#(負(fù)號)* /+(減法)#EQ# #NE# #GT# #GE# #LT# #LE#AND# #OR#v先左后右v先括號內(nèi)，后括號外主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)

15、學(xué)建模課件Module 00-17vLingo內(nèi)部函數(shù)”2、Lingo基本數(shù)學(xué)函數(shù)基本數(shù)學(xué)函數(shù)ABS( X)SIN( X)COS( X)TAN( X)LOG( X)EXP( X)SMAX( list )SMIN( list )SIGN(X)FLOOR ( X)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-183、Lingo集合循環(huán)函數(shù)集合循環(huán)函數(shù)v難點！重點！v循環(huán)操作函數(shù)集合上的元素下標(biāo)：集合函數(shù)名function(setname (set_index_list)|condition:expression_list);集合名集合索引列表過濾條件表達(dá)式FORMAXMI

16、NSUM PROD主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-19v目標(biāo)：噸公里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約束：需求 供應(yīng)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);主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Modul

17、e 00-20v默認(rèn)：？4、Lingo變量定界函數(shù)變量定界函數(shù)BND(L, X,N)L，UFREE( X)自由變量BIN( X)0-1變量GIN( X)整數(shù)變量例1v約束：非負(fù)for(supply: free(X); free(Y); );.0,1,.,6,1,2ijs tcij主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-214、Lingo條件判斷函數(shù)條件判斷函數(shù)v分段函數(shù)IF(logical_condition, true_result, false_result)vOKv例1主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-22目標(biāo)與約束段

18、目標(biāo)與約束段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 ( i , j ) : c ( i , j ) * ( ( x ( j ) - a ( i ) )

19、 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集合段集合段數(shù)據(jù)段數(shù)據(jù)段初始段初始段location供應(yīng)約束供應(yīng)約束需求需求需求點的位置需求點的位置供需量供需量供應(yīng)供應(yīng)初始點初始點目標(biāo)目標(biāo)需求約束需求約束連接連接主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-23v2）改建兩個新料場，需要確定新料場位置(xj,yj)和運(yùn)量ci

20、j ，在其它條件不變下使總噸公里數(shù)最小。決策變量：決策變量：ci j，(xj,yj)16維維非線性規(guī)劃模型非線性規(guī)劃模型例1 選址問題2622 1/2112161min()() .,1,.,6,1,20,1,.,6,1,2ijjijijiijijijjiijcxaybs tcdicejcij目標(biāo)：噸公里約束：需求 供應(yīng)主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-24LINGO集合段集合段數(shù)據(jù)段數(shù)據(jù)段初始段初始段目標(biāo)與目標(biāo)與約束段約束段 最優(yōu)：最優(yōu)：89.8835(噸公里噸公里 ) LP：移到數(shù)據(jù)段移到數(shù)據(jù)段location局部局部主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模

21、課件數(shù)學(xué)建模課件Module 00-25邊界邊界主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 00-26例2 連續(xù)投資問題v設(shè)有一筆資金M=10萬，未來5年內(nèi)可以投資4個項目。其中項目1要到第二年末才可回收資金，本利為115%；項目2只能在第三年初投資，到第五年末回收本利125%,但不超過3萬；項目3在第二年初投資，第五年末回收本利140%，但不超過4萬；項目4每年初投資，年末回收本利106%；試確定5年內(nèi)如何安排投資？年一二三四五限額回收時間本利項目1第二年末115%項目23萬第五年末125%項目34萬第五年末140%項目4年末106%主講人：孫云龍主講人：孫云龍數(shù)學(xué)建模課件數(shù)學(xué)建模課件Module 0