摘 要： 在航天測控系統(tǒng)中，天線系統(tǒng)由電動機、變速箱等機電設備組成。為了使天線快速、安全、準確地對準目標，必須對其進行精確控制。通常情況下，由于天線系統(tǒng)是一個時變、非線性、復雜的動力學系統(tǒng)，其參數(shù)難以精確測量，在設計控制器時難以達到最優(yōu)。利用系統(tǒng)辨識方法，通過天線系統(tǒng)的輸入、輸出數(shù)據(jù)，利用Matlab軟件，對天線模型進行離線辨識，并對辨識出的模型進行了驗證，證明模型能夠準確地表現(xiàn)天線系統(tǒng)，滿足應用要求。

關鍵詞： 系統(tǒng)辨識； 天線； 雷達； Matlab

中圖分類號： TN82?34； TP273 文獻標識碼： A 文章編號： 1004?373X（2016）02?0104?03

Study and realization of identification for aerospace TTC antenna system

DU Peiming， ZHANG Yongsheng

（Unit 63717 of PLA， Taiyuan 030031， China）

Abstract： In aerospace TTC （tracing， telemetry and command） system， the antenna system is composed of electric motor， gearbox and other electromechanical devices. In order to make the antenna point to the target quickly， safely and accurately， it is necessary to control the antenna precisely. In general， because the antenna is a nonlinear， time-variant and complex dynamics system， it is difficult to measure the parameters accurately， and the controller design is hard to reach the optimum. The system identification method， input and output data of the antenna system， and Matlab software are used to conduct off-line identification for the antenna model. The identified model was verified to prove that the model can accurately express the antenna system， and meet the application requirements.

Keywords： system identification； antenna； radar； Matlab

0 引 言

在航天外彈道測量和遙測系統(tǒng)中，天伺饋分系統(tǒng)承擔著信號收發(fā)的任務，所以必須保證天線時刻準確跟蹤飛行目標，跟蹤性能的優(yōu)劣決定于伺服系統(tǒng)性能。在天伺饋分系統(tǒng)的機械結構確定后，其性能由控制方法以及設計參數(shù)決定。通常情況下，控制方法采用三環(huán)（電流環(huán)、速度環(huán)和位置環(huán)）PID校正方法，所以可變的部分只有設計參數(shù)。如果能夠獲取天線參數(shù)，則可以根據(jù)數(shù)學算法推導出所需的控制參數(shù)。但是，一般情況下，天線參數(shù)難以準確測量，采用近似或者試湊法確定控制參數(shù)，導致伺服系統(tǒng)性能難以達到最佳。

1 系統(tǒng)辨識原理

系統(tǒng)辨識建模是一種試驗統(tǒng)計的方法[1]，所辨識的模型是與實際過程等價的一種近似的描述。

1.1 系統(tǒng)辨識定義

系統(tǒng)辨識是指在輸入、輸出數(shù)據(jù)的基礎上，從一組候選的模型中，確定一個與所研究的系統(tǒng)等價的模型[2]。通常情況下，系統(tǒng)參數(shù)是已知的，給定輸入將產(chǎn)生期望的輸出。系統(tǒng)辨識是一個逆過程，系統(tǒng)參數(shù)未知，通過輸入、輸出數(shù)據(jù)估計系統(tǒng)參數(shù)。

1.2 系統(tǒng)辨識原理

系統(tǒng)辨識原理框圖如圖1所示。通過實驗過程獲取系統(tǒng)的輸入和輸出數(shù)據(jù)，利用獲取的輸入和輸出數(shù)據(jù)進行系統(tǒng)辨識，辨識過程是一個不斷迭代的過程。將輸入數(shù)據(jù)輸入到識別出的模型中，模型將產(chǎn)生輸出，并與實際系統(tǒng)輸出數(shù)據(jù)進行比較，比較結果通過一定的算法進行處理，也就是對模型進行評價[3]，計算出模型參數(shù)修正結果，并對模型參數(shù)進行修改，此過程不斷循環(huán)，直到模型誤差滿足要求為止。對模型的評價準則算法一般采用均方差算法。

