CN113536707A - Aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression - Google Patents

Aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression Download PDF

Info

Publication number
CN113536707A
CN113536707A CN202110828849.8A CN202110828849A CN113536707A CN 113536707 A CN113536707 A CN 113536707A CN 202110828849 A CN202110828849 A CN 202110828849A CN 113536707 A CN113536707 A CN 113536707A
Authority
CN
China
Prior art keywords
aircraft
error
angle
radome
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202110828849.8A
Other languages
Chinese (zh)
Other versions
CN113536707B (en
Inventor
陆科林
符启恩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Southeast University
Original Assignee
Southeast University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Southeast University filed Critical Southeast University
Priority to CN202110828849.8A priority Critical patent/CN113536707B/en
Publication of CN113536707A publication Critical patent/CN113536707A/en
Application granted granted Critical
Publication of CN113536707B publication Critical patent/CN113536707B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention discloses an aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression, which is characterized in that an aircraft dynamics model is established for an aircraft guidance system which is only subjected to line-of-sight angle measurement and is influenced by radome errors, each state in a guidance process is estimated by using an interactive multi-model filtering method, a mapping relation between an aircraft visual angle and a radome error angle is established based on the Gaussian process regression model, an analytic expression form of the radome error slope is obtained based on a derivative of the Gaussian process model, and the analytic expression form is finally used for the compensation of the radome errors of the aircraft seeker, so that the guidance performance is effectively improved.

Description

Aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression
Technical Field
The invention relates to the field of aircraft guidance, in particular to an aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression.
Background
In aircraft guidance missions, radomes are often used to protect aircraft seeker antennas from air currents while reducing aerodynamic drag on aircraft flight. However, the radome refracts electromagnetic wave signals entering the seeker, so that the line-of-sight angle measurement obtained by the seeker deviates, parasitic loops are generated in the aircraft guidance system, and the guidance system is unstable. The method for reducing the influence of errors brought by the antenna housing is mainly divided into three categories, the first category is a hardware-based compensation method, namely, the process processing is carried out by the method of inner profile surface grinding and the like when the antenna housing is manufactured; the second type is that the error angle and the error slope of the antenna housing are directly measured and directly compensated during guidance; the last type is estimation and compensation of the error of the antenna housing on the algorithm level, and the traditional method comprises compensation by using a jitter signal and a low-pass filter, online measurement compensation based on Kalman filtering and multi-model filtering methods, compensation of the error angle of the antenna housing based on a neural network and the like. In recent years, the research on the antenna housing error compensation algorithm is further promoted by the development of technologies such as machine learning and adaptive filtering.
A Loop-Shaping method based on the Radome error compensation method is disclosed in "Loop-Shaping Approach to Mitigate antenna errors in Home issues", which is disclosed in Journal of guidelines, Control, and Dynamics (Klein, I.and Russak, I.,2017.Loop-Shaping Approach to mit antenna errors in homes, Journal of guidelines, Control, and Dynamics,40(7), pp.1789-1795.) to reduce the influence of Radome errors and improve the stability margin of the system by adding a phase lead compensation Loop between the Guidance law and the flight Control system in the Guidance Loop. Time-Varying radar Slope Estimation for Passive antenna cancellation-Ship mismatches, IEEE 58 Conference on Decision and Control (CDC) (Ra, w.s., Ahn, s., Lee, y.and whisg, i.h.,2019, Decumber. Time-Varying radar Slope Estimation for Passive antenna cancellation-Ship mismatches. in 2019 IEEE 58 Conference on Decision and Control (CDC) (pp.4940-4945) discloses a method for compensating for antenna cover errors based on a dither signal and a Time Varying kalman filter, wherein a bandpass filter is used to extract the influence of the dither signal on the line angle, and the Estimation and compensation of errors are combined with filtering. An antenna housing aiming error compensation method based on an EKF technology is disclosed in the systematic simulation science newspaper (Zhou di yearn, Li Junlong, Yuanyuqi, a radar seeker antenna housing slope error real-time estimation method [ J ]. modern defense technology 2020,48(05):1-9.) and an antenna housing error compensation method based on an extended Kalman filter is disclosed, the antenna housing error slope is established as a system state, and estimation is carried out by using an extended Kalman filtering algorithm, so that further compensation is carried out. However, the method based on the traditional filtering and control theory needs to accurately model the system model, the error of the model can generate great influence on the guidance effect, and an accurate dynamic model is often difficult to obtain in practice.
In recent years, a radome error estimation and compensation method based on data has been developed. Adaptive Scale Factor Compensation for satellites with satellite vectors, and arXiv (Gaudet, b.,2020.Adaptive Scale Factor Compensation for satellites with satellite vectors Predictive Coding. arXiv Predictive Coding: 2009.00975.) disclose a Predictive Coding method based on a Predictive Coding method, which predicts a radome error angle using a recurrent neural network, adaptively corrects an observation signal, and compensates in an aircraft guidance system. However, the above method has the disadvantage that it requires a large amount of computation and is difficult to implement well in an aircraft guidance mission with a very high real-time requirement.
The Gaussian process is a non-parameterized machine learning model, and compared with a neural network, the posterior covariance obtained in the Gaussian process model prediction can be used as the measurement of the model accuracy, and the model is flexibly applied to model application and has natural advantages. Meanwhile, the data required for training the Gaussian process model is relatively less. A nonlinear filtering method based on Gaussian process is disclosed in IEEE Transactions on Automatic Control (J.Pruher and O.Straka, "Gaussian process motion transform," IEEE Transactions on Automatic Control, vol.63, No.9, pp.2844-2854,2017.), fitting a system dynamics model and a measurement model by using the Gaussian process, performing a statistical moment conversion step in the filtering process based on an identified posterior model of the Gaussian process, and combining a traditional Bayesian filtering framework to be applied to tracking a moving target. However, the method and other control and filtering methods based on the Gaussian process are limited to algorithm improvement, and the combination of the method and the control and filtering method applied to the field of aircraft guidance is not seen.
In the existing error estimation and compensation method for the aircraft guidance antenna housing, a hardware-based method is limited by a process level and needs to be balanced with the protection effect of the antenna housing, and the measurement-based method has high requirements on a sensor for accurately measuring the error of the antenna housing; the methods for compensating at the algorithm level are mostly limited by the accuracy of the system model, and the methods based on data are relatively few and difficult to meet the real-time requirement. At present, no radome error compensation method capable of well applying an unparameterized Bayesian machine learning method such as a Gaussian process model exists.
Disclosure of Invention
The purpose of the invention is as follows: the invention aims to provide an aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression, which is accurate and effective and can carry out online estimation and real-time compensation on seeker radome errors.
The technical scheme is as follows: the invention relates to an aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression, which comprises the following steps:
(1) establishing an aircraft guidance dynamic model and a measurement model influenced by a seeker antenna housing error;
(2) estimating each state in the guidance process of the aircraft based on an interactive multi-model filtering method;
(3) based on a Gaussian process regression method, establishing a mapping relation between an aircraft view angle and an antenna housing error angle, and obtaining an antenna housing error slope;
(4) and compensating the error slope of the antenna housing in the process of manufacturing and guiding based on the established Gaussian process regression model.
Further, the step (1) includes the steps of:
(11) establishing an aircraft guidance dynamics model
The aircraft guidance task is guided by using a proportional guidance method, and a first-order flight control system is configured. The target line-of-sight angle measured by the seeker of the aircraft can generate certain deviation under the influence of the seeker radome, namely a radome error angle, and the expression of the error angle is
Figure BDA0003174750030000035
Figure BDA0003174750030000031
Wherein λ isrIs the radome error angle, thetasIs the angle of view, λ is the target line of view angle, θMIs the projectile attitude angle.
Defining the radome error slope as
Figure BDA0003174750030000032
A geometric diagram of an aircraft seeker with radome error is shown in fig. 1. Can establish a system dynamics model as
Figure BDA0003174750030000033
Wherein
Figure BDA0003174750030000034
Is the system state, R is the distance between the aircraft and the target, γMIs the aircraft flight path angle, AMIn order to be the actual guidance instruction of the aircraft,
Figure BDA0003174750030000041
for the kinetic equation, it was constructed as follows
Figure BDA0003174750030000042
Figure BDA0003174750030000043
Figure BDA0003174750030000044
Figure BDA0003174750030000045
Wherein VMFor aircraft speed, N is the proportional guidance factor, τ is the autopilot time constant, TαIs the time constant of the rate of change of heading,
Figure BDA0003174750030000046
the estimated value of the error slope of the antenna housing is obtained. The aircraft guidance loop including the dynamics model is shown in fig. 2.
(12) Establishing an aircraft guidance measurement model
The aircraft only has a target line-of-sight angle measured value, and a measurement model is established by considering the influence of an error angle of the antenna housing
Figure BDA0003174750030000047
Wherein v iskN (0, R) is measurement noise, h (x)k;ρθ,k) To be a measurement equation
Figure BDA0003174750030000048
Further, the step (2) comprises the steps of:
(21) discretized system dynamics model
Discretizing the aircraft guidance dynamic model established in the step (1) based on a four-order Runge Kutta method, namely discretizing the aircraft guidance dynamic model
xk+1=φ(xk;Δt,ρθ,k)+wk
Wherein wkN (0, Q) is the process noise, representing the discretization error, and Q is its covariance matrix.
(22) Setting local filtering model
Building multiple local filtering models using a set of preset radome error slope values, i.e.
Figure BDA0003174750030000049
For each error slope parameter value
Figure BDA00031747500300000410
Establishing corresponding dynamic and measurement model
Figure BDA0003174750030000051
Figure BDA0003174750030000052
And carrying out filtering estimation on each model by using an unscented Kalman filtering algorithm.
(23) Hybrid local estimation results
Calculating a mixed estimation mean and covariance as estimation results of the multi-model filtering method based on the local model filtering results obtained in the step (22), that is
Figure BDA0003174750030000053
Figure BDA0003174750030000054
Wherein
Figure BDA0003174750030000055
And Pk|kIs the k stepThe estimated mean and covariance of the time,
Figure BDA0003174750030000056
and
Figure BDA0003174750030000057
for the estimated value obtained by the ith local filter,
Figure BDA0003174750030000058
and the model probability corresponding to the ith local model. A schematic diagram of the interactive multi-model filtering method is shown in fig. 3.
(24) Calculating estimated view angle and radome error
According to the multi-model filtering estimation result obtained in the step (23)
Figure BDA0003174750030000059
Computing an estimated view angle
Figure BDA00031747500300000510
As follows
Figure BDA00031747500300000511
Calculating and estimating antenna housing error angle
Figure BDA00031747500300000512
As follows
Figure BDA00031747500300000513
Further, the step (3) includes the steps of:
(31) establishing a Gaussian process regression model from a visual angle to an antenna housing error angle
Using vr,i~N(0,rr,i) Representing the estimation error between the estimated antenna housing error angle and the actual antenna housing error angle in the ith step to obtain the relation between the estimated antenna housing error angle and the estimated visual angle
Figure BDA00031747500300000514
Figure BDA00031747500300000515
Considering arbitrary perspective input
Figure BDA00031747500300000516
The corresponding required predicted error angle of the antenna housing is
Figure BDA00031747500300000517
Establishing
Figure BDA00031747500300000518
And λr,jGaussian process prior distribution between
Figure BDA0003174750030000061
Wherein
Figure BDA0003174750030000062
Is a prior mean, K is a covariance matrix, formed by a covariance function K (x)1,x2) The components of the composition are as follows,
Figure BDA0003174750030000063
based on training data of the estimated radome error angle and the estimated view angle, the posterior distribution of the radome error angle to be predicted can be obtained
Figure BDA0003174750030000064
Wherein
Figure BDA0003174750030000065
Figure BDA0003174750030000066
(32) Calculating radome error slope
The obtained posterior distribution of the Gaussian process is subjected to derivation to obtain any input
Figure BDA0003174750030000067
The error slope of the radome is
Figure BDA0003174750030000068
Wherein
Figure BDA0003174750030000069
Further, the step (4) comprises the steps of:
(41) calculating a corrected line-of-sight angular rate
Constructing a corrected line-of-sight angular rate according to the estimated radome error slope obtained in the step (3)
Figure BDA00031747500300000610
Based on the relationship between the measured line-of-sight angular rate and the actual line-of-sight angular rate
Figure BDA00031747500300000611
Obtaining a corrected line-of-sight angular rate of
Figure BDA00031747500300000612
The aircraft guidance loop after line-of-sight angular rate correction is shown in FIG. 4.
(42) Calculating and correcting actual guidance instruction of aircraft
Based on the estimated value of each state of the guidance system obtained in the step (2)
Figure BDA00031747500300000613
Calculating to obtain a corrected aircraft actual guidance instruction in discrete time
Figure BDA00031747500300000614
Wherein
Figure BDA0003174750030000071
And using the corrected guidance instruction to complete the compensation of the error slope of the antenna housing. An overall schematic of the estimation and compensation of the radome error slope is shown in fig. 5.
Has the advantages that: compared with a method for designing and polishing from a hardware level, the method does not need to rely on a process level and does not need to consider the balance with the radome protection effect, compared with a method for directly measuring the radome error, the method does not need to be provided with such a sensor with strict requirements, and is easier to realize.
Drawings
FIG. 1 is a geometric block diagram of an aircraft seeker;
FIG. 2 is a diagram of an uncompensated aircraft guidance loop using proportional guidance;
FIG. 3 is a schematic diagram of an interactive multi-model filtering algorithm used in the present invention;
FIG. 4 is a diagram of an aircraft guidance loop after line-of-sight angular rate correction;
FIG. 5 is a general schematic diagram of a method for estimating and compensating an error slope of an antenna radome in accordance with the present invention;
FIG. 6 is a schematic diagram of a mean square error estimation of the line-of-sight angle of an aircraft guidance system;
FIG. 7 is a schematic diagram of an aircraft guidance system aircraft relative distance estimation mean square error with a target;
FIG. 8 is a schematic diagram of a mean square error estimation of a flight path angle for an aircraft guidance system;
FIG. 9 is a schematic diagram of an estimated mean square error of an actual guidance instruction for an aircraft guidance system;
FIG. 10 is a schematic diagram of a target view true trajectory and an estimated trajectory for an aircraft guidance system;
FIG. 11 is a schematic diagram of an actual trajectory and an estimated trajectory of an error angle of an antenna housing of an aircraft guidance system;
fig. 12 is a schematic diagram of an error slope true trajectory of an antenna cover and an estimated trajectory based on gaussian process regression;
fig. 13 is a gaussian process regression-based radome error slope estimation mean square error;
FIG. 14 is a graph comparing the final guidance miss distance of the proposed method with other compensation methods;
Detailed Description
The technical scheme of the invention is further described in the following by combining the attached drawings and the detailed description.
Initial values considering the aircraft and target states are as follows
Figure BDA0003174750030000081
Figure BDA0003174750030000082
Wherein (X)M,YM) And (X)T,YT) Is the initial position of the aircraft with respect to the target,
Figure BDA0003174750030000083
and
Figure BDA0003174750030000084
for an initial velocity of the aircraft with the target, the target is assumed to be a fixed target. Obtaining the initial state of the guidance system according to the establishment of the dynamic model in the step (2) as follows
Figure BDA0003174750030000085
The corresponding model parameters are as follows
{VM,N,ρθ,τ,Tα}={500m/s,4,0.025°/°,0.1s,1s}
In the interactive multi-model filter, the discretization time interval is set to Δ t 0.001s, and the estimated mean and variance of the initial state are
Figure BDA0003174750030000086
Figure BDA0003174750030000087
Where n-3 is the number of guessed models, the transition probability matrix for the Markov chain is as follows
Figure BDA0003174750030000088
In addition, process noise wkHas a covariance matrix Q of
Figure BDA0003174750030000089
Observation noise vkVariance of (R) 1.74532×10-12rad2Is as follows.
Based on the above settings, the states in the aircraft guidance process are estimated according to the steps (1) to (4) and the method shown in fig. 5, and the radome error slope is estimated and compensated. In the gaussian process regression, the prediction of the radome error angle and the estimation of the radome error slope are performed in two ways, namely prediction by using all historical data (full history) and prediction by using sliding window historical data (sliding window), respectively. Fig. 6-9 show the mean square error of each state estimation of the guidance system, fig. 10-11 show the real track and the estimated track of the target view angle and the antenna housing error angle, fig. 12 shows the real value and the estimated value of the antenna housing error slope, fig. 13 shows the mean square error of the antenna housing error slope estimation, it can be seen that better effect can be obtained by using sliding window data to perform gaussian process prediction, and fig. 14 shows the final guidance miss distance after the gaussian process estimation and compensation method provided by the invention is adopted, and the comparison of the guidance miss distance based on the traditional multi-model filter compensation (IMM) and the extended kalman filter compensation (EKF) is carried out under the condition of no compensation. The result shown in the attached drawing shows that the method provided by the invention can effectively improve the estimation precision of each state of the guidance system, can effectively estimate the error slope of the antenna housing, and can obtain a better guidance effect compared with other compensation methods.

