簡(jiǎn)介這是一套Matlab信號(hào)處理工具箱離線資源面向已經(jīng)安裝Matlab但缺少信號(hào)處理組件、以及希望單獨(dú)補(bǔ)全該工具箱的工程師、教師和科研人員。內(nèi)容覆蓋濾波器設(shè)計(jì)、頻譜分析、信號(hào)重采樣、Levinson遞推、序列周期估計(jì)等常見任務(wù)既能服務(wù)于課堂實(shí)驗(yàn)和畢業(yè)設(shè)計(jì)也能支撐工程項(xiàng)目中的快速原型驗(yàn)證。壓縮包共405個(gè)文件體量約2.51MB擴(kuò)展名類型較集中316個(gè)m腳本便于閱讀和二次修改69個(gè)p文件用于算法封裝12個(gè)dll動(dòng)態(tài)庫(kù)負(fù)責(zé)底層加速另有少量mat數(shù)據(jù)、頭文件和說明文檔輔助調(diào)用與結(jié)果驗(yàn)證。已有7995人瀏覽學(xué)習(xí)說明該集合經(jīng)過較多同行驗(yàn)證具備一定可用性。下載后把整個(gè)目錄加入Matlab路徑即可直接調(diào)用相關(guān)功能函數(shù)省去自行查找匹配版本的繁瑣過程同時(shí)也能從m腳本中學(xué)習(xí)實(shí)現(xiàn)思路對(duì)初、中級(jí)Matlab使用者尤為友好。 每次拿到傳感器采集回來的原始數(shù)據(jù)第一件事絕對(duì)不應(yīng)該是扔進(jìn)模型里訓(xùn)練也不是直接畫個(gè)波形圖看一眼就完事。數(shù)據(jù)在進(jìn)業(yè)務(wù)邏輯之前必須經(jīng)過信號(hào)處理這一關(guān)濾波、去趨勢(shì)、頻譜分析、特征提取一套流程走完你才有底氣說“這個(gè)信號(hào)是干凈的”。這份signalprocessingtoolbox信號(hào)處理工具箱就是干這個(gè)用的它圍繞信號(hào)處理最常見的幾類操作把生成測(cè)試信號(hào)、濾波器設(shè)計(jì)、FFT頻譜分析、批量特征提取這些腳本整理成了一個(gè)可以直接拿來用的工具包。做嵌入式、搞硬件調(diào)試、寫算法驗(yàn)證的同學(xué)還有正在做振動(dòng)分析或傳感器數(shù)據(jù)預(yù)處理的人都可以直接拿它當(dāng)?shù)讓右蕾囀〉裘看维F(xiàn)翻文檔寫函數(shù)的重復(fù)勞動(dòng)。1. 項(xiàng)目定位與整體設(shè)計(jì)思路1.1 這個(gè)工具箱到底裝了什么很多人第一次看到signalprocessingtoolbox會(huì)以為它指的是MATLAB官方那個(gè)Signal Processing Toolbox。其實(shí)這個(gè)壓縮包里裝的是一套按我自己的使用習(xí)慣整理出來的腳本集合相當(dāng)于把官方工具箱里常用的函數(shù)、自己寫的封裝腳本、還有若干示例數(shù)據(jù)文件打包在一起。這么做的好處是明顯的官方工具箱功能全但調(diào)用鏈長(zhǎng)每次做一個(gè)濾波分析要翻好幾個(gè)函數(shù)頁(yè)而自己的工具包可以把“讀數(shù)據(jù)—濾波—畫頻譜—提特征”串成一條流水線一行命令出結(jié)果。這個(gè)包里的核心內(nèi)容按功能分成這么幾塊信號(hào)生成器生成正弦、方波、白噪聲、調(diào)幅信號(hào)等測(cè)試源濾波器集合低通、高通、帶通、帶阻附加自動(dòng)截止頻率計(jì)算頻譜分析模塊FFT、功率譜密度、語(yǔ)譜圖還有一組針對(duì)工程現(xiàn)場(chǎng)的數(shù)據(jù)導(dǎo)入函數(shù)支持CSV、TXT和常見的二進(jìn)制采樣文件。整個(gè)包沒有做復(fù)雜的GUI界面全部是函數(shù)和腳本因?yàn)閷?shí)際干活的時(shí)候命令行方式反而更靈活尤其是在批處理幾十個(gè)文件的時(shí)候GUI會(huì)把人點(diǎn)瘋。1.2 為什么把散落的腳本整理成工具箱這個(gè)習(xí)慣是從一次項(xiàng)目事故之后養(yǎng)成的。當(dāng)時(shí)做一個(gè)電機(jī)振動(dòng)監(jiān)測(cè)的預(yù)研數(shù)據(jù)采了一整天晚上準(zhǔn)備分析的時(shí)候發(fā)現(xiàn)下午寫的濾波腳本和上午用的函數(shù)居然不兼容參數(shù)格式對(duì)不上得重新改。從那之后我就堅(jiān)持把所有信號(hào)處理相關(guān)的函數(shù)統(tǒng)一收口到一個(gè)工具箱目錄下統(tǒng)一命名、統(tǒng)一接口、統(tǒng)一數(shù)據(jù)單位。另一個(gè)原因是可復(fù)現(xiàn)性。工程上一個(gè)結(jié)論要被認(rèn)可必須能反復(fù)跑出同樣的結(jié)果。如果每次分析都臨時(shí)寫一段腳本中間換個(gè)參數(shù)、換個(gè)文件結(jié)果可能就飄了。工具包把處理流程固定下來輸入文件路徑和幾個(gè)關(guān)鍵參數(shù)輸出就是標(biāo)準(zhǔn)化的頻譜圖和特征表出問題也好溯源。這也是這個(gè)工具箱設(shè)計(jì)上最核心的思路把重復(fù)勞動(dòng)收斂把可變參數(shù)顯性化。1.3 適用人群誰(shuí)最需要這份資源如果你屬于下面任何一類人這份工具箱的整理思路都值得參考做嵌入式開發(fā)的同學(xué)手里有ADC采集數(shù)據(jù)或者通過CAN總線拿到的傳感器原始值需要先做濾波再送給控制算法搞硬件調(diào)試的工程師手上有示波器導(dǎo)出的采樣數(shù)據(jù)想快速看看頻率成分做算法驗(yàn)證的研究生需要對(duì)比不同濾波器的效果或者提取信號(hào)特征作為機(jī)器學(xué)習(xí)輸入剛接觸信號(hào)處理的學(xué)生想通過一套完整的示例理解采樣定理、濾波器設(shè)計(jì)、FFT這些概念到底怎么落地這個(gè)工具箱的價(jià)值不在于代碼本身有多高深而在于它把“分析一段信號(hào)的標(biāo)準(zhǔn)動(dòng)作”沉淀成了可以復(fù)用的資產(chǎn)而且邊界清楚接入新項(xiàng)目的時(shí)候不用傷筋動(dòng)骨。2. 核心功能拆解信號(hào)處理工具箱里最值錢的幾個(gè)模塊2.1 信號(hào)生成與導(dǎo)入先有一份干凈的測(cè)試數(shù)據(jù)信號(hào)處理這條鏈路上第一步往往是“沒有數(shù)據(jù)”。很多初學(xué)者直接拿實(shí)際采集的信號(hào)就來分析噪聲成分未知、干擾未知分析結(jié)果出了問題都不知道該怪誰(shuí)。正確的做法是先構(gòu)造已知的測(cè)試信號(hào)比如一個(gè)50Hz的正弦疊加300Hz的干擾和白噪聲用這個(gè)信號(hào)驗(yàn)證你的濾波器參數(shù)對(duì)不對(duì)、頻譜分析準(zhǔn)不準(zhǔn)再切換到真實(shí)數(shù)據(jù)。工具箱里的信號(hào)生成器封裝了一個(gè)很實(shí)用的函數(shù)function [t, x] siggen(fs, dur, components) % SIGGEN 生成復(fù)合測(cè)試信號(hào) % fs 采樣率單位Hz % dur 時(shí)長(zhǎng)單位秒 % components 元胞數(shù)組如 {50, 1; 300, 0.5} 表示50Hz幅值1 300Hz幅值0.5 t (0:1/fs:dur-1/fs).; x zeros(size(t)); for k 1:size(components, 1) f components{k, 1}; a components{k, 2}; x x a * sin(2*pi*f*t); end x x 0.1 * randn(size(t)); % 加一點(diǎn)白噪聲模擬真實(shí)環(huán)境 end這段代碼的邏輯很直白但有一個(gè)細(xì)節(jié)值得注意白噪聲幅值系數(shù)0.1是刻意加進(jìn)去的。如果你噪聲加得太大濾波器效果好不好一眼看不出來加得太小又起不到驗(yàn)證魯棒性的作用。實(shí)際工程里信噪比設(shè)置多少要看目標(biāo)場(chǎng)景工具箱里默認(rèn)讓測(cè)試信號(hào)比噪聲高20dB左右這個(gè)經(jīng)驗(yàn)值大部分情況下都?jí)蛴?。?shù)據(jù)導(dǎo)入這塊我踩過不少坑。示波器導(dǎo)出的CSV文件經(jīng)常是“第一列時(shí)間第二列通道1第三列通道2”注釋行還特別多直接用csvread往往會(huì)報(bào)錯(cuò)。工具箱里統(tǒng)一用readtable加參數(shù)配置的方式讀取通過函數(shù)參數(shù)指定表頭行數(shù)、時(shí)間列、數(shù)據(jù)列兼容性好了很多。2.2 濾波器設(shè)計(jì)從頻譜里撈出你要的頻率成分濾波器是整個(gè)工具箱里含金量最高的一塊。很多人一上來就用butter隨便寫個(gè)截止頻率結(jié)果濾波后信號(hào)相位亂了或者高頻干擾沒濾干凈。設(shè)計(jì)濾波器本質(zhì)上是在做三件事定類型、定階數(shù)、定截止頻率。濾波器的四個(gè)基本類型對(duì)應(yīng)不同的工程場(chǎng)景低通用于去除高頻噪聲高通用于去除直流漂移和低頻趨勢(shì)項(xiàng)帶通用于只保留某個(gè)頻段比如振動(dòng)分析里的1倍頻、2倍頻成分帶阻用于剔除特定干擾比如50Hz工頻噪聲。選擇哪種類型取決于你的信號(hào)特征和噪聲分布。截止頻率的設(shè)定是新手最容易出問題的地方。工具箱里實(shí)現(xiàn)了一個(gè)輔助函數(shù)根據(jù)采樣率自動(dòng)計(jì)算歸一化截止頻率function Wn normcutoff(fc, fs, ftype) % 根據(jù)截止頻率和采樣率計(jì)算歸一化頻率 % ftype: low, high, band if strcmp(ftype, band) Wn fc / (fs/2); % 這里fc是兩個(gè)元素 [fc1 fc2] else Wn fc / (fs/2); end Wn max(Wn, 0); Wn min(Wn, 0.9999); % 防止越界 end這里必須解釋一個(gè)基礎(chǔ)概念但凡是涉及數(shù)字濾波器的設(shè)計(jì)所有頻率參數(shù)都必須相對(duì)于奈奎斯特頻率采樣率的一半做歸一化。如果采樣率是1000Hz想濾掉100Hz以上的成分歸一化截止頻率就是100 / (1000/2) 0.2。這個(gè)計(jì)算一步都不能省不然濾波器設(shè)計(jì)函數(shù)會(huì)給你一個(gè)完全錯(cuò)誤的系數(shù)。工具箱里默認(rèn)使用butter設(shè)計(jì)巴特沃斯濾波器因?yàn)樗谕◣?nèi)最平坦沒有紋波適合大多數(shù)信號(hào)調(diào)理場(chǎng)景。如果對(duì)相位有嚴(yán)格要求比如后續(xù)要做波形相關(guān)性分析就得改用filtfilt做零相位濾波。filtfilt的原理是把信號(hào)正著濾一遍再反著濾一遍相位延遲互相抵消代價(jià)是計(jì)算量翻倍但換來的零相位偏移在工程上常常是值得的。2.3 頻譜分析與特征提取別只畫圖要能讀到信息FFT大概是信號(hào)處理里被用爛又用錯(cuò)最多的工具。很多人把波形數(shù)據(jù)扔進(jìn)fft畫出幅值譜就完事了。實(shí)際上FFT的正確用法牽扯到三個(gè)關(guān)鍵參數(shù)采樣點(diǎn)數(shù)N、采樣率fs、頻率分辨率df fs / N。頻率分辨率是什么意思就是頻譜圖上相鄰兩根譜線之間的距離。如果你采樣1秒fs是1000Hz那么N 1000df 1Hz頻譜上每根譜線代表1Hz的間隔。如果你的有用信號(hào)是50.5Hz而頻譜分辨率只有1Hz那么50.5Hz的成分會(huì)泄漏到50Hz和51Hz兩根譜線上幅值還會(huì)被拉低。解決兩個(gè)辦法增加采樣時(shí)長(zhǎng)提高分辨率或者做零填充提高譜線密度。工具箱里的頻譜分析函數(shù)默認(rèn)會(huì)做零填充到2的冪次同時(shí)計(jì)算并返回幅值、頻率和對(duì)應(yīng)的功率譜密度f(wàn)unction [f, mag] spec_analyze(x, fs) % 計(jì)算單邊幅值譜 N length(x); X fft(x, 2^nextpow2(N)); % 補(bǔ)零到2的冪次 X X(1:floor(length(X)/2)1); mag abs(X) * 2 / N; % 單邊譜幅值還原 f (0:floor(length(X)/2)) * fs / length(X); end注意這里有個(gè)細(xì)節(jié)單邊譜的幅值必須乘以2再除以N因?yàn)镕FT的結(jié)果是雙邊對(duì)稱的單邊只取一半能量要按比例還原回去。如果忘記乘2你的信號(hào)幅值會(huì)顯示成實(shí)際值的一半這種低級(jí)錯(cuò)誤在工程報(bào)告里出現(xiàn)得相當(dāng)頻繁。用這個(gè)函數(shù)可以順手提取頻譜峰值位置、能量集中頻段、以及諧波畸變率這些特征值可以作為后續(xù)機(jī)器學(xué)習(xí)模型的輸入也可以用于設(shè)備健康狀態(tài)的判斷。3. 實(shí)操過程拿一段真實(shí)信號(hào)走一遍全流程3.1 環(huán)境準(zhǔn)備與數(shù)據(jù)導(dǎo)入動(dòng)手之前先把工具包路徑加進(jìn)來。我這邊的習(xí)慣是在項(xiàng)目啟動(dòng)腳本里統(tǒng)一配置addpath(genpath(signalprocessingtoolbox));這個(gè)命令會(huì)把工具箱目錄下所有子目錄都添加到MATLAB搜索路徑里包括函數(shù)和示例數(shù)據(jù)。如果你在Octave里跑同樣支持只是個(gè)別圖形窗口函數(shù)寫法略有差異。接下來模擬一段現(xiàn)場(chǎng)數(shù)據(jù)。假設(shè)我們有一個(gè)轉(zhuǎn)速傳感器的輸出信號(hào)采樣率2000Hz采集時(shí)長(zhǎng)2秒信號(hào)主體是80Hz的旋轉(zhuǎn)頻率但疊加了40Hz的電源諧波干擾和隨機(jī)振動(dòng)噪聲。先調(diào)用信號(hào)生成器構(gòu)造測(cè)試信號(hào)fs 2000; dur 2; components {80, 1.0; 240, 0.4; 40, 0.6}; [t, x] siggen(fs, dur, components);如果用的是現(xiàn)場(chǎng)采集的真實(shí)數(shù)據(jù)把x替換成readtable讀進(jìn)來的列就行。比如從CSV讀入data readtable(sensor_data.csv, VariableNamingRule, preserve); x data.Ch1; % 假設(shè)通道1是我們要分析的信號(hào) fs 2000; % 注意CSV里通常沒有采樣率得從采集配置里拿這里必須提醒一個(gè)現(xiàn)場(chǎng)高頻踩坑點(diǎn)CSV導(dǎo)出的時(shí)間列經(jīng)常是字符串格式別直接用str2num轉(zhuǎn)完就完事。很多采集軟件用的是“相對(duì)時(shí)間”而非絕對(duì)時(shí)間戳?xí)r間列根本不可靠直接用采樣率生成時(shí)間軸更穩(wěn)。如果一定要用時(shí)間列算采樣率建議先diff一下時(shí)間向量看看采樣間隔是否均勻不均勻的先插值重采樣。3.2 濾波器參數(shù)的設(shè)計(jì)與計(jì)算拿到信號(hào)后先畫個(gè)原始波形和粗略頻譜看噪聲集中在哪里。以我們構(gòu)造的信號(hào)為例有用成分在80Hz干擾在40Hz和240Hz目標(biāo)是把240Hz這個(gè)高頻干擾和40Hz這個(gè)低頻干擾都處理掉同時(shí)保留80Hz。最直接的做法是設(shè)計(jì)一個(gè)帶通濾波器通帶設(shè)在60到120Hzfc_low 60; % 高通截止頻率 fc_high 120; % 低通截止頻率 Wn [fc_low/(fs/2), fc_high/(fs/2)]; [b, a] butter(4, Wn, bandpass); y filtfilt(b, a, x);這里為什么選4階階數(shù)越高過渡帶越窄濾波效果越陡峭但帶來的副作用是相位延遲更大、數(shù)值穩(wěn)定性變差。對(duì)于大部分傳感器信號(hào)處理2到4階巴特沃斯已經(jīng)完全夠用。4階的意思實(shí)際是內(nèi)部等效8階因?yàn)閎utter返回的b和a是二階節(jié)級(jí)聯(lián)的系數(shù)filtfilt在零相位模式下等效階數(shù)還要翻倍效果上已經(jīng)是16階的滾降特性再高容易出現(xiàn)數(shù)值問題。濾波完成后對(duì)比一下濾波前后的頻譜。你會(huì)看到40Hz和240Hz的成分被壓到接近底噪水平80Hz成分的幅值保持在原樣附近。有一點(diǎn)要特別注意filtfilt是零相位濾波信號(hào)在時(shí)域上不會(huì)偏移但如果改用傳統(tǒng)的filter濾波器引入的相位延遲會(huì)讓你的時(shí)間和信號(hào)對(duì)不上做時(shí)間對(duì)齊分析比如兩路信號(hào)的互相關(guān)時(shí)會(huì)直接導(dǎo)致錯(cuò)誤結(jié)論。3.3 頻譜分析與特征提取濾波之后進(jìn)行頻譜分析[f, mag] spec_analyze(y, fs); figure; plot(f, mag); xlabel(頻率 (Hz)); ylabel(幅值);運(yùn)行下來在80Hz處會(huì)看到一個(gè)明顯的單峰。如果用findpeaks提取峰值[pks, locs] findpeaks(mag, MinPeakHeight, 0.3, MinPeakDistance, 20); peak_freq f(locs);MinPeakDistance這個(gè)參數(shù)很關(guān)鍵單位不是Hz而是樣本點(diǎn)。我們的頻率分辨率是fs / N 2000 / 4096 ≈ 0.49Hz如果兩個(gè)峰值相距很近比如80Hz和85Hz對(duì)應(yīng)樣本點(diǎn)距離約10個(gè)點(diǎn)那MinPeakDistance設(shè)置成20就能把兩個(gè)峰分開。如果設(shè)置太小同一片諧波會(huì)被識(shí)別成一堆假峰設(shè)置太大真正的多峰會(huì)被漏掉。特征提取這步我最常用的幾個(gè)指標(biāo)是峰值頻率、峰值幅值、以及頻段能量占比。比如把80Hz附近±2Hz范圍內(nèi)的能量加起來除以整個(gè)通帶能量就能得到“旋轉(zhuǎn)頻率能量占比”這個(gè)特征。這個(gè)比值對(duì)轉(zhuǎn)速波動(dòng)和設(shè)備健康狀況非常敏感是設(shè)備故障診斷里一個(gè)實(shí)用的指標(biāo)。工具箱里順手加了一個(gè)批量處理腳本接收一個(gè)文件夾路徑循環(huán)讀取所有CSV自動(dòng)做同樣的濾波和特征提取最后匯總成一張表格files dir(data/*.csv); result table(); for k 1:length(files) data readtable(fullfile(files(k).folder, files(k).name)); x data.Ch1; y filtfilt(b, a, x); [freq, mag] spec_analyze(y, fs); [pks, locs] findpeaks(mag, MinPeakHeight, 0.3, MinPeakDistance, 20); result [result; table({files(k).name}, f(locs(1)), pks(1), VariableNames, {文件, 峰值頻率Hz, 峰值幅值})]; end writetable(result, feature_summary.csv);就這么一段腳本原來手動(dòng)處理一個(gè)文件大概5分鐘現(xiàn)在幾十個(gè)文件一口氣跑完效率提升是非常直觀的。這也是整個(gè)工具箱最值錢的地方——不是單個(gè)函數(shù)多高級(jí)而是把重復(fù)動(dòng)作壓縮成了一鍵操作。3.4 一鍵批量處理把單次分析變成可復(fù)用工具批量處理腳本最后如果直接收尾那還能更進(jìn)一步。實(shí)際工程里每次項(xiàng)目的采樣率不同、有用頻段不同總不能每次改腳本里的硬編碼參數(shù)。我的做法是抽一個(gè)主入口函數(shù)出來用結(jié)構(gòu)體傳參cfg.fs 2000; cfg.filter_type bandpass; cfg.fc [60, 120]; cfg.filter_order 4; cfg.min_peak_height 0.3; cfg.min_peak_distance 20; cfg.input_dir data/; cfg.output_file feature_summary.csv; batch_analyze(cfg);這樣換項(xiàng)目的時(shí)候只需要改配置結(jié)構(gòu)體函數(shù)體一行不用動(dòng)。配置項(xiàng)集中管理還有個(gè)隱藏好處報(bào)告里要寫清數(shù)據(jù)處理參數(shù)的時(shí)候直接把cfg打印出來就是一份完整的參數(shù)清單審閱的人一眼能看清你做了什么處理對(duì)可復(fù)現(xiàn)性要求高的場(chǎng)景非常有用。4. 常見問題與排查技巧實(shí)錄4.1 autosar can signal如何連接rte做嵌入式控制器的朋友經(jīng)常會(huì)遇到一個(gè)有點(diǎn)類似但更偏軟件架構(gòu)的問題AUTOSAR架構(gòu)下CAN信號(hào)怎么和RTE連接。這雖然不是純信號(hào)處理問題但和信號(hào)鏈路密切相關(guān)我在實(shí)際項(xiàng)目中踩過不少。核心流程是先在CAN通信矩陣DBC或ARXML文件里定義好Signal然后通過DaVinci Configurator或EB tresos這類工具生成RTE接口最后在SWC軟件組件的端口里把RTE Event和數(shù)據(jù)元素映射到對(duì)應(yīng)的CAN Signal上。這里面最容易卡住的是信號(hào)字節(jié)序和縮放因子對(duì)不上的問題。CAN Signal在DBC里定義的是motorola字節(jié)序而RTE生成代碼時(shí)默認(rèn)按intel字節(jié)序處理數(shù)據(jù)就會(huì)完全錯(cuò)亂。排查方法很簡(jiǎn)單抓一幀CAN報(bào)文對(duì)比工具生成的解析值如果出現(xiàn)明顯的數(shù)值量級(jí)異常優(yōu)先檢查字節(jié)序和縮放因子。另外還有一點(diǎn)RTE Event觸發(fā)方式要選對(duì)——數(shù)據(jù)是周期性更新還是事件觸發(fā)直接決定了RTE端口能不能及時(shí)拿到新數(shù)據(jù)選錯(cuò)的話信號(hào)會(huì)一直顯示舊值。4.2 signal tap ii報(bào)錯(cuò)invalid jtag configurationFPGA調(diào)試時(shí)用SignalTap II邏輯分析儀有段時(shí)間一連接就報(bào)invalid jtag configuration。這個(gè)問題的根因通常不在SignalTap本身而在JTAG鏈的配置上。尤其是多片F(xiàn)PGA或FPGA和CPU混接的板子上JTAG鏈上有多個(gè)設(shè)備Quartus如果沒有識(shí)別到正確的鏈結(jié)構(gòu)就會(huì)報(bào)這個(gè)錯(cuò)。排查思路很明確先打開Quartus的Hardware Manager看JTAG鏈掃描結(jié)果確認(rèn)鏈上有幾個(gè)設(shè)備、每個(gè)設(shè)備的IR長(zhǎng)度和IDCODE是否和預(yù)期一致。最常見的坑是板子上FPGA的JTAG引腳被復(fù)用成普通IO了或者JTAG鏈信號(hào)經(jīng)過的電平轉(zhuǎn)換芯片沒有正確使能。之前遇到一次報(bào)錯(cuò)就是因?yàn)镴TAG的TMS引腳被一個(gè)下拉電阻拉低BSDL鏈路上設(shè)備被跳過把電阻改成上拉就好了。另外SignalTap實(shí)例的采樣時(shí)鐘如果沒連接到全局時(shí)鐘網(wǎng)絡(luò)也會(huì)引起配置數(shù)據(jù)加載后無(wú)法啟動(dòng)采樣的詭異現(xiàn)場(chǎng)這個(gè)和JTAG配置錯(cuò)誤是兩個(gè)方向別混在一起排查。4.3 runtime error received signal 11 的排查思路這個(gè)報(bào)錯(cuò)常見于C/C環(huán)境received signal 11其實(shí)就是段錯(cuò)誤訪問了非法內(nèi)存地址。在信號(hào)處理類的程序里段錯(cuò)誤的高發(fā)原因和無(wú)符號(hào)整數(shù)索引跑到負(fù)值、數(shù)組越界、濾波器狀態(tài)緩沖區(qū)沒初始化這三類問題高度相關(guān)。我在項(xiàng)目里遇到過一次原因是FFT輸入緩沖區(qū)大小是動(dòng)態(tài)分配的但某次傳入的數(shù)據(jù)長(zhǎng)度不是2的冪次FFT庫(kù)內(nèi)部索引越界。報(bào)錯(cuò)信息只顯示signal 11找半天才發(fā)現(xiàn)是緩沖區(qū)長(zhǎng)度問題。如果早一步在入口處加個(gè)斷言assert(n 64 (n (n-1)) 0)問題當(dāng)場(chǎng)就能暴露。另外一個(gè)排查技巧是在Linux下用valgrind跑一遍它能直接告訴你非法訪問發(fā)生在哪一行代碼、訪問了哪個(gè)地址定位速度比逐行打斷點(diǎn)快得多。Windows下可以用Application Verifier配合VS調(diào)試器也能拿到類似的越界信息。這類錯(cuò)誤在數(shù)值計(jì)算和信號(hào)處理代碼里尤其隱蔽因?yàn)楦↑c(diǎn)數(shù)計(jì)算不報(bào)錯(cuò)數(shù)組越界可能不是立刻崩潰而是“碰巧”寫進(jìn)了一個(gè)暫時(shí)不影響結(jié)果的內(nèi)存位置等到后續(xù)邏輯用到那塊數(shù)據(jù)才爆雷。所以寫完濾波器或FFT相關(guān)代碼建議第一時(shí)間喂一遍邊界條件長(zhǎng)度0、長(zhǎng)度1、長(zhǎng)度65535、數(shù)據(jù)全0、數(shù)據(jù)全NaN能有效提前暴露問題。4.4 常見問題速查表把上面這些排查經(jīng)驗(yàn)整理成一張速查表方便現(xiàn)場(chǎng)直接對(duì)照問題現(xiàn)象可能原因排查方向CAN信號(hào)解析值不對(duì)字節(jié)序或縮放因子配置錯(cuò)誤檢查DBC/ARXML對(duì)比原始報(bào)文和解析值RTE端口始終拿不到新數(shù)據(jù)觸發(fā)方式選錯(cuò)或周期不匹配檢查RTE Event類型和發(fā)送周期SignalTap報(bào)invalid jtag configurationJTAG鏈設(shè)備識(shí)別問題掃描JTAG鏈檢查TMS上下拉和電平轉(zhuǎn)換頻譜圖上信號(hào)幅度減半FFT單邊譜忘記乘2檢查幅值還原公式單邊譜要乘2再除以N濾波后波形時(shí)間對(duì)不上用了filter而非filtfilt需要零相位偏移的時(shí)候必須用filtfilt程序報(bào)received signal 11數(shù)組越界或緩沖區(qū)未初始化上valgrind或Application Verifier定位頻率分辨率不夠采樣時(shí)長(zhǎng)太短或N太小加長(zhǎng)采樣時(shí)間或補(bǔ)零到更長(zhǎng)FFT點(diǎn)數(shù)濾波器過渡帶太緩濾波器階數(shù)過低適當(dāng)提高階數(shù)注意數(shù)值穩(wěn)定性最后再分享一個(gè)我自用的習(xí)慣每次拿到一批新數(shù)據(jù)我會(huì)先跑一遍不加任何濾波的原始頻譜分析把頻譜圖存成PNG歸檔然后再跑濾波后的版本。這樣不管后面參數(shù)怎么調(diào)整原始數(shù)據(jù)長(zhǎng)什么模樣始終有據(jù)可查也方便復(fù)盤“我到底通過濾波改變了什么”。信號(hào)處理這套東西很多時(shí)候問題不是出在算法不夠高級(jí)而是出在過程不可控、結(jié)果不可復(fù)現(xiàn)一個(gè)把流程固定下來的工具箱能幫你省下大量返工的精力。本文還有配套的精品資源點(diǎn)擊獲取