一種數(shù)據(jù)鏈傳輸延遲建模及其補償方法

馮忠華 王新龍 王 彬

(北京航空航天大學(xué) 宇航學(xué)院，北京 100191)

數(shù)據(jù)鏈作為各作戰(zhàn)單元之間進行信息傳輸?shù)耐ㄐ沛溌?，是實現(xiàn)協(xié)同制導(dǎo)的關(guān)鍵環(huán)節(jié)［1－4］.通過數(shù)據(jù)鏈，將各作戰(zhàn)單元有機結(jié)合為一體，實現(xiàn)了各作戰(zhàn)平臺之間信息共享、統(tǒng)一指揮，協(xié)同攻擊及遠程精確打擊，從而能夠大幅度提高作戰(zhàn)效能.然而，由于數(shù)據(jù)的發(fā)送、傳輸和處理需要一定的時間，因此，當(dāng)目標(biāo)信息由一個作戰(zhàn)單元傳送給另一個作戰(zhàn)單元時，由于數(shù)據(jù)傳輸延遲，使接收到的目標(biāo)信息產(chǎn)生誤差［5－7］.美國先進 AIM-120D 超視距空空導(dǎo)彈與數(shù)據(jù)鏈相結(jié)合使用，使得預(yù)警機探測，引導(dǎo)載機發(fā)射，并協(xié)同友機制導(dǎo)的攻擊方式(即所謂的第3方制導(dǎo))已成為可能［2］.

鑒于此，本文對協(xié)同制導(dǎo)過程中的數(shù)據(jù)鏈傳輸延遲誤差進行了建模，并提出了一種基于“當(dāng)前”統(tǒng)計模型的數(shù)據(jù)傳輸延遲誤差補償?shù)淖赃m應(yīng)卡爾曼濾波方法.

1 數(shù)據(jù)傳輸延遲分析及建模

1.1 數(shù)據(jù)傳輸延遲分析

消息傳輸時間延遲是指從發(fā)送站產(chǎn)生消息到接收站接收到消息這兩個時刻之間的時間間隔.消息在一個作戰(zhàn)單元生成后首先存儲到緩沖區(qū)，若系統(tǒng)空閑，并且作戰(zhàn)單元中的可用服務(wù)時隙到達時就接受服務(wù)，否則按先到先服務(wù)的順序排隊等待接受服務(wù).消息服務(wù)完后通過天線發(fā)送出去，以電磁信號的形式在介質(zhì)中進行傳播.另一個作戰(zhàn)單元的接收站接收消息，這樣就完成了數(shù)據(jù)傳輸.

由數(shù)據(jù)傳輸過程可知，數(shù)據(jù)在傳輸過程中產(chǎn)生的時延td主要由排隊等待時間Wq、系統(tǒng)服務(wù)時間Ws和消息傳播時間Wp構(gòu)成［5］，即

目前，各作戰(zhàn)單元之間的通信距離一般在550 km范圍內(nèi)，因此，數(shù)據(jù)的最大傳播時間Wp大約為1.83ms.如果目標(biāo)飛行速度為500m/s，這時目標(biāo)在數(shù)據(jù)傳播時間內(nèi)的運動距離為0.915m，此時目標(biāo)位置誤差遠小于雷達的測量精度，因此，傳播時間Wp一般不作考慮.

消息服務(wù)時間Ws是指作戰(zhàn)單元開始對消息進行服務(wù)直到完成服務(wù)所需的時間，如果消息僅由一個數(shù)據(jù)包組成，則消息服務(wù)時間可表示為分配給作戰(zhàn)單元中兩時隙之間的時間間隔，即

式中，n為分配給作戰(zhàn)單元可用服務(wù)時隙中兩時隙之間的時隙數(shù);τ為時隙大小.

由于消息只有在作戰(zhàn)單元所分配的時隙中才能傳輸，當(dāng)消息到達時，它不僅要等比它先到發(fā)送隊列的消息傳送完畢，而且還要等作戰(zhàn)單元中可用服務(wù)時隙的到來.因此，消息等待時延Wq是影響數(shù)據(jù)傳輸延遲的主要因素.

