張吉偉 時杰飛 詹永衛(wèi) 張超

【摘要】? ? 本文針對溫度變化影響頻譜分析儀中的多級放大、濾波等模塊的頻響，從而導(dǎo)致其幅度測量精度下降的問題，提出了一種基于最小二乘法的分段線性擬合頻響補償權(quán)值曲線的方法;并根據(jù)頻譜分析儀的實際工作溫度變化曲線模型，提出了一種與開機時間相關(guān)的溫度線性空間的分段劃分方法;同時給出了此方法在頻譜分析儀中的應(yīng)用補償流程。通過仿真分析，該方法有效的使頻響補償值跟蹤溫度的變化，提高了幅度測量精度，解決了分析儀在開機預(yù)熱階段的幅度偏差大以及測量連續(xù)性問題。

【關(guān)鍵詞】? ? 最小二乘法? ? 分段擬合? ? 溫度? ? 頻響補償

A method of compensation about spectrum analyzers frequency response effected by temperature drift based on least squares

Zhang Jiwei? ?Shi Jiefei? ?Zhan Yongwei? ?Zhang Chao

（The 41st? Research? Insitute of? CETC，? Qingdao? 266555，China ）

Abstract： This paper proposes a method of piecewise linear fitting of the frequency response compensation weight curve of a spectrum analyzer based on the least squares， which focuses on the influence of temperature changes on the frequency response of multi-stage amplification and filtering modules in spectrum analyzer. A segmentation method of temperature linear space related to start-up time is proposed according to the actual working temperature curve model of the spectrum analyzer. At the same time， the application compensation process of this method is given. Through simulation analysis， this method effectively makes the frequency response compensation value track the change of temperature， improves the accuracy of amplitude measurement， and solves the problem of large amplitude deviation in the preheating phase of the analyzer and the continuity of high-precision amplitude measurement.

Key words： least squares method? ? Linear fitting method? ? temperature? ? frequency response

引言：

在頻譜分析儀中，幅度精度是一個非常重要的技術(shù)指標，其在儀器測試中體現(xiàn)在信號的頻率響（簡稱頻響）應(yīng)測量誤差上[1]，目前國外的頻譜儀，像是德科技的N9040B系列頻譜分析儀到50GHz頻段的頻率響應(yīng)在20到30℃的溫度范圍內(nèi)誤差在±3.2dB以內(nèi);R&S公司的FSVA3000信號和頻譜分析儀44GHz頻段的頻響在20到30℃條件下誤差在±1.97dB以內(nèi)，在0到50℃條件下誤差在±4dB以內(nèi)。信號經(jīng)過各級通路在不同的溫度環(huán)境下其頻率響應(yīng)誤差不同，即產(chǎn)生溫度漂移（簡稱溫漂）[2-3]。

最小二乘法廣泛應(yīng)用于系統(tǒng)優(yōu)化、最優(yōu)化求解等領(lǐng)域的成熟的統(tǒng)計數(shù)據(jù)分析方法[4-6]。利用最小二乘法訓練學習溫漂的補償權(quán)值，擬合補償權(quán)值曲線，使得頻譜分析儀獲得高的幅度測量精度。

一、信號和頻譜分析儀頻響補償

頻譜分析儀頻率響應(yīng)補償方法如圖1所示，信號源產(chǎn)生不同頻率的信號，經(jīng)過功分器功分兩路分別給頻譜分析儀和功率計，通過計算測量功率計與頻譜分析儀的差值即為當前頻率下的補償值。

在大多數(shù)的頻譜分析儀中，僅僅對頻響補償一次，在所有的測試條件下用這一次的補償值。頻譜分析儀在剛開機階段，其內(nèi)部溫度變化較大，頻率響應(yīng)也變化較大。為了解決開機的溫度變化影響，目前常用的補償方法是待頻譜分析儀開機預(yù)熱一段時間，使其工作在溫度相對穩(wěn)定的條件下再進行頻響補償，如圖2所示。但是，頻譜儀內(nèi)部模塊較多、電路復(fù)雜，所需要的預(yù)熱時間較長，通常需要20分鐘以上，且在補償之后，每次開機預(yù)熱期間頻響誤差大，是不進行幅度測試的;另外頻譜分析儀的工作環(huán)境復(fù)雜多樣，在其補償時差別較大時，其幅度測量誤差也會增大。因此，采用開機預(yù)熱穩(wěn)定后測試存在很多弊端的。

為了克服因溫度環(huán)境變化較大帶來大的幅度測量誤差，一種方法是在頻譜儀內(nèi)部設(shè)計標準的校準源信號。由于頻譜儀的頻段很寬，校準信號數(shù)量有限，僅在校準的頻點幅度準確度高，并且內(nèi)部校準源也會受溫度的影響，同時也增加了成本。另外定時校準時需要停止測量，對用戶使用體驗不佳，在連續(xù)測試的需求下這種校準模式也是不可用的。

二、基于最小二乘的溫度補償方法

2.1 最小二乘法的溫度補償模型

受溫漂的影響，單次補償以及引入校準信號，無法實現(xiàn)連續(xù)測試條件下的高精度幅度測量。因此，獲取寬頻帶頻譜分析儀的頻響溫度特性，并自動進行實時補償，有效的解決溫漂導(dǎo)致大的幅度測量誤差。利用最小二乘法將權(quán)系數(shù)曲線分段擬合，在每一段內(nèi)近似成線性，可獲得高精度擬合曲線。

在每個受溫度影響的模擬前端模塊電路上引入溫度傳感器，實時監(jiān)測各個模塊的溫度，假設(shè)有M個模塊受溫度的影響會影響到頻響參數(shù)，如圖3所示

