← 所有文章

FIELD NOTES / 感測器與量測

取樣頻率不足如何辨識混疊

感測器與量測作者:站長預估 6 分鐘閱讀

以fs=100、f=70計算30 Hz混疊,說明相位、Nyquist邊界與ADC前抗混疊濾波的驗收流程。

本文目錄

一 取樣率先決定能看見什麼

取樣不是把類比訊號任意切成數字,而是決定離散資料能否分辨原波形。若取樣頻率fs低於訊號最高頻率的兩倍,會發生混疊:不同連續訊號在取樣點上留下相同或近似數列,後端再怎麼運算也無法知道原本是哪一個。工程上要先列出可能的頻帶、取樣時鐘、時間戳與允許誤差,不能只看畫面上是否偶爾出現數值。

本篇用正弦訊號說明計算。設定fs=100 Hz、輸入f=70 Hz,取樣角度每點增加2π乘0.70;離散序列可由相位相配的30 Hz訊號表示,因為70=100-30;零相位70Hz正弦等於負的零相位30Hz正弦。這個30 Hz就是混疊頻率。若工程師把數列直接標成70 Hz,頻譜判讀與後續控制都會錯。正常結果應報告「在目前取樣率下觀測到30 Hz等效成分」,而不是宣稱已恢復真實70 Hz。

每筆資料保存sample_index、採樣時間、fs版本與濾波設定。不要用牆上時鐘差值反推精確取樣週期,先確認時鐘是否單調與是否有丟樣。測試可送入已知頻率並改變fs,觀察估計頻率何時折返;若已知輸入超過Nyquist邊界,應標示無法由慢速取樣唯一辨識;未知輸入不能單靠頻譜自動知道已混疊並停止把結果當成原訊號。

先問「想保留哪個最高頻率」再決定fs。若最高頻率是40 Hz,fs=100 Hz雖然滿足理想兩倍條件,仍可能沒有給濾波器過渡帶留下空間;設計常需更高取樣率或更低的通帶上限。需求、硬體能力與儲存成本要一起評估,不能把Nyquist條件當成充分的工程裕度。

二 混疊計算與相位符號

對輸入頻率f與取樣率fs,可先找出整數k,使f_alias等於絕對值f-kfs後落在0到fs/2。fs=100、f=70時取k=1,得到30 Hz;f=130也得到30 Hz。這個算式只說明離散序列的等效頻率,不代表類比端真的變成30 Hz。當訊號含多個頻率,每一條都可能折返,必須分別計算並標記來源。

相位也會改變符號意義。取樣70 Hz的正弦可寫成sin(2π70 n/100+φ),因為sin(2πn-2π30n/100+φ)等於sin(-2π30n/100+φ),它在離散時間呈現反向旋轉的30 Hz成分。若只比較振幅而忽略相位,可能把相序、方向或同步關係判斷反。報表應保留參考時間點與相位定義,不能只寫一個頻率數字。

f=fs/2是邊界案例。理想數學取樣在每半週期落點,某些相位會全落在零點,某些相位會呈交替正負;因此不能說邊界能完美重建原訊號。實際時鐘誤差、雜訊與相位漂移會讓結果更不穩。驗收時把fs/2附近分成低於、等於、高於三組測試,並記錄是否出現折返與不確定標記。

若資料已有抗混疊濾波,仍須記錄截止頻率、衰減、相位延遲與穩定時間。濾波器不能在ADC取樣後消除已經折返的成分;抗混疊濾波必須放在ADC前,並讓通帶、過渡帶與取樣率留出設計裕度。未知道路徑時,不要用軟體平均冒充ADC前濾波。

相位符號最好用實際序列做小表驗證:列出n=0、1、2的角度與數值,再和30 Hz正向、反向表示比較。這能抓出週期、角度單位與正負號錯誤。報告應說明頻率是連續時間輸入、離散等效頻率,還是頻譜峰值,三者名稱不能混用。

三 驗證頻帶與失敗案例

