黃 斌，尹 光，張 昊

（1 解放軍92941 部隊(duì)，遼寧 葫蘆島 125000；2 南京長(zhǎng)峰航天電子科技有限公司，南京 210018）

0 引言

隨著雷達(dá)系統(tǒng)和數(shù)字信號(hào)處理技術(shù)的發(fā)展，雷達(dá)系統(tǒng)的功能更加復(fù)雜，系統(tǒng)集成度越來(lái)越高。模擬器件固有的幅相誤差，不可避免地會(huì)帶來(lái)相控陣?yán)走_(dá)各收發(fā)通道間的幅相特性不一致，從而導(dǎo)致相控陣?yán)走_(dá)性能下降等問(wèn)題。為此，相關(guān)人員對(duì)雷達(dá)收發(fā)通道幅相一致性校準(zhǔn)問(wèn)題進(jìn)行了研究。如：文獻(xiàn)［2］中通過(guò)構(gòu)造以陣列幅相誤差為未知量的代價(jià)函數(shù)，利用迭代的方法求取未知量的估計(jì)值，但該方法初始迭代值的選擇會(huì)對(duì)迭代收斂速度和結(jié)果造成了一定影響。文獻(xiàn)［3-5］采用最大似然原理構(gòu)建非線性的代價(jià)函數(shù)，得到了一種幅相誤差的自校正算法。雖然該方法估計(jì)精度高，但計(jì)算量較大。文獻(xiàn)［6-8］中提出根據(jù)通道幅相數(shù)據(jù)的希望權(quán)值，對(duì)通道幅相誤差進(jìn)行估計(jì)，按照不同的頻點(diǎn)進(jìn)行通道校準(zhǔn)。文獻(xiàn)中采用頻域校準(zhǔn)方法，需要將一定帶寬的接收通道劃分成帶寬較窄的多個(gè)子通道，同時(shí)對(duì)各個(gè)子通道接收到的數(shù)據(jù)進(jìn)行FFT 變換，在頻域提取對(duì)應(yīng)的幅度信息和相位信息。針對(duì)具有一定帶寬的信號(hào)（接收通道），該方法仍無(wú)法精確計(jì)算其通道之間的幅度差和相位差，計(jì)算時(shí)間明顯加長(zhǎng)。由此可以看出，目前陣列幅相誤差的物理測(cè)量或自校正算法仍存在不足，亟待后續(xù)改進(jìn)。

本文研究從時(shí)域信號(hào)出發(fā)，在時(shí)域完成延時(shí)和初相誤差的估計(jì)；采用過(guò)門(mén)限檢測(cè)技術(shù)，完成信號(hào)幅度校準(zhǔn)計(jì)算。相位提取時(shí)，利用寬帶信號(hào)的瞬時(shí)窄帶特性，構(gòu)造初相和觀測(cè)信號(hào)模型，采用最小二乘法實(shí)現(xiàn)相位特征參數(shù)的提取，并對(duì)其進(jìn)行優(yōu)化、轉(zhuǎn)化為二位矩陣計(jì)算；利用FFT 計(jì)算包絡(luò)自相關(guān)。因此，該方法不僅適用于任意樣式的雷達(dá)發(fā)射信號(hào)，而且由于所采用的方法均對(duì)噪聲有一定的適應(yīng)性，能夠?qū)π旁氡鹊陀?0 dB 以下的系統(tǒng)進(jìn)行校準(zhǔn)，適用于大多數(shù)的雷達(dá)發(fā)射系統(tǒng)。此外，算法還具有易于硬件實(shí)現(xiàn)、計(jì)算復(fù)雜度低、計(jì)算速度快等優(yōu)勢(shì)。

1 通道幅相校準(zhǔn)

