摘 要： 針對無人機機載武器低成本精確打擊的需求，提出利用機載光電吊艙獲取制導(dǎo)信息、采用視線制導(dǎo)法進行制導(dǎo)律設(shè)計的方案。通過分析載機、導(dǎo)彈和目標(biāo)的運動關(guān)系，建立了系統(tǒng)狀態(tài)空間模型，基于反步法和滑?？刂圃O(shè)計了改進視線制導(dǎo)律。對制導(dǎo)律中無法由載機獲取的參數(shù)進行了討論，采用跟蹤微分器提取視線角并估計視線角速度信息。數(shù)值仿真結(jié)果表明，在視線角存在噪聲且導(dǎo)彈矢徑拉偏的情況下，改進視線制導(dǎo)律仍可保證對地面目標(biāo)的精準(zhǔn)打擊。與傳統(tǒng)視線制導(dǎo)律相比，載機的運動方向不受限制，可在發(fā)射導(dǎo)彈后遠離目標(biāo)以免受到威脅，符合現(xiàn)實需求。

關(guān)鍵詞： 視線制導(dǎo); 反步法; 滑?？刂? 機載空地導(dǎo)彈; 跟蹤微分器

中圖分類號： TJ 765.3

文獻標(biāo)志碼： A

DOI：10.12305/j.issn.1001-506X.2024.08.25

Design of improved line-of-sight guidance law based on aircraft

visual information

ZHANG Wenwen¹， ZHANG Cheng^1，2，*， ZHENG Chenming¹， CHENG Runbei¹， CHEN Tianle³

（1. School of Aerospace Engineering， Beijing Institute of Technology， Beijing 100081， China; 2. Key

Laboratory of Dynamics and Control of Flight Vehicle， Ministry of Education， Beijing 100081， China;

3. Shanghai Institute of Aerospace System Engineering， Shanghai 201109， China）

Abstract： For the demand of low-cost precision strike of unmanned airborne weapons， the proposal of using airborne optoelectronic pods to obtain guidance information and using line-of-sight （LOS） guidance method for guidance law design is proposed. By analyzing the motion relationship of the carrier aircraft， missile， and target， a system state space model is established， and an improved LOS guidance law is designed based on the backstepping method and sliding mode control. The parameters in the guidance law that cannot be obtained by the aircraft are discussed. The tracking differentiator is used to extract the LOS angle and estimate the LOS velocity information. Numerical simulation results show that the improved LOS guidance law can still ensure the accurate strike on the ground target in the presence of LOS angle noise and missile vector pull-off. Compared with the traditional LOS guidance law， the movement direction of the aircraft is not limited， so it can stay away from the target after missile launching to avoid being threatened， which meets the practical demand.

Keywords： line-of-sight （LOS） guidance; backstepping method; sliding mode control; airborne air-to-ground missile; tracking differentiator

0 引 言

近年來，無人機在軍事行動中得到了廣泛的應(yīng)用，其功能不再局限于單一的偵察通信，而是具備了直接參與目標(biāo)精確打擊的能力，已成為現(xiàn)代智能化戰(zhàn)場中至關(guān)重要的武器裝備^［1-3］。低成本精確打擊是無人機機載武器發(fā)展中持久關(guān)注的主題^［4］，為了在保證打擊精度的前提下降低成本，本文提出基于載機光電吊艙獲取制導(dǎo)信息的方案。通過充分利用無人機偵查目標(biāo)的機載光電吊艙，無需在導(dǎo)彈上安裝導(dǎo)引頭或在載機平臺上增設(shè)雷達制導(dǎo)站等設(shè)備，實現(xiàn)低成本制導(dǎo)的目的。考慮到光電吊艙一般僅能獲取導(dǎo)彈和目標(biāo)的視線角等基本信息，因此本文采用視線制導(dǎo)方法來實現(xiàn)上述方案。

視線制導(dǎo)通常也稱三點法制導(dǎo)^［5］，其基本原理是在制導(dǎo)過程中令導(dǎo)彈始終保持在目標(biāo)跟蹤裝置與目標(biāo)構(gòu)成的視線上。傳統(tǒng)的視線制導(dǎo)法通常基于導(dǎo)彈與視線之間的偏離量生成制導(dǎo)指令。例如，王狂飆等^［6］將導(dǎo)彈偏離視線的距離乘以比例系數(shù)，并通過超前-滯后網(wǎng)絡(luò)進行校正，從而得到了制導(dǎo)指令的形式。針對傳統(tǒng)視線制導(dǎo)法中存在的動態(tài)誤差問題，Zarchan^［7］提出一種前饋補償方法，即在制導(dǎo)指令中加入視線移動的加速度，將其作為補償項。Ratnoo^［8］將此方法應(yīng)用于正弦機動目標(biāo)的攔截，發(fā)現(xiàn)相比于比例導(dǎo)引法，該方法在脫靶量和法向過載方面的性能指標(biāo)表現(xiàn)更為優(yōu)越。

