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