基于最小二乘法的有源相控陣幅相校準(zhǔn)流程中，首先選定參考通道，采用希爾伯特變換，將采集的信號(hào)變換為解析信號(hào)；采用包絡(luò)檢波模塊提取信號(hào)包絡(luò)，利用過(guò)門(mén)限檢測(cè)電路計(jì)算幅度校準(zhǔn)系數(shù)，完成信號(hào)幅度校準(zhǔn)；對(duì)幅度校準(zhǔn)后的信號(hào)，利用包絡(luò)相關(guān)法計(jì)算各待校準(zhǔn)通道相對(duì)于參考通道的傳輸延時(shí)，并通過(guò)拋物線插值計(jì)算精延時(shí)，根據(jù)得到的傳輸延時(shí)，進(jìn)行待校準(zhǔn)通道的信號(hào)延時(shí)校準(zhǔn)；在完成各待校準(zhǔn)通道的延時(shí)校準(zhǔn)后，選取一定的信號(hào)樣本，采用最小二乘法計(jì)算各待校準(zhǔn)通道的初相，從而確定待校準(zhǔn)通道與參考通道之間的相位差，根據(jù)計(jì)算得到的相位差進(jìn)行各待校準(zhǔn)通道的相位差異校準(zhǔn)，使得各信號(hào)通道的幅相特性保持一致。

1.1 通道幅度校準(zhǔn)

在選定參考通道后，對(duì)各待校準(zhǔn)通道的幅度（增益）差異進(jìn)行校準(zhǔn)。將各通道模擬信號(hào)轉(zhuǎn)變?yōu)閿?shù)字信號(hào)，再利用Hilbert 變換，將其轉(zhuǎn)換成解析信號(hào)、即復(fù)信號(hào)。通過(guò)對(duì)各通道的解析信號(hào)進(jìn)行包絡(luò)檢波，從而計(jì)算得到復(fù)信號(hào)的幅度值，即：

其中，為復(fù)信號(hào)的實(shí)部；為復(fù)信號(hào)的虛部；為復(fù)信號(hào)的幅度值。

利用2 個(gè)通道過(guò)門(mén)限信號(hào)包絡(luò)幅度樣本點(diǎn)的平均值，計(jì)算幅度校準(zhǔn)系數(shù)。門(mén)限值通常設(shè)定為當(dāng)前接脈沖信號(hào)幅度最大值的0.2 倍。假定參考通道過(guò)門(mén)限信號(hào)包絡(luò)幅度的平均值為＿，待校準(zhǔn)通道過(guò)門(mén)限的信號(hào)包絡(luò)幅度的平均值為＿，則幅度校準(zhǔn)系數(shù)為：

根據(jù)求得的幅度校準(zhǔn)系數(shù)＿，對(duì)待校準(zhǔn)通道的信號(hào)包絡(luò)幅度和信號(hào)幅度進(jìn)行校準(zhǔn)。使得待校準(zhǔn)通道的信號(hào)包絡(luò)幅度和信號(hào)幅度與參考通道保持一致。在FPGA 中，采用一個(gè)復(fù)數(shù)乘法器即可完成校準(zhǔn)。

1.2 通道延時(shí)校準(zhǔn)

完成包絡(luò)幅度校準(zhǔn)后，待校準(zhǔn)通道信號(hào)包絡(luò)進(jìn)行自相關(guān)粗延時(shí)計(jì)算。為了避免快速傅里葉變換計(jì)算卷積帶來(lái)的疊加模糊問(wèn)題，對(duì)計(jì)算輸入的信號(hào)包絡(luò)進(jìn)行擴(kuò)展延拓、即補(bǔ)零操作，使得2 個(gè)序列的長(zhǎng)度相同。參考通道延拓后的序列可直接進(jìn)行快速傅里葉變換，而對(duì)于待校準(zhǔn)通道，則先將延拓后的信號(hào)包絡(luò)序列進(jìn)行翻轉(zhuǎn)后，再進(jìn)行快速傅里葉變換。該算法實(shí)現(xiàn)過(guò)程可詳述如下。

（1）將參考通道的傅里葉變換結(jié)果與待校準(zhǔn)通道的傅里葉變換結(jié)果依次相乘，并對(duì)相乘后的序列進(jìn)行傅里葉逆變換，得到變換后的序列，其幅度序列為（）。找到傅里葉逆變換后幅度序列幅度最大值對(duì)應(yīng)的位置，并將其投影到下標(biāo)投影序列上，從而可以計(jì)算得到待校準(zhǔn)通道相對(duì)于參考通道的粗延時(shí)。最大值對(duì)應(yīng)的位置為，下標(biāo)投影計(jì)算得到的位置對(duì)應(yīng)的投影序列值為（），則粗延時(shí)的值為：

其中，為信號(hào)AD 的采樣率。

