岳瑞華 徐中英 楊學(xué)猛

(第二炮兵工程大學(xué)自動控制工程系，陜西 西安 710025)

0 引言

數(shù)據(jù)采集器的采樣速率是指采集通道在單位時間內(nèi)所采集的數(shù)據(jù)個數(shù)，通道間延遲時間是指在不同通道間切換時所需時間。

數(shù)據(jù)采集器得到的是離散化數(shù)據(jù)，通過對離散化數(shù)據(jù)的處理可得到需要的測量結(jié)果。數(shù)據(jù)采集器的采樣速率和通道間延遲時間［1－2］直接決定離散數(shù)據(jù)和時域的關(guān)系。因此，對數(shù)據(jù)采集器的采樣速率和通道間延遲時間的評價是正確選擇和使用數(shù)據(jù)采集器的基本條件。

1 采樣速率的評價方法

一般評價系統(tǒng)使用標(biāo)準(zhǔn)周期計(jì)數(shù)法對采樣速率進(jìn)行評價［3－4］，但某些測量系統(tǒng)中配置的是標(biāo)準(zhǔn)信號源，如標(biāo)準(zhǔn)信號源Fluke5520A，標(biāo)準(zhǔn)信號源能夠提供標(biāo)準(zhǔn)正弦信號。因此，本文采用標(biāo)準(zhǔn)正弦波信號和自相關(guān)方法評價數(shù)據(jù)采集器的采樣速率。其基本思想為:標(biāo)準(zhǔn)信號源送出標(biāo)準(zhǔn)正弦信號，頻率為f=100 Hz;標(biāo)準(zhǔn)頻率信號1 年絕對不確定度為(2.5 ±5)×10－6Hz［5－6］，啟動數(shù)據(jù)采集器采集N個數(shù)據(jù)。利用自相關(guān)技術(shù)［7－10］消除采集數(shù)據(jù)中的噪聲和干擾，然后，尋找出N個數(shù)據(jù)中所包含的輸入信號整周期的個數(shù)nc和nc個整周期所包含的采集數(shù)據(jù)個數(shù)Nc，則采樣速率fs為:

式中:nc為輸入信號整周期的個數(shù);Nc為nc個整周期所包含的采集數(shù)據(jù)個數(shù);f為標(biāo)準(zhǔn)信號的頻率，Hz;KT為在一個信號周期中包含的采樣個數(shù)，也稱為信號的數(shù)字化周期。

由于標(biāo)準(zhǔn)信號的頻率f已知，由式(1)可知，評價工作的重點(diǎn)是KT的獲取。

式(1)與一般的標(biāo)準(zhǔn)周期計(jì)數(shù)法的計(jì)算公式相同，本文不再贅述。不同的是，一般的標(biāo)準(zhǔn)周期計(jì)數(shù)法直接對應(yīng)采集數(shù)據(jù)，式(1)卻對應(yīng)采集數(shù)據(jù)的自相關(guān)離散序列，原因如下。

原始的采集數(shù)據(jù)不可避免地包含噪聲和干擾，在自相關(guān)處理后，可以有效地抑制噪聲和干擾而不改變數(shù)據(jù)序列的周期。判斷自相關(guān)后的數(shù)據(jù)序列的可信度比直接判斷采集數(shù)據(jù)序列高。自相關(guān)分析原理如下。

所謂自相關(guān)，就是信號與自身進(jìn)行相關(guān)運(yùn)算得到其自相關(guān)函數(shù)，信號的自相關(guān)函數(shù)具有保持信號周期不變的特性。為了說明這一點(diǎn)，假設(shè)信號x(t)具有周期T，即:x(t)=x(t+nT)。x(t)的傅里葉級數(shù)表達(dá)式為:

這時信號x(t)的自相關(guān)函數(shù)Rxx(t)可用傅里葉級數(shù)表示為:

由式(3)可知，周期信號的自相關(guān)函數(shù)也是同一基頻的周期函數(shù)。Rxx(τ)包含x(t)的全部頻率分量，各頻率的振幅等于相應(yīng)的傅里葉系數(shù)的均方值，振幅變化依賴于振幅本身，且是非線性變化。但x(t)中不同頻率分量的相位在Rxx(τ)中歸零，采用過零點(diǎn)判斷信號數(shù)字化周期時，只關(guān)心過零點(diǎn)時數(shù)據(jù)符號的變化。因此，Rxx(τ)中不同頻率分量的相位歸零對于采用過零點(diǎn)判斷信號數(shù)字化周期是非常有利的。