Claims (5)

1. An aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression; the method is characterized in that: the method comprises the following steps:
(1) establishing a guidance dynamic model and a measurement model of the aircraft influenced by the error of the radome of the seeker;
(2) estimating each state in the guidance process of the aircraft based on an interactive multi-model filtering method;
(3) based on a Gaussian process regression method, establishing a mapping relation between an aircraft view angle and an antenna housing error angle, and obtaining an antenna housing error slope;
(4) and compensating the error slope of the antenna housing in the process of manufacturing and guiding based on the established Gaussian process regression model.
2. The method for estimating and compensating the error slope of the radome of the aircraft seeker based on the gaussian process regression as recited in claim 1, wherein: the step (1) comprises the following steps:
(11) establishing an aircraft guidance dynamics model
The aircraft guidance task is guided by using a proportional guidance method, a first-order flight control system is configured, a target line-of-sight angle measured by a guidance head of the aircraft is influenced by a guidance head antenna housing to generate a certain deviation, namely an antenna housing error angle, and the expression of the error angle is an antenna housing error angle
Figure FDA0003174750020000011
Figure FDA0003174750020000012
Wherein λ isrIs the radome error angle, thetasIs the angle of view, λ is the target line of view angle, θMIs the projectile attitude angle. Defining the radome error slope as
Figure FDA0003174750020000013
Can establish an aircraft guidance dynamic model as
Figure FDA0003174750020000014
Wherein
Figure FDA0003174750020000015
Is the system state, R is the distance between the aircraft and the target, γMIs the aircraft flight path angle, AMIn order to be the actual guidance instruction of the aircraft,
Figure FDA0003174750020000016
as a kinetic squareThe process is constructed as follows
Figure FDA0003174750020000017
Figure FDA0003174750020000018
Figure FDA0003174750020000019
Figure FDA0003174750020000021
Wherein VMFor aircraft speed, N is the proportional guidance factor, τ is the autopilot time constant, TαIs the time constant of the rate of change of heading,
Figure FDA0003174750020000022
Figure FDA0003174750020000023
the estimated value of the error slope of the antenna housing is obtained.
(12) Establishing an aircraft guidance measurement model
The aircraft only has a target line-of-sight angle measured value, and the aircraft guidance measurement model is established by considering the influence of an antenna housing error angle as follows:
Figure FDA0003174750020000024
wherein v iskN (0, R) is measurement noise, h (x)k;ρθ,k) To be a measurement equation
Figure FDA0003174750020000025
3. The method for estimating and compensating the error slope of the radome of the aircraft seeker based on the gaussian process regression as recited in claim 1, wherein: the step (2) comprises the following steps:
(21) discretized system dynamics model
Discretizing the aircraft guidance dynamic model established in the step (1) based on a four-order Runge Kutta method, namely discretizing the aircraft guidance dynamic model
xk+1=φ(xk;Δt,ρθ,k)+wk
Wherein wkN (0, Q) is the process noise, representing the discretization error, and Q is its covariance matrix.
(22) Setting local filtering model
Building multiple local filtering models using a set of preset radome error slope values, i.e.
Figure FDA0003174750020000026
For each error slope parameter value
Figure FDA0003174750020000027
Establishing corresponding dynamic and measurement model
Figure FDA0003174750020000028
Figure FDA0003174750020000029
And carrying out filtering estimation on each model by using an unscented Kalman filtering algorithm.
(23) Hybrid local estimation results
Calculating a mixed estimation mean and covariance as estimation results of the multi-model filtering method based on the local model filtering results obtained in the step (22), that is
Figure FDA0003174750020000031
Figure FDA0003174750020000032
Wherein
Figure FDA0003174750020000033
And Pk|kAs the estimated mean and covariance at step k,
Figure FDA0003174750020000034
and
Figure FDA0003174750020000035
for the estimated value obtained by the ith local filter,
Figure FDA0003174750020000036
and the model probability corresponding to the ith local model.
(24) Calculating estimated view angle and radome error
According to the multi-model filtering estimation result obtained in the step (23)
Figure FDA0003174750020000037
Computing an estimated view angle
Figure FDA0003174750020000038
As follows
Figure FDA0003174750020000039
Calculating and estimating antenna housing error angle
Figure FDA00031747500200000310
As follows
Figure FDA00031747500200000311
4. The method for estimating and compensating the error slope of the radome of the aircraft seeker based on the gaussian process regression as recited in claim 1, wherein: the step (3) comprises the following steps:
(31) establishing a Gaussian process regression model from a visual angle to an antenna housing error angle
Using vr,i~N(0,rr,i) Representing the estimation error between the estimated antenna housing error angle and the actual antenna housing error angle in the ith step to obtain the relation between the estimated antenna housing error angle and the estimated visual angle
Figure FDA00031747500200000312
Figure FDA00031747500200000313
Considering arbitrary perspective input
Figure FDA00031747500200000314
The corresponding required predicted error angle of the antenna housing is
Figure FDA00031747500200000315
Establishing
Figure FDA00031747500200000316
And λr,jGaussian process prior distribution between
Figure FDA00031747500200000317
Wherein
Figure FDA00031747500200000318
Is a prior mean, K is a covariance matrix, formed by a covariance function K (x)1,x2) The components of the composition are as follows,
Figure FDA00031747500200000319
based on training data of the estimated radome error angle and the estimated view angle, the posterior distribution of the radome error angle to be predicted can be obtained
Figure FDA0003174750020000041
Wherein
Figure FDA0003174750020000042
Figure FDA0003174750020000043
(32) Calculating radome error slope
The obtained posterior distribution of the Gaussian process is subjected to derivation to obtain any input
Figure FDA0003174750020000044
The error slope of the radome is
Figure FDA0003174750020000045
Wherein
Figure 1
5. The method for estimating and compensating the error slope of the radome of the aircraft seeker based on the gaussian process regression as recited in claim 1, wherein: the step (4) comprises the following steps:
(41) calculating a corrected line-of-sight angular rate
Constructing a corrected line-of-sight angular rate according to the estimated radome error slope obtained in the step (3)
Figure FDA0003174750020000047
Based on the relationship between the measured line-of-sight angular rate and the actual line-of-sight angular rate
Figure FDA0003174750020000048
Obtaining a corrected line-of-sight angular rate of
Figure FDA0003174750020000049
(42) Calculating and correcting actual guidance instruction of aircraft
Based on the estimated value of each state of the guidance system obtained in the step (2)
Figure FDA00031747500200000410
Calculating to obtain a corrected aircraft actual guidance instruction in discrete time
Figure FDA00031747500200000411
Wherein
Figure FDA00031747500200000412
And using the corrected guidance instruction to complete the compensation of the error slope of the antenna housing.
CN202110828849.8A 2021-07-22 2021-07-22 Aircraft seeker antenna cover error slope estimation and compensation method based on Gaussian process regression Active CN113536707B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110828849.8A CN113536707B (en) 2021-07-22 2021-07-22 Aircraft seeker antenna cover error slope estimation and compensation method based on Gaussian process regression

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110828849.8A CN113536707B (en) 2021-07-22 2021-07-22 Aircraft seeker antenna cover error slope estimation and compensation method based on Gaussian process regression