（2）利用插值計(jì)算精延時(shí)計(jì)算方法如下：根據(jù)相關(guān)系列的幅度序列，計(jì)算最大幅度值（）、最大幅度值左側(cè)的幅度值（1）以及右側(cè)幅度值（1）。精延時(shí)計(jì)算公式如下：

其中，t為計(jì)算得到的精延時(shí)。

將信號(hào)粗延時(shí)與精延時(shí)相加，就會(huì)得到待校準(zhǔn)通道相對(duì)于參考通道的精確延時(shí)，即：

（3）延時(shí)校準(zhǔn)。根據(jù)計(jì)算得到的精確延時(shí)，對(duì)待校準(zhǔn)通道進(jìn)行信號(hào)延時(shí)校準(zhǔn)。校準(zhǔn)實(shí)現(xiàn)如下。

計(jì)算待校準(zhǔn)信號(hào)為：

進(jìn)行延時(shí)校準(zhǔn)后的信號(hào)為：

其中，為待校準(zhǔn)通道信號(hào)進(jìn)行幅度校準(zhǔn)后的數(shù)字信號(hào)。

1.3 相位差異校準(zhǔn)

在完成待校準(zhǔn)通道的信號(hào)延時(shí)校準(zhǔn)后，采用最小二乘法計(jì)算參考通道與待校準(zhǔn)通道的初相，從而得到待校準(zhǔn)通道與參考通道之間的相位差。構(gòu)建通道信號(hào)與其相位的數(shù)學(xué)模型為：

其中，表示通道信號(hào)的原始數(shù)據(jù)；為該通道信號(hào)的頻率值；［（1）（2）］，這里，（1）為通道信號(hào)初相對(duì)應(yīng)的余弦值，（2）為通道信號(hào)初相對(duì)應(yīng)的正弦值，即為該通道信號(hào)初相的正弦值和余弦值組成的二維向量。對(duì)此研究，擬展開(kāi)研究論述如下。

（1）計(jì)算通道信號(hào)初相的正弦值和余弦值。根據(jù)通道信號(hào)與其相位的數(shù)學(xué)模型公式（8），可以按照公式（9）計(jì)算該通道信號(hào)初相的正弦值和余弦值，即：

（2）計(jì)算通道信號(hào)的初相。根據(jù)通道信號(hào)的初相的正弦值和余弦值，可以計(jì)算得到通道信號(hào)S 對(duì)應(yīng)的初相，即：

利用解模糊算法，計(jì)算得到對(duì)應(yīng)時(shí)刻通道信號(hào)的相位值。算法如下：

由于通道信號(hào)已經(jīng)過(guò)翻轉(zhuǎn)處理，因此需要對(duì)初相值進(jìn)行校準(zhǔn)。校準(zhǔn)后的通道信號(hào)的初相為360。

根據(jù)上述流程，分別計(jì)算得到參考通道信號(hào)的初相值與待校準(zhǔn)通道信號(hào)的初相值，即計(jì)算待校準(zhǔn)通道信號(hào)與參考通道信號(hào)的相位差為。

（3）相位差異校準(zhǔn)。在得到2 個(gè)通道的相位差后，利用公式來(lái)進(jìn)行待校準(zhǔn)通道信號(hào)相位差異校準(zhǔn)。此處需用到的公式為：

其中，＿為完成延時(shí)校準(zhǔn)后的信號(hào)，＿＿為完成幅度、延時(shí)和相位校準(zhǔn)后的待校準(zhǔn)通道輸出信號(hào)。

2 仿真實(shí)現(xiàn)

為了驗(yàn)證本文校準(zhǔn)方法的性能，對(duì)多通道射頻信號(hào)的幅相延時(shí)特性差異校準(zhǔn)進(jìn)行仿真分析。由于多通道的校準(zhǔn)與雙通道類似，均以其中一個(gè)通道為參考，剩余通道以此為基準(zhǔn)校準(zhǔn)。因此，本文僅對(duì)選取其中2 個(gè)信號(hào)通道的校準(zhǔn)過(guò)程和結(jié)果進(jìn)行詳細(xì)討論。

仿真實(shí)驗(yàn)中，相關(guān)參數(shù)設(shè)置如下：射頻校準(zhǔn)信號(hào)的中心頻率為2 GHz，脈寬為5 us，脈沖重復(fù)周期為100 us。采集數(shù)據(jù)板的中頻頻率為250 MHz，A／D 采樣頻率為1 GHz，校準(zhǔn)信號(hào)樣本點(diǎn)數(shù)為8 192 點(diǎn)，采集起始距脈沖前沿800 ns，發(fā)射信號(hào)信噪比15 dB。