在每個分段內(nèi)近頻響變化與溫度似成線性，關(guān)于擬合權(quán)系數(shù)a?？色@得方程組如下式所示：

系數(shù)矩陣T為：

系數(shù)矩陣T，以及補償值p為頻響補償時的訓練數(shù)據(jù)。在進行數(shù)據(jù)補償過程中，獲取多組訓練數(shù)據(jù)，即進行多次頻響補償，使得N>M。

在最小二乘意義下，使得誤差的模的平方和

取得最小值，方程的唯一解為：

2.2 溫度區(qū)間的分段劃分

由于整機和模塊的散熱設(shè)計，每此個模塊升溫后的工作溫度是相對穩(wěn)定的。由于信號和頻譜分析儀從環(huán)境溫度到已預(yù)熱階段，溫度傳感器1的升溫變化區(qū)間[T11，TN1]到溫度傳感器M的升溫變化區(qū)間[TM1，TNM]的升溫趨勢相似，如圖4所示，所以溫度線性空間的分段劃分可以近似在時間段上劃分[8-10]。

雖然同一臺頻譜儀內(nèi)部的不同傳感器的溫度變化趨勢相似，但是不同的環(huán)境溫度[T11，TN1]與[TM1，TNM]的劃分區(qū)間的對應(yīng)關(guān)系不同。

為了解決不同環(huán)境起始溫度帶來的分段差異性，使的分段劃分更合理，對信號分析儀進行溫度變化曲線統(tǒng)計學習，將溫度補償系數(shù)分為多組，使其更加逼近線性。同時學習到穩(wěn)定階段的時間，對穩(wěn)定階段學習訓練盡量多的數(shù)據(jù)，從而提高穩(wěn)定階段的精度。

2.3 溫漂補償流程

基于最小二乘法的信號和頻譜分析儀頻率響應(yīng)溫漂補償方法分為三個階段：

1.溫度曲線學習階段，頻譜儀放在環(huán)境試驗箱里，設(shè)置不同的環(huán)境溫度，學習不同的環(huán)境溫度下的各模塊溫度傳感器的變化曲線，并劃分不同的線性段;

2.訓練補償階段，頻譜儀在環(huán)境試驗箱里，在不同的環(huán)境溫度下觀測溫度值與頻響補償值，根據(jù)階段1將補償數(shù)據(jù)分組，利用最小二乘法擬合溫度與頻響曲線;

3.用戶使用測量階段，實時獲取各模塊的溫度值，并匹配所處于的線性段，得到當前狀態(tài)下的頻響補償值。

三、建模仿真分析

根據(jù)對頻譜儀的頻響溫漂特性以及溫度試驗箱的實測數(shù)據(jù)，頻譜儀一般的工作環(huán)境溫度為0℃到50℃范圍內(nèi)，同一環(huán)境溫度不同的預(yù)熱時間段信號的頻響差異達±6dB，不同的環(huán)境溫度同一信號的頻響差異達±4dB。據(jù)此進行建模仿真。

假設(shè)在信號分析儀中有三個模塊的頻響受溫度影響。在環(huán)境溫度為18℃的環(huán)境下的溫度變化曲線如圖5所示。將曲線劃分為4個線性擬合區(qū)間，其中分段1、2、3為預(yù)熱階段，分段4為穩(wěn)定階段，即在此環(huán)境溫度下18min之后，頻譜分析儀的溫度擬合系數(shù)到穩(wěn)定階段。

對環(huán)境溫度為18℃和34℃的條件下的溫度與頻響進行擬合并與仿真補償值對比。如圖6所示，常用的單次補償方法受溫度影響誤差很大，在預(yù)熱階段更明顯;基于最小二乘法的頻響溫漂擬合方法可以有效的使補償值跟蹤溫度變化，且誤差降到很小。

四、結(jié)束語

針對頻譜分析儀的頻響受環(huán)境溫度以及工作狀態(tài)溫度的影響大，提出了引入溫度傳感器，通過最小二乘法分段線性擬合的方法學習訓練頻響補償值，并根據(jù)針對頻譜分析儀的溫度變化特性，提出了溫度分段劃分方法，降低了在多元分段線性擬合的線性區(qū)間的排列組合的多樣性，通過仿真分析，該算法可以有效的使頻響補償值跟蹤溫度的變化，且誤差也小，解決了頻譜分析儀在預(yù)熱階段的測試問題，且該算法不用間歇校準，保證了頻譜儀使用的連續(xù)性。

參? 考? 文? 獻

[1] 李凡，莊濤，徐宏光. 頻譜分析儀測量電平的不確定度分析與評定 [J]. 品牌與標準化. 2016 （03）

[2] 李劍雄.頻譜分析儀與測量技術(shù)基礎(chǔ) [M]. 北京人民郵電出版社. 2011：13-16.

[3] 卞劍. 頻譜分析儀自動校準中的問題及解決方式 [J]. 計測技術(shù). 2016 （01）

[4] 孟紅波. 一種基于最小二乘法和RBF神經(jīng)網(wǎng)絡(luò)的多步誤差補償方法? [J]. 測試技術(shù)學報. 2015 （06）

[5] 于湘濤. 基于小波最小二乘支持向量機的加速度計溫度建模和補償? [J]. 中國慣性技術(shù)學報. 2011 （01）

[6] 楊賓峰. 基于最小二乘的地磁場測量誤差補償技術(shù) [J]. 空軍工程大學學報. 2017 （06）

[9] 郭星星. 基于循環(huán)神經(jīng)網(wǎng)絡(luò)的多核芯片熱管理技術(shù)研究 [D]. 電子科技大學. 2019