簡介面向電子電路仿真與數(shù)據(jù)分析場景LTspice2Matlab工具包專門解決LTspice與Matlab之間的數(shù)據(jù)導(dǎo)入問題適用于開關(guān)電源設(shè)計、模擬電路分析等需要將仿真波形轉(zhuǎn)換為可編程處理數(shù)據(jù)的工程師與科研人員。資源包共153個文件主要包含raw仿真原始數(shù)據(jù)、asc電路原理圖、plt繪圖配置、m腳本及md說明文檔整體大小僅2.37MB結(jié)構(gòu)緊湊便于查閱。其中m腳本提供封裝好的導(dǎo)入函數(shù)配合直流掃描、交流掃描、瞬態(tài)掃描等多個示例電路可在Matlab中直接運(yùn)行并完成波形繪制、快速傅里葉變換、濾波器設(shè)計等后續(xù)任務(wù)省去手動解析文本數(shù)據(jù)的繁瑣過程。借助這套工具還能將仿真結(jié)果與實際硬件測試數(shù)據(jù)對比輔助電路模型校準(zhǔn)與性能預(yù)測。已有112人學(xué)習(xí)下載適合具備一定仿真基礎(chǔ)的中高級用戶能有效縮短從仿真到數(shù)據(jù)分析的鏈路提升電路驗證與優(yōu)化效率。 做電路仿真的人十有八九都遇到過同一個尷尬LTspice波形窗口里看得好好的數(shù)據(jù)想拿進(jìn)Matlab做FFT、參數(shù)辨識或者和實測數(shù)據(jù)對比時卻不知道怎么把它干凈地導(dǎo)出來。LTspice本身有導(dǎo)出文本的功能但每次手動導(dǎo)出、整理格式非常費(fèi)勁仿真跑幾十組參數(shù)的話光導(dǎo)出數(shù)據(jù)就想罵人。今天想聊的就是專門解決這個痛點的工具LTspice2Matlab。它在MathWorks官網(wǎng)文件交換區(qū)就能找到作用只有一個——把LTspice生成的.raw二進(jìn)制數(shù)據(jù)文件直接讀進(jìn)Matlab工作區(qū)返回時間向量和變量矩陣幾行代碼就能完成數(shù)據(jù)對接。適合經(jīng)常用LTspice做模擬電路、電源、濾波器仿真又需要在Matlab里做后續(xù)分析的人。1. 為什么在Matlab里看LTspice數(shù)據(jù)——打破仿真與后處理的隔閡1.1 我為什么被逼到找LTspice2Matlab一個真實場景先講個我自己的經(jīng)歷。之前做一個電源環(huán)路補(bǔ)償設(shè)計LTspice里搭好了Buck電路跑完交流掃描波特圖在LTspice窗口里看得很漂亮。但問題是我要把仿真結(jié)果和實驗室里實測的環(huán)路增益曲線放在同一張圖里對比還得算出相位裕度、穿越頻率這些指標(biāo)。LTspice的波形查看器能看但不能做復(fù)雜的數(shù)值提取更不方便和測試數(shù)據(jù)做歸一化對比。那時候我試過幾種土辦法截屏后手動描點、用Export data as text導(dǎo)成CSV再拖進(jìn)Matlab、甚至想過把波形圖放大后用鼠標(biāo)一個一個讀數(shù)。說實話效率低到離譜。一個波特圖幾百個點手動描點描到懷疑人生而且一旦改了電路參數(shù)就得重新來一遍。LTspice2Matlab就是在這種情況下進(jìn)入視野的它直接把.raw文件里的所有仿真數(shù)據(jù)整體搬運(yùn)到Matlab工作區(qū)變量名、時間軸、數(shù)值全給你整整齊齊地擺好后面你想怎么處理怎么處理。1.2 為什么不用手動導(dǎo)出CSV對比一下效率有人可能會問LTspice本身就能導(dǎo)出文本數(shù)據(jù)我把它存成CSV再導(dǎo)入Matlab不就行了嗎確實行但只適合一次性、小規(guī)模的數(shù)據(jù)。我把兩種方式放在一起對比維度手動導(dǎo)出CSV再導(dǎo)入LTspice2Matlab直讀單次操作耗時1-2分鐘左右?guī)酌攵辔募幚砻總€文件都要手動操作腳本循環(huán)自動處理變量信息需要手動記順序自動保留變量名列表重復(fù)修改參數(shù)后重新走一遍全部流程重新調(diào)用一次函數(shù)即可數(shù)據(jù)完整性可能漏掉部分變量.raw內(nèi)全部數(shù)據(jù)都在尤其在批量掃參的場景下dir(*.raw)加上循環(huán)就能把整個文件夾的仿真結(jié)果一次性導(dǎo)入手動導(dǎo)出根本做不到這種自動化程度。這也是我后來一直用LTspice2Matlab的原因——它把“仿真完成”到“數(shù)據(jù)可用”這段路直接鋪平了。2. 搭建環(huán)境LTspice、Matlab和LTspice2Matlab的安裝與配置2.1 LTspice側(cè)確保.raw文件是二進(jìn)制格式LTspice2Matlab讀取的是LTspice仿真后生成的.raw文件這個文件默認(rèn)是二進(jìn)制格式里面保存了所有仿真變量在每個時間點的數(shù)值。LTspice仿真完成后.raw文件會自動出現(xiàn)在原理圖所在目錄下文件名與原理圖一致后綴為.raw。這里有一個很關(guān)鍵的設(shè)置項。LTspice的控制面板里有一項輸出格式選項如果你不小心勾選了“保存為ASCII文本”之類的選項生成的.raw文件就不是二進(jìn)制LTspice2Matlab讀取時會直接報錯。具體路徑在Control Panel - Operation務(wù)必保證二進(jìn)制raw格式處于開啟狀態(tài)。我見過好幾個同事導(dǎo)入失敗最后發(fā)現(xiàn)都是因為之前在別的地方調(diào)過這個選項。另外要注意的是LTspice有多個版本LTspice IV和LTspice XVII的raw格式在歷史上有些微差異。如果你用的是很老的LTspice IV建議去文件交換區(qū)下載最新版的LTspice2Matlab腳本新版本對老格式兼容性更好。2.2 在Matlab里接入LTspice2Matlab下載、放路徑、調(diào)用LTspice2Matlab本質(zhì)上是一個.m函數(shù)文件不涉及編譯安裝部署非常簡單。第一步在MathWorks官網(wǎng)的File Exchange搜索“LTspice2Matlab”下載LTspice2Matlab.m文件。如果你用公司內(nèi)網(wǎng)訪問不了也可以找同事要一份現(xiàn)成的腳本這個文件是開源的傳播很普遍。第二步把LTspice2Matlab.m放到Matlab可以找到的位置。我習(xí)慣建一個專門放自用腳本的目錄比如D:\MyMatlabTools然后在Matlab里執(zhí)行addpath(D:\MyMatlabTools);為了避免每次啟動Matlab都要重新添加路徑可以在預(yù)設(shè)項 - MATLAB - 常規(guī) - 自定義路徑里把這個目錄加進(jìn)去這樣以后啟動就自動加載了。第三步驗證是否安裝成功。隨便跑一個LTspice仿真到仿真目錄下執(zhí)行data LTspice2Matlab(your_sim.raw);如果不報錯說明環(huán)境已經(jīng)通了。這一步建議先拿最簡單的電路驗證不要一上來就處理大型掃描數(shù)據(jù)否則出了問題很難判斷是環(huán)境問題還是數(shù)據(jù)問題。3. 實操第一個完整的LTspice-to-Matlab導(dǎo)入任務(wù)3.1 先跑一個最簡單的RC電路紙上談兵沒意思直接拿一個RC電路完整走一遍流程。在LTspice里搭一個最簡單的低通濾波器一個電壓源V1一個電阻R1一個電容C1節(jié)點命名分別是IN和OUT。給V1設(shè)置一個階躍源或者脈沖源仿真類型選暫態(tài)分析.tran 5m也就是仿真5毫秒。仿真跑完后在原理圖同目錄下會生成一個RC_test.raw文件。現(xiàn)在打開Matlab把當(dāng)前工作目錄切到原理圖所在文件夾然后執(zhí)行data LTspice2Matlab(RC_test.raw);第一次看到data這個結(jié)構(gòu)體的時候你可能會覺得有點亂但其實字段很清晰。我用得最多的三個字段是data.time_vec時間向量也就是橫軸數(shù)據(jù)。注意這個時間向量不一定是等間隔的LTspice在波形變化劇烈時會自動加密步長。data.variable_name_list變量名列表按順序?qū)?yīng)每一列數(shù)據(jù)里面能看到V(in)、V(out)、I(R1)這些名字。data.variable_mat實際的數(shù)值矩陣每一列對應(yīng)一個變量。3.2 核心命令與數(shù)據(jù)結(jié)構(gòu)解析來看一段完整的讀取代碼% 讀取LTspice仿真數(shù)據(jù) data LTspice2Matlab(RC_test.raw); % 提取時間軸 t data.time_vec; % 查看變量列表確認(rèn)輸入輸出電壓對應(yīng)的列 disp(data.variable_name_list); % 提取輸出電壓數(shù)據(jù) vout data.variable_mat(:, 2); % 假設(shè)第2列是V(out) % 繪制波形 plot(t * 1e3, vout, LineWidth, 1.5); xlabel(Time (ms)); ylabel(Voltage (V)); title(RC Step Response); grid on;這里有個容易踩坑的地方variable_mat的列順序和variable_name_list里的順序是對應(yīng)的但每次仿真變量順序可能不一樣。我吃過一次虧直接用第2列當(dāng)輸出電壓結(jié)果那次仿真的第2列恰好是輸入電壓圖畫出來完全不對。所以正式處理之前一定先打印variable_name_list確認(rèn)一下或者用代碼自動查找目標(biāo)變量的索引varIdx find(strcmp(data.variable_name_list, V(out))); vout data.variable_mat(:, varIdx);3.3 導(dǎo)入之后讓數(shù)據(jù)真正“能用”畫圖、FFT、提取指標(biāo)數(shù)據(jù)導(dǎo)入Matlab只是第一步大部分人的真實需求是在導(dǎo)入后做定量分析。最常見的一個場景是做FFT頻譜分析但這里有個細(xì)節(jié)很多人沒有注意到LTspice暫態(tài)仿真的時間步長是自適應(yīng)變化的不是均勻采樣直接在原始時間軸上做FFT會得到不準(zhǔn)確的結(jié)果。穩(wěn)妥的做法是先做重采樣把數(shù)據(jù)插值到等間隔的時間網(wǎng)格上再做FFT% 重采樣到均勻時間軸fs1MHz fs 1e6; t_uniform (0 : 1/fs : t(end)); vout_uniform interp1(t, vout, t_uniform, linear); % 做FFT Y fft(vout_uniform); f (0 : length(Y)-1) * fs / length(Y); plot(f(1:length(f)/2), abs(Y(1:length(Y)/2))); xlabel(Frequency (Hz)); ylabel(Magnitude);除了FFT另一個高頻需求是從時域波形里提取指標(biāo)比如上升時間、過沖、穩(wěn)定時間。這些在LTspice里也能測但只能通過光標(biāo)一點點讀比較麻煩。導(dǎo)入Matlab后可以用幾行代碼自動算出來尤其適合批量處理不同參數(shù)下的大量波形。4. 數(shù)據(jù)導(dǎo)入之后幾個高頻場景的封裝方法4.1 批量處理多個.raw文件仿真不可能只跑一次。我在做參數(shù)優(yōu)化的時候經(jīng)常一個晚上跑二三十個不同電容、電阻參數(shù)的仿真。如果每個文件都手動讀取效率太低。這時候直接用dir配合循環(huán)處理files dir(*.raw); for k 1:length(files) data LTspice2Matlab(files(k).name); % 根據(jù)文件名提取參數(shù) params extractParams(files(k).name); % 提取輸出變量并計算指標(biāo) idx find(strcmp(data.variable_name_list, V(out))); vout data.variable_mat(:, idx); % 計算并存儲 result(k).name files(k).name; result(k).overshoot calcOvershoot(vout); end批量仿真還有一個更省事的方式用LTspice的命令行模式。在命令行直接執(zhí)行PATH_TO_LTSPICE\XVIIx86.exe -b netlist.ascLTspice會靜默仿真并在后臺生成.raw文件配合Matlab的system命令可以做到全自動化。我現(xiàn)在的流程基本是Matlab腳本自動修改原理圖參數(shù)→調(diào)用LTspice命令行仿真→用LTspice2Matlab讀回數(shù)據(jù)→自動分析并畫圖。4.2 處理.step參數(shù)掃描的數(shù)據(jù)LTspice的.step指令可以一次仿真跑多個參數(shù)比如.step param R 1k 10k 1k會跑10組電阻。這種仿真生成的.raw文件里包含多組數(shù)據(jù)處理起來相對麻煩。我的個人建議是如果后續(xù)要批量導(dǎo)入Matlab做定量分析盡量少用.step而是把不同的參數(shù)組合拆成單獨(dú)的原理圖文件用批量仿真方式生成多個.raw文件。這樣每個.raw對應(yīng)一組參數(shù)數(shù)據(jù)結(jié)構(gòu)清晰不用在Matlab里再去拆分。.step適合在LTspice里快速掃參看趨勢但如果要精確提取每組參數(shù)對應(yīng)的指標(biāo)單個文件導(dǎo)入會更省心。如果你確實需要處理.step生成的raw文件LTspice2Matlab也能讀但返回的數(shù)據(jù)結(jié)構(gòu)會包含所有step的數(shù)據(jù)需要根據(jù)變量名或者數(shù)據(jù)長度自行切分。我之前試過幾次感覺切分邏輯有點繞后來還是改成了批量生成獨(dú)立文件的方式效率反而更高。4.3 對齊時間軸非均勻采樣的處理LTspice的自適應(yīng)步長對仿真精度很有幫助但對后期數(shù)據(jù)分析卻是個障礙。特別是當(dāng)你需要把多個不同仿真條件的結(jié)果放在同一張圖里對比或者要把不同仿真的時間軸對齊時必須統(tǒng)一采樣網(wǎng)格。我的標(biāo)準(zhǔn)做法是所有從LTspice導(dǎo)入的數(shù)據(jù)統(tǒng)一用interp1插值到相同的時間軸。比如所有對比實驗都重采樣到0.1微秒間隔這樣后面的計算和繪圖都不用再去處理時間軸的不一致性。% 定義統(tǒng)一時間軸 t_target (0 : 1e-7 : 5e-3); vout_resampled interp1(t, vout, t_target, pchip);這里我用pchip而不是linear因為三次Hermite插值在波形轉(zhuǎn)折處的處理更平滑不會像線性插值那樣引入明顯的棱角。不過如果你的波形本身很平滑兩者差別不大。5. 常見問題與排查技巧實錄5.1 五類高頻報錯速查表用LTspice2Matlab的過程中我從新手到熟練踩了不少坑總結(jié)成下面這張速查表報錯現(xiàn)象原因解決辦法找不到文件當(dāng)前目錄不對或文件名拼寫錯誤用cd切換目錄檢查后綴是否為.rawInvalid file identifier.raw文件被保存成ASCII文本格式在LTspice控制面板改回二進(jìn)制格式重新仿真變量列表為空仿真沒有保存目標(biāo)變量在LTspice中使用.save指令指定保存變量讀取結(jié)果明顯不對變量列順序和預(yù)期不一致先用variable_name_list確認(rèn)列順序中文路徑報錯Matlab在中文路徑下讀文件異常把原理圖和.raw文件放到英文路徑第二條值得多說一句。LTspice的“ASCII raw”選項是個老坑一旦勾上后續(xù)所有仿真都會輸出ASCII格式文件體積大且LTspice2Matlab讀不了。如果你發(fā)現(xiàn)raw文件比平時大很多先懷疑這個設(shè)置。5.2 避坑小貼士文件占用、命名規(guī)范、變量順序除了報錯還有一些比較隱蔽的細(xì)節(jié)第一個是數(shù)據(jù)文件的占用問題。LTspice仿真結(jié)束后如果波形窗口還開著或者正在查看數(shù)據(jù).raw文件可能被程序鎖定。這時候Matlab去讀取會報錯。保險起見等LTspice完全釋放文件后再執(zhí)行導(dǎo)入。我自己通常先關(guān)掉波形窗口再切到Matlab執(zhí)行讀取。第二個是命名規(guī)范。.raw文件的文件名最好不要包含中文、空格、特殊符號。LTspice本身支持中文路徑但Matlab讀取時很容易因為編碼問題報錯。把工程文件統(tǒng)一用英文命名是省心省力的好習(xí)慣。第三個是變量順序。這一點前面提過但值得反復(fù)強(qiáng)調(diào)永遠(yuǎn)不要靠記憶猜變量列數(shù)。每換一個電路或原理圖變量順序很可能就變了。用strcmp自動查索引是必須寫進(jìn)腳本的習(xí)慣。如果你要做長期的數(shù)據(jù)分析建議在LTspice2Matlab外面再包一層自己的封裝函數(shù)。比如我寫了一個loadLTspiceData(fileName, signalName)內(nèi)部自動調(diào)LTspice2Matlab然后根據(jù)信號名自動找到對應(yīng)的列返回時間軸和信號數(shù)據(jù)。這樣主腳本里不需要每次都處理變量順序的邏輯代碼可讀性也高。用了一兩年LTspice2Matlab之后我最深的體會是真正難的不是導(dǎo)入這個動作而是想清楚數(shù)據(jù)導(dǎo)入之后要怎么組織。如果你只是偶爾看一次波形LTspice自帶窗口完全夠用一旦涉及批量掃參、與實測數(shù)據(jù)對比、自動化分析這套流程會讓你省下大量時間。先用最簡單的RC電路把流程跑通確認(rèn)環(huán)境沒問題再把復(fù)雜度逐步加進(jìn)去這個工具就能變成你仿真流程里一個很順手的環(huán)節(jié)。本文還有配套的精品資源點擊獲取