在一般情況下，信號不僅由基波及其諧波組成，還不可避免地包含噪聲和干擾，在數(shù)字化過程中，還要加入量化噪聲等。因此，一般可假定記錄x(t)是由周期信號s(t)和隨機(jī)噪聲n(t)組成，即:x(t)=s(t)+n(t)。噪聲n(t)的均值為零，與s(t)不相關(guān)。x(t)的自相關(guān)函數(shù)為:

式中:Rss(τ)為周期信號的自相關(guān)函數(shù);Rnn(τ)為噪聲信號的自相關(guān)函數(shù);Rsn(τ)為噪聲信號與周期信號的互相關(guān)函數(shù);Rns(τ)為周期信號與噪聲信號的互相關(guān)函數(shù)。

由式(4)可知，由于s(t)和n(t)不相關(guān)，式(4)中的后兩項(xiàng)即Rsn(τ)+Rns(τ)等于0。n(t)為非周期的噪聲，所以Rnn(τ)隨τ的增大而減小。因而當(dāng)τ=τ1時，有:

當(dāng) τ＞ τ1時，可用 Rxx(τ)近似地表示 Rss(τ)。所以當(dāng)τ足夠大時，Rxx(τ)所表現(xiàn)出的周期性就是x(t)中所包含的周期信號s(t)的周期。

采用自相關(guān)理論處理數(shù)字化數(shù)據(jù)，得到的自相關(guān)序列保持了原數(shù)據(jù)的周期，且抑制了噪聲對數(shù)據(jù)的污染。由于實(shí)際測量數(shù)據(jù)中諧波和干擾幅度一般小于有用周期信號幅度，自相關(guān)運(yùn)算的振幅響應(yīng)是非線性函數(shù)，即:原振幅為A，則自相關(guān)函數(shù)振幅為，振幅改變量為倍。相對而言，自相關(guān)運(yùn)算加強(qiáng)了有用周期信號的幅度，換言之，其抑制了諧波和干擾的作用，對于信號周期確定的準(zhǔn)確性十分有利。

式中:x(nΔt)為第n個數(shù)采值，簡記為x(n)。

由于數(shù)采點(diǎn)為有限值，L不可能為無窮，由相關(guān)理論分析可知，信號中包含的噪聲是隨機(jī)的，具有獨(dú)立性，且與周期性信號不相關(guān)，所以假定在一個完整周期內(nèi)進(jìn)行自相關(guān)，Rsn(τ)和 Rsn(τ)也可視為0，Rnn(τ)也得到了充分抑制。根據(jù)信號的周期性可知，在一個周期內(nèi)的自相關(guān)與在若干個(包括無窮)周期內(nèi)的自相關(guān)是一樣的。考慮到經(jīng)自相關(guān)運(yùn)算后得到的自相關(guān)序列應(yīng)至少體現(xiàn)一個完整周期，所以L的長度至少應(yīng)包括兩個完整的信號周期。由于自相關(guān)后的相位歸零，為保證自相關(guān)序列中至少有兩個同方向過零點(diǎn)，L的長度應(yīng)擴(kuò)大到包含3個信號周期。

2 運(yùn)算長度的周期性影響

在一般情況下，本文提出的采樣速率評價法只能估計(jì)信號的周期而無法準(zhǔn)確預(yù)知，即自相關(guān)運(yùn)算的長度L無法保證恰好是信號周期的整數(shù)倍。同時，為便于說明運(yùn)算長度對自相關(guān)序列的周期性影響，假設(shè)信號為 x(t)=sinΩt，積分時間為 T'，則:

式中:等號右側(cè)的第2項(xiàng)幅度可近似地認(rèn)為與ΩT'成反比。當(dāng)采樣周期足夠長時，非整數(shù)周期的影響可以忽略不計(jì)，則與數(shù)采長度 T'無關(guān)。實(shí)際應(yīng)用中，ΩT'＞6π，即大于3個周期即可。一般情況下，數(shù)采長度取信號周期上限的5～10倍并不困難。

式(10)即為最終的自相關(guān)序列計(jì)算公式。通過對Rxx(τ)的邏輯判斷(過零點(diǎn)判斷)，可以方便地得到KT，進(jìn)而求得數(shù)據(jù)采集器的采樣速率。

3 仿真試驗(yàn)及結(jié)論

下面對自相關(guān)方法求取KT進(jìn)行簡單的仿真試驗(yàn)，以驗(yàn)證上述論證。