在系統(tǒng)通道校準(zhǔn)過(guò)程中，由信號(hào)源提供中頻測(cè)試信號(hào)，測(cè)試信號(hào)通過(guò)功分網(wǎng)絡(luò)分成多路，分別注入到相控陣?yán)走_(dá)系統(tǒng)的各個(gè)發(fā)射通道中，完成發(fā)射通道對(duì)應(yīng)的變頻、放大等發(fā)射處理。在射頻輸出端，采用多通道數(shù)采器對(duì)相控陣?yán)走_(dá)發(fā)射輸出端的多路射頻信號(hào)進(jìn)行采集，選取其中一路通道作為參考通道、即通道1，其余通道作為待校準(zhǔn)通道。數(shù)采器對(duì)經(jīng)過(guò)各通道的信號(hào)進(jìn)行采集，并將輸入的射頻信號(hào)進(jìn)行下變頻，得到250 MHz 的中頻信號(hào)，再通過(guò)A／D 模塊采集信號(hào)中頻、數(shù)據(jù)，數(shù)據(jù)采樣率為1 GHz。利用本文的幅相校準(zhǔn)方法，對(duì)中頻信號(hào)進(jìn)行Hilbert 變換后，采用相關(guān)法計(jì)算待校準(zhǔn)通道信號(hào)相對(duì)于參考通道信號(hào)的幅度差異及延時(shí)差異，并進(jìn)行校準(zhǔn)；再利用最小二乘法，分別估計(jì)2 個(gè)通道的相位差，并對(duì)待校準(zhǔn)通道的信號(hào)進(jìn)行相位補(bǔ)償。對(duì)于多個(gè)射頻通道，實(shí)現(xiàn)方式采用兩兩校準(zhǔn)，完成整個(gè)系統(tǒng)差異性補(bǔ)償。

2.1 采集校準(zhǔn)參考信號(hào)

經(jīng)過(guò)不同通道的信號(hào)時(shí)域圖及信號(hào)包絡(luò)圖如圖1、圖2 所示。由圖1、圖2 可知，不同發(fā)射通道之間，由于器件特性存在差異，經(jīng)過(guò)通道后輸出的射頻信號(hào)在幅度增益、鏈路延時(shí)和信號(hào)初相之間存在一定的差異。

圖1 經(jīng)過(guò)不同通道后的信號(hào)時(shí)域圖Fig.1 Time-domain signals after different channels

圖2 經(jīng)過(guò)不同通道后的信號(hào)包絡(luò)圖Fig.2 Envelope diagram of signals after different channels

2.2 幅度延時(shí)校準(zhǔn)

計(jì)算參考通道信號(hào)與待校準(zhǔn)通道信號(hào)的幅度差異和延時(shí)差異。其中，幅度差異由包絡(luò)均值對(duì)比得到，而高精度延時(shí)差異值則采用包絡(luò)相關(guān)和多項(xiàng)式插值法計(jì)算得到，并對(duì)幅度延時(shí)進(jìn)行補(bǔ)償。與傳統(tǒng)相關(guān)法相比，采用包絡(luò)相關(guān)法可以有效避免信號(hào)周期起伏帶來(lái)的峰值穩(wěn)定性差異，使該求解方法能夠適用于較低的信噪比環(huán)境；采用拋物線插值可以提高延時(shí)求解的精度。對(duì)于傳統(tǒng)相關(guān)法，求解精度為1（為采樣率）。本次仿真中，傳統(tǒng)相關(guān)法求解精度只能達(dá)到1 ns，而采用插值后，求解精度能到達(dá)到0.1 ns。對(duì)于上述校準(zhǔn)信號(hào)，計(jì)算得到延時(shí)差異為1.416 428 077 101 513e-08 s，即14.164 280 771 ns，其延時(shí)精度為0.1 ns。根據(jù)計(jì)算得到的延時(shí)差對(duì)通道2的信號(hào)進(jìn)行補(bǔ)償，結(jié)果如圖3 和圖4 所示。

