簡(jiǎn)介這份MATLAB計(jì)算程序面向信號(hào)處理、ADC/DAC測(cè)試與電路設(shè)計(jì)方向的學(xué)習(xí)者和工程師用于評(píng)估數(shù)字信號(hào)質(zhì)量與系統(tǒng)性能。程序圍繞SNR、SNDR、THD、ENOB、SFDR五項(xiàng)核心指標(biāo)展開通過讀取數(shù)據(jù)文件、預(yù)處理、FFT頻譜分析分離信號(hào)、噪聲與失真分量并計(jì)算各參數(shù)功率比幫助讀者理解指標(biāo)背后的理論含義與實(shí)現(xiàn)方式。資源包共2個(gè)文件包含1個(gè)m腳本與1個(gè)txt數(shù)據(jù)文件壓縮包約4KB腳本負(fù)責(zé)計(jì)算流程數(shù)據(jù)文件提供待分析樣本便于直接運(yùn)行驗(yàn)證。目前已有4557人學(xué)習(xí)下載適合需要快速搭建指標(biāo)計(jì)算框架、對(duì)照理論公式排查結(jié)果偏差的讀者參考也可作為信號(hào)鏈性能評(píng)估與優(yōu)化的小型工具腳本使用。1. 五個(gè)指標(biāo)算錯(cuò)一個(gè)整條信號(hào)鏈的評(píng)估就翻車了做數(shù)據(jù)采集或者 ADC 評(píng)估的人幾乎都遇到過這種場(chǎng)景辛辛苦苦采了一段正弦波丟進(jìn) FFT想看看 SNR 到底多少結(jié)果頻譜上底噪忽高忽低算出來的 ENOB 一會(huì)兒 10.2 位一會(huì)兒 11.8 位完全沒法寫進(jìn)報(bào)告。問題往往不在硬件而在計(jì)算程序本身——窗函數(shù)選錯(cuò)、相干采樣沒做、直流和高次諧波沒剔除、歸一化基準(zhǔn)搞混任何一個(gè)環(huán)節(jié)出問題SNR、SNDR、THD、ENOB、SFDR 這五個(gè)指標(biāo)就會(huì)集體失真。這套 MATLAB 計(jì)算程序要解決的正是把「從一段時(shí)域采樣到五個(gè)頻域指標(biāo)」的整條鏈路固化下來讓每一次評(píng)估都可復(fù)現(xiàn)、可對(duì)比。搞懂它背后的機(jī)制比會(huì)調(diào)用幾個(gè)函數(shù)重要得多因?yàn)橹笜?biāo)之間的耦合關(guān)系才是判斷一個(gè)轉(zhuǎn)換器真實(shí)性能的關(guān)鍵。2. 從時(shí)域采樣到頻域指標(biāo)FFT 之后到底發(fā)生了什么2.1 為什么必須做相干采樣泄漏是怎么毀掉 SNR 的FFT 的本質(zhì)是把有限長(zhǎng)序列當(dāng)成一個(gè)周期信號(hào)的完整周期來展開。如果采樣點(diǎn)數(shù)和信號(hào)頻率不滿足整數(shù)周期關(guān)系首尾就接不上頻譜上會(huì)出現(xiàn)能量泄漏——主瓣能量向兩側(cè)擴(kuò)散原本干凈的底噪被抬高SNR 直接虛低。這就是為什么相干采樣是第一步讓信號(hào)頻率、采樣率、采樣點(diǎn)數(shù)三者滿足fin M/N * fsM 取與 N 互質(zhì)的整數(shù)。常見的實(shí)現(xiàn)思路是先確定采樣率 fs 和點(diǎn)數(shù) N一般取 2 的冪比如 8192、16384、65536再反推一個(gè)滿足相干條件的輸入頻率。程序里通常會(huì)有這樣一段N 2^13; % FFT 點(diǎn)數(shù)取 2 的冪便于基 2 運(yùn)算 fs 100e6; % 采樣率 100 MSPS M 31; % 與 N 互質(zhì)的整數(shù)決定相干頻率 fin M/N * fs; % 相干采樣頻率保證整數(shù)周期這里 M 必須與 N 互質(zhì)否則信號(hào)會(huì)在多個(gè)周期內(nèi)重復(fù)等效于降低了有效采樣點(diǎn)數(shù)。選 31 這種質(zhì)數(shù)是最省事的做法。如果實(shí)在做不到相干采樣就必須加窗但加窗會(huì)引入主瓣展寬和旁瓣衰減后續(xù)算 SNR 時(shí)要把窗函數(shù)的等效噪聲帶寬ENBW補(bǔ)償回去否則又是一筆糊涂賬。2.2 加窗與 ENBW 補(bǔ)償非相干采樣下的后悔藥工程現(xiàn)場(chǎng)很多時(shí)候沒法精確控制輸入頻率尤其是評(píng)估板上已經(jīng)焊死的信號(hào)源。這時(shí)候只能加窗。窗函數(shù)的作用是讓截?cái)嗵幍耐蛔冏兤交瑝褐婆园?。但代價(jià)是主瓣變寬頻譜上每個(gè) bin 代表的噪聲帶寬不再是 1而是 ENBW。以 Hann 窗為例ENBW 約為 1.5 個(gè) bin。也就是說你把所有噪聲 bin 的功率加起來之后要除以 ENBW 才是真實(shí)的噪聲功率。這一步漏掉SNR 會(huì)偏高約 1.76 dB。程序里一般這樣處理win hanning(N); % 生成 Hann 窗 xw x .* win; % 加窗 X fft(xw)/sum(win); % 歸一化補(bǔ)償窗增益 P abs(X).^2; % 功率譜 enbw sum(win.^2)/sum(win)^2 * N; % 等效噪聲帶寬bin注意fft之后除以sum(win)而不是 N這是為了補(bǔ)償窗函數(shù)帶來的幅度衰減。歸一化基準(zhǔn)搞錯(cuò)是所有指標(biāo)集體偏移的常見原因。加窗是權(quán)宜之計(jì)相干采樣才是正解這一點(diǎn)在寫評(píng)估報(bào)告時(shí)心里要有數(shù)。2.3 直流、諧波、噪聲的分離邏輯拿到功率譜之后五個(gè)指標(biāo)的計(jì)算本質(zhì)上是對(duì)不同 bin 集合做功率求和。直流分量bin 0 附近要剔除因?yàn)樗淮斫涣餍阅苄盘?hào)主瓣fin 對(duì)應(yīng)的 bin 及其相鄰泄漏 bin要單獨(dú)拎出來作為信號(hào)功率諧波位置2fin、3fin…的 bin 要?dú)w到諧波功率剩下的才是噪聲。這里有個(gè)容易踩的坑信號(hào)主瓣到底占幾個(gè) bin。相干采樣下理想情況只占 1 個(gè) bin但實(shí)際會(huì)有輕微泄漏通常取主瓣中心左右各 1 個(gè) bin 共 3 個(gè) bin 作為信號(hào)功率。諧波 bin 的選取同理要按諧波頻率對(duì)應(yīng)的 bin 索引去定位而不是簡(jiǎn)單按倍數(shù)關(guān)系取整。程序里一般會(huì)寫一個(gè)輔助函數(shù)來定位某個(gè)頻率對(duì)應(yīng)的 binfunction k freq2bin(f, fs, N) k round(f/fs*N) 1; % MATLAB 索引從 1 開始 end這個(gè)1是 MATLAB 數(shù)組下標(biāo)從 1 開始的經(jīng)典坑漏掉它整個(gè)頻譜會(huì)偏移一個(gè) bin諧波定位全錯(cuò)。血淚經(jīng)驗(yàn)每次換語言移植這段邏輯第一件事就是檢查索引基準(zhǔn)。3. 五個(gè)指標(biāo)的計(jì)算機(jī)制誰是誰的分子分母3.1 SNR 與 SNDR一字之差差在諧波算不算噪聲SNR 的定義是信號(hào)功率與噪聲功率之比這里的噪聲不包括諧波。SNDR也叫 SINAD則是信號(hào)功率與「噪聲 諧波」總功率之比。兩者差多少直接反映諧波失真有多嚴(yán)重。計(jì)算邏輯上先把信號(hào)功率 Ps 求出來再把除直流和信號(hào)外的所有 bin 功率求和得到總噪聲諧波功率 Pn_total把諧波 bin 功率求和得到 Ph。那么SNR 10*log10(Ps / (Pn_total - Ph))SNDR 10*log10(Ps / Pn_total)Ps sum(P(sig_bins)); % 信號(hào)功率 Ph sum(P(harm_bins)); % 諧波功率 Pn_total sum(P) - Ps - P(dc_bins); % 噪聲諧波 SNR 10*log10(Ps / (Pn_total - Ph)); SNDR 10*log10(Ps / Pn_total);很多人第一次算會(huì)發(fā)現(xiàn) SNR 比 SNDR 高出一大截誤以為程序有 bug其實(shí)這正是諧波貢獻(xiàn)的差距。如果兩者幾乎相等說明系統(tǒng)諧波失真很低是好事。反過來如果 SNR 和 SNDR 都很低那問題在底噪不在諧波。3.2 THD諧波功率到底要取到第幾次THD 是總諧波功率與信號(hào)功率之比通常用百分比或 dB 表示。關(guān)鍵問題是取到第幾次諧波。理論上越高越好但實(shí)際受限于 Nyquist 帶寬——高于 fs/2 的諧波會(huì)折疊回來混進(jìn)噪聲里取多了反而把噪聲當(dāng)諧波。常見的做法是取前 5 到 10 次諧波且只取那些落在 Nyquist 以內(nèi)的。程序里一般會(huì)設(shè)一個(gè)nHarm參數(shù)nHarm 5; % 取前 5 次諧波 harm_bins []; for h 2:nHarm fh h * fin; if fh fs/2 % 只取 Nyquist 以內(nèi) harm_bins [harm_bins, freq2bin(fh, fs, N)]; end end THD 10*log10(sum(P(harm_bins)) / Ps);注意諧波 bin 定位后通常也要像信號(hào)主瓣一樣取左右相鄰 bin防止泄漏導(dǎo)致諧波功率被低估。取幾次諧波沒有絕對(duì)標(biāo)準(zhǔn)但同一份報(bào)告里必須一致否則橫向?qū)Ρ葲]有意義。3.3 ENOB 與 SFDR一個(gè)看整體一個(gè)看最壞點(diǎn)ENOB有效位數(shù)是把 SNDR 換算成等效的理想 ADC 位數(shù)公式是ENOB (SNDR - 1.76) / 6.02。這個(gè) 1.76 和 6.02 來自理想量化噪聲模型滿量程正弦輸入下量化信噪比約為 6.02N 1.76 dB。所以 ENOB 本質(zhì)上是 SNDR 的另一種表達(dá)單位從 dB 換成了 bit。SFDR無雜散動(dòng)態(tài)范圍則是信號(hào)功率與最大雜散不一定是諧波也可能是時(shí)鐘饋通、交調(diào)產(chǎn)物功率之比。它看的是頻譜上最壞的那個(gè)點(diǎn)而不是整體噪聲。ENOB (SNDR - 1.76) / 6.02; P_spur max(P(setdiff(1:N/2, [dc_bins, sig_bins]))); % 最大雜散 SFDR 10*log10(Ps / P_spur);ENOB 和 SFDR 經(jīng)常一起看ENOB 高但 SFDR 低說明底噪干凈但有個(gè)別雜散ENOB 低但 SFDR 高說明整體噪聲大但沒有突出雜散。這兩個(gè)指標(biāo)的組合比單看任何一個(gè)都更能定位問題。4. 避坑與排查五個(gè)指標(biāo)算不對(duì)時(shí)先查這幾處4.1 現(xiàn)象SNR 算出來是負(fù)數(shù)或個(gè)位數(shù)原因通常是歸一化基準(zhǔn)錯(cuò)了。比如fft之后沒有除以點(diǎn)數(shù)或窗增益導(dǎo)致信號(hào)功率和噪聲功率的絕對(duì)數(shù)值都失真比值自然離譜。另一種可能是信號(hào) bin 定位偏了把噪聲當(dāng)成了信號(hào)。解決先畫出功率譜肉眼確認(rèn)信號(hào)主瓣位置是否正確。再檢查fft后的歸一化系數(shù)相干采樣下除以 N加窗時(shí)除以sum(win)。把 Ps 和 Pn 的絕對(duì)值打印出來和理論值對(duì)一下。4.2 現(xiàn)象ENOB 超過 ADC 標(biāo)稱位數(shù)原因一般是噪聲 bin 求和時(shí)把直流或信號(hào)泄漏 bin 也排除了導(dǎo)致噪聲功率偏小SNDR 虛高。也可能是諧波沒算進(jìn) SNDR誤用了 SNR 去算 ENOB。解決確認(rèn) SNDR 用的是「噪聲 諧波」總功率而不是 SNR。檢查直流剔除范圍是否過寬把信號(hào)主瓣附近的 bin 也誤刪了。ENOB 超過標(biāo)稱位數(shù)在物理上不可能一定是計(jì)算口徑問題。4.3 現(xiàn)象THD 每次算出來都不一樣原因是非相干采樣導(dǎo)致泄漏隨機(jī)諧波 bin 定位不穩(wěn)定?；蛘咧C波次數(shù)取得太多把折疊回來的高頻噪聲也算進(jìn)去了。解決優(yōu)先保證相干采樣。做不到就加窗并做 ENBW 補(bǔ)償。諧波次數(shù)固定下來且只取 Nyquist 以內(nèi)的。每次評(píng)估用同樣的 N、fs、nHarm保證可復(fù)現(xiàn)。4.4 現(xiàn)象SFDR 比 SNR 還低原因通常是最大雜散定位時(shí)把信號(hào)主瓣的泄漏 bin 也算進(jìn)去了導(dǎo)致「雜散」功率虛高。也可能是直流附近的低頻雜散沒排除。解決雜散搜索范圍要排除直流 bin、信號(hào)主瓣 bin 及其相鄰泄漏 bin。如果確實(shí)存在靠近信號(hào)的雜散那 SFDR 低是真實(shí)的需要從時(shí)鐘或電源上找原因。4.5 現(xiàn)象換一段數(shù)據(jù)五個(gè)指標(biāo)集體偏移原因多半是采樣率或點(diǎn)數(shù)變了但程序里的相干頻率、諧波次數(shù)沒跟著調(diào)整。或者窗函數(shù)換了但 ENBW 沒重算。解決把 N、fs、fin、nHarm、窗類型作為一組參數(shù)統(tǒng)一管理換數(shù)據(jù)時(shí)整組更新。程序里加一個(gè)自檢確認(rèn)fin/fs*N接近整數(shù)否則給出警告。5. 設(shè)計(jì)取舍為什么這套程序不追求「全自動(dòng)」寫這類計(jì)算程序最容易掉進(jìn)的坑是追求「一鍵出結(jié)果」——自動(dòng)找信號(hào)、自動(dòng)定諧波、自動(dòng)判雜散。但實(shí)際評(píng)估中信號(hào)頻率、諧波次數(shù)、雜散范圍這些參數(shù)往往需要根據(jù)被測(cè)對(duì)象和測(cè)試條件人為判斷。全自動(dòng)的代價(jià)是黑匣子化一旦結(jié)果異常你根本不知道它把哪些 bin 算進(jìn)了噪聲。這套程序的取舍是把參數(shù)暴露出來把中間量功率譜、信號(hào) bin、諧波 bin、噪聲 bin都保留讓使用者能自己核對(duì)。代價(jià)是上手要多填幾個(gè)參數(shù)收益是每一個(gè)數(shù)字都能追溯到具體的 bin。我一般會(huì)這么理解評(píng)估程序的價(jià)值不在于算得快而在于算得可解釋。取舍點(diǎn)全自動(dòng)方案本程序方案信號(hào)定位峰值搜索按相干頻率計(jì)算 bin諧波次數(shù)自適應(yīng)手動(dòng)指定 nHarm雜散范圍全譜搜索排除直流與信號(hào)主瓣可復(fù)現(xiàn)性依賴算法穩(wěn)定性依賴參數(shù)一致性排查難度高黑匣子低中間量可見如果你要驗(yàn)證自己的實(shí)現(xiàn)對(duì)不對(duì)一個(gè)實(shí)用的方法是用理想正弦波加已知噪聲生成一段數(shù)據(jù)理論 SNR 可以手算然后看程序輸出是否吻合。再人為加入二次諧波看 THD 是否按預(yù)期變化。這種「已知答案反推」的驗(yàn)證比拿真實(shí)數(shù)據(jù)反復(fù)調(diào)參靠譜得多。從那以后我每次拿到一段新的采樣數(shù)據(jù)都強(qiáng)制先畫功率譜、確認(rèn)相干條件、核對(duì)信號(hào) bin 位置再去看那五個(gè)指標(biāo)。希望幫到你。本文還有配套的精品資源點(diǎn)擊獲取