利用計(jì)算機(jī)產(chǎn)生x(n)數(shù)據(jù)采集序列，對這些序列進(jìn)行自相關(guān)運(yùn)算，得到的結(jié)果如下。

①非整周期運(yùn)算對結(jié)果的影響

當(dāng)x(n)=sin(8π/1024n)時，相當(dāng)于每周期包括256個數(shù)據(jù)采集點(diǎn)，采樣長度為1024，自相關(guān)運(yùn)算長度分別取 1024、924、824、724，得到的自相關(guān)序列采用過零點(diǎn)判斷，數(shù)字化周期均為256。

由以上分析可知，非整周期運(yùn)算對數(shù)字化周期的影響可以忽略。

②噪聲對運(yùn)算結(jié)果的影響

當(dāng)x(n)=sin(8π/1024n)+rand(－0.5，+0.5)時，(rand()表示服從均勻分布的隨機(jī)數(shù))，原始數(shù)據(jù)根本無法采用過零點(diǎn)判斷數(shù)字化周期。采用自相關(guān)運(yùn)算后，當(dāng)運(yùn)算長度為1024、924、824時，數(shù)字化周期為256;當(dāng)運(yùn)算長度為724時，數(shù)字化周期為254或255或256。由以上分析可知，隨著運(yùn)算長度的減小，抑制噪聲能力有所下降，因此，應(yīng)保證足夠的運(yùn)算長度。

③諧波(干擾)對運(yùn)算結(jié)果的影響

當(dāng)自相關(guān)運(yùn)算長度分別為1024、924、824、724時，數(shù)字化周期分別為257、256、255、255。由以上分析可知，自相關(guān)算法具有足夠的抗干擾能力。

將激勵信號加載到單通道和雙通道后，分別得到采樣時間(采樣速率的倒數(shù))。單通道測量時采樣速率快，即采樣時間短;雙通道測量時采樣速率慢，其原因是通道轉(zhuǎn)換需要一定時間。因此，兩者之差即為通道間延遲時間。

4 結(jié)束語

本文為有效抑制采樣速率評價過程中采樣數(shù)據(jù)受噪聲和干擾的影響，提出了一種采用標(biāo)準(zhǔn)信號源提供標(biāo)準(zhǔn)正弦信號和采用自相關(guān)濾波技術(shù)對采樣速率進(jìn)行評價的方法。

該方法能夠有效地抑制采樣數(shù)據(jù)受噪聲和干擾的影響，方法簡單、計(jì)算量小，評價數(shù)據(jù)的一致性好，仿真驗(yàn)證了該方法的有效性。在某測試設(shè)備計(jì)量檢定系統(tǒng)中，該方法利用評價得出的采樣速率和數(shù)字化周期進(jìn)行交流電壓測量和相位測量，滿足了計(jì)量檢定系統(tǒng)的要求。

［1］趙冬青，任勇峰.高速數(shù)字采集系統(tǒng)通道間信號傳輸延遲時間分析［J］.測試技術(shù)學(xué)報，2011，25(2):178 －182.

［2］梁志國，朱濟(jì)杰.數(shù)據(jù)采集系統(tǒng)通道間延遲時間差的精確評價［J］.儀器儀表學(xué)報，1999，20(6):619 －623.

［3］梁志國，孫璟宇.采集速率評價結(jié)果的不確定度［J］.華北工學(xué)院測試技術(shù)學(xué)報，2000，14(1):35 －38.

［4］國家技術(shù)監(jiān)督局 JJG 1048-95.數(shù)據(jù)采集系統(tǒng)校準(zhǔn)規(guī)范［S］.1995.

［5］中國人民解放軍總裝備部 GB 3756-99.測量不確定度的表示及評定［S］.1999.

［6］李金海.誤差理論與測量不確定度評定［M］.北京:中國計(jì)量出版社，2003.

［7］楊德釗，宋凝芳，林志立，等.基于相關(guān)及相位差法的高精度頻率估計(jì)算法［J］.北京航空航天大學(xué)學(xué)報:自然科學(xué)版，2011，37(8):1030－1033.

［8］岳瑞華，許化龍.數(shù)采信號周期確定自相關(guān)方法研究［J］.電子測量與儀器學(xué)報，2002，16(Z):1248 －1251.

［9］岳瑞華，許化龍，劉根旺.基于相關(guān)分析的交流電壓測量［J］.航空計(jì)測技術(shù)，2004，24(1):17 －18.

［10］郭永彩，廉飛宇，林曉鋼.數(shù)字信號處理［M］.重慶:重慶大學(xué)出版社，2009.