陳長鵬,汪禮勝
(武漢理工大學(xué)物理科學(xué)與技術(shù)系,湖北武漢430070)
當頻率為ω0的激光束垂直入射到以速度 v運動的光柵時,從光柵出射的光束各衍射級的頻率可以合寫成[1]
式(1)中,n是衍射級次,nωd是多普勒頻移.ωd與光柵常量 d及速度v的關(guān)系為[1]
因此,如果能測出多普勒頻移 nωd,就可以確定運動光柵的速度v.假如光柵是固定在一運動的物體上,通過測量衍射光的瞬時多普勒頻移就可以定出物體的瞬時運動速度.
要測量瞬時多普勒頻移,首先要分析雙光柵的衍射光和相應(yīng)的光電流.圖1是激光束依次通過動光柵、靜光柵后的衍射光束分布圖,其中光束1和光束2是相互平行的.由于 He-Ne激光器發(fā)出的是有一定寬度的光束(束腰直徑約為1.2 mm),所以只要兩光柵的距離比較小,則兩光束在空間會有重疊的部分,故兩光束能疊加從而形成光拍[1].光拍信號進入硅光電池后,輸出的光電流為[1]
其中ξ為硅光電池的光電轉(zhuǎn)換系數(shù),E10和 E20分別是兩光束光場的振幅,φ1和φ2依次是兩光束光電場的初相位,I0是光電流的直流部分.硅光電池輸出的光電流是很小的,為了便于處理應(yīng)進行放大.
圖1 ±1級衍射光波的多普勒頻移
用硅光電池作檢測元件時,是把它作為電流源來使用的(如圖2的虛線框中所示).但式(3)所表示的光電流is通常很小,必須利用運算放大器進行放大.圖2是 T形網(wǎng)絡(luò)電路組成的微電流放大和 I/V轉(zhuǎn)換電路,在理想的情況下有[2]
將 u=-isR1代入_(4)式可求得輸出電壓為
把式(3)代入式(5)可得
式中
由式(2)和式(6)可以看出:要測出瞬時速度,需要解決瞬時頻率ωd的測量問題.
圖2 微電流放大電路
目前,瞬時頻率的估計一般采用3種方法來解決,分別是:相位差分法、相位建模法和時頻分布法[3].使用這3種方法來處理變頻信號,首先是進行數(shù)據(jù)的采集,然后再送到微機系統(tǒng)用專門的應(yīng)用軟件進行處理得出瞬時頻率值.由于這些估計瞬時頻率的方法依賴特定的應(yīng)用軟件,故缺乏實時性.
在式(6)中,K是不變的常量,U0是信號的直流部分.因此,可采用相應(yīng)的運算電路把式(6)轉(zhuǎn)化為余弦函數(shù)的形式[2],即
式(7)中φ=ωdt+(φ2-φ1).余弦函數(shù)cosφ的冪級數(shù)展開式是
忽略8次冪以上的高次項,用圖3所示的電路完成冪級數(shù)前4項的運算.
在圖3中,冪級數(shù)運算電路由4個乘法器和3個運算放大器組成.容易得到冪級數(shù)運算電路的最終輸出是
因此,當 U5=1 V,U1=2 V,U41U4=4!V,U21U22U3=6!V時,該電路能完成余弦函數(shù)的近似計算,也就是
圖3 冪級數(shù)運算電路
為了解出瞬時頻率,需要用到反余弦函數(shù)電路.圖4是反函數(shù)運算電路,其中N4和N5是2個運算放大器,標有cos uo1的方框代表的是圖3所示的冪級數(shù)運算電路.冪級數(shù)運算電路的輸入端連至A點,輸出端與運算放大器N5的負極性輸入端相連,也就是將原函數(shù)運算電路接入反饋回路以達成實現(xiàn)反函數(shù)運算的目的.
圖4 反余弦函數(shù)運算電路
對于理想的運算放大器,有
根據(jù)圖4及前面的討論,易知
由式(13)和式(14),可得到
因此,該電路實現(xiàn)了對輸入信號 u1反余弦函數(shù)的運算.
由式(7)和式(13)可以看出:uo1對時間t的導(dǎo)數(shù)就是瞬時頻率.這一任務(wù)由圖5中的微分電路來完成.
在圖5中,微分電路的輸出是
而 uo1≈arccos u1,u1=cosφ,所以
由式(15)可以得到瞬時頻率為
根據(jù)式(2)和式(16),瞬時速度為
圖5 微分和計算電路
微分電路的輸出接到單片機系統(tǒng)的A/D轉(zhuǎn)換上.在編程求瞬時頻率和瞬時速度時,先把R0C0和的值存放在1個通用寄存器上.A/D轉(zhuǎn)換的結(jié)果除以 R0C0的值就是瞬時角頻率,瞬時角頻率除以就是動光柵的瞬時速度.經(jīng)D/A轉(zhuǎn)換后,由數(shù)碼管顯示實時頻率和速度的瞬時值.這一測量系統(tǒng),可用于物體運動范圍限制在很小長度內(nèi)的微振動速度的監(jiān)測.
圖6是微振動速度測量系統(tǒng)的原理圖.音叉在換能器的作用下做微振動,動光柵粘貼在音叉上.激光束通過動光柵、靜光柵衍射,光電池對準光束1和光束2的方向(如圖1所示);光電池輸出的光電流通過微電流放大器,轉(zhuǎn)換成電壓(uo)輸出;經(jīng)過整形電路,把 uo轉(zhuǎn)換成標準的余弦信號的輸出;然后,經(jīng)過反余弦函數(shù)運算電路、微分電路后,得到瞬時頻率;最后,經(jīng)過單片機的計算求出瞬時速度的值.
圖6 微振動速度測量系統(tǒng)
表1是實測數(shù)據(jù).其中,換能器的激振頻率為50 Hz,激光打在動光柵的點距音叉的自由端0.03 m,動光柵的振幅為2 mm.
表1 微振動速度實測數(shù)據(jù)
上面這組數(shù)據(jù)是取自20 m s內(nèi)的測量結(jié)果,也就是1個音叉振動周期內(nèi)的數(shù)據(jù).不難看出:上面數(shù)據(jù)的變化周期是10 m s,這是由于沒有辨向電路導(dǎo)致的結(jié)果.
在1/4周期(5 m s)內(nèi),同等條件下,對動光柵進行動力學(xué)分析再結(jié)合數(shù)值計算,并且忽略動光柵的質(zhì)量,得到的動光柵瞬時速度的理論值如表2所示.
表2 1/4周期內(nèi)動光柵瞬時速度的理論值 m/s
把實測的數(shù)據(jù)和理論值進行對比,不難發(fā)現(xiàn)除了速度為0的3個數(shù)據(jù)的對偏差超過5%以外,其余的各個數(shù)據(jù)的相對偏差都未超過5%.之所以在速度為0的附近偏差較大,是因為采用了微分電路,速度為0對應(yīng)頻率為0.其余各組數(shù)據(jù)的偏差,主要是用級數(shù)來近似地代替余弦函數(shù)所導(dǎo)致.
雖然在速度為0的附近測量數(shù)據(jù)的偏差較大,但是其余的各個數(shù)據(jù)相對偏差都未超過5%,說明該方法是測量微運動瞬時速度的有效方法.怎樣解決速度為0的附近測量數(shù)據(jù)的誤差較大,這是一個值得探討和有待解決的問題.
[1] 龍作友,戴亞文,楊應(yīng)平,等.大學(xué)物理實驗[M].武漢:武漢理工大學(xué)出版社,2008.
[2] 張國雄.測控電路[M].北京:機械工業(yè)出版社,2001:144-165.
[3] 陳水波.利用雙光柵的多普勒頻移測速度[J].物理實驗,2007,27(7):6-9.
[4] 李正正,蔡虹,洪小剛,等.雙光束激光多普勒測速系統(tǒng)[J].物理實驗,2005,25(3):44-46.
[5] 馮松立,陳高平.瞬時頻率估計的相位建模法及Matlab的實現(xiàn)[J].中國測試技術(shù),2003,29(3)38,41.