秦 風(fēng) 高 原吳 雙

(中國(guó)工程物理研究院應(yīng)用電子學(xué)研究所 綿陽(yáng) 621900)

(中國(guó)工程物理研究院復(fù)雜電磁環(huán)境科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室 綿陽(yáng) 621900)

1 引言

在電磁環(huán)境測(cè)試技術(shù)領(lǐng)域，經(jīng)常采用十幾米甚至幾十米的同軸電纜傳輸測(cè)試信號(hào)，對(duì)于瞬時(shí)高頻或?qū)拵}沖信號(hào)，經(jīng)長(zhǎng)距離同軸電纜傳輸時(shí)易出現(xiàn)畸變，且畸變程度隨著信號(hào)頻率、帶寬以及傳輸距離的增大而急劇增大[1–3]。此外，當(dāng)同軸電纜意外受到擠壓、拉伸或折疊時(shí)，其信號(hào)傳輸特性會(huì)發(fā)生改變，進(jìn)一步加劇信號(hào)的畸變程度，導(dǎo)致測(cè)試儀器記錄的測(cè)量信號(hào)嚴(yán)重偏離實(shí)際被測(cè)信號(hào)[4–7]。為此，需對(duì)測(cè)試儀器記錄的測(cè)量信號(hào)進(jìn)行畸變補(bǔ)償，以獲得準(zhǔn)確的被測(cè)信號(hào)，提高測(cè)試精度。

同軸電纜信號(hào)傳輸畸變補(bǔ)償?shù)谋举|(zhì)是通過(guò)測(cè)試儀器記錄的測(cè)量信號(hào)(即輸出信號(hào))來(lái)完成實(shí)際被測(cè)信號(hào)(輸入信號(hào))的重構(gòu)；目前，主要有3類方法：衰減補(bǔ)償法[8,9]、均衡網(wǎng)絡(luò)法[10,11]和反卷積法[12,13]。衰減補(bǔ)償法利用同軸電纜的衰減常數(shù)重構(gòu)輸入信號(hào)，方法簡(jiǎn)單、應(yīng)用廣泛；該方法中通常假定衰減常數(shù)與頻率平方根成正比，這種近似精度有限。均衡網(wǎng)絡(luò)法根據(jù)同軸電纜的傳輸特性設(shè)計(jì)補(bǔ)償網(wǎng)絡(luò)，以均衡同軸電纜對(duì)不同頻率信號(hào)的衰減；而當(dāng)同軸電纜改變時(shí)，補(bǔ)償網(wǎng)絡(luò)需重新設(shè)計(jì)，重復(fù)使用性低。此外，對(duì)于高頻信號(hào)，補(bǔ)償網(wǎng)絡(luò)的分布參數(shù)會(huì)引入新的畸變因素，影響補(bǔ)償效果。反卷積法通過(guò)同軸電纜的脈沖響應(yīng)函數(shù)與輸出信號(hào)的反卷積運(yùn)算完成輸入信號(hào)的重構(gòu)，雖然計(jì)算模型能夠較好呈現(xiàn)電纜的實(shí)際工作狀態(tài)，但在計(jì)算過(guò)程中易出現(xiàn)病態(tài)矩陣，導(dǎo)致重構(gòu)出的輸入信號(hào)嚴(yán)重偏離實(shí)際信號(hào)。而在圖像重建[14]和沖擊力分析[15]領(lǐng)域，有一種應(yīng)用廣泛的逆問(wèn)題分析方法——基于貝葉斯推理的非負(fù)Tikhonov正則化方法。該方法基于系統(tǒng)真實(shí)脈沖響應(yīng)，在貝葉斯推斷框架下進(jìn)行未知輸入信號(hào)的重構(gòu)，使得計(jì)算模型貼近系統(tǒng)實(shí)際工作狀態(tài)；并且，采用Tikhonov正則化方法處理不適定問(wèn)題，能夠有效解決病態(tài)矩陣計(jì)算過(guò)程中的不穩(wěn)定問(wèn)題[16,17]。該方法為同軸電纜信號(hào)補(bǔ)償提供了一種新的思路，然而該方法對(duì)未知輸入信號(hào)有非負(fù)限制，因而不能直接應(yīng)用于補(bǔ)償含有負(fù)分量的電磁信號(hào)。

