簡(jiǎn)介YALMIP 官方 R20230609 版資源包是面向 MATLAB 環(huán)境下運(yùn)籌優(yōu)化、自動(dòng)控制及算法研究者的完整建模工具箱可用來(lái)快速構(gòu)建與求解線性規(guī)劃、二次規(guī)劃、混合整數(shù)規(guī)劃、非線性規(guī)劃以及動(dòng)態(tài)最優(yōu)控制問(wèn)題。壓縮包共 1307 個(gè)文件以 1295 個(gè) .m 主程序與輔助函數(shù)為主體涵蓋核心建模語(yǔ)法、約束構(gòu)造、變量類型、求解器接口、預(yù)處理與后處理機(jī)制其余 .txt、.mat、.yml、.c 文件分別提供說(shuō)明文檔、測(cè)試數(shù)據(jù)、運(yùn)行配置與底層實(shí)現(xiàn)便于按需查閱和擴(kuò)展。整個(gè)資源僅 1.08MB輕量且結(jié)構(gòu)清晰。已有 1286 人學(xué)習(xí)/下載。讀者可從中獲得完整的 YALMIP 工具箱主體、常用第三方求解器適配接口以及示例測(cè)試腳本既能從簡(jiǎn)單 LP/QP 開(kāi)始逐步上手也可深入探討動(dòng)態(tài)優(yōu)化、靈敏度分析、隨機(jī)優(yōu)化與多目標(biāo)優(yōu)化等高級(jí)功能適合不同層次的 MATLAB 優(yōu)化開(kāi)發(fā)者使用。1. YALMIP 資源下載包不是裝完即用而是接好一條求解鏈路第一次接觸 YALMIP 資源下載包的人八成是從論文復(fù)現(xiàn)里看到這個(gè)名詞一段 MATLAB 代碼用 sdpvar 聲明變量、用 optimize 下發(fā)求解全程沒(méi)手寫(xiě)矩陣分解。YALMIP 是 MATLAB 環(huán)境下的優(yōu)化建模工具箱把線性規(guī)劃、半定規(guī)劃、混合整數(shù)規(guī)劃統(tǒng)一收攏成幾條聲明與約束語(yǔ)句是問(wèn)題描述與底層求解器之間的翻譯官。不少搜索者以為拿到下載包就是裝完即用實(shí)際上它更像一套要自己接線的設(shè)備。真正值錢(qián)的不只是壓縮包里那幾個(gè)文件而是搞清目錄結(jié)構(gòu)、路徑設(shè)置、求解器接口之后能用同一套語(yǔ)法穩(wěn)定切換不同后端。這篇實(shí)操筆記按一線順序講清下載包放哪、怎么接進(jìn) MATLAB、怎么驗(yàn)證以及五個(gè)高頻翻車(chē)點(diǎn)。適合人群很明確復(fù)現(xiàn)論文卡在安裝環(huán)節(jié)的研究生、想引入這套建?？蚣茏鲈万?yàn)證的工程師以及被下載了卻用不了勸退的入門(mén)者。讀完你會(huì)得到一條可復(fù)現(xiàn)路徑而不是玄學(xué)安裝步驟。2. 拆開(kāi)下載包看架構(gòu)YALMIP 為什么能把建模與求解解耦先說(shuō)一個(gè)反直覺(jué)的結(jié)論YALMIP 資源下載包里沒(méi)有任何代碼負(fù)責(zé)真正的數(shù)值計(jì)算。它做的是建模、轉(zhuǎn)換和調(diào)度重活全部交給下游求解器。把這個(gè)邊界劃清楚后面遇到的大部分安裝報(bào)錯(cuò)都能自己定位。下載包這個(gè)關(guān)鍵詞容易讓人誤以為它是一個(gè)完整求解環(huán)境實(shí)際它更像一個(gè)接線板插頭是你寫(xiě)的優(yōu)化模型電器是求解器。理解這個(gè)架構(gòu)之后你再看它的目錄就不會(huì)覺(jué)得亂。2.1 直接調(diào)求解器的痛為什么要有一層翻譯官假設(shè)要手工求解一個(gè)半定規(guī)劃傳統(tǒng)流程是先把問(wèn)題化成標(biāo)準(zhǔn)形式構(gòu)造目標(biāo)矩陣、約束矩陣和右端項(xiàng)再把這些數(shù)據(jù)結(jié)構(gòu)填進(jìn)求解器接口。三五個(gè)約束還能應(yīng)付到了十幾個(gè)約束、變量本身還是矩陣結(jié)構(gòu)的場(chǎng)景手填下標(biāo)基本就是災(zāi)難。更麻煩的是換求解器格式要全部重寫(xiě)調(diào)試一整天起步。早期做控制理論研究的人相當(dāng)一部分時(shí)間就耗在這種重復(fù)勞動(dòng)上這也是 YALMIP 這類建模工具出現(xiàn)的直接原因——它把這一層體力活標(biāo)準(zhǔn)化了。YALMIP 的具體做法是引入符號(hào)決策變量。sdpvar 聲明的變量在求解前不綁定數(shù)值你可以對(duì)它們做線性組合、求跡、取對(duì)數(shù)可以寫(xiě)矩陣不等式Y(jié)ALMIP 把這些操作維護(hù)成一棵表達(dá)式樹(shù)等 optimize 真正下發(fā)的時(shí)刻再按目標(biāo)求解器的標(biāo)準(zhǔn)格式做轉(zhuǎn)換。于是同一個(gè)模型在不同求解器之間切換只需要改一行配置。你寫(xiě)的模型不再屬于某個(gè)求解器而是屬于問(wèn)題本身。這個(gè)先建模、后下發(fā)的架構(gòu)決定了包內(nèi)必然存在兩條線一條是符號(hào)引擎負(fù)責(zé)變量、約束對(duì)象和表達(dá)式求值另一條是接口層負(fù)責(zé)把統(tǒng)一模型翻譯成特定求解器的輸入格式。很多人安裝時(shí)只檢查主函數(shù)能不能用忽略了接口層與本地環(huán)境的匹配這是后面幾乎所有安裝類報(bào)錯(cuò)的共同根源。2.2 下載包內(nèi)部的關(guān)鍵角色類目錄、接口層與輔助腳本把壓縮包解開(kāi)第一眼會(huì)看到一堆以 開(kāi)頭的目錄。在 MATLAB 里這是類文件夾是整個(gè)符號(hào)語(yǔ)法能工作的根基。這些角色可以用一張表說(shuō)清楚目錄或文件在建模鏈路里的角色根目錄頂層函數(shù)yalmip.m 等入口命令、緩存清理與版本信息sdpvar、constraint 等類目錄符號(hào)變量與約束對(duì)象的類型定義和運(yùn)算符重載求解器接口目錄把統(tǒng)一模型轉(zhuǎn)換成各求解器認(rèn)識(shí)的標(biāo)準(zhǔn)輸入示例與輔助腳本演示用法幫助理解參數(shù)語(yǔ)義類文件夾的名字必須保留原名不能為了好看去改。類方法能否被正常調(diào)用完全依賴 MATLAB 的類目錄機(jī)制一旦移動(dòng)或改名構(gòu)造函數(shù)立刻失效報(bào)錯(cuò)還經(jīng)常是未定義函數(shù)或變量這種與直覺(jué)相反的信息。這就是為什么教程反復(fù)強(qiáng)調(diào)用 genpath 遞歸加路徑只加根目錄時(shí) which yalmip 是正常的但 sdpvar 用不了讓人百思不得其解。目錄解壓的位置同樣有講究。有人圖省事直接把包解進(jìn) MATLAB 安裝目錄內(nèi)部看似路徑不用配了實(shí)則埋雷。MATLAB 大版本升級(jí)時(shí)會(huì)對(duì)內(nèi)置目錄做掃描和重置第三方文件混在中間很容易出現(xiàn)路徑優(yōu)先級(jí)上的詭異沖突。我一般單獨(dú)建一個(gè)第三方工具箱目錄解壓后整棵目錄樹(shù)有據(jù)可查備份回滾都是一條命令的事。2.3 多數(shù)人缺的不是 YALMIP 本體而是求解器后端資源下載包這個(gè)搜索詞背后真實(shí)常見(jiàn)的情景是包也解壓了路徑也加了一跑卻提示找不到可用求解器。這不是包壞了而是機(jī)器上沒(méi)有任何求解器。YALMIP 只負(fù)責(zé)把模型翻譯成標(biāo)準(zhǔn)形式真正做迭代計(jì)算的是求解器。它刻意不捆綁求解器為的是保持前端工具的中立性代價(jià)就是新用戶必須先補(bǔ)齊后端。常見(jiàn)的可靠做法是先裝一個(gè)開(kāi)源免費(fèi)求解器把它的目錄也加入 MATLAB 路徑再重跑 yalmiptest 確認(rèn)接口被識(shí)別。選擇依據(jù)是問(wèn)題類型——只解線性規(guī)劃一個(gè)基礎(chǔ)求解器就夠涉及半定規(guī)劃必須選支持半定約束的混合整數(shù)規(guī)劃則需要支持分支定界機(jī)制的求解器。拿名氣去選不如拿問(wèn)題特征去選。等前端和后端組裝到一起這個(gè)下載包才真正變成能出數(shù)值結(jié)果的工具鏈。3. 把資源包接進(jìn) MATLAB路徑配置與最小驗(yàn)證YALMIP 沒(méi)有安裝程序?qū)儆诘湫偷木G色工具箱。安裝的本質(zhì)工作只有一個(gè)讓 MATLAB 的搜索路徑覆蓋下載包的所有目錄。把這句話記住后面所有路徑類報(bào)錯(cuò)都能自己推導(dǎo)出答案不用到處問(wèn)人。3.1 目錄規(guī)劃給第三方工具箱一個(gè)固定位置我通常會(huì)在工作區(qū)建一個(gè)專門(mén)的第三方工具箱目錄所有這類下載來(lái)的工具箱都放同一個(gè)地方不跟業(yè)務(wù)代碼混在一起。好處很明顯備份時(shí)一條命令打包出問(wèn)題時(shí)整體移除換新版本時(shí)目標(biāo)明確不會(huì)出現(xiàn)這個(gè)文件夾到底是不是正在用的這種靈魂拷問(wèn)。mkdir -p ~/tools/third_party unzip yalmip-download.zip -d ~/tools/third_partyWindows 用戶用解壓軟件或 PowerShell 的 Expand-Archive 完成同樣操作。解壓完先看一眼目錄結(jié)構(gòu)如果出現(xiàn) yalmip/yalmip 這樣的雙層嵌套說(shuō)明壓縮包內(nèi)層還有一個(gè)同名頂層目錄先整理成單層根目錄再進(jìn)行下一步。這一步往往決定了后面路徑配置是對(duì)是錯(cuò)別急著往下走。3.2 路徑配置addpath、genpath 與 savepath 的配合% 遞歸添加 YALMIP 根目錄到搜索路徑 addpath(genpath(D:\tools\third_party\yalmip)); % 持久化保存下次啟動(dòng) MATLAB 依然有效 savepath; % 驗(yàn)證入口函數(shù)是否可見(jiàn) which yalmip參數(shù)說(shuō)明genpath 的關(guān)鍵是遞歸它把根目錄下所有子目錄展開(kāi)后一次性加入路徑。YALMIP 的核心功能大量放在 開(kāi)頭的類目錄里只加根目錄會(huì)漏掉這些類目錄sdpvar、constraint 等構(gòu)造函數(shù)全部不可見(jiàn)。savepath 把當(dāng)前路徑寫(xiě)入啟動(dòng)配置文件但在某些環(huán)境下會(huì)因?qū)憴?quán)限失敗報(bào)錯(cuò)時(shí)先檢查當(dāng)前用戶對(duì)該配置文件所在目錄的寫(xiě)權(quán)限受限時(shí)可以在 MATLAB 啟動(dòng)腳本里加 addpath 語(yǔ)句效果等同。which yalmip 能返回完整路徑是安裝通過(guò)的第一個(gè)標(biāo)志但不是全部。3.3 最小驗(yàn)證清緩存、測(cè)接口、看報(bào)告% 清掉 YALMIP 內(nèi)部緩存換路徑或換版本后必須執(zhí)行 yalmip(clear); % 運(yùn)行自帶的測(cè)試套件檢查各求解器接口 yalmiptest邏輯說(shuō)明yalmip(clear) 重置 YALMIP 內(nèi)部狀態(tài)特別是路徑變更之后殘留的舊變量定義不清理會(huì)出現(xiàn)難以復(fù)現(xiàn)的詭異報(bào)錯(cuò)。yalmiptest 會(huì)依次跑若干類最小測(cè)試問(wèn)題并在輸出里報(bào)告每個(gè)問(wèn)題對(duì)應(yīng)的求解器接口狀態(tài)?？唇Y(jié)果時(shí)別盯著失敗兩個(gè)字直接崩潰要分類看大面積失敗大概率是路徑或求解器缺失只有半定規(guī)劃、混合整數(shù)那幾項(xiàng)失敗說(shuō)明當(dāng)前求解器的能力覆蓋不全需要換支持對(duì)應(yīng)問(wèn)題類型的后端。這一步通過(guò)之后安裝環(huán)節(jié)才算真正結(jié)束。4. 用下載包的語(yǔ)法跑通三類問(wèn)題LP、SDP 與 MILP 最小模板路徑接好只是入場(chǎng)券能不能用順手取決于建模語(yǔ)法是否熟練。下面三張模板覆蓋了 YALMIP 最主流的應(yīng)用場(chǎng)景新項(xiàng)目開(kāi)工時(shí)可以直接抄改改變量和約束就變成你自己的模型。4.1 LP線性規(guī)劃的最小模板% 聲明決策變量2 行 1 列的實(shí)數(shù)符號(hào)向量 x sdpvar(2, 1); % 目標(biāo)函數(shù)f x1 2*x2 objective x(1) 2*x(2); % 約束x1 x2 1且 x 各分量非負(fù) constraints [x(1) x(2) 1, x 0]; % 求解最小化 objective返回求解狀態(tài) optimize(constraints, objective); % 提取數(shù)值解 value(x)邏輯說(shuō)明sdpvar(2,1) 生成一個(gè)二行一列的符號(hào)向量默認(rèn)類型是連續(xù)實(shí)數(shù)。constraints 用方括號(hào)把多個(gè)不等式拼接成一個(gè)約束數(shù)組這是 YALMIP 里最常見(jiàn)的組合方式。optimize 默認(rèn)做最小化第一個(gè)參數(shù)傳約束集合第二個(gè)參數(shù)傳目標(biāo)表達(dá)式。value(x) 在求解完成后把符號(hào)變量替換成數(shù)值向量不調(diào)用 valuex 永遠(yuǎn)保持符號(hào)狀態(tài)——不少新手在這里直接打印 x看到一堆符號(hào)表達(dá)式以為求解失敗了實(shí)際上只是沒(méi)取數(shù)值。模型里沒(méi)有二次項(xiàng)、沒(méi)有整數(shù)變量YALMIP 會(huì)自動(dòng)識(shí)別為線性規(guī)劃并走對(duì)應(yīng)的求解路徑。如果你還需要目標(biāo)函數(shù)對(duì)應(yīng)的對(duì)偶變量可以用 dual(constraints) 取回拉格朗日乘子這是做敏感性分析時(shí)很常用的手段。4.2 SDP矩陣變量與半定約束的正確寫(xiě)法% 聲明 2x2 對(duì)稱矩陣符號(hào)變量 X sdpvar(2, 2, symmetric); % 半定約束與等式約束X 0 且跡等于 1 constraints [X 0, trace(X) 1]; % 最小化 -logdet(X)等價(jià)于最大化 log det(X) optimize(constraints, -logdet(X)); % 查看最優(yōu) X value(X)參數(shù)說(shuō)明sdpvar 的第三個(gè)參數(shù)傳 symmetricYALMIP 就把 X 當(dāng)作對(duì)稱矩陣處理內(nèi)部只保留上三角元素避免重復(fù)自由變量拖慢求解。這里的 X 0 在矩陣變量下是半定約束含義是 X 的特征值全部非負(fù)——這是一個(gè)經(jīng)常被誤解的語(yǔ)法點(diǎn)它不等于逐元素大于零。想要逐元素非負(fù)必須顯式寫(xiě) X(:) 0兩者語(yǔ)義完全不同混用會(huì)讓模型面目全非。logdet 是 YALMIP 內(nèi)置算子專門(mén)處理半定規(guī)劃里的對(duì)數(shù)行列式目標(biāo)。它會(huì)自動(dòng)對(duì)模型做重構(gòu)以保證凸性這塊也正是手寫(xiě)標(biāo)準(zhǔn)形式時(shí)最容易出錯(cuò)、最耗時(shí)間的部分。用 YALMIP 寫(xiě)這類問(wèn)題的體驗(yàn)跟手填矩陣相比基本是降維打擊。4.3 MILP離散變量用 binvar% 3 個(gè) 0/1 變量3 個(gè)連續(xù)變量 x binvar(3, 1); y sdpvar(3, 1); % 約束y 落在 [0, 10*x]x0 時(shí)強(qiáng)制 y0 constraints [y 0, y 10*x, sum(y) 5]; % 最大化 y 的總和寫(xiě)成最小化 -sum(y) objective -sum(y); % 求解讓 YALMIP 自動(dòng)選擇支持混合整數(shù)的求解器 optimize(constraints, objective);邏輯說(shuō)明binvar 生成 0/1 離散變量intvar 生成一般整數(shù)變量與 sdpvar 的唯一區(qū)別就在變量類型?；旌险麛?shù)問(wèn)題對(duì)求解器能力要求更高如果當(dāng)前環(huán)境沒(méi)有支持分支定界的求解器optimize 會(huì)直接報(bào)錯(cuò)而不會(huì)退化成連續(xù)求解。這個(gè)寧可報(bào)錯(cuò)也不靜默降級(jí)的行為比某些自動(dòng)忽略整數(shù)約束的做法安全得多。約束里的 y 10*x 是典型的 big-M 寫(xiě)法10 是業(yè)務(wù)上給出的上界。M 值要按實(shí)際問(wèn)題取合理上界不是越大越安全——過(guò)大的 M 會(huì)顯著惡化數(shù)值條件導(dǎo)致分支定界過(guò)程反復(fù)震蕩求解時(shí)間成倍上漲。4.4 參數(shù)下發(fā)sdpsettings 是唯一入口% 自動(dòng)選擇求解器關(guān)閉冗長(zhǎng)輸出開(kāi)啟診斷信息 options sdpsettings(solver, , verbose, 0, debug, 1); % 用兩個(gè)返回值接收sol 是求解狀態(tài)diagnostic 是診斷信息 [sol, diagnostic] optimize(constraints, objective, options); % 狀態(tài)碼為 0 表示求解成功 if sol.problem ~ 0 disp(diagnostic.info); end參數(shù)說(shuō)明sdpsettings 返回一個(gè)結(jié)構(gòu)體集中管理所有求解參數(shù)。solver 字段傳空字符串表示自動(dòng)探測(cè)合適求解器手動(dòng)指定求解器名可以鎖定求解環(huán)境復(fù)現(xiàn)論文結(jié)果時(shí)一般都會(huì)鎖防止不同求解器算法差異造成數(shù)值對(duì)不上。verbose 控制命令行輸出量調(diào)試階段設(shè) 1批量跑實(shí)驗(yàn)設(shè) 0否則日志文件會(huì)迅速膨脹。debug 置 1 時(shí)一旦求解器返回不可行或異常YALMIP 會(huì)額外輸出約束層面的分析信息這是排查不可行模型最直接的工具。注意 options 必須放在 optimize 的第三個(gè)參數(shù)位置漏傳等于所有設(shè)置沒(méi)生效。這個(gè)錯(cuò)誤很隱蔽因?yàn)榇a不報(bào)錯(cuò)只是行為不對(duì)。5. YALMIP 安裝與使用避坑指南五個(gè)高頻翻車(chē)現(xiàn)場(chǎng)下面五條按真實(shí)出現(xiàn)頻率排序前兩條集中在安裝環(huán)節(jié)后三條集中在建模環(huán)節(jié)。每一條我都見(jiàn)過(guò)至少一次其中一條自己也貢獻(xiàn)過(guò)。5.1 which yalmip 能找到但 sdpvar 報(bào)未定義函數(shù)現(xiàn)象完成路徑配置后which yalmip 能正常返回路徑可是一運(yùn)行 sdpvarMATLAB 直接報(bào)未定義函數(shù)或變量。原因路徑只加到了根目錄沒(méi)有遞歸覆蓋以 開(kāi)頭的類目錄。yalmip.m 在根目錄所以可見(jiàn)但 sdpvar 的構(gòu)造函數(shù)定義在 sdpvar 類目錄里這個(gè)目錄沒(méi)進(jìn)搜索路徑。解決改用 genpath 遞歸添加再刷新緩存。addpath(genpath(D:\tools\third_party\yalmip)); savepath; rehash toolboxcache;這個(gè)坑的隱蔽點(diǎn)在于入口函數(shù)可見(jiàn)給人安裝成功的錯(cuò)覺(jué)。判斷安裝是否完成標(biāo)準(zhǔn)是 sdpvar、constraint 這類符號(hào)函數(shù)的可見(jiàn)性而不是 which yalmip 的輸出。5.2 optimize 報(bào)錯(cuò)找不到任何可用的求解器現(xiàn)象第一次跑 optimize命令行提示找不到合適的求解器但 YALMIP 本身工作正常yalmiptest 也能運(yùn)行。此時(shí)黑匣子感覺(jué)最強(qiáng)最容易讓人懷疑是包的問(wèn)題。原因機(jī)器上只有 YALMIP沒(méi)有任何求解器或者裝的求解器與當(dāng)前 MATLAB 版本不兼容在接口探測(cè)階段就被過(guò)濾掉了。解決安裝一個(gè)開(kāi)源免費(fèi)求解器加入路徑后重新運(yùn)行 yalmiptest確認(rèn)對(duì)應(yīng)測(cè)試項(xiàng)從失敗變通過(guò)。如果接口仍不被識(shí)別檢查求解器位數(shù)與 MATLAB 位數(shù)是否一致同時(shí)確認(rèn)它要求的 MATLAB 版本范圍。這里要認(rèn)清一個(gè)事實(shí)YALMIP 下載包不捆綁任何求解器是設(shè)計(jì)使然新用戶必須自己補(bǔ)齊后端。5.3 約束方向沒(méi)寫(xiě)錯(cuò)但結(jié)果和論文差幾個(gè)數(shù)量級(jí)現(xiàn)象公式能和論文對(duì)上跑出來(lái)的最優(yōu)值卻完全不同甚至符號(hào)都不對(duì)反復(fù)核對(duì)也看不出問(wèn)題。原因矩陣變量下的 X 0 實(shí)際是半定約束卻被當(dāng)成了逐元素非負(fù)或者不等號(hào)方向與預(yù)期理解反了。前者直接把約束集合換掉后者會(huì)讓對(duì)偶解完全跑偏。解決用 value 取出解后再調(diào)用 check 看約束殘差。residuals check(constraints);check 會(huì)在當(dāng)前解上逐條計(jì)算殘差殘差數(shù)量級(jí)明顯異常的那條約束就是重點(diǎn)懷疑對(duì)象。同時(shí)確認(rèn)變量聲明方式sdpvar(2,2,symmetric) 下 X 0 是半定約束想表達(dá)逐元素非負(fù)顯式寫(xiě) X(:) 0。這兩個(gè)寫(xiě)法語(yǔ)義完全不同混用的結(jié)果就是模型對(duì)、約束錯(cuò)、數(shù)值自然對(duì)不上。排查這個(gè)問(wèn)題的思路是先懷疑約束語(yǔ)義再懷疑求解器精度順序別反。5.4 MATLAB 升級(jí)后舊腳本開(kāi)始莫名報(bào)錯(cuò)現(xiàn)象同一個(gè)資源包MATLAB 大版本升級(jí)之后前一天能跑的腳本開(kāi)始拋錯(cuò)報(bào)錯(cuò)信息指向一個(gè)以前沒(méi)注意過(guò)的函數(shù)。原因新版本 MATLAB 自帶功能與 YALMIP 部分函數(shù)重名路徑搜索順序?qū)е赂采w或者 YALMIP 版本太舊使用了新環(huán)境已移除的舊版 API。解決先執(zhí)行 yalmip(clear) 清理內(nèi)部緩存再用 which 查看沖突函數(shù)實(shí)際來(lái)自哪個(gè)目錄。如果確認(rèn)是舊 API 問(wèn)題升級(jí) YALMIP 資源包并把代碼里的舊式入口調(diào)用替換為當(dāng)前版本的 optimize 寫(xiě)法。換包之后必須重新 addpath(genpath(...))——只替換文件不更新路徑等于沒(méi)換這是最容易犯的遺漏。5.5 解壓出 yalmip/yalmip 雙層嵌套路徑怎么配都別扭現(xiàn)象解壓后看到兩層同名目錄加入路徑后測(cè)試時(shí)好時(shí)壞有些功能能用有些功能隨機(jī)報(bào)錯(cuò)。原因壓縮包內(nèi)有同名頂層目錄解壓時(shí)與本地目錄合并產(chǎn)生了嵌套或者解壓軟件選擇解壓到當(dāng)前文件夾時(shí)多生成了一層。解決先整理成單一根目錄。mv ~/tools/third_party/yalmip/yalmip ~/tools/third_party/確認(rèn)最終只有一條目錄包含 yalmip.m 主文件然后重新 addpath(genpath(...))并運(yùn)行 yalmip(clear) 重測(cè)。這個(gè)坑的隱蔽點(diǎn)在于 which yalmip 能返回內(nèi)層路徑給人路徑?jīng)]問(wèn)題的假象但其他類目錄的映射已經(jīng)錯(cuò)位報(bào)錯(cuò)隨機(jī)且難復(fù)現(xiàn)看起來(lái)像玄學(xué)本質(zhì)上就是目錄結(jié)構(gòu)問(wèn)題。6. 把資源包用出效率三個(gè)建模習(xí)慣與版本管理安裝只是開(kāi)始真正省時(shí)間的是后面的建模習(xí)慣。第一個(gè)習(xí)慣是調(diào)試不可行模型時(shí)先賦值再看殘差而不是反復(fù)改約束重跑 optimize 碰運(yùn)氣% 給變量賦一個(gè)試探性的可行點(diǎn) assign(x, [1; 1]); % 在該點(diǎn)下逐條檢查約束殘差 check(constraints);assign 把符號(hào)變量臨時(shí)替換成數(shù)值check 在這個(gè)點(diǎn)上計(jì)算每條約束的違反程度。哪條約束把模型逼近不可行一眼就能定位比盲目重跑高效得多。排查不可行問(wèn)題時(shí)這個(gè)組合拳比任何日志都直觀。第二個(gè)習(xí)慣是鎖定求解環(huán)境。復(fù)現(xiàn)論文或?qū)ε慕Y(jié)果時(shí)把 sdpsettings 寫(xiě)進(jìn)腳本開(kāi)頭并指定求解器別依賴默認(rèn)自動(dòng)選擇——兩個(gè)求解器之間的算法差異足以讓對(duì)比實(shí)驗(yàn)得出完全相反的結(jié)論。批量實(shí)驗(yàn)時(shí)把 verbose 設(shè) 0日志從幾百行縮成幾行跑批速度也有肉眼可見(jiàn)的提升。第三個(gè)習(xí)慣是資源包的版本管理。目錄命名帶上日期舊版本整體改名備份而不是覆蓋換版本按固定順序執(zhí)行備份舊目錄、解壓新包、更新路徑、yalmip(clear)、rehash toolboxcache、yalmiptest。順序一亂最常見(jiàn)的結(jié)果就是新舊文件混在同一個(gè)類目錄里報(bào)錯(cuò)隨機(jī)到你懷疑人生。我自己就吃過(guò)這個(gè)虧。有一回為了省事把兩個(gè)版本混在同一個(gè)目錄里某天所有約束開(kāi)始間歇性報(bào)錯(cuò)排查了整整一個(gè)晚上最后發(fā)現(xiàn)是 sdpvar 類目錄被舊文件覆蓋了一半。從那以后我裝任何第三方工具箱都先建獨(dú)立目錄、再記錄版本日期寧可多敲兩行命令也不把時(shí)間浪費(fèi)在路徑和目錄的隨機(jī)性上。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取