圖1中：[uk]為系統(tǒng)輸入；[yk]為理想條件下系統(tǒng)輸出；[ek]為噪聲；[zk]為包括噪聲的系統(tǒng)輸出；[zk]為辨識模型輸出；[z（k）]為實際輸出與模型輸出的誤差。

圖1 系統(tǒng)辨識原理框圖

1.3 系統(tǒng)辨識對輸入的要求

在辨識時間內(nèi)系統(tǒng)必須被輸入信號持續(xù)激勵，輸入信號必須充分激勵系統(tǒng)的所有模態(tài)，輸入信號的頻譜必須覆蓋系統(tǒng)的頻譜，且其統(tǒng)計特性與統(tǒng)計時間起點無關[4]。

輸入信號的幅度要合適，不能過大或過小，過大將使系統(tǒng)出現(xiàn)嚴重非線性問題，過小又無法使系統(tǒng)充分激勵。系統(tǒng)辨識最理想輸入是白噪聲，但白噪聲工程上難以實現(xiàn)，且系統(tǒng)也無法對白噪聲做出合適的響應，所以，M序列信號得到了廣泛應用。M序列是一種最長偽隨機序列，能夠最大限度地接近白噪聲性能，其表達式如下：

[xi=a1xi-1⊕a2xi-2⊕…⊕apxi-p] （1）

合理地選擇系數(shù)[a1，a2，…，ap]，就可以產(chǎn)生周期為[2p-1]的最長M序列。需要注意的是避免出現(xiàn)全零序列，否則序列將進入死循環(huán)。

1.4 采樣周期的選擇

采樣周期的選擇主要決定因素有：

（1） 滿足香能采樣定理[5]，采樣頻率必須大于系統(tǒng)帶寬的2倍，采樣周期不能太小，否則系統(tǒng)高頻性能將無法得到辨識。

（2） 采樣頻率也不能太高，否則將出現(xiàn)的問題有：計算量大幅增加，計算效率低下；容易形成病態(tài)方程組，難以準確計算出模型參數(shù)；容易受到噪聲干擾，進而影響辨識精度。

（3） 要充分考慮被辨識對象的性能，做到采樣信號周期與被辨識對象時間常數(shù)和M序列輸入信號周期相一致。

1.5 系統(tǒng)辨識算法

系統(tǒng)辨識最基本算法是最小二乘算法[6]，其計算公式如下：

[θ=ΦTΦ-1ΦTY] （2）

式中：[θ=a1，…，ana，b0，…，bnbT]，為待估參數(shù)；[Y=y1，y2，…，yL]為實測模型輸出，[Φ=?T1 ?T2 … ?TL]為輸出和輸入聯(lián)合矩陣。

2 系統(tǒng)辨識應用研究

2.1 傳輸函數(shù)推導

在某伺服系統(tǒng)中，采用了三環(huán)控制方法。由于速度環(huán)和電流環(huán)的帶寬較位置環(huán)大很多，伺服系統(tǒng)的性能主要決定于位置環(huán)參數(shù)的設計[7]。在系統(tǒng)辨識時，將速度環(huán)和電流環(huán)與天線作為一個整體。由于從位置環(huán)輸入更方便和可靠，辨識的系統(tǒng)模型包括了位置環(huán)部分，是已知部分。某伺服系統(tǒng)框圖如圖2所示。

圖2 某伺服系統(tǒng)原理框圖

圖2中：[G1s]為位置PID傳輸函數(shù)；[G2s]為包括速度環(huán)、電流環(huán)和負載的傳輸函數(shù)；[Fs]為反饋傳輸函數(shù)，且[Fs=1]。

整個閉環(huán)系統(tǒng)傳遞函數(shù)[8]為：

[Φs=G1sG2s1+G1sG2sFs] （3）

