學(xué)建模實(shí)戰(zhàn)：用LINGO求解運(yùn)輸與投資優(yōu)化問(wèn)題)
1. 項(xiàng)目概述當(dāng)數(shù)學(xué)建模遇上LINGO最近在帶學(xué)生做數(shù)學(xué)建模競(jìng)賽的培訓(xùn)發(fā)現(xiàn)很多同學(xué)一看到“運(yùn)輸方案優(yōu)化”或者“連續(xù)投資決策”這類(lèi)題目就有點(diǎn)發(fā)怵。題目描述往往是一大堆數(shù)據(jù)表格要求你給出成本最低的運(yùn)輸計(jì)劃或者未來(lái)幾年收益最大的投資策略。思路大家其實(shí)都有列方程、設(shè)變量嘛但真到了求解那一步尤其是變量一多、約束一復(fù)雜手工計(jì)算幾乎不可能用通用編程語(yǔ)言比如Python從頭寫(xiě)算法又容易在調(diào)試上耗費(fèi)大量時(shí)間。這時(shí)候一個(gè)專(zhuān)業(yè)的工具就能讓你事半功倍——它就是LINGO。簡(jiǎn)單來(lái)說(shuō)LINGO是一款專(zhuān)門(mén)用于求解線(xiàn)性規(guī)劃、非線(xiàn)性規(guī)劃、整數(shù)規(guī)劃等優(yōu)化問(wèn)題的軟件。你可以把它理解為一個(gè)“數(shù)學(xué)規(guī)劃模型的翻譯器和求解器”。你不需要去深究背后的單純形法、分支定界法是什么原理你只需要用接近數(shù)學(xué)語(yǔ)言的格式把目標(biāo)函數(shù)、決策變量和約束條件“告訴”它它就能快速給你算出最優(yōu)解。對(duì)于“運(yùn)輸方案”和“連續(xù)投資”這類(lèi)典型的運(yùn)籌學(xué)問(wèn)題LINGO幾乎是量身定做的解決方案。本文我就結(jié)合這兩個(gè)經(jīng)典場(chǎng)景帶你從零開(kāi)始掌握用LINGO將數(shù)學(xué)模型轉(zhuǎn)化為實(shí)際解決方案的全過(guò)程并分享一些我多年輔導(dǎo)中總結(jié)的、在官方教程里未必會(huì)寫(xiě)的實(shí)戰(zhàn)技巧和避坑指南。2. 核心思路拆解從問(wèn)題到LINGO模型在打開(kāi)LINGO軟件之前最關(guān)鍵的一步是把現(xiàn)實(shí)問(wèn)題抽象成一個(gè)清晰的數(shù)學(xué)模型。很多新手會(huì)直接跳進(jìn)軟件里敲代碼結(jié)果發(fā)現(xiàn)模型本身有邏輯漏洞求解自然失敗。我們分別看看這兩個(gè)問(wèn)題的建模核心。2.1 運(yùn)輸問(wèn)題成本最小化的網(wǎng)絡(luò)流運(yùn)輸問(wèn)題的背景非常直觀你有若干個(gè)工廠(chǎng)產(chǎn)地每個(gè)工廠(chǎng)有固定的產(chǎn)量供給量你也有若干個(gè)銷(xiāo)售點(diǎn)銷(xiāo)地每個(gè)銷(xiāo)售點(diǎn)有固定的需求量。從每個(gè)工廠(chǎng)到每個(gè)銷(xiāo)售點(diǎn)都有單位運(yùn)輸成本。目標(biāo)是如何安排運(yùn)輸量在滿(mǎn)足所有供給和需求的前提下讓總運(yùn)輸成本最低。這本質(zhì)上是一個(gè)線(xiàn)性規(guī)劃問(wèn)題。它的核心結(jié)構(gòu)非常規(guī)整決策變量通常設(shè)為X_ij表示從工廠(chǎng)i運(yùn)到銷(xiāo)售點(diǎn)j的貨物量。目標(biāo)函數(shù)總成本最小化。即對(duì)所有可能的i, j求單位成本_ij * X_ij的總和。約束條件供給約束對(duì)于每個(gè)工廠(chǎng)i運(yùn)往所有銷(xiāo)售點(diǎn)的貨物總量不能超過(guò)其產(chǎn)量。需求約束對(duì)于每個(gè)銷(xiāo)售點(diǎn)j從所有工廠(chǎng)運(yùn)來(lái)的貨物總量必須恰好滿(mǎn)足其需求量平衡型運(yùn)輸問(wèn)題。非負(fù)約束運(yùn)輸量不能為負(fù)數(shù)。這個(gè)模型的妙處在于它雖然變量多工廠(chǎng)數(shù)×銷(xiāo)售點(diǎn)數(shù)但約束和目標(biāo)函數(shù)都是線(xiàn)性的結(jié)構(gòu)清晰是LINGO最擅長(zhǎng)處理的類(lèi)型之一。在LINGO中我們可以利用集合SETS功能來(lái)優(yōu)雅地定義這種多維變量和約束避免為每一個(gè)變量和約束都寫(xiě)一行代碼。2.2 連續(xù)投資問(wèn)題動(dòng)態(tài)決策的跨期優(yōu)化連續(xù)投資問(wèn)題也叫多階段投資問(wèn)題比靜態(tài)的運(yùn)輸問(wèn)題多了一個(gè)時(shí)間維度。典型描述是你有一筆初始資金需要在未來(lái)若干年比如5年內(nèi)進(jìn)行投資。市場(chǎng)上有多種投資項(xiàng)目每種項(xiàng)目每年有特定的投資期限和收益率。每年你會(huì)有新的資金流入如項(xiàng)目到期回收的本息然后需要重新分配投資。目標(biāo)是在規(guī)劃期結(jié)束時(shí)讓總資金額最大。這本質(zhì)上是一個(gè)帶有時(shí)序關(guān)系的線(xiàn)性或非線(xiàn)性規(guī)劃問(wèn)題。其核心在于“動(dòng)態(tài)”和“決策鏈”決策變量通常設(shè)為X_it表示在第t年投入到項(xiàng)目i的資金量。目標(biāo)函數(shù)最終時(shí)刻如第5年末的總資金最大化。約束條件這是關(guān)鍵和難點(diǎn)資金流平衡約束最核心每年可用于投資的資金總額等于上年末回收的資金加上本年可能的外來(lái)資金。這建立起了年份之間的連接。例如第2年的投資總額等于第1年那些為期1年的項(xiàng)目回收的本息和。項(xiàng)目投資限額約束每個(gè)項(xiàng)目每年可能有投資上限。非負(fù)約束投資額不能為負(fù)數(shù)。邏輯約束如果項(xiàng)目期限大于1年那么一旦投資在到期前資金是被鎖定的中間年份不能再對(duì)其追加或減少投資。這通常需要通過(guò)定義變量或附加約束來(lái)實(shí)現(xiàn)。注意連續(xù)投資問(wèn)題的模型建立比運(yùn)輸問(wèn)題更考驗(yàn)邏輯嚴(yán)密性。一個(gè)常見(jiàn)的錯(cuò)誤是忽略了資金的時(shí)間鎖定效應(yīng)或者寫(xiě)錯(cuò)了資金平衡方程的順序?qū)е履Ｐ颓蠼饨Y(jié)果違背常識(shí)比如第一年就把所有錢(qián)投到一個(gè)長(zhǎng)期項(xiàng)目后面幾年無(wú)錢(qián)可投。在LINGO中實(shí)現(xiàn)時(shí)需要仔細(xì)定義集合尤其是時(shí)間集并利用FOR循環(huán)來(lái)精確表達(dá)這些跨期約束。3. LINGO實(shí)戰(zhàn)運(yùn)輸問(wèn)題建模與求解我們用一個(gè)具體例子來(lái)上手。假設(shè)有3個(gè)工廠(chǎng)A1, A2, A3產(chǎn)量分別是7噸、4噸、9噸4個(gè)銷(xiāo)售點(diǎn)B1, B2, B3, B4需求量分別是3噸、6噸、5噸、6噸。單位運(yùn)價(jià)表元/噸如下從\到B1B2B3B4A1311310A21928A3741053.1 LINGO模型代碼詳解打開(kāi)LINGO在模型窗口輸入以下代碼。我會(huì)逐段解釋MODEL: ! 1. 定義集合; SETS: factory /A1, A2, A3/: capacity; ! 工廠(chǎng)集合屬性為產(chǎn)能capacity; market /B1, B2, B3, B4/: demand; ! 市場(chǎng)集合屬性為需求demand; link(factory, market): cost, volume; ! 派生集合表示工廠(chǎng)到市場(chǎng)的路線(xiàn)屬性為單位成本cost和運(yùn)量volume; ENDSETS ! 2. 輸入數(shù)據(jù); DATA: capacity 7 4 9; ! 工廠(chǎng)產(chǎn)能; demand 3 6 5 6; ! 市場(chǎng)需求; cost 3 11 3 10 ! 成本矩陣按行排列: A1-B1,B2,B3,B4; 1 9 2 8 ! A2-B1,B2,B3,B4; 7 4 10 5; ! A3-B1,B2,B3,B4; ENDDATA ! 3. 定義目標(biāo)函數(shù)最小化總運(yùn)輸成本; MIN SUM(link(i, j): cost(i, j) * volume(i, j)); ! 4. 定義約束條件; ! 供給約束每個(gè)工廠(chǎng)運(yùn)出的總量不超過(guò)其產(chǎn)能; FOR(factory(i): SUM(market(j): volume(i, j)) capacity(i) ); ! 需求約束每個(gè)市場(chǎng)運(yùn)入的總量等于其需求; FOR(market(j): SUM(factory(i): volume(i, j)) demand(j) ); ! 5. 變量范圍約束非負(fù)已在LINGO中默認(rèn)如需特殊限制可在此添加; END關(guān)鍵點(diǎn)解析集合定義這是LINGO建模的精髓。link(factory, market)自動(dòng)生成了所有工廠(chǎng)和市場(chǎng)的組合3x412條路線(xiàn)我們無(wú)需手動(dòng)定義12個(gè)變量。cost和volume是附著在這些路線(xiàn)上的屬性。數(shù)據(jù)輸入DATA部分用于賦值。注意cost矩陣的輸入是按行排列的順序必須與link集合中元素的生成順序一致通常是先遍歷第一個(gè)集合的所有元素再遍歷第二個(gè)。這是新手最容易出錯(cuò)的地方之一。目標(biāo)與約束SUM和FOR是LINGO的循環(huán)求和與循環(huán)定義函數(shù)。FOR(factory(i): ...)意味著為factory集合中的每一個(gè)成員i即A1, A2, A3都生成一個(gè)后面的約束不等式。這樣三行供給約束就通過(guò)一行代碼生成了。3.2 求解與結(jié)果分析輸入代碼后點(diǎn)擊工具欄的“求解”按鈕靶心圖標(biāo)。LINGO會(huì)調(diào)用內(nèi)置求解器進(jìn)行計(jì)算。求解狀態(tài)報(bào)告求解完成后會(huì)彈出對(duì)話(huà)框。你應(yīng)該看到“Global optimal solution found.”找到全局最優(yōu)解這表明求解成功。同時(shí)會(huì)顯示目標(biāo)函數(shù)值最小總成本和迭代次數(shù)。查看結(jié)果關(guān)閉報(bào)告窗口回到模型界面。選擇菜單LINGO - Solution或直接查看狀態(tài)窗口下方的“Solution Report”標(biāo)簽頁(yè)。這里會(huì)列出所有變量的最優(yōu)值。 對(duì)于本例最優(yōu)解volume變量大致如下A1-B3: 5噸 A1-B4: 2噸A2-B1: 3噸 A2-B3: 0噸 A2-B4: 1噸A3-B2: 6噸 A3-B4: 3噸 其他路線(xiàn)運(yùn)量為0 最小總成本為86元。靈敏度分析可選但重要對(duì)于線(xiàn)性規(guī)劃LINGO還可以提供靈敏度分析報(bào)告LINGO - Range。這份報(bào)告會(huì)告訴你每個(gè)約束條件的“影子價(jià)格”和每個(gè)變量目標(biāo)函數(shù)系數(shù)的“允許變化范圍”。例如某個(gè)工廠(chǎng)產(chǎn)能約束的影子價(jià)格是5意味著該工廠(chǎng)產(chǎn)能每增加1單位總成本最多能減少5元。這在商業(yè)決策中極具價(jià)值。實(shí)操心得第一次運(yùn)行時(shí)如果遇到“No feasible solution found”無(wú)可行解大概率是約束條件互相矛盾。比如總產(chǎn)能74920小于總需求365620本例是平衡的。如果不平衡產(chǎn)能需求供給約束應(yīng)改為如果產(chǎn)能需求則問(wèn)題無(wú)解或需求約束需改為并修改目標(biāo)。務(wù)必先手算驗(yàn)證一下總量關(guān)系。4. LINGO進(jìn)階連續(xù)投資問(wèn)題建模我們看一個(gè)簡(jiǎn)化版的5年投資問(wèn)題。假設(shè)有3個(gè)項(xiàng)目可供投資項(xiàng)目A每年初投資年末回收本利年利率10%。項(xiàng)目B每年初投資需要兩年跨兩個(gè)年度到期后回收本利年利率40%。項(xiàng)目C每年初投資需要三年跨三個(gè)年度到期后回收本利年利率80%。 期初有10萬(wàn)元資金每年末的項(xiàng)目回收款可以用于下一年初再投資。如何分配每年的投資使第5年末的本利總額最大建模思路我們需要定義5個(gè)時(shí)間點(diǎn)年初。決策變量Xij表示在第i年年初投入到項(xiàng)目j的資金量j1,2,3對(duì)應(yīng)A,B,C。關(guān)鍵在于理清每年年初的可投資資金流。4.1 LINGO模型代碼MODEL: SETS: year /1..5/:; ! 定義5個(gè)年份; project /A, B, C/:; ! 定義3個(gè)項(xiàng)目; invest(year, project): X; ! 決策變量X(i,j)表示第i年初投給項(xiàng)目j的錢(qián); ENDSETS DATA: ! 利率參數(shù); rate_A 0.10; rate_B 0.40; rate_C 0.80; ! 項(xiàng)目期限年; term_A 1; term_B 2; term_C 3; ENDDATA ! 初始資金; INIT: capital 100000; ENDINIT ! 目標(biāo)函數(shù)最大化第五年末的總資金; MAX F; ! 定義每年的資金流平衡; ! 第1年初可投資資金就是初始資金; SUM(project(j): X(1, j)) capital; ! 第2年初可投資資金 第1年項(xiàng)目A到期回收的本利; SUM(project(j): X(2, j)) (1rate_A) * X(1, A); ! 第3年初可投資資金 第2年項(xiàng)目A到期回收的本利 第1年項(xiàng)目B到期回收的本利; SUM(project(j): X(3, j)) (1rate_A) * X(2, A) (1rate_B) * X(1, B); ! 第4年初可投資資金 第3年項(xiàng)目A到期回收的本利 第2年項(xiàng)目B到期回收的本利 第1年項(xiàng)目C到期回收的本利; SUM(project(j): X(4, j)) (1rate_A) * X(3, A) (1rate_B) * X(2, B) (1rate_C) * X(1, C); ! 第5年初可投資資金 第4年項(xiàng)目A到期回收的本利 第3年項(xiàng)目B到期回收的本利 第2年項(xiàng)目C到期回收的本利; SUM(project(j): X(5, j)) (1rate_A) * X(4, A) (1rate_B) * X(3, B) (1rate_C) * X(2, C); ! 第5年末的總資金F等于第5年所有投資到期回收的本利和; F (1rate_A) * X(5, A) (1rate_B) * X(4, B) (1rate_C) * X(3, C); ! 非負(fù)約束; FOR(invest(i, j): X(i, j) 0); END4.2 模型難點(diǎn)與技巧這個(gè)模型比運(yùn)輸問(wèn)題復(fù)雜體現(xiàn)在資金流方程的構(gòu)建每年初的投資總額必須等于上年末回收的資金總額。回收的資金來(lái)自之前年份投資的、且剛好在本年年初到期的項(xiàng)目。必須仔細(xì)追蹤每個(gè)項(xiàng)目的投資年份和到期年份。我上面用注釋詳細(xì)列出了每筆錢(qián)的來(lái)源建議在紙上畫(huà)一個(gè)時(shí)間線(xiàn)圖來(lái)輔助理解。目標(biāo)函數(shù)的定義最終財(cái)富F不是簡(jiǎn)單地加總所有投資而是第五年末時(shí)所有已到期投資的本利和。即第五年初的投資中只有項(xiàng)目A能在當(dāng)年末收回第四年初投的B項(xiàng)目、第三年初投的C項(xiàng)目也會(huì)在第五年末到期。使用INIT部分對(duì)于初始資金這種在模型開(kāi)始前就確定的常量可以用INIT部分來(lái)初始化使模型更清晰。求解結(jié)果分析求解此模型你會(huì)得到每年在各個(gè)項(xiàng)目上的最優(yōu)投資分配X(i,j)以及最終的最大本利和F。你會(huì)發(fā)現(xiàn)盡管項(xiàng)目C的利率最高但因?yàn)樗i定期長(zhǎng)3年所以并不會(huì)把所有錢(qián)都投給它。最優(yōu)策略是在前期利用C和B的高收益進(jìn)行積累后期則更注重資金的流動(dòng)性多投A以確保最后一年有足夠的到期資金來(lái)貢獻(xiàn)給最終收益F。這個(gè)動(dòng)態(tài)權(quán)衡的過(guò)程正是連續(xù)投資問(wèn)題的精髓。踩坑提醒在寫(xiě)資金平衡方程時(shí)極易出現(xiàn)“時(shí)間錯(cuò)配”。例如把X(1,B)的回收算到第2年初實(shí)際上它需要2年應(yīng)在第3年初回收。一個(gè)有效的檢查方法是對(duì)于任意X(i, j)其回收年份一定是i term_j。在LINGO中你可以嘗試用更復(fù)雜的集合和SUM函數(shù)來(lái)生成這些約束以增強(qiáng)模型的通用性比如年份數(shù)變化時(shí)無(wú)需重寫(xiě)方程但這需要更熟練的集合操作。初學(xué)時(shí)像上面這樣“笨拙”地明確寫(xiě)出每個(gè)方程反而更不容易出錯(cuò)也便于調(diào)試。5. 常見(jiàn)問(wèn)題、調(diào)試技巧與模型優(yōu)化即使模型邏輯正確在LINGO中也可能遇到各種問(wèn)題。以下是一些常見(jiàn)情況及處理辦法。5.1 常見(jiàn)錯(cuò)誤與排查錯(cuò)誤類(lèi)型可能原因排查與解決方法No feasible solution found1. 約束條件相互矛盾。2. 變量上下限設(shè)置過(guò)緊無(wú)解空間。3. 初始值設(shè)置不當(dāng)導(dǎo)致求解器誤判。1.檢查總量平衡如運(yùn)輸問(wèn)題中總供給與總需求關(guān)系。2.簡(jiǎn)化模型先注釋掉部分約束看是否能得到可行解逐步添加以定位矛盾點(diǎn)。3. 檢查BND或變量賦值是否合理。Unbounded solution目標(biāo)函數(shù)可以在不違反約束的情況下無(wú)限優(yōu)化如利潤(rùn)無(wú)限大。這通常是模型構(gòu)建錯(cuò)誤。檢查是否遺漏了關(guān)鍵的資源限制約束。例如在投資問(wèn)題中是否忘記了每年投資總額不能超過(guò)可用資金這一根本約束。求解時(shí)間過(guò)長(zhǎng)或內(nèi)存不足1. 問(wèn)題規(guī)模太大變量/約束過(guò)多。2. 模型是非線(xiàn)性或整數(shù)規(guī)劃本身復(fù)雜。1. 嘗試使用LINGO - Options中的求解器選項(xiàng)如調(diào)整迭代次數(shù)、收斂精度。2. 對(duì)于整數(shù)規(guī)劃可以設(shè)置啟發(fā)式搜索或調(diào)整分支策略。3. 考慮是否能用線(xiàn)性規(guī)劃近似或者分解問(wèn)題。結(jié)果與預(yù)期不符1. 目標(biāo)函數(shù)系數(shù)如成本、利率輸入錯(cuò)誤。2. 約束條件符號(hào),,用錯(cuò)。3. 單位不統(tǒng)一。1.仔細(xì)核對(duì)DATA部分特別是矩陣數(shù)據(jù)的排列順序。2.逐行檢查約束理解其經(jīng)濟(jì)/物理意義。3. 用LINGO - Picture功能可視化模型結(jié)構(gòu)檢查行、列對(duì)應(yīng)關(guān)系。5.2 LINGO編程與調(diào)試技巧多用注釋使用!添加詳細(xì)注釋說(shuō)明每個(gè)集合、變量、約束的含義。這對(duì)于復(fù)雜的、幾天后再回顧的模型至關(guān)重要。模塊化建模對(duì)于大型模型可以分塊編寫(xiě)。先定義集合和數(shù)據(jù)再寫(xiě)目標(biāo)函數(shù)然后分模塊寫(xiě)約束如供給約束塊、需求約束塊、資金平衡塊。每寫(xiě)完一個(gè)模塊可以暫時(shí)注釋掉其他部分單獨(dú)測(cè)試其語(yǔ)法是否正確。利用FOR和SUM簡(jiǎn)化代碼這是LINGO的核心優(yōu)勢(shì)。對(duì)于規(guī)律性強(qiáng)的約束一定要用循環(huán)語(yǔ)句避免手動(dòng)羅列。例如如果我們的投資問(wèn)題有10年、10個(gè)項(xiàng)目手動(dòng)寫(xiě)資金平衡方程將是災(zāi)難。此時(shí)應(yīng)該定義“項(xiàng)目期限”屬性然后利用SUM函數(shù)和條件判斷來(lái)動(dòng)態(tài)生成回收項(xiàng)。! 假設(shè)years是年份集合projects是項(xiàng)目集合term(p)是項(xiàng)目p的期限; FOR(years(t) | t #GT# 1: ! 從第2年開(kāi)始; SUM(projects(p): X(t, p)) SUM(projects(p): (1rate(p)) * X(t - term(p), p)) ! 注意這里需要確保 t-term(p) 1否則會(huì)引用不存在的變量需用IF函數(shù)處理; );查看模型概況求解前點(diǎn)擊LINGO - Generate - Display model可以查看LINGO將你的集合語(yǔ)言展開(kāi)后的完整線(xiàn)性規(guī)劃模型。這是驗(yàn)證模型是否按你設(shè)想的方式構(gòu)建的終極手段。從簡(jiǎn)單案例開(kāi)始不要一開(kāi)始就構(gòu)建復(fù)雜的模型。先用2個(gè)工廠(chǎng)、2個(gè)市場(chǎng)或者2年、2個(gè)項(xiàng)目的小數(shù)據(jù)測(cè)試模型。用手算能驗(yàn)證結(jié)果確認(rèn)模型邏輯無(wú)誤后再替換成真實(shí)的大規(guī)模數(shù)據(jù)。5.3 模型擴(kuò)展與優(yōu)化思路掌握了基礎(chǔ)模型后你可以應(yīng)對(duì)更復(fù)雜的情況運(yùn)輸問(wèn)題的變種產(chǎn)銷(xiāo)不平衡修改供給或需求約束的等號(hào)為不等號(hào)。帶有中轉(zhuǎn)站引入新的節(jié)點(diǎn)集合中轉(zhuǎn)站變量變?yōu)閄_ik(工廠(chǎng)到中轉(zhuǎn)站) 和Y_kj(中轉(zhuǎn)站到市場(chǎng))并增加中轉(zhuǎn)站的流量平衡約束。多商品運(yùn)輸為每種商品增加一個(gè)維度使用三維數(shù)組X_ijk(工廠(chǎng)i, 市場(chǎng)j, 商品k)。投資問(wèn)題的變種風(fēng)險(xiǎn)約束引入風(fēng)險(xiǎn)系數(shù)增加組合風(fēng)險(xiǎn)不超過(guò)某閾值的約束。借貸可能引入“借款”變量利率不同在資金平衡方程中增加借款流入和還款流出項(xiàng)。非線(xiàn)性收益如果收益率與投資額有關(guān)如規(guī)模效應(yīng)則目標(biāo)函數(shù)變?yōu)榉蔷€(xiàn)性需要在LINGO中選用非線(xiàn)性求解器LINGO - Options - Solver中選擇Global或Multistart。最后我個(gè)人在數(shù)學(xué)建模競(jìng)賽培訓(xùn)和實(shí)際項(xiàng)目中最深的體會(huì)是LINGO這類(lèi)工具的強(qiáng)大永遠(yuǎn)建立在你對(duì)問(wèn)題本質(zhì)的深刻理解之上。軟件只是幫你做計(jì)算而建模才是真正的思考?；ㄔ谇逦x變量、梳理約束關(guān)系上的時(shí)間遠(yuǎn)比調(diào)試錯(cuò)誤代碼的時(shí)間更有價(jià)值。當(dāng)你拿到一個(gè)復(fù)雜問(wèn)題不妨先問(wèn)自己決策是什么目標(biāo)是什么限制條件有哪些它們之間如何用數(shù)學(xué)語(yǔ)言關(guān)聯(lián)把這些想明白了翻譯成LINGO代碼就是水到渠成的事。開(kāi)始的時(shí)候可能會(huì)覺(jué)得語(yǔ)法生疏但就像學(xué)任何語(yǔ)言一樣多寫(xiě)幾個(gè)完整的模型你就能越來(lái)越熟練地用它來(lái)表達(dá)你的優(yōu)化思想了。 數(shù)學(xué)建模實(shí)戰(zhàn)：用LINGO求解運(yùn)輸與投資優(yōu)化問(wèn)題
1. 項(xiàng)目概述當(dāng)數(shù)學(xué)建模遇上LINGO最近在帶學(xué)生做數(shù)學(xué)建模競(jìng)賽的培訓(xùn)發(fā)現(xiàn)很多同學(xué)一看到“運(yùn)輸方案優(yōu)化”或者“連續(xù)投資決策”這類(lèi)題目就有點(diǎn)發(fā)怵。題目描述往往是一大堆數(shù)據(jù)表格要求你給出成本最低的運(yùn)輸計(jì)劃或者未來(lái)幾年收益最大的投資策略。思路大家其實(shí)都有列方程、設(shè)變量嘛但真到了求解那一步尤其是變量一多、約束一復(fù)雜手工計(jì)算幾乎不可能用通用編程語(yǔ)言比如Python從頭寫(xiě)算法又容易在調(diào)試上耗費(fèi)大量時(shí)間。這時(shí)候一個(gè)專(zhuān)業(yè)的工具就能讓你事半功倍——它就是LINGO。簡(jiǎn)單來(lái)說(shuō)LINGO是一款專(zhuān)門(mén)用于求解線(xiàn)性規(guī)劃、非線(xiàn)性規(guī)劃、整數(shù)規(guī)劃等優(yōu)化問(wèn)題的軟件。你可以把它理解為一個(gè)“數(shù)學(xué)規(guī)劃模型的翻譯器和求解器”。你不需要去深究背后的單純形法、分支定界法是什么原理你只需要用接近數(shù)學(xué)語(yǔ)言的格式把目標(biāo)函數(shù)、決策變量和約束條件“告訴”它它就能快速給你算出最優(yōu)解。對(duì)于“運(yùn)輸方案”和“連續(xù)投資”這類(lèi)典型的運(yùn)籌學(xué)問(wèn)題LINGO幾乎是量身定做的解決方案。本文我就結(jié)合這兩個(gè)經(jīng)典場(chǎng)景帶你從零開(kāi)始掌握用LINGO將數(shù)學(xué)模型轉(zhuǎn)化為實(shí)際解決方案的全過(guò)程并分享一些我多年輔導(dǎo)中總結(jié)的、在官方教程里未必會(huì)寫(xiě)的實(shí)戰(zhàn)技巧和避坑指南。2. 核心思路拆解從問(wèn)題到LINGO模型在打開(kāi)LINGO軟件之前最關(guān)鍵的一步是把現(xiàn)實(shí)問(wèn)題抽象成一個(gè)清晰的數(shù)學(xué)模型。很多新手會(huì)直接跳進(jìn)軟件里敲代碼結(jié)果發(fā)現(xiàn)模型本身有邏輯漏洞求解自然失敗。我們分別看看這兩個(gè)問(wèn)題的建模核心。2.1 運(yùn)輸問(wèn)題成本最小化的網(wǎng)絡(luò)流運(yùn)輸問(wèn)題的背景非常直觀你有若干個(gè)工廠(chǎng)產(chǎn)地每個(gè)工廠(chǎng)有固定的產(chǎn)量供給量你也有若干個(gè)銷(xiāo)售點(diǎn)銷(xiāo)地每個(gè)銷(xiāo)售點(diǎn)有固定的需求量。從每個(gè)工廠(chǎng)到每個(gè)銷(xiāo)售點(diǎn)都有單位運(yùn)輸成本。目標(biāo)是如何安排運(yùn)輸量在滿(mǎn)足所有供給和需求的前提下讓總運(yùn)輸成本最低。這本質(zhì)上是一個(gè)線(xiàn)性規(guī)劃問(wèn)題。它的核心結(jié)構(gòu)非常規(guī)整決策變量通常設(shè)為X_ij表示從工廠(chǎng)i運(yùn)到銷(xiāo)售點(diǎn)j的貨物量。目標(biāo)函數(shù)總成本最小化。即對(duì)所有可能的i, j求單位成本_ij * X_ij的總和。約束條件供給約束對(duì)于每個(gè)工廠(chǎng)i運(yùn)往所有銷(xiāo)售點(diǎn)的貨物總量不能超過(guò)其產(chǎn)量。需求約束對(duì)于每個(gè)銷(xiāo)售點(diǎn)j從所有工廠(chǎng)運(yùn)來(lái)的貨物總量必須恰好滿(mǎn)足其需求量平衡型運(yùn)輸問(wèn)題。非負(fù)約束運(yùn)輸量不能為負(fù)數(shù)。這個(gè)模型的妙處在于它雖然變量多工廠(chǎng)數(shù)×銷(xiāo)售點(diǎn)數(shù)但約束和目標(biāo)函數(shù)都是線(xiàn)性的結(jié)構(gòu)清晰是LINGO最擅長(zhǎng)處理的類(lèi)型之一。在LINGO中我們可以利用集合SETS功能來(lái)優(yōu)雅地定義這種多維變量和約束避免為每一個(gè)變量和約束都寫(xiě)一行代碼。2.2 連續(xù)投資問(wèn)題動(dòng)態(tài)決策的跨期優(yōu)化連續(xù)投資問(wèn)題也叫多階段投資問(wèn)題比靜態(tài)的運(yùn)輸問(wèn)題多了一個(gè)時(shí)間維度。典型描述是你有一筆初始資金需要在未來(lái)若干年比如5年內(nèi)進(jìn)行投資。市場(chǎng)上有多種投資項(xiàng)目每種項(xiàng)目每年有特定的投資期限和收益率。每年你會(huì)有新的資金流入如項(xiàng)目到期回收的本息然后需要重新分配投資。目標(biāo)是在規(guī)劃期結(jié)束時(shí)讓總資金額最大。這本質(zhì)上是一個(gè)帶有時(shí)序關(guān)系的線(xiàn)性或非線(xiàn)性規(guī)劃問(wèn)題。其核心在于“動(dòng)態(tài)”和“決策鏈”決策變量通常設(shè)為X_it表示在第t年投入到項(xiàng)目i的資金量。目標(biāo)函數(shù)最終時(shí)刻如第5年末的總資金最大化。約束條件這是關(guān)鍵和難點(diǎn)資金流平衡約束最核心每年可用于投資的資金總額等于上年末回收的資金加上本年可能的外來(lái)資金。這建立起了年份之間的連接。例如第2年的投資總額等于第1年那些為期1年的項(xiàng)目回收的本息和。項(xiàng)目投資限額約束每個(gè)項(xiàng)目每年可能有投資上限。非負(fù)約束投資額不能為負(fù)數(shù)。邏輯約束如果項(xiàng)目期限大于1年那么一旦投資在到期前資金是被鎖定的中間年份不能再對(duì)其追加或減少投資。這通常需要通過(guò)定義變量或附加約束來(lái)實(shí)現(xiàn)。注意連續(xù)投資問(wèn)題的模型建立比運(yùn)輸問(wèn)題更考驗(yàn)邏輯嚴(yán)密性。一個(gè)常見(jiàn)的錯(cuò)誤是忽略了資金的時(shí)間鎖定效應(yīng)或者寫(xiě)錯(cuò)了資金平衡方程的順序?qū)е履Ｐ颓蠼饨Y(jié)果違背常識(shí)比如第一年就把所有錢(qián)投到一個(gè)長(zhǎng)期項(xiàng)目后面幾年無(wú)錢(qián)可投。在LINGO中實(shí)現(xiàn)時(shí)需要仔細(xì)定義集合尤其是時(shí)間集并利用FOR循環(huán)來(lái)精確表達(dá)這些跨期約束。3. LINGO實(shí)戰(zhàn)運(yùn)輸問(wèn)題建模與求解我們用一個(gè)具體例子來(lái)上手。假設(shè)有3個(gè)工廠(chǎng)A1, A2, A3產(chǎn)量分別是7噸、4噸、9噸4個(gè)銷(xiāo)售點(diǎn)B1, B2, B3, B4需求量分別是3噸、6噸、5噸、6噸。單位運(yùn)價(jià)表元/噸如下從\到B1B2B3B4A1311310A21928A3741053.1 LINGO模型代碼詳解打開(kāi)LINGO在模型窗口輸入以下代碼。我會(huì)逐段解釋MODEL: ! 1. 定義集合; SETS: factory /A1, A2, A3/: capacity; ! 工廠(chǎng)集合屬性為產(chǎn)能capacity; market /B1, B2, B3, B4/: demand; ! 市場(chǎng)集合屬性為需求demand; link(factory, market): cost, volume; ! 派生集合表示工廠(chǎng)到市場(chǎng)的路線(xiàn)屬性為單位成本cost和運(yùn)量volume; ENDSETS ! 2. 輸入數(shù)據(jù); DATA: capacity 7 4 9; ! 工廠(chǎng)產(chǎn)能; demand 3 6 5 6; ! 市場(chǎng)需求; cost 3 11 3 10 ! 成本矩陣按行排列: A1-B1,B2,B3,B4; 1 9 2 8 ! A2-B1,B2,B3,B4; 7 4 10 5; ! A3-B1,B2,B3,B4; ENDDATA ! 3. 定義目標(biāo)函數(shù)最小化總運(yùn)輸成本; MIN SUM(link(i, j): cost(i, j) * volume(i, j)); ! 4. 定義約束條件; ! 供給約束每個(gè)工廠(chǎng)運(yùn)出的總量不超過(guò)其產(chǎn)能; FOR(factory(i): SUM(market(j): volume(i, j)) capacity(i) ); ! 需求約束每個(gè)市場(chǎng)運(yùn)入的總量等于其需求; FOR(market(j): SUM(factory(i): volume(i, j)) demand(j) ); ! 5. 變量范圍約束非負(fù)已在LINGO中默認(rèn)如需特殊限制可在此添加; END關(guān)鍵點(diǎn)解析集合定義這是LINGO建模的精髓。link(factory, market)自動(dòng)生成了所有工廠(chǎng)和市場(chǎng)的組合3x412條路線(xiàn)我們無(wú)需手動(dòng)定義12個(gè)變量。cost和volume是附著在這些路線(xiàn)上的屬性。數(shù)據(jù)輸入DATA部分用于賦值。注意cost矩陣的輸入是按行排列的順序必須與link集合中元素的生成順序一致通常是先遍歷第一個(gè)集合的所有元素再遍歷第二個(gè)。這是新手最容易出錯(cuò)的地方之一。目標(biāo)與約束SUM和FOR是LINGO的循環(huán)求和與循環(huán)定義函數(shù)。FOR(factory(i): ...)意味著為factory集合中的每一個(gè)成員i即A1, A2, A3都生成一個(gè)后面的約束不等式。這樣三行供給約束就通過(guò)一行代碼生成了。3.2 求解與結(jié)果分析輸入代碼后點(diǎn)擊工具欄的“求解”按鈕靶心圖標(biāo)。LINGO會(huì)調(diào)用內(nèi)置求解器進(jìn)行計(jì)算。求解狀態(tài)報(bào)告求解完成后會(huì)彈出對(duì)話(huà)框。你應(yīng)該看到“Global optimal solution found.”找到全局最優(yōu)解這表明求解成功。同時(shí)會(huì)顯示目標(biāo)函數(shù)值最小總成本和迭代次數(shù)。查看結(jié)果關(guān)閉報(bào)告窗口回到模型界面。選擇菜單LINGO - Solution或直接查看狀態(tài)窗口下方的“Solution Report”標(biāo)簽頁(yè)。這里會(huì)列出所有變量的最優(yōu)值。 對(duì)于本例最優(yōu)解volume變量大致如下A1-B3: 5噸 A1-B4: 2噸A2-B1: 3噸 A2-B3: 0噸 A2-B4: 1噸A3-B2: 6噸 A3-B4: 3噸 其他路線(xiàn)運(yùn)量為0 最小總成本為86元。靈敏度分析可選但重要對(duì)于線(xiàn)性規(guī)劃LINGO還可以提供靈敏度分析報(bào)告LINGO - Range。這份報(bào)告會(huì)告訴你每個(gè)約束條件的“影子價(jià)格”和每個(gè)變量目標(biāo)函數(shù)系數(shù)的“允許變化范圍”。例如某個(gè)工廠(chǎng)產(chǎn)能約束的影子價(jià)格是5意味著該工廠(chǎng)產(chǎn)能每增加1單位總成本最多能減少5元。這在商業(yè)決策中極具價(jià)值。實(shí)操心得第一次運(yùn)行時(shí)如果遇到“No feasible solution found”無(wú)可行解大概率是約束條件互相矛盾。比如總產(chǎn)能74920小于總需求365620本例是平衡的。如果不平衡產(chǎn)能需求供給約束應(yīng)改為如果產(chǎn)能需求則問(wèn)題無(wú)解或需求約束需改為并修改目標(biāo)。務(wù)必先手算驗(yàn)證一下總量關(guān)系。4. LINGO進(jìn)階連續(xù)投資問(wèn)題建模我們看一個(gè)簡(jiǎn)化版的5年投資問(wèn)題。假設(shè)有3個(gè)項(xiàng)目可供投資項(xiàng)目A每年初投資年末回收本利年利率10%。項(xiàng)目B每年初投資需要兩年跨兩個(gè)年度到期后回收本利年利率40%。項(xiàng)目C每年初投資需要三年跨三個(gè)年度到期后回收本利年利率80%。 期初有10萬(wàn)元資金每年末的項(xiàng)目回收款可以用于下一年初再投資。如何分配每年的投資使第5年末的本利總額最大建模思路我們需要定義5個(gè)時(shí)間點(diǎn)年初。決策變量Xij表示在第i年年初投入到項(xiàng)目j的資金量j1,2,3對(duì)應(yīng)A,B,C。關(guān)鍵在于理清每年年初的可投資資金流。4.1 LINGO模型代碼MODEL: SETS: year /1..5/:; ! 定義5個(gè)年份; project /A, B, C/:; ! 定義3個(gè)項(xiàng)目; invest(year, project): X; ! 決策變量X(i,j)表示第i年初投給項(xiàng)目j的錢(qián); ENDSETS DATA: ! 利率參數(shù); rate_A 0.10; rate_B 0.40; rate_C 0.80; ! 項(xiàng)目期限年; term_A 1; term_B 2; term_C 3; ENDDATA ! 初始資金; INIT: capital 100000; ENDINIT ! 目標(biāo)函數(shù)最大化第五年末的總資金; MAX F; ! 定義每年的資金流平衡; ! 第1年初可投資資金就是初始資金; SUM(project(j): X(1, j)) capital; ! 第2年初可投資資金 第1年項(xiàng)目A到期回收的本利; SUM(project(j): X(2, j)) (1rate_A) * X(1, A); ! 第3年初可投資資金 第2年項(xiàng)目A到期回收的本利 第1年項(xiàng)目B到期回收的本利; SUM(project(j): X(3, j)) (1rate_A) * X(2, A) (1rate_B) * X(1, B); ! 第4年初可投資資金 第3年項(xiàng)目A到期回收的本利 第2年項(xiàng)目B到期回收的本利 第1年項(xiàng)目C到期回收的本利; SUM(project(j): X(4, j)) (1rate_A) * X(3, A) (1rate_B) * X(2, B) (1rate_C) * X(1, C); ! 第5年初可投資資金 第4年項(xiàng)目A到期回收的本利 第3年項(xiàng)目B到期回收的本利 第2年項(xiàng)目C到期回收的本利; SUM(project(j): X(5, j)) (1rate_A) * X(4, A) (1rate_B) * X(3, B) (1rate_C) * X(2, C); ! 第5年末的總資金F等于第5年所有投資到期回收的本利和; F (1rate_A) * X(5, A) (1rate_B) * X(4, B) (1rate_C) * X(3, C); ! 非負(fù)約束; FOR(invest(i, j): X(i, j) 0); END4.2 模型難點(diǎn)與技巧這個(gè)模型比運(yùn)輸問(wèn)題復(fù)雜體現(xiàn)在資金流方程的構(gòu)建每年初的投資總額必須等于上年末回收的資金總額。回收的資金來(lái)自之前年份投資的、且剛好在本年年初到期的項(xiàng)目。必須仔細(xì)追蹤每個(gè)項(xiàng)目的投資年份和到期年份。我上面用注釋詳細(xì)列出了每筆錢(qián)的來(lái)源建議在紙上畫(huà)一個(gè)時(shí)間線(xiàn)圖來(lái)輔助理解。目標(biāo)函數(shù)的定義最終財(cái)富F不是簡(jiǎn)單地加總所有投資而是第五年末時(shí)所有已到期投資的本利和。即第五年初的投資中只有項(xiàng)目A能在當(dāng)年末收回第四年初投的B項(xiàng)目、第三年初投的C項(xiàng)目也會(huì)在第五年末到期。使用INIT部分對(duì)于初始資金這種在模型開(kāi)始前就確定的常量可以用INIT部分來(lái)初始化使模型更清晰。求解結(jié)果分析求解此模型你會(huì)得到每年在各個(gè)項(xiàng)目上的最優(yōu)投資分配X(i,j)以及最終的最大本利和F。你會(huì)發(fā)現(xiàn)盡管項(xiàng)目C的利率最高但因?yàn)樗i定期長(zhǎng)3年所以并不會(huì)把所有錢(qián)都投給它。最優(yōu)策略是在前期利用C和B的高收益進(jìn)行積累后期則更注重資金的流動(dòng)性多投A以確保最后一年有足夠的到期資金來(lái)貢獻(xiàn)給最終收益F。這個(gè)動(dòng)態(tài)權(quán)衡的過(guò)程正是連續(xù)投資問(wèn)題的精髓。踩坑提醒在寫(xiě)資金平衡方程時(shí)極易出現(xiàn)“時(shí)間錯(cuò)配”。例如把X(1,B)的回收算到第2年初實(shí)際上它需要2年應(yīng)在第3年初回收。一個(gè)有效的檢查方法是對(duì)于任意X(i, j)其回收年份一定是i term_j。在LINGO中你可以嘗試用更復(fù)雜的集合和SUM函數(shù)來(lái)生成這些約束以增強(qiáng)模型的通用性比如年份數(shù)變化時(shí)無(wú)需重寫(xiě)方程但這需要更熟練的集合操作。初學(xué)時(shí)像上面這樣“笨拙”地明確寫(xiě)出每個(gè)方程反而更不容易出錯(cuò)也便于調(diào)試。5. 常見(jiàn)問(wèn)題、調(diào)試技巧與模型優(yōu)化即使模型邏輯正確在LINGO中也可能遇到各種問(wèn)題。以下是一些常見(jiàn)情況及處理辦法。5.1 常見(jiàn)錯(cuò)誤與排查錯(cuò)誤類(lèi)型可能原因排查與解決方法No feasible solution found1. 約束條件相互矛盾。2. 變量上下限設(shè)置過(guò)緊無(wú)解空間。3. 初始值設(shè)置不當(dāng)導(dǎo)致求解器誤判。1.檢查總量平衡如運(yùn)輸問(wèn)題中總供給與總需求關(guān)系。2.簡(jiǎn)化模型先注釋掉部分約束看是否能得到可行解逐步添加以定位矛盾點(diǎn)。3. 檢查BND或變量賦值是否合理。Unbounded solution目標(biāo)函數(shù)可以在不違反約束的情況下無(wú)限優(yōu)化如利潤(rùn)無(wú)限大。這通常是模型構(gòu)建錯(cuò)誤。檢查是否遺漏了關(guān)鍵的資源限制約束。例如在投資問(wèn)題中是否忘記了每年投資總額不能超過(guò)可用資金這一根本約束。求解時(shí)間過(guò)長(zhǎng)或內(nèi)存不足1. 問(wèn)題規(guī)模太大變量/約束過(guò)多。2. 模型是非線(xiàn)性或整數(shù)規(guī)劃本身復(fù)雜。1. 嘗試使用LINGO - Options中的求解器選項(xiàng)如調(diào)整迭代次數(shù)、收斂精度。2. 對(duì)于整數(shù)規(guī)劃可以設(shè)置啟發(fā)式搜索或調(diào)整分支策略。3. 考慮是否能用線(xiàn)性規(guī)劃近似或者分解問(wèn)題。結(jié)果與預(yù)期不符1. 目標(biāo)函數(shù)系數(shù)如成本、利率輸入錯(cuò)誤。2. 約束條件符號(hào),,用錯(cuò)。3. 單位不統(tǒng)一。1.仔細(xì)核對(duì)DATA部分特別是矩陣數(shù)據(jù)的排列順序。2.逐行檢查約束理解其經(jīng)濟(jì)/物理意義。3. 用LINGO - Picture功能可視化模型結(jié)構(gòu)檢查行、列對(duì)應(yīng)關(guān)系。5.2 LINGO編程與調(diào)試技巧多用注釋使用!添加詳細(xì)注釋說(shuō)明每個(gè)集合、變量、約束的含義。這對(duì)于復(fù)雜的、幾天后再回顧的模型至關(guān)重要。模塊化建模對(duì)于大型模型可以分塊編寫(xiě)。先定義集合和數(shù)據(jù)再寫(xiě)目標(biāo)函數(shù)然后分模塊寫(xiě)約束如供給約束塊、需求約束塊、資金平衡塊。每寫(xiě)完一個(gè)模塊可以暫時(shí)注釋掉其他部分單獨(dú)測(cè)試其語(yǔ)法是否正確。利用FOR和SUM簡(jiǎn)化代碼這是LINGO的核心優(yōu)勢(shì)。對(duì)于規(guī)律性強(qiáng)的約束一定要用循環(huán)語(yǔ)句避免手動(dòng)羅列。例如如果我們的投資問(wèn)題有10年、10個(gè)項(xiàng)目手動(dòng)寫(xiě)資金平衡方程將是災(zāi)難。此時(shí)應(yīng)該定義“項(xiàng)目期限”屬性然后利用SUM函數(shù)和條件判斷來(lái)動(dòng)態(tài)生成回收項(xiàng)。! 假設(shè)years是年份集合projects是項(xiàng)目集合term(p)是項(xiàng)目p的期限; FOR(years(t) | t #GT# 1: ! 從第2年開(kāi)始; SUM(projects(p): X(t, p)) SUM(projects(p): (1rate(p)) * X(t - term(p), p)) ! 注意這里需要確保 t-term(p) 1否則會(huì)引用不存在的變量需用IF函數(shù)處理; );查看模型概況求解前點(diǎn)擊LINGO - Generate - Display model可以查看LINGO將你的集合語(yǔ)言展開(kāi)后的完整線(xiàn)性規(guī)劃模型。這是驗(yàn)證模型是否按你設(shè)想的方式構(gòu)建的終極手段。從簡(jiǎn)單案例開(kāi)始不要一開(kāi)始就構(gòu)建復(fù)雜的模型。先用2個(gè)工廠(chǎng)、2個(gè)市場(chǎng)或者2年、2個(gè)項(xiàng)目的小數(shù)據(jù)測(cè)試模型。用手算能驗(yàn)證結(jié)果確認(rèn)模型邏輯無(wú)誤后再替換成真實(shí)的大規(guī)模數(shù)據(jù)。5.3 模型擴(kuò)展與優(yōu)化思路掌握了基礎(chǔ)模型后你可以應(yīng)對(duì)更復(fù)雜的情況運(yùn)輸問(wèn)題的變種產(chǎn)銷(xiāo)不平衡修改供給或需求約束的等號(hào)為不等號(hào)。帶有中轉(zhuǎn)站引入新的節(jié)點(diǎn)集合中轉(zhuǎn)站變量變?yōu)閄_ik(工廠(chǎng)到中轉(zhuǎn)站) 和Y_kj(中轉(zhuǎn)站到市場(chǎng))并增加中轉(zhuǎn)站的流量平衡約束。多商品運(yùn)輸為每種商品增加一個(gè)維度使用三維數(shù)組X_ijk(工廠(chǎng)i, 市場(chǎng)j, 商品k)。投資問(wèn)題的變種風(fēng)險(xiǎn)約束引入風(fēng)險(xiǎn)系數(shù)增加組合風(fēng)險(xiǎn)不超過(guò)某閾值的約束。借貸可能引入“借款”變量利率不同在資金平衡方程中增加借款流入和還款流出項(xiàng)。非線(xiàn)性收益如果收益率與投資額有關(guān)如規(guī)模效應(yīng)則目標(biāo)函數(shù)變?yōu)榉蔷€(xiàn)性需要在LINGO中選用非線(xiàn)性求解器LINGO - Options - Solver中選擇Global或Multistart。最后我個(gè)人在數(shù)學(xué)建模競(jìng)賽培訓(xùn)和實(shí)際項(xiàng)目中最深的體會(huì)是LINGO這類(lèi)工具的強(qiáng)大永遠(yuǎn)建立在你對(duì)問(wèn)題本質(zhì)的深刻理解之上。軟件只是幫你做計(jì)算而建模才是真正的思考?；ㄔ谇逦x變量、梳理約束關(guān)系上的時(shí)間遠(yuǎn)比調(diào)試錯(cuò)誤代碼的時(shí)間更有價(jià)值。當(dāng)你拿到一個(gè)復(fù)雜問(wèn)題不妨先問(wèn)自己決策是什么目標(biāo)是什么限制條件有哪些它們之間如何用數(shù)學(xué)語(yǔ)言關(guān)聯(lián)把這些想明白了翻譯成LINGO代碼就是水到渠成的事。開(kāi)始的時(shí)候可能會(huì)覺(jué)得語(yǔ)法生疏但就像學(xué)任何語(yǔ)言一樣多寫(xiě)幾個(gè)完整的模型你就能越來(lái)越熟練地用它來(lái)表達(dá)你的優(yōu)化思想了。