先建立輸入頻率表,例如fs=100 Hz時,10 Hz應保留在10 Hz,49 Hz接近上限仍可觀察,51 Hz會折返為49 Hz,70 Hz會折返為30 Hz。每列附上預期等效頻率、相位方向、輸入振幅與測量窗口。正常結果是低頻通過且與預期一致;失敗結果是把70 Hz報成70 Hz,或用短窗口把頻率解析度誤當成取樣率造成的能力。

頻率估計還受窗口長度影響。100 Hz取樣、觀察0.1秒只有10點,精確30Hz對應三個DFT週期,理想時落在頻格;若改33Hz就不是整數週期,頻譜洩漏可能使峰值分散。增加窗口能改善解析度,卻不能恢復被混疊的原頻率。排查時分開記錄取樣率、窗口長度、觸發相位與估計演算法;不要把改變窗口後的外觀改善說成消除混疊。

若要驗證相位,使用同一取樣時鐘同步的參考通道,保存兩通道sample_index與相位差定義。不同時鐘的時間戳若未校正,看到的相位漂移可能來自時鐘偏差。若資料曾經丟樣或重排,應標記序號缺口,不要插值後再宣稱原始相位可追溯。

控制或保護用途須先定義允許頻帶與最壞情境。混疊可能把高頻振動變成低頻假訊號,誘發錯誤濾波或錯誤判斷。這篇只提供取樣與資料設計方法,不指定任何PLC模組、寄存器或設備規格;實際ADC前端、硬體濾波與診斷能力要依目標產品文件確認。

若輸入頻率未知,先做頻帶盤點或用較高率的獨立通道觀察,再決定慢速通道能否安全使用。把未知訊號直接降採樣並只看結果,無法證明沒有高頻成分。遇到頻率隨時間變化,還要以最壞瞬間頻率和掃頻案例驗收。

四 設計流程與限制

流程第一步是列出輸入最高可能頻率,而不是只列正常運轉頻率;第二步選擇fs與抗混疊濾波器,使通帶和過渡帶有明確裕度;第三步在ADC前確認濾波器實際位置;第四步用已知10、49、51、70 Hz案例驗收。每一步都保存設定版本,避免日後只拿到一串無法解釋的數列。

若取樣頻率受通訊或記憶體限制,不能以降低fs換取容量而不改變頻帶要求。可改採分層:快速通道保存高頻診斷,慢速通道供畫面趨勢;兩者要有不同名稱與品質標記。下游看到30 Hz時,要能追查它來自原始頻帶還是混疊,不可混在同一個標籤中。

限制包括未知輸入頻帶、漂移的取樣時鐘、非同步多通道與濾波器延遲。即使數學計算正確,也不能由混疊資料反推出唯一原波形。若安全判斷依賴高頻訊號,應在設計階段增加硬體頻寬與獨立驗證,不能用事後軟體補救。

文件要把濾波器的群延遲列為資料的一部分。抗混疊濾波成功壓制帶外成分,卻可能讓事件時間向後移;控制策略若需要準確時間,須同時驗證振幅衰減與延遲,而不是只看頻譜乾淨。

五 常見問題與官方參考

FAQ1:fs=100 Hz能測70 Hz嗎?答:離散序列會與30 Hz等效,不能僅靠取樣資料辨識真實70 Hz。

FAQ2:f=fs/2是否代表最高可測頻率?答:它是理想邊界,重建受相位與濾波器影響,工程驗收應留裕度。

FAQ3:ADC後平均能消除混疊嗎?答:不能;混疊已進入離散資料,抗混疊濾波應位於ADC前。

FAQ4:增加FFT點數能找回原頻率嗎?答:不能;真正延長有效觀察時間可改善解析,單純補零只加密頻格,不會消除取樣造成的頻率折返。

參考:NI:頻寬、取樣與混疊的量測說明。

參考:Analog Devices MT-002:Nyquist取樣、混疊與ADC前抗混疊濾波的工程說明。

延伸閱讀