通過系統(tǒng)辨識可得[Φs]，[G1s]已知，可計算[G2s]，公式如下：

[G2s=ΦsG1s1-Φs] （4）

2.2 輸入、輸出數(shù)據(jù)選擇

根據(jù)某伺服系統(tǒng)結構和工作原理，輸出數(shù)據(jù)類型只有角度一種，而輸入可以是角度、速度和電流中的一種，通常情況下，系統(tǒng)辨識輸入為階躍數(shù)據(jù)，如果輸入作為電流或速度數(shù)據(jù)，將使整個天線處于運轉狀態(tài)，不利于對天線的控制和數(shù)據(jù)的獲取，同時對天線安全也造成一定的影響。所以系統(tǒng)辨識的輸入、輸出都采用角度數(shù)據(jù)，其曲線如圖3所示。

圖3 模型輸入、輸出數(shù)據(jù)曲線

2.3 采樣率選擇

某伺服系統(tǒng)帶寬最大為2 Hz，根據(jù)采樣定理，采樣頻率必須大于4 Hz，系統(tǒng)辨識數(shù)據(jù)采樣頻率選擇公式如下[9]：

[fc=5～15B] （5）

式中：[B]為伺服系統(tǒng)帶寬，因此選擇采樣頻率20 Hz滿足要求。

3 Matlab仿真和驗證

Matlab軟件提供了可視化系統(tǒng)辨識接口[4]，界面如圖4所示。該軟件提供了所有常用的系統(tǒng)辨識方法和過程，使用非常方便和靈活。將實驗數(shù)據(jù)導入到軟件中，首先完成適當?shù)姆侄?、去除?shù)據(jù)趨勢等預處理，選擇合適的建模方法，在所建模型中選擇最優(yōu)模型作為最終的模型。

圖4 系統(tǒng)辨識操作界面

為了檢驗模型的精確性[10]，用驗證數(shù)據(jù)對所建模型進行了檢驗，設備實際輸出和模型輸出比對圖如圖5所示。可以看出，模型輸出與實際測量數(shù)據(jù)基本一致，誤差較小，最大為0.1°，滿足要求。

4 結 語

系統(tǒng)辨識方法應用于航天測控設備中的模型構建是可行的，所構建的模型與實際設備相符，精度滿足要求。本文采用了離線辨識方法，用所辨識模型優(yōu)化設備參數(shù)設計，達到了預期效果。在線辨識是自適應控制的基礎和前提，但在機電系統(tǒng)中進行在線辨識存在諸多困難，其實時性、準確性和安全性難以保證，需要進一步研究。

圖5 模型輸出和實際輸出比較

參考文獻

[1] 劉興堂.現(xiàn)代辨識工程[M].北京：國防工業(yè)出版社，2006.

[2] 龐中華，崔紅.系統(tǒng)辨識與自適應控制Matlab仿真[M].北京：北京航空航天大學出版社，2013.

[3] 劉宏才.系統(tǒng)辨識與參數(shù)估計[M].北京：冶金工業(yè)出版社，1999.

[4] 齊曉慧，田慶民，董海瑞.基于Matlab系統(tǒng)辨識工具箱的系統(tǒng)建模[J].兵工自動化，2006，25（10）：88?90.

[5] 潘立登，潘仰東.系統(tǒng)辨識與建模[M].北京：化學工業(yè)出版社，2004.

[6] 王秀峰，盧桂章.系統(tǒng)建模與辨識[M].北京：電子工業(yè)出版社，2004.

[7] 易繼鍇，江祥賢.電氣傳動自動控制原理與設計[M].北京：北京工業(yè)大學出版社，1997.

[8] 孔凡才.自動控制系統(tǒng)[M].北京：機械工業(yè)出版社，2003.

[9] 李言俊，張科.系統(tǒng)辨識理論及應用[M].北京：國防工業(yè)出版社，2004.

[10] 高隆昌，楊元.數(shù)學建模基礎理論[M].北京：科學出版社，2007.