視線制導(dǎo)法也可應(yīng)用于目標(biāo)飛行器發(fā)射防御彈，對攻擊彈進行主動防御的“三體對抗”問題。Ratnoo等^［9］在傳統(tǒng)視線制導(dǎo)律中加入了受保護飛行器的加速度項，在攻擊彈分別采用迎擊、尾追和側(cè)向攻擊的方式下，采用視線制導(dǎo)的防御彈均能成功實施攔截，且防御彈的速度和法向過載均低于攻擊彈。Yamasaki等^［10］基于視線制導(dǎo)的思想，提出一種防御導(dǎo)彈的三角制導(dǎo)（triangle guidance， TG）策略。其主要思想是將受保護飛行器、防御彈和攻擊彈所構(gòu)成的碰撞三角形退化成一條直線，即令防御彈保持在受保護飛行器和攻擊彈所構(gòu)成的直線上。受傳統(tǒng)比例導(dǎo)引律零化彈目視線轉(zhuǎn)率的啟發(fā)，Yamasaki等^［11］提出一種機載視線指令制導(dǎo)律（airborne command to line-of-sight guidance， A-CLOSG），通過令防御彈的過載指令正比于期望橫向速度與實際橫向速度之差，從而使防御彈控制在受保護飛行器和攻擊彈的連線上。

為了在三維空間中攔截高速機動目標(biāo)，Ahi等^［12］將導(dǎo)彈偏離視線的角加速度作為控制量，提出了一種由前饋控制和嵌套飽和反饋法構(gòu)成的新型視線制導(dǎo)策略。Ha等^［13］提出一種基于反饋線性化的改進視線制導(dǎo)律設(shè)計方法。該方法結(jié)合了導(dǎo)彈的姿態(tài)角信息，通過構(gòu)建六自由度攔截動力學(xué)模型，將導(dǎo)彈偏離視線的距離作為輸出方程，建立了系統(tǒng)的狀態(tài)方程，從而將攔截問題轉(zhuǎn)化為誤差跟蹤問題?；谏鲜瞿Ｐ停墨I［14-16］分別采用了模糊自適應(yīng)、滑?？刂啤⒎床娇刂频确椒ㄔO(shè)計視線制導(dǎo)律，在攔截高速機動目標(biāo)方面取得了顯著的效果。

目前，對于視線制導(dǎo)律的研究主要集中在制導(dǎo)站固定于地面不動進行對空或?qū)Φ毓舻膱鼍爸?，即使在“三體對抗”這類以目標(biāo)飛行器作為動制導(dǎo)站的場景中，也沒有充分利用制導(dǎo)站的信息來研究制導(dǎo)律，關(guān)于應(yīng)用于無人機機載導(dǎo)彈的視線制導(dǎo)律方面的研究還相對較少。此外，現(xiàn)有視線制導(dǎo)律的研究主要關(guān)注導(dǎo)彈的攔截精度以及需用過載等方面，而對于制導(dǎo)律中參數(shù)的不確定性以及擾動等影響的探討還不夠充分。

反步滑?？刂疲╞ackstepping sliding mode control， BSMC）由于對系統(tǒng)參數(shù)的擾動和不確定性具有較強的魯棒性，在控制系統(tǒng)中得到了廣泛應(yīng)用^［17-22］。相較于傳統(tǒng)的滑模控制，BSMC通過微分迭代可以有效抑制抖振問題^［18-19］。在導(dǎo)彈制導(dǎo)律的設(shè)計中，BSMC常用于考慮將目標(biāo)機動作為干擾的場景^［20-22］。然而，在系統(tǒng)制導(dǎo)參數(shù)獲取過程中，受到噪聲和估計誤差等因素的影響，同樣導(dǎo)致擾動和不確定性。對于這類情形，BSMC是否能夠保持良好的性能仍需深入研究。

本文提出一種基于載機視覺信息的視線制導(dǎo)方法，采用文獻［13］的數(shù)學(xué)模型，并對原模型進行擴展，以整合載機的狀態(tài)信息，將制導(dǎo)律設(shè)計問題轉(zhuǎn)化為三維空間內(nèi)的誤差跟蹤問題。基于擴展后的模型，運用BSMC對視線制導(dǎo)律進行了改進設(shè)計，以應(yīng)對目標(biāo)機動、制導(dǎo)信息受到噪聲及估計誤差干擾等因素引起的擾動和不確定性，確保系統(tǒng)具有良好的魯棒性。針對制導(dǎo)律中所需的視線角速度，由于載機無法直接提供，采用跟蹤微分器進行濾波和估計。最后，在存在視線角噪聲和估計誤差的條件下，對所設(shè)計的改進視線制導(dǎo)律進行了六自由度全彈道仿真，以驗證其在復(fù)雜環(huán)境中的有效性和魯棒性。

1 數(shù)學(xué)模型

1.1 導(dǎo)彈動力學(xué)模型