圖3 經(jīng)過(guò)幅度延時(shí)校準(zhǔn)后的2 個(gè)通道信號(hào)時(shí)域Fig.3 Time-domain signals of two channels after amplitude error and time delay calibration

圖4 經(jīng)過(guò)幅度延時(shí)校準(zhǔn)后的2 個(gè)通道信號(hào)包絡(luò)Fig.4 Envelope diagram of signals after amplitude error and time delay calibration

由圖3～圖4 中可見(jiàn)，經(jīng)過(guò)幅度和延時(shí)校準(zhǔn)后，兩通道信號(hào)的延時(shí)已經(jīng)完全消除，不同發(fā)射通道的輸出信號(hào)在時(shí)域上已經(jīng)能夠完全重合。

2.3 相位差異校準(zhǔn)

經(jīng)過(guò)延時(shí)幅度校準(zhǔn)后，通道1 信號(hào)與通道2 信號(hào)的包絡(luò)完全重合，幅度延時(shí)誤差已經(jīng)完全校準(zhǔn)。由信號(hào)的時(shí)域圖可以看出，信號(hào)并不能夠重合，因?yàn)椴煌l(fā)射通道的鏈路對(duì)信號(hào)初相的影響存在一定的差異，需要求解不同信號(hào)通道差異帶來(lái)的相位誤差，并對(duì)相位誤差進(jìn)行補(bǔ)償。在本次校準(zhǔn)實(shí)驗(yàn)中，采用200 個(gè)樣本點(diǎn)計(jì)算一次相位差，得到的相位差為：-2.413 519 811 395 667e+02°，即通道2 信號(hào)與通道1 信號(hào)由于發(fā)射鏈路帶來(lái)的初相誤差為-241.351 98°。得到相位誤差后，再對(duì)信號(hào)通道的相位差異性進(jìn)行補(bǔ)償，最終的校準(zhǔn)結(jié)果如圖5 所示。

圖5 完成幅相一致性校準(zhǔn)后的2 個(gè)通道信號(hào)Fig.5 Two channel signals with completed amplitude-phase calibration

對(duì)2 個(gè)通道完成幅相誤差校準(zhǔn)后，又利用數(shù)采器對(duì)2 個(gè)通道的信號(hào)進(jìn)行采集，不進(jìn)行任何處理，直接提取校準(zhǔn)通道后的兩通道信號(hào)的幅度、延時(shí)和相位差異，并與未校準(zhǔn)前對(duì)比，對(duì)比結(jié)果見(jiàn)表1。

表1 通道校準(zhǔn)前、后提取幅相差異對(duì)比Tab.1 Comparison of amplitude and phase differences before and after channels calibration

根據(jù)校準(zhǔn)前、后提取的兩通道幅相差異對(duì)比分析可知，采用本文方法進(jìn)行相控陣?yán)走_(dá)發(fā)射通道信號(hào)校準(zhǔn)后，通道差異性大大降低，具有較高的校準(zhǔn)精度。

3 結(jié)束語(yǔ)

本文針對(duì)相控陣?yán)走_(dá)多個(gè)發(fā)射通道的幅相特性差異，以及DBF（Digital Beam Forming）接收整列多通道的幅相差異，避免時(shí)域?qū)拵盘?hào)校準(zhǔn)困難的問(wèn)題，提出了基于最小二乘法的幅相校準(zhǔn)方法。在提取相位時(shí)，利用寬帶信號(hào)的瞬時(shí)窄帶特性，采用最小二乘法實(shí)現(xiàn)相位特征參數(shù)的提取，由于相位提取時(shí)只需要少數(shù)采樣點(diǎn)參與計(jì)算，使其能夠應(yīng)用于任何形式的雷達(dá)信號(hào)。采用一塊獨(dú)立的校準(zhǔn)單元，配置相關(guān)通道數(shù)，即可完成對(duì)應(yīng)通道數(shù)目的發(fā)射鏈路通道校準(zhǔn)。在多個(gè)項(xiàng)目的陣列通道校準(zhǔn)工作中，已經(jīng)驗(yàn)證本文的校準(zhǔn)方法具有較高的信噪比容忍度，在10 dB 以下信噪比仍然具有較高的校準(zhǔn)精度。因此，該方法還可應(yīng)用于接收通道的校準(zhǔn)，目前已在部分DBF 項(xiàng)目中進(jìn)行應(yīng)用。