因此，本文在貝葉斯推理的非負(fù)Tikhonov正則化方法的基礎(chǔ)上，提出一種改進(jìn)的信號(hào)補(bǔ)償方法，本方法可完成任意值信號(hào)的分析重構(gòu)，從而可應(yīng)用于含有負(fù)值電磁脈沖信號(hào)的同軸電纜傳輸畸變補(bǔ)償。在此基礎(chǔ)上，以長(zhǎng)度15 m的受擠壓同軸電纜為對(duì)象，采用此方法完成了對(duì)雙指數(shù)脈沖信號(hào)、調(diào)制方波信號(hào)、雙極脈沖信號(hào)等3種不同樣式脈沖信號(hào)的傳輸畸變補(bǔ)償。結(jié)果表明：本方法補(bǔ)償效果優(yōu)異，相比于傳統(tǒng)的衰減補(bǔ)償法，補(bǔ)償后信號(hào)與輸入信號(hào)間偏差大幅減小。并且，改進(jìn)的信號(hào)補(bǔ)償方法具有較好的魯棒性。

2 同軸電纜信號(hào)傳輸畸變補(bǔ)償方法

2.1 同軸電纜信號(hào)傳輸畸變問(wèn)題

由于同軸電纜的損耗、色散特性，輸出信號(hào)往往會(huì)產(chǎn)生畸變，偏離于輸入信號(hào)；并且隨著信號(hào)頻率、頻帶寬度以及同軸電纜長(zhǎng)度的增加，輸出信號(hào)會(huì)更加偏離輸入信號(hào)。根據(jù)同軸電纜黑箱模型可知，同軸電纜的輸出信號(hào)y(t)是輸入信號(hào)x(t)和電纜脈沖沖擊響應(yīng)h(t)的卷積，即

此處，xi為輸入信號(hào)離散點(diǎn)，yi為輸出信號(hào)離散點(diǎn)，hi為同軸電纜脈沖響應(yīng)離散點(diǎn)，N為信號(hào)離散點(diǎn)數(shù)，即采樣點(diǎn)數(shù)。從式(2)看，當(dāng)輸出信號(hào)y和脈沖響應(yīng)H已知時(shí)，通過(guò)矩陣直接求逆的方法得到輸入信號(hào)x似乎很容易。實(shí)際上，測(cè)量數(shù)據(jù)y通常是含有噪聲的，即y=H x+ω,ω為隨機(jī)噪聲，而矩陣H經(jīng)常是病態(tài)的，導(dǎo)致直接求逆的方法對(duì)噪聲很敏感，很小的噪聲ω也會(huì)導(dǎo)致解出現(xiàn)較大的振蕩，嚴(yán)重時(shí)會(huì)完全覆蓋真實(shí)信息。因此，很難通過(guò)直接求解式(2)逆問(wèn)題的方式獲得準(zhǔn)確的輸入信號(hào)。

2.2 改進(jìn)的基于貝葉斯推理的非負(fù)Tikhonov正則化方法

針對(duì)式(2)所示的不適定逆問(wèn)題，根據(jù)經(jīng)典Tikhonov正則化理論，其求解可轉(zhuǎn)換為含有正則化項(xiàng)的最優(yōu)問(wèn)題

由式(6)可知，正則化參數(shù)δ的選取是求解x的關(guān)鍵，常用方法有廣義交叉驗(yàn)證(Generalized Cross-Validation,GCV)法[18]、L曲線法[19]，但是在δ解值曲線過(guò)于平滑時(shí)，這些方法難以收斂。文獻(xiàn)[20]提出了一種基于貝葉斯推斷的非負(fù)正則化方法，該方法可自動(dòng)選取正則化參數(shù)，同時(shí)給出逆問(wèn)題的正則化解，并且求解過(guò)程的收斂性和解的唯一性得到了嚴(yán)格證明。但是該方法只適用于非負(fù)信號(hào)(即xi≥0)的重構(gòu)分析，難以用于含有負(fù)值電磁脈沖環(huán)境信號(hào)的重構(gòu)或畸變補(bǔ)償；為此，本文提出一種

其中，修正后的未知輸入信號(hào)n滿足ni≥0。

進(jìn)而可以采用基于貝葉斯推理的非負(fù)Tikhonov正則化方法來(lái)求解式(14)，得到未知輸入信號(hào)n的最優(yōu)估計(jì)值n?。根據(jù)文獻(xiàn)[20]提出的貝葉斯統(tǒng)計(jì)推斷方法，未知參數(shù)n的求解可轉(zhuǎn)換為以下最大后驗(yàn)概率估計(jì)的最優(yōu)化問(wèn)題

將式(16)—式(20)代入式(15)中，得到

從而將式(15)所示的最優(yōu)化問(wèn)題等價(jià)于求解如式(22)泛函的極小值