載機、導(dǎo)彈與目標(biāo)在三維空間中的幾何關(guān)系如圖1所示，分別用A、M和T表示。

圖1中，ε和β分別表示由載機光電吊艙探測得到的高低角和方位角，R_m和R_t稱為導(dǎo)彈矢徑和目標(biāo)矢徑，表示載機到導(dǎo)彈/目標(biāo)的距離。為方便表述，本文將載機與導(dǎo)彈/目標(biāo)的連線稱為機彈/機目視線。Oxyz為地面坐標(biāo)系，Mx₂y₂z₂為彈道坐標(biāo)系，V_m為導(dǎo)彈速度矢量，θ和_V分別為彈道傾角和彈道偏角。

設(shè)［x_m，y_m，z_m］^T為導(dǎo)彈在地面坐標(biāo)系下的位置，［a_x2，a_y2，a_z2］^T為導(dǎo)彈在彈道坐標(biāo)系下的加速度，則導(dǎo)彈在地面坐標(biāo)系下的運動方程為

導(dǎo)彈質(zhì)心運動的動力學(xué)方程為

1.2 跟蹤誤差的表達

視線制導(dǎo)跟蹤誤差的表達通常有3種：① 定義為導(dǎo)彈視線角與目標(biāo)視線角之差;② 將導(dǎo)彈與機目視線分別向地面坐標(biāo)系中的側(cè)向平面以及包含機彈視線的縱向平面進行投影，各平面內(nèi)導(dǎo)彈與機目視線投影的垂直距離即為跟蹤誤差;③ 在以載機和目標(biāo)構(gòu)成的目標(biāo)視線坐標(biāo)系下表示，導(dǎo)彈在該坐標(biāo)系下垂直于機目視線的坐標(biāo)即為跟蹤誤差。本文采用第③種定義下的跟蹤誤差，在此定義下，跟蹤誤差的靈敏度在制導(dǎo)后半段不會隨著導(dǎo)彈接近目標(biāo)而降低^［15］。

如圖2所示，定義目標(biāo)視線坐標(biāo)系為Ax_Ly_Lz_L，P_m=［x_m，y_m，z_m］^T和P_a=［x_a，y_a，z_a］^T分別為導(dǎo)彈和載機在地面坐標(biāo)系下的位置，則導(dǎo)彈在目標(biāo)視線坐標(biāo)系下的位置可表示為

P_m－P_a表示導(dǎo)彈在以載機為原點的地面坐標(biāo)系下的位置，可通過載機的極坐標(biāo)系轉(zhuǎn)化得到，即

由前述定義可知，導(dǎo)彈在目標(biāo)視線坐標(biāo)系下垂直于R_t的兩個坐標(biāo)y_L和z_L即為跟蹤誤差。為把跟蹤誤差表示為載機可觀測的量，聯(lián)立式（3）和式（4）可得

式（5）中的導(dǎo)彈/目標(biāo)視線角信息由載機提供，導(dǎo)彈矢徑R_m將在后文說明，可當(dāng)作已知信息處理。因此，由式（5）得到的跟蹤誤差可當(dāng)成已知量，在后續(xù)狀態(tài)空間模型的建立中可以作為輸出量。

1.3 狀態(tài)空間表達式

根據(jù)上述動力學(xué)關(guān)系和跟蹤誤差的表達，選取系統(tǒng)狀態(tài)變量為

X=［x_m，y_m，z_m，x·_m，y·_m，z·_m，_V，θ］^T（6）

控制量：

U=［a_x2，a_y2，a_z2］^T（7）

輸出量：

為方便描述，定義以下矩陣：

聯(lián)立式（1）、式（7）和式（9），可得

聯(lián)立式（2）、式（7）和式（9），可得

聯(lián)立式（1）、式（8）和式（9），可得

Y=L_t（P_m－P_a）（12）

至此可構(gòu)造出系統(tǒng)的狀態(tài)空間表達式

對于實際的系統(tǒng)，導(dǎo)彈的輸入控制量只有a_y2和a_z2兩個法向加速度，式（7）把切向加速度a_x2當(dāng)作控制量處理，是為了避免出現(xiàn)奇異問題。因此，在后續(xù)求得的控制量中，只有a_y2和a_z2起作用，但可保證跟蹤誤差Y快速衰減至零，從而使導(dǎo)彈沿機目視線飛行直至命中目標(biāo)。

2 改進視線制導(dǎo)律的設(shè)計

由狀態(tài)空間表達式（13），對輸出跟蹤誤差Y分別求一階、二階導(dǎo)數(shù)并化簡，可得

2.1 反步滑模視線制導(dǎo)律

根據(jù)反步法的原理，結(jié)合滑?？刂频乃枷?，基于反步滑模控制的改進視線制導(dǎo)律設(shè)計步驟如下。

步驟1 設(shè)計Y·的期望虛擬控制量。

定義系統(tǒng)輸出量的誤差為

e₁=Y_d－Y（16）

