1 引 言
電容式濕度傳感器由于具有測濕范圍寬、響應(yīng)速度快、溫漂小、穩(wěn)定性好、使用方便等優(yōu)點得到了廣泛的應(yīng)用. 濕度傳感器生產(chǎn)企業(yè)在批量生產(chǎn)和降低成本上要解決的關(guān)鍵問題是產(chǎn)品批量檢測問題.本文介紹了一種可用于電容式濕度傳感器測量的諧波分析測量法.
在電容式濕度傳感器測量實際得到的測量信號中存在頻率范圍較寬的低頻和高頻干擾. 本文的諧波分析測量法是基于所測信號具有穩(wěn)定頻率(1kHz)的特點,通過對傅里葉變換得到的頻譜進(jìn)行重構(gòu),有效地將非信號頻率的各種干擾消除,使信噪比得到較大提高,獲得了理想的測量結(jié)果 .
2 離散傅里葉變換
數(shù)字頻譜分析最基本的方法是傅里葉變換. 在分析帶有許多不同頻率的復(fù)雜信號時, 傅里葉變換是一種廣泛應(yīng)用而有效的信號處理方法. 對信號x ( t)的采樣在有限長度的樣本空間記錄上進(jìn)行, 將對x ( t)采樣得到數(shù)據(jù)序列記為x ( n) ( n = 0, 1, ., N- 1) . 設(shè)信號的周期為T, 樣本長度為N, 采樣時間間隔為Δt,即T =NΔt,采樣頻率fs = 1 /Δt,則譜線的頻率間隔Δf = 1 / T, 即在頻域內(nèi)其幅頻曲線是由N條離散譜線獨立組成的x ( n)的離散傅里葉變換為
離散傅里葉變換的意義在于可以對任意連續(xù)的時域信號進(jìn)行抽樣和截斷, 然后進(jìn)行傅里葉變換得到一系列離散型頻譜,該頻譜的包絡(luò)線,就是原來連續(xù)信號真實頻譜的估計值.
3 諧波分析測量法原理
在A /D轉(zhuǎn)換后采集的數(shù)字信號中主要的頻譜成分是正弦波信號源頻率(1kHz)的成分,并包含一定成分的A /D量化的噪聲、工頻干擾和其他高頻干擾. 諧波分析測量方法是在對A /D轉(zhuǎn)換得到的數(shù)字信號采用數(shù)字頻譜分析的方法求出信號所對應(yīng)的諧波分量,保留1kHz的成分,并從信號中將1kHz之外頻率成分去掉,再對得到的1kHz的成分的幅值和相位進(jìn)行數(shù)據(jù)處理就可獲得電容式傳感器的電容值C和品質(zhì)因數(shù)Q. 為消除工頻干擾,對信號進(jìn)行測量時要保證采樣長度達(dá)到工頻周期,也要保證對1kHz成分周期內(nèi)有足夠數(shù)量的采樣數(shù)據(jù). 所以,需要較多的采樣數(shù)據(jù).
首先,對傳感器信號進(jìn)行采集數(shù)據(jù),之后對獲得的數(shù)據(jù)進(jìn)行傅里葉變換, 求出信號的1kHz所對應(yīng)的諧波成分的幅值A(chǔ)1和相位θ1. 其次, 對正弦波信號源進(jìn)行同樣次數(shù)的采樣, 對獲得的數(shù)據(jù)進(jìn)行傅里葉變換,求出正弦波信號的1kHz所對應(yīng)的諧波成分的幅值A(chǔ)2和相位θ2. 通過θ1和θ2的差可得到傳感器的品質(zhì)因數(shù)Q, 通過A1 得到傳感器的電容值C. A2用于對由于正弦波信號源幅值變化及模擬環(huán)節(jié)其他因素產(chǎn)生比例變化對電容值測量結(jié)果的影響進(jìn)行修正.
取N = 128 ×20 = 2 560,以128kHz的采樣頻率分別對信號和正弦波信號源進(jìn)行A /D轉(zhuǎn)換,在一個工頻周期時間內(nèi)得到兩個長為2 560的數(shù)據(jù)序列.此時, 1kHz就是20次諧波, 工頻干擾在基頻上, 即信號數(shù)據(jù)序列xn、正弦波信號源數(shù)據(jù)序列yn ( n = 0,1, ., 2 559) . 采用式( 1)對xn、yn 進(jìn)行離散傅里葉變換,得到X20、Y20. 對k≠20情況的Xk、Yk 由于不是信號的頻譜, 作為噪聲處理(這樣處理相當(dāng)于進(jìn)行了數(shù)字帶通濾波) ,就不再進(jìn)行計算,所以計算量很小. 之后按下述方法對X20、Y20進(jìn)行數(shù)據(jù)處理.設(shè)X20的實部為Xa20、X20的虛部為Xb20 ,即
進(jìn)行這樣的測量后,可以將工頻干擾和其他低頻干擾作為19次以下諧波分量消除,將高頻干擾作為21次以上諧波分量消除.
4 仿真分析和編程實現(xiàn)
為了對諧波分析測量法的測量效果進(jìn)行分析,采用Matlab對合成信號進(jìn)行了仿真 . 在考慮工頻干擾和高頻隨機(jī)干擾成分的情況下,構(gòu)造的仿真信號為x ( t) = sin (2πf) + 0.1sin (40πf) + 0.1sin (200πf) +0.1sin (0.1πf). 其中, sin ( t)為1kHz正弦信號,信號的幅值為1; sin (40πf) 、sin (200πf)分別為20、100次諧波的高頻干擾, 強(qiáng)度為0.1; sin (0.1πf)為工頻干擾,強(qiáng)度為0.1.
對此模擬信號按前面所述的方法, 以128kHz的速度進(jìn)行數(shù)據(jù)采集之后, 采用Matlab進(jìn)行128 ×20點DFT變換處理、重構(gòu)信號,將重構(gòu)信號與無干擾信號進(jìn)行對比. 仿真結(jié)果如圖1, 圖1a是仿真信號的波形, 圖1b是采用諧波分析法處理后的波形.對上述仿真結(jié)果進(jìn)行分析可見, 經(jīng)過諧波分析處理后信號的信噪比得到大幅度提高. 仿真結(jié)果表明, 采用諧波分析法對信號處理可以有效地消除高次諧波和低頻干擾的影響, 可以提高信噪比約12dB.
從前面諧波分析測量法的原理中可以看出,本測量法只需計算一個頻譜,進(jìn)行x ( k) 與e- j2πnN ( k =20 )的N次復(fù)數(shù)乘法運(yùn)算和N - 1次的復(fù)數(shù)加法運(yùn)算. 因此計算量很小,在DSP編程時一個N 次循環(huán)就可以完成計算.
5 結(jié) 語
將電容式濕度傳感器測量的諧波分析測量法用于電容式濕度傳感器測試系統(tǒng)中,與其他電容式濕度傳感器的測量方法(如,復(fù)數(shù)電壓測量法)相比具有較強(qiáng)的消除高頻干擾和低頻干擾能力,計算量比復(fù)數(shù)電壓測量法要大,但對于DSP系統(tǒng)來說已經(jīng)不是問題. 此方法還可以用于通過加載固定頻率正弦波進(jìn)行測量的測試儀器中,將會增強(qiáng)此類測試儀器的測量精確度和抗干擾能力.