1.2 數(shù)據(jù)傳輸延遲建模

時分多址(TDMA，Time Division Multiple Access)作為一種有效的多址接入方式在數(shù)據(jù)鏈傳輸系統(tǒng)中得到了廣泛應(yīng)用.設(shè)在一個TDMA系統(tǒng)中有N個作戰(zhàn)單元，消息傳輸時間被分成連續(xù)的時幀，每個時幀又由連續(xù)的時隙組成.令τ為時隙大小，M為每幀中的時隙數(shù)，則時幀長度TF=Mτ.作戰(zhàn)單元中傳輸?shù)南⑹怯梢粋€或多個固定長度的數(shù)據(jù)包組成，并且每個數(shù)據(jù)包僅需一個時隙即可發(fā)送完畢.如果作戰(zhàn)單元i在一個時幀中分配的時隙數(shù)ni=1，則其每隔M個時隙傳送一個數(shù)據(jù)包.

N(t)表示t時刻作戰(zhàn)單元中的消息隊長，其概率分布可表示為

概率母函數(shù)為

若消息到達時系統(tǒng)空閑，則消息中第一個到達數(shù)據(jù)包的服務(wù)時間服從參數(shù)為TF/M到TF+TF/M的隨機分布，設(shè)其分布為B^(t)，一階矩和二階矩分別為b^1，b^2;而消息中隨后到達的數(shù)據(jù)包服務(wù)時間為TF，設(shè)其服務(wù)時間分布為B(t)，一階矩和二階矩分別為b1，b2.假設(shè)作戰(zhàn)單元中第n個消息由Fn個數(shù)據(jù)包組成，且{Fn，n≥1}是獨立同分布的，其一階矩和二階矩分別為Fn1和Fn2.并且第n個消息的服務(wù)時間為Sn(t)，則{Sn(t)，n≥1}也是獨立同分布的.

假設(shè)作戰(zhàn)單元中消息的生成服從參數(shù)為λ的Poisson過程，根據(jù)消息服務(wù)時間分布以及消息先到先服務(wù)的規(guī)則，可以將消息傳輸過程抽象為一個排隊系統(tǒng)，進而利用排隊理論來分析數(shù)據(jù)傳輸時間延遲.因此，當(dāng) λb1＜1 時，有［8］

式中，B*(s)，B^*(s)是消息服務(wù)時間分布的拉普拉斯—斯蒂爾切斯變換;P(z)是作戰(zhàn)單元中消息數(shù)為n的概率母函數(shù);P0是上式在z=1處得到的.

進一步，根據(jù)文獻［8］可得

則，消息傳輸延遲為

式(8)給出了作戰(zhàn)單元在每個時幀內(nèi)僅分配一個時隙時產(chǎn)生的消息傳輸延遲，將該結(jié)論進行推廣，即作戰(zhàn)單元在每個時幀內(nèi)分配多個時隙，并且分配的時隙位置在時幀中是均勻分布.設(shè)作戰(zhàn)單元i(1≤i≤N)在每個時幀中分配的時隙數(shù)為 ni，則將式(8)中的用代替，可得

根據(jù)TDMA協(xié)議［5］，假設(shè)每條消息僅由一個數(shù)據(jù)包組成，即Fn1=Fn2=1.再根據(jù)式(9)，可得到消息到達率λ和服務(wù)率μ與傳輸延遲時間的變化情況(其中μ為服務(wù)時間的倒數(shù))，如圖1和圖2所示.可以看出，傳輸延遲時間與消息的服務(wù)率μ和到達率λ有關(guān):當(dāng)λ一定時，延遲時間隨著μ增大而減小;當(dāng)μ一定時，延遲時間隨著λ增大而增大;當(dāng)λ=μ時，傳輸延遲時間急劇增大.

圖1 延遲時間隨消息服務(wù)率μ的變化情況

圖2 延遲時間隨消息到達率λ的變化情況

2 數(shù)據(jù)鏈延遲補償方法

2.1 機動目標(biāo)運動模型

對于機動目標(biāo)來說，目標(biāo)下一時刻的位置是未知的，因此需要建立更合理、更準(zhǔn)確的機動目標(biāo)的運動模型來描述目標(biāo)的運動.由于當(dāng)目標(biāo)以某一加速度機動時，下一時刻目標(biāo)加速度的取值只能在“當(dāng)前”加速度的鄰域內(nèi).因此，通常采用“當(dāng)前”統(tǒng)計模型來描述機動目標(biāo)的運動［9］.

由于目標(biāo)的隨機機動加速度在時間軸上符合一階時間相關(guān)過程(目標(biāo)在3個方向上的運動是相互獨立的，因此僅對x方向進行分析):

根據(jù)式(10)，可得直角坐標(biāo)系下表示的目標(biāo)運動的狀態(tài)方程為

式中，Q0是與α和采樣周期T有關(guān)的常量矩陣.

展開式(11)，取其第3個分量，可得

對式(13)兩邊再取均值，則有

式中，zk為到當(dāng)前時刻為止目標(biāo)的所有量測值;a^(k)為目標(biāo)“當(dāng)前”加速度的預(yù)測值.

接著，利用目標(biāo)機動加速度與相鄰采樣時刻目標(biāo)位置變化之間的關(guān)系對Q(k)陣進行自適應(yīng)調(diào)整［10］.由于目標(biāo)在k－1時刻與k時刻的位移Δd與其加速度的變化量Δa之間的關(guān)系為

而Δd又可表示為

因此，由式(15)、式(16)可得

另外，目標(biāo)當(dāng)前機動加速度a(k)可表示為當(dāng)前加速度的均值與兩采樣時刻間隔之間加速度的變化Δa之和，即

并且，由于當(dāng)前機動加速度a(k)服從修正的瑞利分布［10］，則可得到

取目標(biāo)在地球直角坐標(biāo)系下的位置x為量測量，則機動目標(biāo)的量測方程為

式中，H(k)=［1 0 0］是觀測矩陣;V(k)是均值為0;方差為R(k)的觀測噪聲.

2.2 基于自適應(yīng)卡爾曼濾波的數(shù)據(jù)傳輸延遲誤差補償方法

采用機動目標(biāo)預(yù)測算法，根據(jù)k+td時刻得到的目標(biāo)信息Z(k)，通過自適應(yīng)卡爾曼濾波跟蹤算法可得到k時刻目標(biāo)信息的估計值X(k).

此時，目標(biāo)量測信息的預(yù)測值為

而目標(biāo)在信息傳輸延遲時間間隔td內(nèi)移動的距離可表示為

從而，可得到k+td時刻延遲補償后的目標(biāo)信息為

3 仿真驗證與分析

目標(biāo)的仿真條件設(shè)定為:目標(biāo)初始位置為［100;6;80］km，初始速度為300m/s，初始航跡傾角為0°，航跡偏角為150°;目標(biāo)采用的機動方式為:目標(biāo)先勻速直線飛行，10 s后以8 g的加速度右轉(zhuǎn)彎機動，轉(zhuǎn)過90°后勻速直線飛行，10 s后再以8g的加速度左轉(zhuǎn)彎機動，轉(zhuǎn)過90°后勻速直線飛行，10 s后再以2 g的加速度爬升飛行，當(dāng)俯仰角為45°時再以2 g的加速度改平飛行，最后目標(biāo)在水平面內(nèi)平飛.仿真中目標(biāo)的機動頻率α=0.05，采樣周期T=1 s，仿真時間80 s(其中前30 s為預(yù)警機探測目標(biāo)信息，后50 s為友機探測目標(biāo)信息).

數(shù)據(jù)傳輸過程中目標(biāo)信息的傳輸延遲除了受信息的排隊等待時間、系統(tǒng)服務(wù)時間、消息傳播時間等因素的影響外，還會受到作戰(zhàn)單元對目標(biāo)信息的處理以及算法運行時間等誤差因素的影響，使得數(shù)據(jù)傳輸延遲時間具有隨機性，并服從均勻分布.因此，總的延遲時間td可設(shè)定為0.5 s到1 s的隨機數(shù)［11］.

采用機動目標(biāo)3個方向的位置、速度和加速度作為目標(biāo)的狀態(tài)變量，對目標(biāo)的狀態(tài)進行估計，再根據(jù)式(22)對目標(biāo)的測量信息進行延遲補償，計算得到延遲補償前、后目標(biāo)的量測信息.將補償前、后的目標(biāo)測量值與目標(biāo)真實的測量值之差作為目標(biāo)信息的測量誤差，仿真結(jié)果見圖3～圖5.

圖3 目標(biāo)距離測量誤差

圖4 目標(biāo)俯仰角測量誤差

圖5 目標(biāo)偏航角測量誤差

由圖3～圖5可看出，數(shù)據(jù)延遲補償前目標(biāo)測量誤差很大，經(jīng)延遲補償后誤差明顯降低，位置測量誤差最大減少307.268 9m，俯仰角測量誤差最大減少0.0926°，偏航角測量誤差最大減少0.1506°.可見，延遲補償后的目標(biāo)信息更接近于目標(biāo)的真實運動信息.

進一步，采用蒙特卡羅方法對數(shù)據(jù)鏈傳輸延遲補償前后導(dǎo)引頭對目標(biāo)的截獲情況進行分析.

蒙特卡羅仿真條件為:導(dǎo)彈初始位置為［0;5000;0］m，初始速度為300m/s，初始彈道傾角10°，彈道偏角 － 30°;友機初始位置為［1 000;6000;2000］m，初始速度為［0;0;－200］m/s;目標(biāo)位置測量誤差為100m;導(dǎo)引頭天線3 dB波束寬度為4°，蒙特卡羅仿真次數(shù)為200，截獲時彈目距離小于20 km.

根據(jù)目標(biāo)截獲時的導(dǎo)彈信息、目標(biāo)的真實信息以及數(shù)據(jù)傳輸延遲補償前、后的目標(biāo)信息可以得到補償前、后的目標(biāo)截獲概率，如圖6所示.

圖6 目標(biāo)在導(dǎo)引頭視場中的分布

圖6中，虛線代表導(dǎo)引頭視場的范圍，目標(biāo)落在導(dǎo)引頭視場范圍內(nèi)時表示可以被截獲，否則不能被截獲.根據(jù)200次蒙特卡羅仿真實驗結(jié)果，延遲補償前目標(biāo)在導(dǎo)引頭視場中的分布范圍較大，其中有154次落在導(dǎo)引頭視場范圍內(nèi)，這種情況下目標(biāo)截獲概率為77%;而延遲補償后目標(biāo)在導(dǎo)引頭視場中的分布比較集中，有198次落在導(dǎo)引頭視場范圍內(nèi)，此時目標(biāo)截獲概率為99%.可見，目標(biāo)測量信息經(jīng)延遲補償后能夠大大提高中末制導(dǎo)交班時導(dǎo)彈對目標(biāo)的截獲概率.

4 結(jié)論

通過對數(shù)據(jù)傳輸延遲過程進行分析，建立了數(shù)據(jù)傳輸延遲模型，并利用目標(biāo)機動“當(dāng)前”統(tǒng)計模型與自適應(yīng)卡爾曼濾波器相結(jié)合對數(shù)據(jù)鏈傳輸延遲誤差進行了補償，可以得到以下結(jié)論:

1)數(shù)據(jù)在傳輸過程中產(chǎn)生的時延主要由排隊等待時間、系統(tǒng)服務(wù)時間和消息傳播時間構(gòu)成.其中，系統(tǒng)服務(wù)時間和排隊等待時間是影響數(shù)據(jù)傳輸延遲的主要因素.

2)利用機動目標(biāo)加速度與相鄰采樣時刻目標(biāo)位置估計變化之間的關(guān)系，可實現(xiàn)對系統(tǒng)噪聲方差陣的自適應(yīng)調(diào)整.

3)根據(jù)目標(biāo)測量信息的預(yù)測值以及數(shù)據(jù)傳輸延遲間隔，利用自適應(yīng)卡爾曼濾波可遞推得到延遲誤差補償后的目標(biāo)信息，實現(xiàn)對目標(biāo)信息傳輸延遲誤差的補償.

References)

［1］李林森，于海勛，韓志剛，等.論協(xié)同空戰(zhàn)及其關(guān)鍵技術(shù)［J］.電光與控制，2000(1):29－34 Li Linsen，Yu Haixun，Han Zhigang，et al.The discussion of cooperative air combat and its key technology［J］.Electronics Optics＆ Control，2000(1):29 －34(in Chinese)

［2］童明成，張鳳霞，賀建良.網(wǎng)絡(luò)環(huán)境下空中多平臺聯(lián)合攻擊作戰(zhàn)方式［J］.電光與控制，2009，16(5):19 －22 Tong Mingcheng，Zhang Fengxia，He Jianliang.Aerial multi-platform joint-attack warfare modes in network environment［J］.Electronics Optics ＆ Control，2009，16(5):19 －22(in Chinese)

［3］劉志華，季海波.數(shù)據(jù)鏈技術(shù)在戰(zhàn)術(shù)武器中的應(yīng)用方案研究［J］.航天控制，2009，27(3):61 －64 Liu Zhihua，Ji Haibo.Study of application project of data link technologies in tactical weapons［J］.Aerospace Control，2009，27(3):61－64(in Chinese)

［4］張斌，黃俊，寇英信，等.雙機合作下的遠程空空導(dǎo)彈攻擊技術(shù)［J］.火力與指揮控制，2009，34(5):92 －94 Zhang Bin，Huang Jun，Kou Yingxin，et al.Research on technology of long-range air-to-air missile attacking by two cooperative fighters［J］.Fire Control ＆ Command Control，2009，34(5):92－94(in Chinese)

［5］任培，周經(jīng)倫，羅鵬程，等.JTIDS信息傳輸時間延遲分析［J］.計算機工程，2008，34(21):23 －25 Ren Pei，Zhou Jinglun，Luo Pengcheng，et al.Time delay analysis of JTIDS message transmission［J］.Computer Engineering，2008，34(21):23 －25(in Chinese)

［6］ Mitchell G E，Gander J F.Impact of time delays on the JTIDS network［C］//IEEE National Aerospace and Electronics Conference.Dayton，OH，USA:IEEE，1990:440 －445

［7］ Rubin I.Message delays in FDMA and TDMA communication channels［J］.IEEE Transactions and Communication，1979，27(5):769－777

［8］ Lam SS.Delay analysis of a time division multiple access(TDMA)channel［J］.IEEE Transactions on Communication，1977，25(12):1489－1494

［9］樊宗臣，吳曉燕.機動目標(biāo)跟蹤的自適應(yīng)卡爾曼濾波算法實現(xiàn)［J］.戰(zhàn)術(shù)導(dǎo)彈技術(shù)，2009(1):79－81 Fan Zongchen，Wu Xiaoyan.Study of adaptive kalman filtering algorithm in maneuvering target track［J］.Tactical Missile Technology，2009(1):79 －81(in Chinese)

［10］胡振濤，劉先省.基于“當(dāng)前”統(tǒng)計模型的一種改進機動目標(biāo)跟蹤算法［J］.山東大學(xué)學(xué)報，2005，35(3):111 －114 Hu Zhentao，Liu Xianxing.An improved maneuvering target tracking algorithm based on the“current”statistical model［J］.Journal of Shangdong University(Engineering Science)，2005，35(3):111 －114(in Chinese)

［11］吳龍，劉正敏，姜長生，等.基于卡爾曼濾波算法的數(shù)據(jù)鏈延時補償［J］.航空兵器，2010，(5):14－17 Wu Long，Liu Zhengmin，Jiang Changsheng，et al.Data link delay-time compensation based on Kalman filtering algorithm［J］.Aero Weaponry，2010，(5):14 －17(in Chinese)