式中：Y_d表示期望的跟蹤誤差輸出量。

對式（16）求導(dǎo)有

e·₁=Y·_d－Y·（17）

構(gòu)造第一個Lyapunov函數(shù)為

V₁=12e^T₁e₁（18）

則V₁對時間的導(dǎo)數(shù)為

V·₁=e^T₁e·₁=e^T₁（Y·_d－Y·）（19）

為使系統(tǒng)達到穩(wěn)定，需使V·₁負(fù)定。將Y·視為虛擬控制量，令其期望值為η₁，設(shè)計如下：

η₁=Y·_d+K₁e₁（20）

式中：K₁=k₁

k₂，且k₁，k₂gt;0。

步驟2 設(shè)計Y¨的期望虛擬控制量。

定義Y·的誤差為

e₂=η₁－Y·（21）

把式（20）代入式（21）可得

e₂=Y·_d+K₁e₁－Y·=e·₁+K₁e₁（22）

選取線性滑模面

s=e₂=K₁e₁+e·₁（23）

對線性滑模面求導(dǎo)可得

s·=e·₂=η·₁－Y¨（24）

構(gòu)造第二個Lyapunov函數(shù)為

V₂=V₁+12s^Ts（25）

則V₂對時間的導(dǎo)數(shù)為

V·₂=e^T₁e·₁+s^Ts·（26）

把式（23）和式（24）代入式（26），可得

V·₂=e^T₁（s－K₁e₁）+s^T（η·₁－Y¨）=

－e^T₁K₁e₁+s^T（e₁+η·₁－Y¨）（27）

將Y¨視為虛擬控制量，令其期望值為η₂。為使系統(tǒng)達到穩(wěn)定，需使V·₂負(fù)定，η₂設(shè)計如下：

η₂=e₁+η·₁+K₂s+ε· sgn（s）（28）

步驟3 設(shè)計真實控制量U_bs。

令Y¨=η₂，即

Y¨=η₂=e₁+η·₁+K₂s+ε·sgn（s）（29）

把式（20）和式（23）代入式（29），有

Y¨=e₁+Y¨_d+K₁e·₁+

K₂（K₁e₁+e·₁）+ε·sgn（K₁e₁+e·₁）（30）

再把式（16）和式（17）代入式（30），整理得

Y¨=Y¨_d+（K₁+K₂）（Y·_d－Y·）+（K₂K₁+I₂）（Y_d－Y）（31）

式中：I₂表示2階單位矩陣。

對于視線制導(dǎo)而言，其輸出量及其一階和二階導(dǎo)數(shù)的期望值均為零，即Y_d=Y·_d=Y¨d=0。據(jù)此對式（31）進行化簡，并把式（15）代入。用U_bs表示利用反步滑模法求得的控制量，并考慮到滑?？刂乒逃械亩墩駟栴}，采用飽和函數(shù)sat（s）代替符號函數(shù)sgn（s），于是有

U_bs=L^T_mL^T_t［－（K₁+K₂）Y·－（K₂K₁+I₂）Y－

ε·sat（K₁Y+Y·）－（C¨+2C·A）X－

G+L¨_tP_a+2L·_tP·_a+L_tP¨_a］（32）

式中：Δ為邊界層;sat（x）為飽和函數(shù)，表達式為

最終得到基于BSMC的視線制導(dǎo)（line-of-sight guidance based on BSMC， BS-LOSG）律，即改進視線制導(dǎo)律為

2.2 穩(wěn)定性分析

將式（29）代入式（27），可得

由式（23）可知s與輸出量Y同維度，為方便表述，設(shè)s=［s₁，s₂］^T，則

V·₂=－e^T₁K₁e₁－s^TK₂s－ε₁|s₁|－ε₂|s₂|（35）

式中：K₁，K₂選取為正定矩陣，且ε₁，ε₂gt;0，故V·₂為負(fù)定的。

V₂正定而V·₂負(fù)定，由Lyapunov穩(wěn)定性原理可知，e₁和s能夠漸進穩(wěn)定地趨于零。這說明所設(shè)計的虛擬控制量Y¨可使Y·漸進穩(wěn)定地趨于期望值η₁，且輸出量Y漸進穩(wěn)定地趨于期望值Y_d，從而說明所設(shè)計的BSMC律的穩(wěn)定性可以得到保證。

2.3 制導(dǎo)律的實施

由于機載光電吊艙僅提供導(dǎo)彈和目標(biāo)的視線角度作為制導(dǎo)信息，而上述制導(dǎo)律的表達式中還涉及其他信息，需要進行綜合分析，以解決制導(dǎo)律在實施過程中的現(xiàn)實問題。

（1） 分析可視為已知量處理的項。導(dǎo)彈矢徑R_m及其變化率R·_m無法通過載機進行實時測量，但可事先預(yù)估其隨時間變化的曲線，并在仿真中進行拉偏處理來模擬其估計值，因此可將該項視為已知量。此外，對于載機而言，其位置P_a、速度P·_a及加速度P¨_a可作為已知信息，直接用于制導(dǎo)律的計算。