式(22)即為帶非負(fù)限制的增廣Tikhonov正則化泛函，而求解該泛函的極小值，相當(dāng)于求解其各項(xiàng)偏導(dǎo)數(shù)為0，即?x即為基于改進(jìn)補(bǔ)償方法計(jì)算得到的輸入信號(hào)估計(jì)值(即補(bǔ)償信號(hào))。

2.3 補(bǔ)償效果評(píng)價(jià)方法

本文采用補(bǔ)償信號(hào)或輸出信號(hào)與輸入信號(hào)之間的相對(duì)偏差RE來(lái)評(píng)價(jià)信號(hào)畸變補(bǔ)償方法的補(bǔ)償效果，具體計(jì)算方式為

3 實(shí)驗(yàn)與討論

3.1 實(shí)驗(yàn)設(shè)置

為驗(yàn)證改進(jìn)補(bǔ)償方法的有效性和適用性，以受擠壓的15 m同軸電纜為對(duì)象(圖2)，開(kāi)展同軸電纜傳輸畸變補(bǔ)償實(shí)驗(yàn)。實(shí)驗(yàn)設(shè)置如圖3所示，主要由任意波形發(fā)生器、功分器、0.5 m同軸電纜和受擠壓的15 m同軸電纜以及數(shù)字示波器組成。其中，通過(guò)受擠壓的15 m同軸電纜兩個(gè)端口的信號(hào)分別為輸入信號(hào)x(t)和輸出信號(hào)y(t)，輸入信號(hào)和輸出信號(hào)由數(shù)字示波器采集得到；且為了方便，采用長(zhǎng)度0.5 m、損耗極低的穩(wěn)幅穩(wěn)相同軸電纜實(shí)現(xiàn)受擠壓的15 m同軸電纜輸入信號(hào)x(t)的采集。此外，由于同軸電纜的脈沖響應(yīng)難以直接測(cè)量；因此，在本文中，首先通過(guò)矢量網(wǎng)絡(luò)分析儀測(cè)量同軸電纜的S21參數(shù)(見(jiàn)圖4(a))，其中：測(cè)量頻率范圍為10 kHz～5 GHz、頻率間隔為250 k Hz；然后進(jìn)行快速傅里葉逆變換(Inverse Fast Fourier Transform,IFFT)，獲得脈沖沖擊響應(yīng)h(t)(圖4(b))。從圖4可以看出，受試同軸電纜的S21曲線抖動(dòng)明顯，其原因是擠壓變形改變了同軸電纜的特性阻抗，從而造成了S21參數(shù)的顯著變化；并且，與正常同軸電纜的脈沖沖擊響應(yīng)相比，受試同軸電纜的脈沖沖擊響應(yīng)有明顯的毛刺和振蕩，這將會(huì)對(duì)傳輸信號(hào)造成更為嚴(yán)重的畸變。需要注意的是，在測(cè)量輸入輸出信號(hào)和S21參數(shù)過(guò)程中，數(shù)字示波器、矢量網(wǎng)絡(luò)分析儀與同軸電纜之間始終保持阻抗匹配，以消除阻抗失配對(duì)實(shí)驗(yàn)結(jié)果的影響。

圖1 改進(jìn)方法的算法流程

圖2 受擠壓的15 m被測(cè)同軸電纜

圖3 實(shí)驗(yàn)測(cè)試設(shè)置示意圖

圖4 同軸電纜脈沖沖擊響應(yīng)計(jì)算

3.2 改進(jìn)補(bǔ)償方法對(duì)不同類型信號(hào)的補(bǔ)償

基于上述改進(jìn)后的非負(fù)Tikhonov正則化方法，首先對(duì)上升時(shí)間2.5 ns、脈沖半高寬23 ns的類核電磁脈沖雙指數(shù)脈沖信號(hào)(輸入信號(hào))在15 m受擠壓同軸電纜中傳輸、畸變后的信號(hào)(輸出信號(hào))進(jìn)行補(bǔ)償。同軸電纜的輸入信號(hào)、輸出信號(hào)以及采用改進(jìn)方法計(jì)算得到的補(bǔ)償信號(hào)如圖5所示，其中：采樣時(shí)間間隔為0.2 ns，補(bǔ)償過(guò)程中采用的超參數(shù)α1,α2,α3,β1,β2,β3和收斂判別參數(shù)ε分別設(shè)置為3,1,1,2×10–10,150,1×10–4和1×10–5，偏置信號(hào)d(t)采用數(shù)值為1的常量信號(hào)。15 m受擠壓同軸電纜對(duì)雙指數(shù)脈沖信號(hào)傳輸影響較大，從圖5可以明顯看出，輸出信號(hào)與輸入信號(hào)之間存在明顯差異，輸出信號(hào)與輸入信號(hào)之間的相對(duì)偏差達(dá)到16.5%；而采用改進(jìn)補(bǔ)償方法補(bǔ)償后的信號(hào)與輸入信號(hào)吻合較好、基本一致，補(bǔ)償信號(hào)與輸入信號(hào)間的相對(duì)偏差僅為1.5%(表1)；表明該改進(jìn)后的補(bǔ)償方法能夠很好地補(bǔ)償輸出信號(hào)。