Publications (2)

Publication Number Publication Date
CN113536707A true CN113536707A (en) 2021-10-22
CN113536707B CN113536707B (en) 2024-03-01

Family

ID=78120355

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110828849.8A Active CN113536707B (en) 2021-07-22 2021-07-22 Aircraft seeker antenna cover error slope estimation and compensation method based on Gaussian process regression

Country Status (1)

Country Link
CN (1) CN113536707B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3940767A (en) * 1955-01-21 1976-02-24 Hughes Aircraft Company Electronic radome-error compensation system
CN106507895B (en) * 2010-11-29 2014-06-04 中国空空导弹研究院 A kind of seeker antenna cover collimating fault compensation method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3940767A (en) * 1955-01-21 1976-02-24 Hughes Aircraft Company Electronic radome-error compensation system
CN106507895B (en) * 2010-11-29 2014-06-04 中国空空导弹研究院 A kind of seeker antenna cover collimating fault compensation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KELIN LU等: "Gaussian process-based Bayesian non-linear filtering for online target tracking", IET RADAR SONAR NAVIG, vol. 14, no. 3, pages 448 - 458, XP006089200, DOI: 10.1049/iet-rsn.2019.0495 *
许海深: "导弹导引头天线罩误差斜率补偿研究", 中国优秀硕士学位论文全文数据库(电子期刊), no. 1, pages 032 - 190 *