（2） 導(dǎo)彈彈道角的近似。鑒于本文所設(shè)定的情景，制導(dǎo)信息僅能從載機處獲取，彈上缺乏能夠測量彈道傾角θ和彈道偏角_V的裝置。因此，需要使用近似值來替代這些參數(shù)。綜合考慮視線制導(dǎo)的特點以及制導(dǎo)信息獲取的難易程度，采用導(dǎo)彈的高低角ε_m和方位角β_m作為導(dǎo)彈彈道傾角θ和彈道偏角_V的近似代替。

（3） 目標(biāo)視線角加速度的簡化。制導(dǎo)律表達式中涉及到目標(biāo)視線角加速度信息，從理論上可以通過對測得的目標(biāo)視線角進行二階導(dǎo)數(shù)計算來獲取。然而，在實際測量中存在噪聲和其他干擾因素，直接對其進行微分運算可能會放大噪聲的影響?？紤]到本文所研究的場景為地面目標(biāo)且載機的機動性不高，故目標(biāo)視線角加速度的值較小。因此，為了簡化制導(dǎo)律，可以忽略此項，后續(xù)的仿真驗證了這種簡化處理的合理性。

（4） 導(dǎo)彈/目標(biāo)視線角速度的提取。同理，視線角速度不采用直接對視線角求微分的方法獲取，可采用卡爾曼濾波及其推廣形式^［23-25］或自抗擾控制中的跟蹤微分器^［26-28］進行提取。本文選用最速跟蹤微分器進行提取，如圖3所示，跟蹤微分器的輸入信號為v（t），輸出信號v₁（t）和v₂（t）分別為v（t）的跟蹤信號和微分信號。

最速跟蹤微分器可在最大加速度的限制下保證對輸入信號的最速跟蹤，其具體形式如下：

x₁（k+1）=x₁（k）+h·x₂（k）

x₂（k+1）=x₂（k）+h·fhan（x₁（k）－ν（t），x₂（k），r，h₀）

（36）

式中：x₁，x₂分別為輸入信號ν（t）的跟蹤信號和微分信號;h為積分步長;fhan（·）為一種最速綜合函數(shù)，可較好地安排輸入信號的過渡過程;r為速度因子，決定跟蹤快慢;h₀為濾波因子，決定噪聲濾波效應(yīng)。

經(jīng)過上述分析，制導(dǎo)律中的各項參數(shù)，除待設(shè)計參數(shù)外均可直接或間接獲取，其關(guān)系如圖4所示。其中，X由式（4）和式（6）計算得到;Y，Y·由式（5）和式（8）計算得到;L_m，L_t，L·_t，L¨_t由式（9）計算得到。

3 仿真驗證

為驗證改進視線制導(dǎo)律的有效性，本文在六自由度導(dǎo)彈模型中進行仿真驗證。設(shè)載機、導(dǎo)彈和目標(biāo)三者的初始狀態(tài)如表1所示，其中載機的速度大小恒定，其運動方向?qū)⒃诤罄m(xù)設(shè)定;導(dǎo)彈的初始速度方向沿著機目視線，其火箭發(fā)動機主要參數(shù)為總沖量I=3 000 N·s，工作時間t_p=0.7 s，推進劑質(zhì)量m_r=2.8 kg;目標(biāo)的運動狀態(tài)在后文中分情況討論。

考慮到本文主要研究內(nèi)容為制導(dǎo)律，因此對于仿真中的控制部分采用簡化形式。如圖5所示，自動駕駛儀采用開環(huán)結(jié)構(gòu)。圖5中，a_yc為法向過載指令，δ_ec為俯仰舵偏角指令；a_y和δ_e表示實際法向過載與實際俯仰舵偏角。

首先由制導(dǎo)律計算得出過載指令a_yc，進而求出舵偏角指令δ_ec，再通過舵機得出實際舵偏角δ_e，作用于彈體后產(chǎn)生實際過載a_y。其中，舵機可等效為二階振蕩環(huán)節(jié)，其傳遞函數(shù)為

δ_eδ_ec=ω²_ns²+2ξω_ns+ω²_n（37）

式中：ω_n為無阻尼固有頻率，仿真中取為30π rad/s;ξ為阻尼比，仿真中取為0.7。

3.1 不同目標(biāo)運動狀態(tài)下的改進視線制導(dǎo)律驗證

采用所提出的改進視線制導(dǎo)律對靜止、勻速直線運動和圓弧機動目標(biāo)進行仿真驗證。設(shè)目標(biāo)勻速直線運動時的速度為20 m/s，航向傾角θ_t為0°，航向偏角_t為45°。圓弧機動的形式采用文獻［29］中的表達式

式中：θ·_t，·_t分別為目標(biāo)的航向傾角變化率和航向偏角變化率。

在仿真過程中，假設(shè)載機的航向傾角θ_a為0°，航向偏角_a為90°。設(shè)載機測得的視線角中含有均值為0、方差為0.01°的高斯白噪聲，導(dǎo)彈矢徑估計值分別拉偏±10%。

經(jīng)過整定的跟蹤微分器參數(shù)設(shè)計為：h=0.001，r=180，h₀=0.06。制導(dǎo)律的相關(guān)設(shè)計參數(shù)為：k₁=k₂=1，k₃=k₄=1，ε₁=ε₂=1，Δ=1。

仿真結(jié)果如圖6～圖9所示。

圖6為導(dǎo)彈矢徑估計值拉偏+10%時，三維空間中載機、導(dǎo)彈和目標(biāo)的運動軌跡圖。表2為相應(yīng)的脫靶量，可以看出，在導(dǎo)彈矢徑估計值分別拉偏±10%的條件下，本文所推導(dǎo)的改進視線制導(dǎo)律均能較好地命中靜止和運動目標(biāo)，且脫靶量大部分都在0.5 m以內(nèi)。

圖7為圖6相應(yīng)的目標(biāo)視線角與導(dǎo)彈視線角之差隨時間變化的曲線?？梢钥闯?，視線角偏差能夠快速衰減至零，符合視線制導(dǎo)規(guī)律。

為了驗證跟蹤微分器的效果，以采用改進視線制導(dǎo)律攻擊靜止目標(biāo)為例，圖8為原始帶有噪聲的目標(biāo)方位角信號及其經(jīng)過跟蹤微分器提取后的角與角速度信號。相比于原始信號，經(jīng)過處理的視線角信號更加平滑，角速度信號也有效地抑制了噪聲放大效應(yīng)，可用于制導(dǎo)指令的計算。

圖9為導(dǎo)彈法向過載指令與實際過載指令的對比圖，考慮到導(dǎo)彈的實際過載限制，設(shè)定過載指令限幅為±5g。由于在角速度的提取過程中存在一定的抖振，因此過載指令曲線也呈現(xiàn)出一定的不光滑性，但其對制導(dǎo)效果產(chǎn)生的影響較小。從圖9可以觀察到，實際過載能夠較好地跟蹤過載指令。

3.2 不同載機運動方向下改進視線制導(dǎo)律驗證

為了驗證本文改進視線制導(dǎo)律在面對載機運動方向變化時的魯棒性，將其與傳統(tǒng)視線制導(dǎo)律進行對比。傳統(tǒng)視線制導(dǎo)律根據(jù)導(dǎo)彈偏離機目視線的距離產(chǎn)生制導(dǎo)指令，參考文獻［30］，其制導(dǎo)律表達式為

a_yc=－k_py_q－k_dy·_q－k_i∫^t₀y_qdt（39）

式中：k_p，k_i，k_d為待設(shè)計參數(shù);y_q為導(dǎo)彈偏離機目視線的距離，其在俯仰和偏航平面的表達式分別為R_m（ε_m－ε_t）和R_mcos ε_m（β_t－β_m）。

仿真參數(shù)和改進視線制導(dǎo)律的設(shè)計參數(shù)與第3.1節(jié)相同，傳統(tǒng)視線制導(dǎo)律的參數(shù)設(shè)計為k_p=6，k_i=3，k_d=4。為體現(xiàn)載機運動方向變化的影響，在仿真過程中不引入視線角噪聲且對導(dǎo)彈矢徑不進行拉偏，視線角速度仍采用前述跟蹤微分器進行估計。假設(shè)載機保持固定高度飛行，運動方向的改變通過其航向偏角的變化實現(xiàn)，取值范圍從－135°開始，逐次遞增45°，直至180°，仿真結(jié)果如下。

設(shè)定脫靶量指標(biāo)為1 m，圖10為兩種視線制導(dǎo)律在目標(biāo)靜止情況下的彈道曲線，其中實線表示滿足指標(biāo)要求，而雙劃線則表示不滿足，表3為其對應(yīng)的脫靶量數(shù)值。通過觀察可以發(fā)現(xiàn)，對于靜止目標(biāo)，改進視線制導(dǎo)律不受載機運動方向的影響，能夠在各個方向上精準(zhǔn)命中目標(biāo)。而在傳統(tǒng)視線制導(dǎo)律下，載機運動方向偏離目標(biāo)時的命中精度較低。

目標(biāo)做勻速直線運動時，彈道曲線如圖11所示，對應(yīng)的脫靶量數(shù)值如表4所示。目標(biāo)的航向偏角為45°，在傳統(tǒng)視線制導(dǎo)律下，只有當(dāng)載機的航向偏角為0°和45°時才能滿足脫靶量指標(biāo)，表明載機運動方向應(yīng)與目標(biāo)運動方向近似一致。然而，通過應(yīng)用本文提出的改進視線制導(dǎo)律，即使在載機和目標(biāo)的運動方向不一致甚至相反時，仍能夠達到脫靶量指標(biāo)。

圖12所示為兩種制導(dǎo)律分別在不同的載機航向偏角下攻擊相同的圓弧機動目標(biāo)，表5為其對應(yīng)的脫靶量數(shù)值。與前兩種目標(biāo)相似，在目標(biāo)圓弧機動情況下，傳統(tǒng)視線制導(dǎo)律只有當(dāng)載機朝著目標(biāo)運動時才能滿足脫靶量指標(biāo)，而改進視線制導(dǎo)律則不受此限制。

總體而言，利用本文提出的改進視線制導(dǎo)律，可以有效解決傳統(tǒng)視線制導(dǎo)律在載機遠離目標(biāo)運動時命中精度較低甚至脫靶的問題，解除了制導(dǎo)律在實施過程中對載機運動方向的限制。

4 結(jié) 論

針對機載武器低成本精確打擊的需求，本文提出在不額外增設(shè)制導(dǎo)站的前提下充分利用載機原有的視線角信息進行空地導(dǎo)彈改進視線制導(dǎo)律的設(shè)計?；谌S空間模型建立了系統(tǒng)狀態(tài)空間方程，將導(dǎo)彈加速度作為控制量、導(dǎo)彈偏離機目視線的位置作為輸出量，通過設(shè)計控制量使輸出量漸進穩(wěn)定地趨于零，從而實現(xiàn)視線制導(dǎo)。在此基礎(chǔ)上，基于反步法和滑?？刂圃O(shè)計了反步滑模視線制導(dǎo)律，并考慮了制導(dǎo)律在應(yīng)用過程中的實施問題。最后，在六自由度彈道中進行了仿真驗證。在同時存在噪聲干擾及導(dǎo)彈矢徑拉偏的條件下，改進視線制導(dǎo)律可以命中靜止、勻速直線運動及圓弧機動的地面目標(biāo)。與傳統(tǒng)視線制導(dǎo)律相比，本文的改進視線制導(dǎo)律不對載機運動方向設(shè)限，在制導(dǎo)過程中載機無需向著目標(biāo)飛行，提高了載機的安全性，具有實際的工程研究意義。

參考文獻

［1］任新聯(lián)， 李剛. 美國無人機載空地導(dǎo)彈發(fā)展現(xiàn)狀及趨勢［J］. 飛航導(dǎo)彈， 2021（4）： 59-62， 67.

REN X L， LI G. Development status and trends of UAV-borne air-to-ground missiles in US［J］. Aerodynamic Missile Journal， 2021（4）： 59-62， 67.

［2］CHAVEZ K. Learning on the fly： drones in the Russian-Ukrainian war［J］. Arms Control Today， 2023， 53（1）： 6-11.

［3］YANG K， XU C， XU Y. Design and implementation of trajectory planning algorithm for UAV borne air-to-surface missile［J］. Journal of Physics： Conference Series， 2021， 1721（1）： 012038.

［4］杜澤弘， 樓俏， 賀敏， 等. 無人機機載武器發(fā)展現(xiàn)狀及展望［J］. 軍事文摘， 2020（5）： 28-32.

DU Z H， LOU Q， HE M， et al. Current status and outlook of UAV airborne weapons development［J］. Military Digest， 2020（5）： 28-32.

［5］QI Z K， LIN D F. Design of guidance and control systems for tactical missiles［M］. Beijing： Beijing Institute of Technology Press， 2019.

［6］王狂飆， 張?zhí)鞓? 激光駕束制導(dǎo)的仿真研究［J］. 南京理工大學(xué)學(xué)報， 2001， 25（2）： 125-128.

WANG K B， ZHANG T Q. The simulation study on laser beam riding guidance［J］. Journal of Nanjing University of Science and Technology， 2001， 25（2）： 125-128.

［7］ZARCHAN P. Tactical and strategic missile guidance［M］. 6th ed. Reston： AIAA， 2012.

［8］RATNOO A. Three-point guidance for intercepting weaving targets［J］. Journal of Guidance， Control， and Dynamics， 2016， 39（8）： 1879-1884.

［9］RATNOO A， SHIMA T. Line-of-sight interceptor guidance for defending an aircraft［J］. Journal of Guidance， Control， and Dynamics， 2011， 34（2）： 522-532.

［10］YAMASAKI T， BALAKRISHNAN S N. Triangle intercept guidance for aerial defense［C］∥Proc.of the AIAA Guidance， Navigation， and Control Conference， 2010.

［11］YAMASAKI T， BALAKRISHNAN S N. Intercept guidance for cooperative aircraft defense against a guided missile［J］. IFAC Proceedings Volumes， 2010， 43（15）： 118-123.

［12］AHI B， HAERI M. Novel command to line of sight guidance with practical limitations［J］. Asian Journal of Control， 2022， 24（3）： 1426-1436.

［13］HA I J， CHONG S. Design of a CLOS guidance law via feedback linearization［J］. IEEE Trans.on Aerospace and Electronic Systems， 1992， 28（1）： 51-63.