圖5 雙指數(shù)脈沖信號(hào)在受擠壓15 m同軸電纜中的傳輸畸變補(bǔ)償

進(jìn)一步地，采用此改進(jìn)后的補(bǔ)償方法，分別對(duì)調(diào)制方波信號(hào)(類窄帶高功率微波信號(hào)、脈寬15 ns、中心頻率2 GHz、信號(hào)頻率范圍1940～2060 MHz)、雙極脈沖信號(hào)(類寬帶電磁脈沖信號(hào)、脈寬3 ns、信號(hào)頻率范圍33～705 MHz)經(jīng)15 m受擠壓同軸電纜傳輸后的輸出信號(hào)進(jìn)行補(bǔ)償計(jì)算，結(jié)果分別如圖6(a)、圖6(b)所示；補(bǔ)償過(guò)程中，調(diào)制方波信號(hào)、雙極脈沖信號(hào)的采樣時(shí)間間隔分別為0.025 ns,0.2 ns，所采用的超參數(shù)和收斂判別參數(shù)與上述雙指數(shù)脈沖信號(hào)相同，所采用的偏置信號(hào)分別設(shè)置為數(shù)值為2.0,2.0的常量信號(hào)。從圖6可以看出，對(duì)于頻率較高的調(diào)制方波信號(hào)、雙極脈沖信號(hào)，輸出信號(hào)與輸入信號(hào)之間也存在明顯的差異，其相對(duì)偏差分別為24.4%,25.3%(表1)，大于頻率較低的類核電磁脈沖雙指數(shù)脈沖信號(hào)，更是難以直接通過(guò)輸出信號(hào)準(zhǔn)確描述輸入信號(hào)。而采用改進(jìn)補(bǔ)償方法對(duì)輸出信號(hào)進(jìn)行補(bǔ)償，補(bǔ)償信號(hào)與輸入信號(hào)幾乎相同，相對(duì)偏差大幅減小(分別為2.4%,2.3%)，此實(shí)驗(yàn)結(jié)果進(jìn)一步證實(shí)該改進(jìn)補(bǔ)償方法具有良好的適用性。然而，隨著信號(hào)頻率的增大，該改進(jìn)補(bǔ)償方法的補(bǔ)償精度還是有稍許的下降。

表1 輸出信號(hào)、采用改進(jìn)補(bǔ)償方法補(bǔ)償后信號(hào)與輸入信號(hào)間相對(duì)偏差(%)

圖6 改進(jìn)補(bǔ)償法對(duì)于不同類型信號(hào)的畸變補(bǔ)償

3.3 改進(jìn)補(bǔ)償方法與衰減補(bǔ)償法的補(bǔ)償結(jié)果對(duì)比

同樣地，以上述3種類型的信號(hào)作為同軸電纜的輸入信號(hào)，對(duì)改進(jìn)補(bǔ)償方法和常用的衰減補(bǔ)償法的補(bǔ)償效果進(jìn)行對(duì)比研究。同軸電纜的輸入信號(hào)、采用改進(jìn)補(bǔ)償方法與衰減補(bǔ)償法計(jì)算得到的補(bǔ)償后信號(hào)如圖7所示；從中可以看出，改進(jìn)補(bǔ)償方法補(bǔ)償后信號(hào)與輸入信號(hào)基本一致，衰減補(bǔ)償法雖然也能大幅減小輸出信號(hào)與輸入信號(hào)的差別，但相比于改進(jìn)補(bǔ)償方法，補(bǔ)償效果還是存在明顯差距。

圖7 改進(jìn)補(bǔ)償法與衰減補(bǔ)償法的補(bǔ)償效果對(duì)比

同時(shí)，針對(duì)3種不同類型的信號(hào)，還分別計(jì)算了衰減法補(bǔ)償信號(hào)與輸入信號(hào)、改進(jìn)方法補(bǔ)償信號(hào)與輸入信號(hào)之間的相對(duì)偏差，結(jié)果如表2所示。當(dāng)采用衰減補(bǔ)償法開(kāi)展信號(hào)補(bǔ)償時(shí)，此3種不同類型信號(hào)的補(bǔ)償信號(hào)與輸入信號(hào)的相對(duì)偏差從9.0%～13.3%；相比之下，采用改進(jìn)補(bǔ)償法補(bǔ)償后信號(hào)與輸入信號(hào)的最大相對(duì)偏差不超過(guò)2.4%，遠(yuǎn)低于衰減補(bǔ)償法。實(shí)驗(yàn)結(jié)果表明：相比于常用的衰減補(bǔ)償法，提出的改進(jìn)補(bǔ)償方法具有更好的補(bǔ)償效果。由于受試同軸電纜S21參數(shù)抖動(dòng)明顯，補(bǔ)償過(guò)程中近似計(jì)算得到的衰減常數(shù)與實(shí)際衰減值偏差較大，進(jìn)而降低了衰減補(bǔ)償法補(bǔ)償信號(hào)的精度；而對(duì)于改進(jìn)補(bǔ)償方法，在補(bǔ)償過(guò)程中，未對(duì)受試同軸電纜的S21參數(shù)進(jìn)行近似處理，同時(shí)還考慮了測(cè)量數(shù)據(jù)(同軸電纜輸出信號(hào))中的噪聲影響，并且整個(gè)計(jì)算過(guò)程在貝葉斯非負(fù)正則化處理下穩(wěn)定收斂，確保了信號(hào)畸變補(bǔ)償精度。

表2 不同補(bǔ)償方法補(bǔ)償后信號(hào)與輸入信號(hào)間相對(duì)偏差(%)

3.4 不同噪聲下補(bǔ)償結(jié)果對(duì)比

為了研究改進(jìn)補(bǔ)償方法的魯棒性，本文在3種不同輸出信號(hào)y(t)中添加高斯白噪聲(白噪聲根據(jù)給定的信噪比(Signal Noise Ratio,SNR)隨機(jī)生成)，并分別采用改進(jìn)補(bǔ)償方法和衰減補(bǔ)償法進(jìn)行信號(hào)畸變補(bǔ)償計(jì)算。不同信噪比情況下，兩種補(bǔ)償方法補(bǔ)償后信號(hào)與輸入信號(hào)之間的相對(duì)偏差計(jì)算結(jié)果分別如表3、表4所示?？梢钥闯觯捎诟倪M(jìn)補(bǔ)償方法在補(bǔ)償過(guò)程中考慮了輸出信號(hào)的噪聲水平，因此在相同信噪比下，改進(jìn)補(bǔ)償方法的補(bǔ)償效果大幅優(yōu)于衰減補(bǔ)償法；并且在信噪比大于30 d B時(shí)，改進(jìn)補(bǔ)償方法可保持良好的穩(wěn)定性，表明其具有較強(qiáng)的魯棒性。當(dāng)信噪比小于30 d B時(shí)，由于輸出信號(hào)受噪聲影響較大，測(cè)量的信號(hào)質(zhì)量明顯下降，導(dǎo)致補(bǔ)償精度隨之降低。

表3 不同信噪比下改進(jìn)補(bǔ)償方法補(bǔ)償后信號(hào)與輸入信號(hào)間相對(duì)偏差(%)

表4 不同信噪比下衰減補(bǔ)償法補(bǔ)償后信號(hào)與輸入信號(hào)間相對(duì)偏差(%)

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

本文針對(duì)同軸電纜信號(hào)傳輸畸變問(wèn)題，基于貝葉斯推理的非負(fù)Tikhonov正則化方法，提出一種改進(jìn)的同軸電纜信號(hào)傳輸畸變補(bǔ)償方法。并以受擠壓的15 m同軸電纜為對(duì)象，完成了3種不同類型信號(hào)的同軸電纜傳輸畸變補(bǔ)償實(shí)驗(yàn)，并將改進(jìn)方法的補(bǔ)償結(jié)果與傳統(tǒng)衰減補(bǔ)償法的補(bǔ)償結(jié)果進(jìn)行了比對(duì)。結(jié)果表明，改進(jìn)補(bǔ)償方法具有優(yōu)異的補(bǔ)償效果、較強(qiáng)的魯棒性及良好的普適性。所提改進(jìn)補(bǔ)償方法僅需利用系統(tǒng)的脈沖響應(yīng)函數(shù)，整個(gè)計(jì)算過(guò)程無(wú)特別限制，本方法可推廣應(yīng)用于其他系統(tǒng)(如天線、濾波器等)的信號(hào)傳輸畸變補(bǔ)償或輸入信號(hào)重構(gòu)。