寫MATLAB程序這件事我前后寫了小十年。從最早用它做課程作業(yè)里的矩陣運(yùn)算到后來處理現(xiàn)場采集的傳感器數(shù)據(jù)、搭建仿真模型、甚至給別人做工具包可以說MATLAB既是入門最容易的工具也是最能暴露一個(gè)人編程習(xí)慣的鏡子。這個(gè)系列走到第六篇我想認(rèn)真聊聊“程序設(shè)計(jì)”本身——不是某個(gè)工具箱的某個(gè)函數(shù)怎么用而是當(dāng)你面對(duì)一個(gè)需要長期維護(hù)、反復(fù)修改、甚至交給別人接著開發(fā)的MATLAB項(xiàng)目時(shí)那些真正決定成敗的細(xì)節(jié)。這篇文章適合兩類人一類是已經(jīng)能用MATLAB跑通腳本、但總覺得代碼“能跑但很亂”的進(jìn)階初學(xué)者另一類是在其他語言里寫慣了嚴(yán)謹(jǐn)程序、剛轉(zhuǎn)到MATLAB、被它的靈活弄得不知所措的開發(fā)者。我會(huì)從工程級(jí)代碼的組織方式、控制流與向量化的取舍、調(diào)試與性能優(yōu)化的實(shí)際手法、再到把代碼變成可復(fù)用工具的封裝思路一層層講清楚。里面所有經(jīng)驗(yàn)都來自實(shí)際項(xiàng)目中的試錯(cuò)和總結(jié)不是教科書里的標(biāo)準(zhǔn)答案但每一招我都親自驗(yàn)證過。1. 工程級(jí)代碼的起點(diǎn)先把編程風(fēng)格立起來1.1 認(rèn)清MATLAB的兩種程序形態(tài)初學(xué)者經(jīng)常搞混一件事腳本和函數(shù)到底有什么區(qū)別。簡單說腳本就是一段按順序執(zhí)行的命令集合它操作的是工作區(qū)里的現(xiàn)有變量而函數(shù)有自己獨(dú)立的工作區(qū)通過輸入?yún)?shù)和輸出參數(shù)與外界交換數(shù)據(jù)。我見過太多人把幾百行代碼全塞在一個(gè)腳本里跑完一遍工作區(qū)堆了幾十個(gè)變量第二天想改個(gè)參數(shù)得從頭捋一遍才知道哪個(gè)變量是干嘛的。這就是典型的“腳本思維”。工程級(jí)代碼的第一步就是學(xué)會(huì)把功能拆成函數(shù)。舉個(gè)實(shí)際例子。假設(shè)你要處理一批實(shí)驗(yàn)數(shù)據(jù)讀取文件、去噪、提取特征、畫圖。用腳本寫全部擠在一起用函數(shù)拆結(jié)構(gòu)就非常清楚data loadData(experiment_01.csv); % 讀取 clean removeNoise(data, median, 5); % 去噪 features extractFeatures(clean); % 提取特征 plotResults(clean, features); % 畫圖每一行都是一個(gè)函數(shù)調(diào)用你不需要關(guān)注每個(gè)函數(shù)內(nèi)部怎么實(shí)現(xiàn)只要確認(rèn)輸入輸出符合約定就行。這就是模塊化的意義代碼的可讀性不再靠“讀”而是靠“信任函數(shù)名”。1.2 函數(shù)編寫的關(guān)鍵約定寫函數(shù)的時(shí)候有幾個(gè)習(xí)慣我從第三年才開始堅(jiān)持但一旦堅(jiān)持下來返工率直線下降每個(gè)函數(shù)只做一件事。如果一個(gè)函數(shù)既做濾波又做特征提取還順帶畫圖那它就是一個(gè)定時(shí)炸彈。日后任何一處需要修改你都得小心翼翼生怕動(dòng)一個(gè)功能影響另一個(gè)。參數(shù)校驗(yàn)要放在函數(shù)開頭。別嫌麻煩。validateattributes或assert幾行代碼能省掉你后面排查半天“為什么這次數(shù)據(jù)跑出來全是NaN”的時(shí)間。優(yōu)先使用varargin做可選參數(shù)。配合inputParser你可以寫出非常優(yōu)雅的參數(shù)接口function result processData(data, varargin) p inputParser; addRequired(p, data, isnumeric); addParameter(p, Window, 5, (x) isnumeric(x) x 0); addParameter(p, Method, median, (x) ismember(x, {median,mean})); parse(p, data, varargin{:}); window p.Results.Window; method p.Results.Method; % ... end調(diào)用的時(shí)候可以寫成processData(data, Window, 10, Method, mean)看代碼的人一眼就知道每個(gè)參數(shù)在干什么。這種命名參數(shù)的方式比單純按位置傳參要友好得多。1.3 命名規(guī)范與區(qū)塊化注釋MATLAB的編輯器自帶%%區(qū)塊功能這是一個(gè)被嚴(yán)重低估的生產(chǎn)力工具。用%%把代碼分成邏輯區(qū)塊每個(gè)區(qū)塊可以單獨(dú)運(yùn)行調(diào)試時(shí)不用從頭跑一遍。我習(xí)慣在每個(gè)區(qū)塊開頭寫一行注釋說明這段在干什么。配合section標(biāo)題整個(gè)腳本就像一本書的目錄。命名方面我個(gè)人遵循一套簡單約定函數(shù)名用小寫開頭的駝峰式如extractFeatures變量名用有意義的英文單詞組合如rawData、filteredSignal常量用全大寫如SAMPLE_RATE。這套規(guī)范堅(jiān)持了五年最大的好處是代碼不需要額外文檔——名字自己就是文檔。提示不要在變量名里用拼音縮寫。當(dāng)時(shí)寫的時(shí)候覺得方便三個(gè)月后你自己都看不懂“sjcl”是“數(shù)據(jù)處理”還是“時(shí)間序列”。英文單詞即使拼錯(cuò)也比拼音好猜。2. 控制流與向量化告別又慢又蠢的 for 循環(huán)2.1 條件分支和循環(huán)的進(jìn)階寫法MATLAB里的if-elseif-else和switch-case表面上是常識(shí)但實(shí)際項(xiàng)目里分支條件往往不是簡單的大于小于而是復(fù)合邏輯。這時(shí)候務(wù)必把最可能命中的分支放在最前面——不是玄學(xué)這是利用分支預(yù)測減少判斷次數(shù)的通用優(yōu)化思路。循環(huán)方面for和while各有適用場景。迭代次數(shù)明確、可預(yù)知的用for依賴條件終止、次數(shù)不確定的用while。但無論是哪種只要能用向量化解決的問題就盡量不要用循環(huán)解決。下面這個(gè)對(duì)比我做了無數(shù)次每次都能把人說服n 1e6; data randn(n, 1); % 循環(huán)寫法 tic; result1 zeros(n, 1); for i 1:n result1(i) data(i)^2 3*data(i) 1; end toc; % 向量化寫法 tic; result2 data.^2 3*data 1; toc;在我目前用的機(jī)器上循環(huán)版本大概耗時(shí)0.05秒向量化版本0.003秒差距接近一個(gè)數(shù)量級(jí)。當(dāng)數(shù)據(jù)量再擴(kuò)大十倍、百倍這個(gè)差距會(huì)變得更加明顯。MATLAB的底層矩陣運(yùn)算庫是高度優(yōu)化的你寫data.^2 3*data 1時(shí)它調(diào)用的是一連串經(jīng)過手冊(cè)級(jí)優(yōu)化的C/Fortran例程而你自己寫循環(huán)則是讓MATLAB解釋器逐行逐元素地跑——兩者根本不在一個(gè)量級(jí)。2.2 向量化的核心邏輯逐元素操作與邏輯索引向量化的本質(zhì)是把“對(duì)每個(gè)元素做操作”的思維替換成“對(duì)整個(gè)數(shù)組做操作”的思維。一開始不習(xí)慣很正常我花了很久才真正轉(zhuǎn)變過來。但掌握兩個(gè)工具后多數(shù)場景都能直接化解第一個(gè)是逐元素操作符。.*、./、.^這三個(gè)帶點(diǎn)的運(yùn)算符專門處理數(shù)組逐元素運(yùn)算。寫多項(xiàng)式擬合、歸一化、逐點(diǎn)運(yùn)算它們都是主力。很多人容易漏掉點(diǎn)號(hào)寫成data^2就會(huì)報(bào)錯(cuò)——因?yàn)閊是矩陣乘冪和逐元素的.^是兩回事。這個(gè)坑每個(gè)MATLAB用戶都踩過。第二個(gè)是邏輯索引。這是MATLAB最獨(dú)特也最強(qiáng)大的特性之一。你可以直接用邏輯數(shù)組做索引data randn(1000, 1); outliers data(abs(data) 2); % 提取所有絕對(duì)值大于2的元素 data(abs(data) 2) 0; % 把異常值置零 condition data 0 data 1; selected data(condition); % 滿足條件的結(jié)果這種寫法不僅代碼短而且執(zhí)行效率遠(yuǎn)高于在循環(huán)里加if判斷。核心原因在于邏輯索引直接映射到內(nèi)存層面的批量操作而循環(huán)內(nèi)的條件判斷會(huì)讓解釋器頻繁分支跳轉(zhuǎn)。2.3 什么時(shí)候必須用循環(huán)向量化不是銀彈。有些場景用循環(huán)反而更合理迭代過程依賴前一步結(jié)果比如遞推公式、狀態(tài)轉(zhuǎn)移這種天然就是串行的硬要向量化只能寫出極其晦澀的代碼。循環(huán)體內(nèi)有復(fù)雜的條件分支每個(gè)分支的處理邏輯完全不同強(qiáng)行向量化需要構(gòu)造多層掩碼矩陣可讀性急劇下降。需要在循環(huán)中動(dòng)態(tài)改變數(shù)組大小——當(dāng)然這種情況建議先預(yù)分配見后面的性能優(yōu)化部分。我個(gè)人的判斷標(biāo)準(zhǔn)很簡單如果量級(jí)在百萬級(jí)以上且沒有數(shù)據(jù)依賴優(yōu)先向量化如果量級(jí)小、邏輯復(fù)雜、或者涉及串行遞推直接用循環(huán)不用糾結(jié)。代碼清晰和運(yùn)行速度之間需要你自己權(quán)衡不必走極端。3. 調(diào)試實(shí)戰(zhàn)讓代碼自己“說話”3.1 斷點(diǎn)調(diào)試的正確姿勢很多人調(diào)試MATLAB代碼還是靠到處插disp打印變量。這在腳本時(shí)代沒問題但工程代碼一旦龐大打印信息會(huì)淹沒在輸出窗口里而且你得反復(fù)刪除調(diào)試代碼。正確做法是用斷點(diǎn)調(diào)試。在編輯器左側(cè)點(diǎn)擊行號(hào)旁的空格就設(shè)了一個(gè)斷點(diǎn)。程序運(yùn)行到這里會(huì)暫停這時(shí)你可以在“編輯器”選項(xiàng)卡里查看工作區(qū)懸停變量看當(dāng)前值甚至直接在命令行窗口輸入代碼來測試下一步的計(jì)算。配合step單步和step in進(jìn)入函數(shù)內(nèi)部你可以像放慢電影一樣逐幀觀察程序的每一個(gè)動(dòng)作。實(shí)際項(xiàng)目中我最常用的組合是先在主函數(shù)入口設(shè)斷點(diǎn)確認(rèn)輸入?yún)?shù)正確然后在我懷疑出問題的模塊前設(shè)斷點(diǎn)運(yùn)行后檢查中間變量再配合step in進(jìn)入函數(shù)內(nèi)部一步步看到底是哪一行產(chǎn)生了錯(cuò)誤結(jié)果。整個(gè)過程不需要改動(dòng)任何代碼排查完直接繼續(xù)運(yùn)行調(diào)試痕跡為零。3.2 條件斷點(diǎn)省掉重復(fù)的F5有一種情況特別煩人循環(huán)跑幾千次到第158次才出錯(cuò)。你不可能每次都手動(dòng)跑到那一次再斷下來。條件斷點(diǎn)就是為了解決這個(gè)問題的。在斷點(diǎn)圖標(biāo)上右鍵選擇“設(shè)置條件”輸入類似i 158或者更復(fù)雜的條件abs(value - 0.618) 1e-6滿足條件時(shí)程序才會(huì)暫停。這個(gè)功能能讓排查效率提升一個(gè)量級(jí)。我在處理長循環(huán)里的偶發(fā)錯(cuò)誤時(shí)幾乎全靠它。注意條件斷點(diǎn)里的表達(dá)式每次迭代都會(huì)求值所以不要寫太復(fù)雜的內(nèi)容否則你自己會(huì)變成性能瓶頸。3.3 性能瓶頸定位從猜想到實(shí)測程序“有點(diǎn)慢”是模糊的說法要優(yōu)化它得先知道時(shí)間花在哪。MATLAB自帶的profile工具就是干這個(gè)的。運(yùn)行前敲profile on myScript % 或者 myFunction(...) profile viewer它會(huì)生成一份詳細(xì)報(bào)告列出每個(gè)函數(shù)的調(diào)用次數(shù)、總耗時(shí)、自身耗時(shí)。我見過太多人憑感覺優(yōu)化代碼——把某個(gè)函數(shù)重寫了一遍結(jié)果才發(fā)現(xiàn)瓶頸根本不在這里。在實(shí)際項(xiàng)目里我一般遵循這樣的流程先profile跑一遍看耗時(shí)排名揪出前幾名再針對(duì)性地分析是算法問題、數(shù)據(jù)結(jié)構(gòu)問題、還是循環(huán)問題修改后再跑一次profile對(duì)比。優(yōu)化的前提是量化量化之后才不會(huì)做無用功。3.4 異常處理與“試錯(cuò)”邏輯工程級(jí)代碼必須考慮異常情況文件不存在、數(shù)據(jù)格式不對(duì)、除零、矩陣維度不匹配。裸奔的代碼遇到這些情況要么直接報(bào)錯(cuò)崩潰要么給出一個(gè)完全看不懂的提示。try-catch可以讓你優(yōu)雅地處理這些問題try data readtable(data.xlsx); catch ME warning(讀取文件失敗%s, ME.message); data []; % 這里可以執(zhí)行備用方案 endME是MException對(duì)象里面包含錯(cuò)誤類型、錯(cuò)誤消息、錯(cuò)誤堆棧。用warning而不是error可以讓程序繼續(xù)運(yùn)行但提醒你出了問題。這樣在批量處理多個(gè)文件時(shí)一個(gè)文件出錯(cuò)不會(huì)中斷整個(gè)任務(wù)最后你還能匯總所有警告統(tǒng)一查看。還應(yīng)該養(yǎng)成一個(gè)習(xí)慣對(duì)外提供的函數(shù)入口要給出友好報(bào)錯(cuò)。在函數(shù)最前面檢查輸入一旦不滿足要求就用error拋出一個(gè)帶說明的錯(cuò)誤。這個(gè)習(xí)慣能讓調(diào)用者第一時(shí)間明白問題出在哪而不是被一個(gè)莫名其妙的維度錯(cuò)誤帶跑偏。4. 性能優(yōu)化從“能跑”到“跑得漂亮”4.1 預(yù)分配最簡單但最有效的優(yōu)化MATLAB的數(shù)組如果頻繁擴(kuò)展會(huì)反復(fù)重新分配內(nèi)存?？催@個(gè)例子tic; result []; % 空數(shù)組開始 for i 1:100000 result(i) i^2; % 每次擴(kuò)展都可能重新分配內(nèi)存 end toc; tic; result zeros(100000, 1); % 預(yù)先分配好空間 for i 1:100000 result(i) i^2; end toc;在我的機(jī)器上預(yù)分配版本幾乎快一個(gè)數(shù)量級(jí)。原因是當(dāng)數(shù)組長度變化時(shí)MATLAB需要在內(nèi)存中尋找更大的連續(xù)區(qū)塊并拷貝舊數(shù)據(jù)這個(gè)操作的時(shí)間成本隨數(shù)組長度增長而增長。如果你預(yù)先知道數(shù)組最終大小直接zeros或NaN初始化循環(huán)只負(fù)責(zé)填值不負(fù)責(zé)擴(kuò)容。這個(gè)優(yōu)化簡單到不需要?jiǎng)幽X但在實(shí)際工作中效果極其明顯。凡是遇到循環(huán)里不斷往數(shù)組里追加數(shù)據(jù)的我第一反應(yīng)就是找循環(huán)外面能不能預(yù)分配。4.2 算法優(yōu)化向量化與“分而治之”很多人以為性能優(yōu)化就是寫更快的代碼其實(shí)真正重要的是選更快的算法。同樣的功能用不同的數(shù)據(jù)結(jié)構(gòu)和算法性能差距可以超過100倍。舉一個(gè)常見場景判斷一個(gè)元素是否屬于集合。用ismember解決很快但如果你在一個(gè)大循環(huán)里反復(fù)調(diào)用ismember耗時(shí)可能還是會(huì)超出預(yù)期。更好的方式是先把集合轉(zhuǎn)成containers.Map查詢是常數(shù)時(shí)間keys {alpha, beta, gamma, delta}; map containers.Map(keys, 1:length(keys)); % 在循環(huán)里快速查詢 for i 1:numel(testList) if isKey(map, testList{i}) % 命中邏輯 end end這個(gè)例子的核心思想是用哈希結(jié)構(gòu)替代線性查找把查詢復(fù)雜度從O(n)降到O(1)。很多MATLAB性能問題追根溯源都是算法復(fù)雜度的問題而不是代碼寫法的問題。當(dāng)你覺得“MATLAB好慢”的時(shí)候先問自己算法是不是有更優(yōu)的選擇4.3 內(nèi)存連續(xù)性與列優(yōu)先順序MATLAB存儲(chǔ)矩陣是列優(yōu)先的。這意味著遍歷一個(gè)矩陣時(shí)按列訪問比按行訪問要快得多。很多從C/C轉(zhuǎn)過來的朋友容易忽略這一點(diǎn)。看一個(gè)經(jīng)典對(duì)比A randn(5000, 5000); % 按列求和快速 tic; for j 1:size(A, 2) colSum sum(A(:, j)); end toc; % 按行求和慢速 tic; for i 1:size(A, 1) rowSum sum(A(i, :)); end toc;兩者邏輯相似但按列訪問比按行訪問快得多。大數(shù)據(jù)量下甚至可能差好幾倍。所以寫循環(huán)遍歷矩陣時(shí)盡量讓循環(huán)的外層對(duì)應(yīng)行、內(nèi)層對(duì)應(yīng)列也就是“先列后行”的順序這樣內(nèi)層訪問的是連續(xù)內(nèi)存。此外在數(shù)據(jù)導(dǎo)入時(shí)盡量讓最終使用的矩陣是連續(xù)的數(shù)值矩陣避免頻繁使用cell數(shù)組存取數(shù)值數(shù)據(jù)。cell數(shù)組提供了靈活性但它的存儲(chǔ)開銷和訪問速度都不如普通數(shù)值數(shù)組。如果在性能敏感的代碼里發(fā)現(xiàn)cell數(shù)組是瓶頸可以嘗試重構(gòu)為數(shù)值矩陣加索引的方式。4.4 大數(shù)據(jù)場景下的注意事項(xiàng)當(dāng)數(shù)據(jù)量大到內(nèi)存吃緊時(shí)普通的優(yōu)化思路就不夠了。這時(shí)有幾個(gè)原則需要特別留意避免產(chǎn)生大型中間變量。比如A B * C; D A E;A作為中間變量占了一次完整的內(nèi)存拷貝。如果B*C的結(jié)果只是用來和E相加可以嘗試約減運(yùn)算符。MATLAB對(duì)表達(dá)式的中間結(jié)果做了不少優(yōu)化但你自己也要注意保持表達(dá)式簡潔。用single精度代替double精度。很多科學(xué)計(jì)算場景其實(shí)用單精度就足夠了。內(nèi)存占用直接減半帶寬壓力也減半。但要注意單精度的數(shù)值誤差會(huì)增大涉及敏感運(yùn)算時(shí)務(wù)必驗(yàn)證。分塊處理代替全量加載。當(dāng)數(shù)據(jù)文件有幾十GB時(shí)不可能一次性讀進(jìn)內(nèi)存。用datastore或自定義分塊讀取邏輯每次只處理一個(gè)塊處理完即釋放。這是工程級(jí)的“流式處理”思路對(duì)MATLAB同樣適用。5. 把代碼變成工具面向?qū)ο笈cApp設(shè)計(jì)5.1 面向?qū)ο蟪绦蛟O(shè)計(jì)的實(shí)用價(jià)值很多人一聽說“面向?qū)ο蟆本陀X得是其他語言的專利MATLAB只是“數(shù)學(xué)計(jì)算器”。其實(shí)MATLAB的classdef機(jī)制相當(dāng)成熟而且在復(fù)雜項(xiàng)目中價(jià)值非常大。什么場景適合用面向?qū)ο笪业呐袛鄻?biāo)準(zhǔn)是當(dāng)你有多個(gè)變量需要捆綁在一起處理并且有一組操作總是圍繞它們展開時(shí)。比如你處理的是一個(gè)采集信號(hào)它包含數(shù)據(jù)、采樣率、時(shí)間戳、采集通道名稱、單位等屬性同時(shí)需要對(duì)它做濾波、去基線、切片、畫圖等操作。用結(jié)構(gòu)體固然可以但每次都要手動(dòng)帶上參數(shù)而且無法保證“數(shù)據(jù)”和“采樣率”來自同一個(gè)信號(hào)——類可以把這一切封裝起來。classdef SignalData properties Data SampleRate Timestamp ChannelNames end methods function obj SignalData(data, fs) obj.Data data; obj.SampleRate fs; obj.Timestamp datetime(now); end function filtered removeBaseline(obj, windowSize) filtered obj.Data - movmean(obj.Data, windowSize); end function plotTimeDomain(obj) t (0:length(obj.Data)-1) / obj.SampleRate; plot(t, obj.Data); xlabel(時(shí)間/s); ylabel(幅度); end end end使用的時(shí)候sig SignalData(workData, 1000); % 創(chuàng)建對(duì)象 cleanData removeBaseline(sig, 50); % 調(diào)用方法自動(dòng)帶上對(duì)象屬性 sig.plotTimeDomain(); % 直接畫圖這樣的代碼不僅整潔而且降低了調(diào)用方的認(rèn)知負(fù)擔(dān)——不需要把一堆分散的參數(shù)到處傳遞。對(duì)于項(xiàng)目越來越大、協(xié)作人數(shù)越來越多的情況類定義本身就是一份活文檔。5.2 什么時(shí)候不要過度設(shè)計(jì)但我也遇到過不少人把MATLAB寫成了Java為每個(gè)小功能建一個(gè)類接口抽象了三層。這是另一條彎路。MATLAB的核心優(yōu)勢是快速開發(fā)、快速驗(yàn)證過度設(shè)計(jì)會(huì)讓簡單問題復(fù)雜化。我個(gè)人的原則是腳本階段允許“臟亂快”一旦代碼要整個(gè)項(xiàng)目復(fù)用、要交給別人維護(hù)才重構(gòu)為函數(shù)或類。這個(gè)臨界點(diǎn)需要你自己根據(jù)項(xiàng)目規(guī)模判斷?？傊嫦?qū)ο笫菫榻鉀Q協(xié)作和擴(kuò)展問題的不是為了炫耀技巧。5.3 用App Designer把程序變成產(chǎn)品如果做了工具給別人用命令行界面不夠友好可以試試MATLAB自帶的App Designer。它相當(dāng)于一個(gè)可視化的界面設(shè)計(jì)器可以拖拽按鈕、文本框、坐標(biāo)區(qū)、儀表盤等控件然后為每個(gè)控件編寫回調(diào)函數(shù)。我之前做過一個(gè)數(shù)據(jù)查看工具用了不到半天時(shí)間就把一個(gè)命令行腳本包裝成了帶界面的小程序左邊是文件列表中間是曲線圖右邊是參數(shù)面板底部是日志輸出。App Designer在生成代碼時(shí)自動(dòng)產(chǎn)生一個(gè)類定義基于matlab.apps.AppBase你在回調(diào)里訪問控件屬性即可例如% 按鈕回調(diào)示例 properties (Access private) Data end methods (Access private) function onLoadButtonPushed(app, event) [file, path] uigetfile(*.csv); app.Data readtable(fullfile(path, file)); plot(app.UIAxes, app.Data.Time, app.Data.Value); end end把App打包成.exe可執(zhí)行文件配合MATLAB Compiler你的算法、界面、邏輯全部可以交付給不裝MATLAB的用戶使用。這一步完成之后“程序”就變成了真正的“工具”。6. 避坑實(shí)錄我從無數(shù)次實(shí)驗(yàn)里撿出來的經(jīng)驗(yàn)6.1 常見問題速查表我把多年來遇到的高頻問題整理成了一張表格幾乎每個(gè)項(xiàng)目都會(huì)用到問題場景典型報(bào)錯(cuò)解決思路矩陣維度不匹配Matrix dimensions must agree檢查size注意行向量與列向量的區(qū)別必要時(shí)用reshape或轉(zhuǎn)置索引越界Index exceeds the number of array elements檢查循環(huán)邊界考慮使用numel或size動(dòng)態(tài)獲取變量被覆蓋工作區(qū)變量與函數(shù)名沖突避免用函數(shù)名作變量名比如不要命名變量sum或mean字符串比較失敗if str abc總是不成立用strcmp不要用比較字符串文件路徑錯(cuò)誤File not found使用fullfile拼接路徑注意相對(duì)路徑與當(dāng)前工作目錄不一致循環(huán)慢到崩潰等了好久不出結(jié)果查profile檢查是否有預(yù)分配思考能否向量化程序報(bào)錯(cuò)但不知錯(cuò)哪堆棧深錯(cuò)誤信息晦澀用斷點(diǎn)調(diào)試配合dbstack定位出錯(cuò)行浮點(diǎn)比較失敗if (a/b) 0.333不成立浮點(diǎn)數(shù)用abs(x-y) tol方式比較不要用第1條和第8條是我見過最多的兩個(gè)錯(cuò)誤。矩陣維度問題大多源于對(duì)數(shù)據(jù)形狀不夠敏感浮點(diǎn)比較問題則是對(duì)計(jì)算機(jī)浮點(diǎn)表示缺乏認(rèn)識(shí)。理解這兩點(diǎn)能省掉大量排錯(cuò)時(shí)間。6.2 新手最容易踩的五個(gè)坑用i或j做循環(huán)變量然后又把復(fù)數(shù)單位弄丟了。MATLAB里i和j默認(rèn)是虛數(shù)單位。一旦你寫for i 1:10循環(huán)結(jié)束后i就被改成了數(shù)值10。如果后面某個(gè)表達(dá)式需要用到虛數(shù)單位你就得到莫名其妙的結(jié)果。建議循環(huán)變量用idx、k這類名稱徹底避開沖突。腳本里直接硬編碼文件路徑。今天在自己機(jī)器上能跑明天換臺(tái)電腦就崩潰。把輸入路徑放在腳本開頭集中管理甚至用input交互式獲取才是工程級(jí)的做法。調(diào)試代碼忘刪。disp、keyboard、plot臨時(shí)驗(yàn)證代碼留在函數(shù)里導(dǎo)致輸出亂成一團(tuán)。我建議臨時(shí)驗(yàn)證代碼統(tǒng)一用一個(gè)變量控制開關(guān)比如DEBUG true;發(fā)布前改成false。忽視矩陣方向。讀入數(shù)據(jù)時(shí)行是樣本、列是特征還是反過來的這個(gè)約定如果前后不統(tǒng)一程序跑出來全是錯(cuò)的還特別難查。建議在數(shù)據(jù)讀取后立刻驗(yàn)證size并寫注釋。循環(huán)里動(dòng)態(tài)命名變量。比如eval([var, num2str(i), data])。這類代碼讀起來極其痛苦調(diào)試更是噩夢。應(yīng)該用cell數(shù)組或結(jié)構(gòu)體比如var{i} data通過索引訪問優(yōu)雅又安全。6.3 我自己的“代碼體檢”清單每次交付一段比較重要的代碼之前我會(huì)從頭到尾過一遍這個(gè)清單所有函數(shù)都有明確輸入輸出注釋說明了參數(shù)含義和期望的數(shù)據(jù)類型。所有循環(huán)都有預(yù)分配或者有足夠的理由不預(yù)分配。所有文件訪問都通過fullfile拼接路徑且考慮了文件不存在的情況。所有比較浮點(diǎn)數(shù)的地方都用容差方式而不是。所有可能在批處理中出錯(cuò)的函數(shù)都包了try-catch且給出可讀的警告信息。所有臨時(shí)驗(yàn)證代碼都被移除或關(guān)閉。腳本和函數(shù)的工作區(qū)清晰沒有殘留無用變量。這套清單并不復(fù)雜但每次按它走一遍代碼質(zhì)量都能上一個(gè)臺(tái)階。很多“莫名奇妙跑不通”的問題其實(shí)早在寫代碼的階段就可以通過習(xí)慣避免。6.4 一個(gè)真實(shí)的排查案例最后分享一個(gè)我印象深刻的排查經(jīng)歷。某次處理一批實(shí)測振動(dòng)數(shù)據(jù)程序跑得挺快但結(jié)果曲線總有毛刺。一開始我懷疑是傳感器本身的問題后來把數(shù)據(jù)導(dǎo)入到一個(gè)簡單的繪制腳本里發(fā)現(xiàn)原始數(shù)據(jù)就沒有毛刺說明問題出在我自己的處理流程里。于是我從頭逐步執(zhí)行每一個(gè)處理步驟對(duì)比每一步前后的數(shù)據(jù)形狀和范圍。最后發(fā)現(xiàn)是在去噪步驟里我寫的窗口參數(shù)直接用了默認(rèn)值movmean(data, 5)但這個(gè)默認(rèn)值對(duì)高頻采樣的數(shù)據(jù)來說窗口太小完全沒起到去噪作用。把窗口大小改成根據(jù)采樣率自適應(yīng)計(jì)算后問題迎刃而解。這個(gè)案例給我最大的啟發(fā)是程序出問題時(shí)先確認(rèn)數(shù)據(jù)流再懷疑算法最后才懷疑硬件和工具。多數(shù)時(shí)候問題都在自己的代碼邏輯里用斷點(diǎn)加對(duì)比的方法能快速定位。寫MATLAB這些年我最深的體會(huì)是這個(gè)工具真正強(qiáng)大的地方不在于它會(huì)多少內(nèi)置函數(shù)而在于你能不能把一段只服務(wù)于當(dāng)前任務(wù)的“草稿代碼”錘煉成結(jié)構(gòu)清晰、性能達(dá)標(biāo)、邊界健壯的“工程代碼”。第六篇講的這些內(nèi)容不會(huì)讓你一次就寫出完美程序但至少下次再面對(duì)一個(gè)復(fù)雜任務(wù)時(shí)你會(huì)知道從哪里下手組織代碼、用什么手段定位問題、在哪一步去追求性能。照著這些思路去寫代碼自己會(huì)變得越來越干凈。