［14］SADEGHINASAB E， KOOFIGAR H R， ATAEI M. Design of robust command to line-of-sight guidance law： a fuzzy adaptive approach［J］. Journal of Engineering Science and Technology， 2016， 11（11）： 1528-1542.

［15］ELHALWAGY Y Z， TARBOUCHI M. Fuzzy logic sliding mode control for command guidance law design［J］. ISA Trans-actions， 2004， 43（2）： 231-242.

［16］ZAIDI H， WU P L， BELLAHCENE A. Missile guidance law design via backstepping technique［J］. International Journal of Engineering and Applied Sciences， 2016， 3（4）： 85-90.

［17］MA’ARIF A， VERA M A M， MAHMOUD M S， et al. Backstepping sliding mode control for inverted pendulum system with disturbance and parameter uncertainty［J］. Journal of Robotics and Control， 2022， 3（1）： 86-92.

［18］JIA Z Y， YU J Q， MEI Y S， et al. Integral backstepping sliding mode control for quadrotor helicopter under external uncertain disturbances［J］. Aerospace Science and Technology， 2017， 68： 299-307.

［19］HUANG S， HUANG J F， CAI Z Q， et al. Adaptive backstepping sliding mode control for quadrotor UAV［J］. Scientific Programming， 2021， 2021： 3997648.

［20］SHI Q L， WANG H， CHENG H. Multiple constraints-based adaptive three-dimensional back-stepping sliding mode guidance law against a maneuvering target［J］. Aerospace， 2022， 9（12）： 796-817.

［21］XIA Y Q， LU K F， ZHU Z， et al. Adaptive back-stepping sliding mode attitude control of missile systems［J］. International Journal of Robust and Nonlinear Control， 2013， 23（15）： 1699-1717.

［22］ZHU C Q， GUO Z Y. Design of head-pursuit guidance law based on backstepping sliding mode control［J］. International Journal of Aerospace Engineering， 2019， 2019： 8214042.

［23］MALEY J M. Line of sight rate estimation for guided projectiles with strapdown seekers［C］∥Proc.of the AIAA Guidance， Navigation， and Control Conference， 2015： 749-772.

［24］GONTUMUKKALA S S T， GODAVARTHI Y S V， GONUGUNTA B R R T， et al. Kalman filter and proportional navigation based missile guidance system［C］∥Proc.of the IEEE 8th International Conference on Advanced Computing and Communication Systems， 2022： 1731-1736.

［25］CHEN T L， ZHANG C， SUN T， et al. Research on guidance information extraction technology of strap-down seeker based on adaptive cubature Kalman filter［J］. Journal of Physics： Conference Series， 2023， 2558（1）： 012018.

［26］韓京清， 王偉. 非線性跟蹤─微分器［J］. 系統(tǒng)科學(xué)與數(shù)學(xué)， 1994， 14（2）： 177-183.

HAN J Q， WANG W. Nonlinear tracking-differentiator［J］. Journal of Systems Science and Mathematical Sciences， 1994， 14（2）： 177-183.

［27］WANG P， ZHANG K. Research on line-of-sight rate extraction of strapdown seeker［C］∥Proc.of the IEEE 33rd Chinese Control Conference， 2014： 859-863.

［28］WEN Q Q， WANG M K， LI X， et al. Learning-based design optimization of second-order tracking differentiator with application to missile guidance law［J］. Aerospace Science and Technology， 2023， 137： 108302.

［29］周運強. 空地導(dǎo)彈先進制導(dǎo)律研究［D］. 北京： 北京理工大學(xué)， 2016.

ZHOU Y Q. Research on advanced guidance law of air-to-ground missiles［D］. Beijing： Beijing Institute of Technology， 2016.

［30］吳光輝， 方東洋， 陳意芬， 等. 基于卡爾曼濾波器的激光駕束制導(dǎo)導(dǎo)彈制導(dǎo)方法研究［J］. 彈箭與制導(dǎo)學(xué)報， 2022， 42（3）： 42-46.

WU G H， FANG D Y， CHEN Y F， et al. Research on guidance method of laser beam riding guided missile based on Kalman filter［J］. Journal of Projectiles， Rockets， Missiles and Guidance， 2022， 42（3）： 42-46.

作者簡介

張文穩(wěn)（2000—），男，碩士研究生，主要研究方向為飛行器制導(dǎo)與控制。

張 成（1973—），男，研究員，博士研究生導(dǎo)師，博士，主要研究方向為飛行器制導(dǎo)與控制、飛行器視覺與人工智能。

鄭晨明（1993—），男，博士研究生，主要研究方向為飛行器制導(dǎo)與控制。

程潤北（2000—），男，碩士研究生，主要研究方向為飛行器制導(dǎo)與控制。

陳天樂（1999—），男，助理工程師，碩士，主要研究方向為飛行器制導(dǎo)與控制。