Also Published As

Publication number Publication date
CN113536707B (en) 2024-03-01

Similar Documents

Publication Publication Date Title
CN108036676B (en) A kind of autonomous reentry guidance method of full directive based on three-dimensional resolution Value of Reentry Vehicle
Lee et al. Polynomial guidance laws considering terminal impact angle and acceleration constraints
Wang et al. Huber-based unscented filtering and its application to vision-based relative navigation
Shin et al. Adaptive support vector regression for UAV flight control
CN110347170B (en) Reusable carrier reentry segment robust fault-tolerant guidance control system and working method
CN105222648A (en) A kind of linear pseudo-spectrum broad sense mark control miss distance Guidance and control method
CN108073742B (en) Method for estimating flight state of intercepted missile tail section based on improved particle filter algorithm
CN111913175A (en) Water surface target tracking method with compensation mechanism under transient failure of sensor
Zhao et al. ADRC based integrated guidance and control scheme for the interception of maneuvering targets with desired los angle
CN109062241A (en) Autonomous full directive reentry guidance method based on linear pseudo- spectrum model PREDICTIVE CONTROL
CN114740497A (en) Unmanned aerial vehicle deception method based on UKF multi-source fusion detection
Liao et al. Guidance laws with input saturation and nonlinear robust H∞ observers
Yang et al. Online midcourse guidance method for boost phase interception via adaptive convex programming
Guo et al. Qualitative indicator-based guidance scheme for bank-to-turn missiles against couplings and maneuvering targets
CN113536707A (en) Aircraft seeker radome error slope estimation and compensation method based on Gaussian process regression
Yuqi et al. Time-varying parameters estimation with adaptive neural network EKF for missile-dual control system
CN113110428A (en) Carrier-based aircraft landing fixed time trajectory tracking method based on limited backstepping control
CN114660587A (en) Jump and glide trajectory target tracking method and system based on Jerk model
Léchevin et al. Backstepping guidance for missiles modeled as uncertain time-varying first-order systems
Gong et al. Unscented particle smoother and its application to transfer alignment of airborne distributed POS
CN112595319A (en) Model adaptive compensation return trajectory estimation algorithm
Yang et al. Design of Nonsingular Terminal Sliding Mode Guidance Law with Disturbance Observer
Feng et al. Trajectory tracking for hypersonic glide vehicles based on improved sine-AIMM
Wan et al. A Midcourse Guidance Method Combined with Trajectory Prediction for Antinear‐Space‐Gliding Vehicles
He et al. Global finite-time PLOS-OCCs guidance and adaptive integral sliding mode path following control for unmanned surface vehicles with ocean currents and input saturation

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant