姚述福，瞿　洋，儲飛黃

(電子工程學(xué)院，合肥 230037)

?

基于基帶噪聲互相關(guān)法的同源多路信號相位差測量

姚述福，瞿洋，儲飛黃

(電子工程學(xué)院，合肥 230037)

摘要：根據(jù)同一激勵源信號通過功分器得到的同源多路信號間表現(xiàn)為純時延的特性,本測量方法對僅存在時延差異的同源多路噪聲調(diào)頻信號進(jìn)行自相關(guān)運算和互相關(guān)運算，再根據(jù)相關(guān)運算得到的不同譜峰位置計算出信號間的相對時延，也就是信號間的相位差。信號間相位差的計算精度與采樣頻率、采樣時間、調(diào)頻指數(shù)有關(guān)。隨著上述3個參數(shù)值的增大，相位差的測量精度增加。

關(guān)鍵詞：噪聲調(diào)頻；互相關(guān)；采樣頻率；相位差；調(diào)頻指數(shù)

0引言

空間功率合成技術(shù)突破了功率器件和傳統(tǒng)單機(jī)的功率限制，是當(dāng)前獲得超大功率最有效的方法。在通信干擾領(lǐng)域，基于空間功率合成技術(shù)的超短波通信干擾站等干擾系統(tǒng)已相繼投入使用?？臻g功率合成技術(shù)利用同源多路信號分別經(jīng)過功放放大后在空間遠(yuǎn)場某個方位進(jìn)行疊加后形成超大功率，而影響其合成功率大小的最主要因素是同源多路信號間的相位差?，F(xiàn)有的測量多路信號間相位差的方法主要有IQ正交求解[1]、數(shù)字相關(guān)法[2]、傅里葉變換(FFT)法[3]和相位差分法[4]等，上述方法各有利弊，但主要研究對象都是單頻信號，相對簡單。本文針對裝備中常用的噪聲調(diào)頻干擾信號，提出一種基于基帶噪聲互相關(guān)法的同源多路信號相位差測量方法，該方法利用噪聲信號互相關(guān)峰更加尖銳的優(yōu)點，大大提高了相位差測量的精度。該方法可以用于精確測量路間相位差，也可以為進(jìn)一步測量空間功率合成的等效合成功率打下技術(shù)基礎(chǔ)。

1算法理論描述

對于平穩(wěn)的隨機(jī)信號x(t)和y(t)，其互相關(guān)函數(shù)可以表示為：

Rxy(τ)=E[x(t)y(t-τ)]=

(1)

當(dāng)2路甚至多路信號之間表現(xiàn)為純延時的特性，即當(dāng)y (t)=x (t-D)時，兩者的互相關(guān)函數(shù)變?yōu)椋?/p>

Rxy(τ)=E[x(t)x(t-D-τ)]=Rxx(τ-D)

(2)

式中：Rxx為x(t)的自相關(guān)函數(shù)，此時兩信號的互相關(guān)結(jié)果就成為自相關(guān)函數(shù)的時延。

而根據(jù)信號的自相關(guān)函數(shù)的性質(zhì)，對于任何τ≠0， 都有Rxx(τ)≤Rxx(0)，即τ=0時Rxx(τ)為其最大值。因此易知Rxy(τ)在τ=D時取最大值，這樣就可以從Rxy(τ)的峰值點位置對應(yīng)的τ值測出延時D， 如圖1所示。

圖1　互相關(guān)法測量時延

而對于白噪聲來說，它的自相關(guān)函數(shù)為沖激函數(shù)，因此白噪聲與它的時延信號的互相關(guān)函數(shù)就是沖激函數(shù)的時延。而沖激函數(shù)的峰值非常尖銳，能夠更加準(zhǔn)確地定位其時延，如圖2所示。

圖2　白噪聲與其時延的互相關(guān)

裝備中常用的噪聲調(diào)頻干擾信號是將一段帶限的白噪聲調(diào)制在頻率上，實現(xiàn)寬帶或窄帶的噪聲干擾。而對于采用空間功率合成干擾技術(shù)的裝備來說，其干擾激勵器會產(chǎn)生同源多路干擾激勵信號，這樣的多路激勵信號間的相位差對于最終合成干擾功率有著決定性的影響，因而需要重點研究求解。考慮到噪聲互相關(guān)峰的尖銳性對準(zhǔn)確求解其相位差的優(yōu)點，這里采取先將噪聲調(diào)頻干擾信號解調(diào)，得到基帶噪聲信號，再將多路基帶噪聲信號進(jìn)行互相關(guān)運算，最終得到多路信號的相位差。

2測量方法仿真分析

為了驗證算法的可行性，進(jìn)行了MATLAB仿真。以2個通道為例，對一段隨機(jī)產(chǎn)生的高斯白噪聲用低通濾波器進(jìn)行限帶，然后對上述帶限高斯白噪聲進(jìn)行頻率調(diào)制，得到噪聲調(diào)頻信號。對仿真產(chǎn)生的噪聲調(diào)頻信號先加入不同的高斯白噪聲和時延，來模擬不同通道含噪聲和時延的實際信號，然后分別對上述不同通道信號進(jìn)行非相干解調(diào)，最后對各通道解調(diào)出的基帶噪聲信號進(jìn)行互相關(guān)運算，得到通道間的時延。

針對窄帶噪聲調(diào)頻和寬帶噪聲調(diào)頻2種情況分別做了MATLAB仿真。

主要仿真參數(shù)：fc=100MHz，fs=800MHz，采樣點數(shù)N=80 000，調(diào)制噪聲帶寬300kHz，調(diào)制噪聲功率10dBW，已調(diào)信號信噪比30dB。

(1) 窄帶調(diào)頻，有效調(diào)頻指數(shù)mf=0.5。圖3、圖4、圖5分別表示生成噪聲調(diào)頻信號、對其進(jìn)行時延處理、對信號進(jìn)行相關(guān)得到的時域圖或頻譜圖。根據(jù)多次仿真結(jié)果，此時時延誤差小于等于1個采樣點的概率為40.25%。

(2) 寬帶調(diào)頻，有效調(diào)頻指數(shù)mf=20。圖6、圖7、圖8分別表示生成的噪聲調(diào)頻信號、對其進(jìn)行時延處理、對信號進(jìn)行相關(guān)得到的時域圖或頻譜圖。根據(jù)多次仿真結(jié)果，此時時延誤差小于等于1個采樣點的概率為58.63%。

3測量仿真的精度分析

噪聲調(diào)頻信號相位參數(shù)的測量主要考慮在不同的信噪比和調(diào)頻指數(shù)的情況下，相位參數(shù)計算方法的精度問題。

為了分析測量精度，對噪聲調(diào)頻信號間的相位差測量進(jìn)行多次仿真。仿真條件如前所示：此時信號頻率為f0=100MHz，采樣頻率為400MHz，采樣時間為0.1ms，采樣點數(shù)為40 000，調(diào)制噪聲的帶寬為300kHz，調(diào)制噪聲的功率為10dBW，調(diào)頻指數(shù)為5，信號信噪比為30dB時，假設(shè)2路信號間的時延為100個采樣點，使用相關(guān)法計算得到2路信號間的時延。在此情況下仿真5 000次，得到計算的時延誤差點數(shù)的分布如圖9所示。

圖3　噪聲波形與頻譜

圖4　噪聲波形調(diào)制與解調(diào)

圖5　基帶噪聲互相關(guān)結(jié)果

圖6　噪聲波形與頻譜

從圖9圖可以看到，測量誤差小于等于1采樣點的概率為80%左右。不同的參數(shù)會影響到測量誤差，下面分析不同參數(shù)對時延計算誤差的影響。

由于在實際的測量當(dāng)中，對于2路信號間的相位時延是在納秒級的精度，以信號頻率為100 MHz為例，只有當(dāng)采樣點的精度達(dá)到1 ns時，對于相位的精度才能達(dá)到±18°。因此，在對噪聲調(diào)頻采用相關(guān)法測量信號間相位差時需要很高的采樣精度。而為了達(dá)到1 ns的采樣精度，需要采樣頻率達(dá)到1 GHz，下面在采樣頻率為1 GHz時進(jìn)行分析。此時，影響其精度的因素有信噪比、調(diào)頻指數(shù)、采樣時長(采樣點數(shù))。

表1分析了在采樣頻率為1 GHz、信噪比為30 dB的情況下，不同的采樣時間和調(diào)頻指數(shù)情況下，時延誤差小于1 ns的概率。

圖7　噪聲波形調(diào)制與解調(diào)

圖8　基帶噪聲互相關(guān)結(jié)果

采樣時間調(diào)頻指數(shù)0.1ms0.2ms0.5ms1ms0.540.8%63.2%62.5%89.5%142.6%49.4%91.5%99.7%248.8%67.2%98.5%99.8%551.4%72.4%95%99.8%1056.0%75.4%93%99.7%

圖9　時延誤差點數(shù)的分布

從表1可以得到，調(diào)頻指數(shù)的增加可以降低誤差，但是在增大到一定的數(shù)值時對減小誤差的作用就不再有幫助。而增大采樣時間對減小誤差的作用則是非常明顯，在采樣時間達(dá)到1 ms時，對大于等于1的調(diào)頻指數(shù)來說誤差大于1 ns的概率小于1%，達(dá)到很高的精度。

為了得到更高的相位精度，需要對采樣頻率進(jìn)行進(jìn)一步的提升，以100 MHz的信號頻率為例，為了得到±1.8°的精度，需要采樣頻率達(dá)到10 GHz，同時進(jìn)行1 ms時間的采樣，此時對測量設(shè)備的要求也就更高了。

4結(jié)束語

從以上的分析來看，在本測量方法中，使用同一采樣頻率的情況下，窄帶調(diào)頻(mf?1)的測量精度比寬帶調(diào)頻(mf?1)稍差，隨著采樣時間的增加，算法的計算精度迅速升高?？偟膩碚f，本文算法的精度值較高，但對設(shè)備提出了一定要求。

參考文獻(xiàn)

[1]戈穩(wěn).雷達(dá)接收機(jī)技術(shù)[M].北京：電子工業(yè)出版社，2005.

[2]路艷潔，席志紅，王姜鉑.FFT法與數(shù)字相關(guān)法在相位測量上的比較[J].信息技術(shù)，2007(12):105-108.

[3]張海濤，涂亞慶.基于FFT的一種計及負(fù)頻率影響的相位差測量新方法[J].計量學(xué)報，2008,29(2):168- 171.

[4]李輝，王巖飛.正弦信號的直接FFT參數(shù)估計與相位差分法對比研究[J].電子與信息學(xué)報，2010,32(3):544-547.

Phase Difference Measurement of Homologous Multi-channel Signals Based on Baseband Noise Cross-correlation Method

YAO Shu-fu,QU Yang,CHU Fei-huang

(Electronic Engineering Institute,Hefei 230037，China)

Abstract:According to the characteristics of pure time delay revealed among homologous multi-channel signals generated by the same excitation source signal getting through the power divider,the measurement method performs autocorrelation and cross-correlation calculation to homologous multi-channel noise frequency modulation (FM) signals with only time-delay difference,then the relative time delay among signals is calculated according to the different spectrum peak position deduced by correlation calculation,i.e.phase difference among signals.The calculation accuracy of the phase difference among signals is affected by the sampling frequency,sampling time,FM index.With the value of above three parameters increases,the measurement accuracy of phase difference improves.

Key words:noise frequency modulation;cross-correlation;sampling frequency;phase difference;frequency modulation index

DOI:10.16426/j.cnki.jcdzdk.2016.01.026

中圖分類號：TN98

文獻(xiàn)標(biāo)識碼：A

文章編號：CN32-1413(2016)01-0115-05

收稿日期:2015-05-10