■п О
Розглянуто розв'язання термталь-ног задачi управлтня та синтезований параметризований закон управлтня в аналтичному виглядi, який залежить вид змтного параметра глибини прогнозу. Дослиджено особливостi впли-ву величини параметра управлтня на точтсть досягнення ктцевог точки, Ваш рекомендаци з вибору параметра для твелювання помилки терцшних вимiрювань. Синтез управлтня здшс-нюеться методом переслидування ведучог точки за тформащею, отри-маною ттегруванням вимiрювань фактичного прискорення i метить помилку, характерну для акселеро-метрiв
Ключовi слова: точтсть термi-нального управлтня, помилка терщ-альних вимiрювань, бесплатформена
терщальна навкацшна система
□-□
Рассмотрено решение терминальной задачи управления и синтезирован параметризованный закон управления в аналитическом виде, обеспечивающий движение объекта в заданную точку и содержащий варьируемый параметр глубины прогноза. Исследованы особенности влияния величины параметра управления на точность достижения конечной точки, даны рекомендации по выбору параметра для нивелирования ошибки инерциальных измерений. Синтез управления осуществляется методом преследования ведущей точки по информации, полученной интегрированием измерений фактического ускорения и содержащей ошибку, характерную для акселерометров
Ключевые слова: точность терминального управления, ошибка инерци-альных измерений, бесплатформнная инерциальная навигационная система
■о о
UDC 629.78
|DOI: 10.15587/1729-4061.2018.123271]
DEVELOPMENT OF THE ALGORITHM FOR AIRCRAFT CONTROL AT INACCURATE MEASUREMENT OF THE STATE VECTOR AND VARIABLE ACCURACY PARAMETER
N. Khatsko
PhD, Senior Lecturer Department of Computer Modeling of Processes and Systems National Technical University "Kharkiv Polytechnic Institute" Kyrpychova str., 2, Kharkiv, Ukraine, 61002 E-mail: n.khatzko@gmail.com H. Makogon PhD, Senior Lecturer* E-mail: helmkg@ukr.net О. I sa kov Lecturer* E-mail: alexalex999@ukr.net S. Krepchenko Lecturer*
E-mail: s.krepchenko@gmail.com I. Cherepnov
PhD, Associate Professor** E-mail: voenpred314@ukr.net N. Kiriyenko PhD, Associate Professor** E-mail: hfinpomt@ukr.net V. Zadorozhnya PhD, Associate Professor** E-mail: zadorvika@ukr.net A. Polyansky Doctor of Technical Sciences, Professor** E-mail: khadi.pas@gmail.com *Department of Armored vehicles and military equipment Military institute of Tank Troops of National Technical University
"Kharkiv Polytechnic Institute" Poltavskyi Shliakh, str., 198, Kharkiv, Ukraine, 61098 **Department of Life Safety and Law Kharkiv Petro Vasylenko National Technical University of Agriculture Alchevskykh str., 44, Kharkiv, Ukraine, 61002
1. Introduction
In the era of significant expansion of air traffic, which is characterized by a wide spread of unmanned aerial
vehicles (UAV), automatic methods of flight control are becoming increasingly important. Automatic control of an aerial vehicle (AV) allows solving practical problems in various areas.
©
One of the elements of the systems of automatic control of AV motion is an on-board navigation system which provides obtaining information about the position, speed and orientation of a vehicle in space. Integrated inertial-satellite navigation systems (IISNS), which have become widespread due to advances of satellite radio navigation [1], have proven themselves as an effective way to achieve necessary precision of air navigation [2]. Operation of the systems is based on aggregation of information, received from inertial sensors (gyroscopes and accelerometers), and satellite information on current coordinates and speed of an AV. Thus, at regular and sufficiently frequent information flow from both navigation systems, it is possible, as long as you want, to determine navigational parameters with small errors, the value of which are virtually in the vicinity of zero. However, in practice, the order of arrival of information from each system has its own specific features.
Thus, a signal from the satellite navigation system can be of inadequate quality. It may be associated with magnetic storms, low inclination of satellites' orbits, re-reflection of a radio signal. Another challenge may be a search for a satellite signal by a receiver after disconnection, because initially, a receiver does not know where a satellite is, whether it approaches or retreats and what the shift of the frequency of its signal is. In addition, the contact with the found satellite continues only when a complete set of information from it is received. Otherwise, the search for a suitable satellite continues.
Even at fully deployed satellite orbital groupings, over 24 hours in some areas of the globe, there can be time intervals of more than 15 minutes when satellite information is not suitable for high-precision navigation due to special reciprocal location of satellites in space (geometry factor). Under these conditions, control of AV is carried out only by information of the inertial navigation subsystem with its inherent errors of determining [3]. In most cases, the onboard strapdown inertial navigation system (SINS) is used. The enumerated circumstances determine the relevance of the considered problem in the context of modern flight control technology.
2. Literature review and problem statement
During implementation of flight tasks, there often arise the tasks of AV flight to a specified point in space with assigned parameters of motion within the allotted time. Such tasks are called terminal problems. Of course, the best solution to a terminal task in the algorithmic sense is to obtain explicit analytical dependences for control functions [4]. In this case, it is possible at any moment to synthesize the values of control actions for AV reaching the destination.
Modern control technology is based on the theory of feedback. Control of current parameters of the state of an AV and formation of motion correcting control increases efficiency of an automatic flight control system (AFCS). Typically, the control process includes the following stages [5]:
- development of a reference trajectory of AV motion to a specified destination, that is, obtaining the required state of AV after a specified period of time under ideal flight conditions;
- obtaining information about the current state of AV from navigation systems on the pitch of AFCS operation;
- analysis of obtained information and development of control actions for continuation of the AV motion.
Thus, so that AFCS can operate correctly, the data exchange between a navigation system and an automatic control system is required. In modern AFCS, the main tool for obtaining information about the current vector of the AV state is often IISNS, uniting gauges of various physical nature - satellite and inertial navigation systems. However, global satellite navigation systems are subjected to various threats (ionospheric disturbances, multi-path scattering, unfavorable location of satellites and so on), any of which may lead to partial or complete loss of a signal [6].
To improve the accuracy of determining navigation and angular orientation parameters, IISNS use different methods of information aggregation: separate, weakly connected, rigidly connected and deeply integrated. Possibilities of obtaining higher precision of estimation of navigation parameters are studied in [7].
However, regardless of the data integration degree, at a delay of obtaining of reliable information from the on-board satellite navigation equipment, SINS operate in the autonomous mode. At autonomous operation of SINS, the known shortcomings fully affects accuracy of measurements of navigation parameters and control, generated on their basis. That is why autonomous operation of SINS is often accompanied by data correction from the other, non-inertial measurement sources [8, 9].
The accuracy of obtaining navigation parameters of SINS is directly influenced by accuracy of its sensors of primary information. While aviation used to apply only expensive devices, lately, there has been a strong tendency to make inertial systems cheaper. Technology of micro-electro-mechanical systems (MEMS) has been developing rapidly for 20 last years.
Accuracy of MEMS-based sensors used to be low, and they could be used only in cars and household appliances. However, lately, the accuracy of MEMS-based accelerometers increased and reached the navigation accuracy, MEMS-gyroscopes also have a high precision level [10]. The world's leading companies organized production of navigation modules with MEMS sensors for small and medium aircraft aviation.
As already described, during its operation, SINS accumulates measurement errors of inertial sensors [11]. The nature of accumulation of SINS errors significantly depends on maneuvers, performed by AV in the process of SINS operation, so relationship of errors of SINS with errors of sensors is extremely complicated. Under these conditions, numerical modeling is most often applied for the study of peculiarities of AV control.
At the same time, to clarify some peculiarities of interaction of a control system and an inertial navigation system, it seems necessary to perform analytical studies of the impact of errors of SINS in general and of inertial sensors in particular on accuracy of control of three-dimensional AV motion.
3. The purposes and objectives of the research
The aim of present research is to obtain analytic expressions for the synthesis of control actions, taking into account the errors of inertial navigation system and to estimate accuracy of AFCS, receiving input data of the inertial module.
To accomplish the set goal, the following tasks were to be solved:
- to obtain the law of control in the analytical form, including the generalized error of PFINS, for synthesis for control actions of the closed AFCS;
- to study the influence of the generalized error of SINS on accuracy of solution of the terminal problem and to consider the ways of accuracy improvement by changing the values of the prediction depth parameter in the law of control;
- to verify analytical conclusions by the method of computer simulation.
4. Obtaining the law of control for a closed AFCS
General formulation of the problem is as follows. Let us assume that at the initial moment, the motion, position and speed of AV are known, they are contained in the state vector. It is required to transfer an AV to a desired point of space with the help of control actions within the specified time at a required ultimate speed. It is also known, that information about the current value of the state vector of an AV is available with an error, correspondent to the accepted model of the inertial system (in a general case) or of the inertial sensor (in a simplified one-dimensional case).
A fundamental course of solving a problem will be considered in a one-dimensional case, so as not to clutter presentation of analytical results. Let us describe an undisturbed controlled motion of AV by equation
r(t) = u(t),
(1)
where r(t) is the coordinate of an object, u(t) is the sought-for controlling acceleration. At the initial (t=0) and final (t=T) moment, the state of a control object is characterized by a coordinate and speed of the motion:
r (0) = ^ r (0) = ^ r (T )= rT, r (T ) = rT.
(2) (3)
bo = bi = r0,
3(rT - r0) rT + 2r0
b2 =
b =
T2
T
2(r0 - rT ) , rT + r0
T3
T2
Knowing the equation of software trajectory (4), we easily find program control from (1) in the form of
u(t ) = 6b3t + 2b2.
After substituting coefficients (5) into this expression, we will obtain
u(t) =
+6 rT—T- - 2 ^
T2
12 + 6 ÎÂ T 3 T 2
rT + 2r„
t +
T
(6)
An object must reach the final state within time T.
Taking into account obtained information about the state vector, it is necessary to synthesize the law of control u(t, r, r), providing satisfaction of conditions (3) due to (1) and (2).
The basic methodology of solving this problem is the methodology for solving inverse problems of dynamics, where it is necessary to determine the program control u*(t) by the known trajectory of motion [4]. We will briefly describe obtaining of program control in the analytical form for control actions at non-disturbed motion.
Let us assign the program trajectory as the extreme of functional
T T
J r2(t) ■ dt = J u2(t) ■ dt.
0 0
In this case, the trajectory of motion takes the form
r (t) = b3 ■ t3 + b2 ■ t2 + b1 ■ t + b0, (4)
and coefficients bi, i = 0,3 are determined from conditions of agreement of solution (4) with boundary conditions (2), (3), and takes the following form:
Motion of AV, caused by influence (6), taking into account the assumed ideal character of motion, will be subsequently called reference. Control (6) will be called program control and the trajectory (4) will be called a program or reference trajectory.
Program control (6) corresponds to control in the open control scheme and is not often used in practice, because usually control is adjusted taking feedback into consideration. To implement the motion along the program trajectory under conditions of possible, parametric and dynamic disturbances, it is proposed to use the algorithm of synthesis by the method of pursuing the leading point.
Let us carry out synthesis of control by the current value of the state vector of AV. Let us assume that measurements of coordinate r(t) and speed r(t) are available at any moment. We will solve the problem of transferring the state vector of object (r,r) from point (r,r) to point (r *(t + t), r(t + t)) within some time t. New variable t will be called the prediction depth.
It should be noted that according to the problem statement, the destination point is selected on the reference trajectory (4). It is ahead of the current point by time t and by point of time T, in accordance with the procedure of plotting of trajectory (4), it comes to the desired final position, providing a solution to the original problem of terminal control.
We will introduce additional disturbance in motion (1) and limit it by module | f |< f *. Then the mathematical model of the problem of synthesis takes the form
r(t') = u(t') + f (t'), under boundary conditions:
r (0) = r(t), r (0) = r (t),
r (t) = r * (t + t), r (t) = r*(t + t), (7)
where t' e [0, t] is the variable of local time.
We will join the initial and the final points (7) by a cubic polynomial of the variable of local time, and obtain
r (t') = a3t'3 + a2t'2 + a11' + a0,
where
a0 = r (t ), a1 = r(t ),
a = 3
r * (t + t) -r(t ) r* (t + t) + 2r(t )
a = 2 r(t) -r* (t + t) + r* (t + t) + Tr(t)
Then we will obtain dependence for control by the local trajectory
u(t ') = 6a3t ' + 2a2.
(8)
We will limit the time interval of using control (8) only by one moment t' = 0, as it will be necessary to develop another local trajectory at the next moment.
As the law of control of synthesis of control, we will finally obtain
«(t, f, f) = 6 •r *(t + t)-f(t) - 2 •r*(t + t) + 2-f(t), (9)
in which both factual current values of assessment of the state vector and its predicted reference value, which is the known function of time, are used.
Fig. 1 shows that motion of the point deviated and point r (t) actually does not lie on the reference trajectory. Deviation in position occurred due to disturbing influence. It is necessary to develop control, which will return the object to the reference trajectory. Local return trajectory is shown in the Figure by blue dots.
Fig. 1. Deviation of disturbed motion from the reference trajectory, where r (t) is the measured state vector at time t, r (t) is the reference value of state vector at time t, r (t + t) is the reference value of state vector after time t (prediction depth)
We will analyze the accuracy of the obtained control of AV in a closed system. To analyze effectiveness of the selected scheme of synthesis, we will explore the closed control system in terms of possibility of implementation of the reference trajectory.
We will assume that the actual value of variable r (t), corresponding to solution of the disturbed equation of motion is different from reference value r*(t), assigned by (4), by magnitude
dr (t) = r (t) - r '(t).
Suppose that measurements of the state vector of SINS are ideal, in particular r(t) = r(t), r(t) = r(t). Then control (9) can be represented as
u (t, r, r) = u (t) + du (dr, dr), where control error:
du(t) = 48 r(t)-48 r(t). t t
Let us write equation for error as dr (t ) = du + f (t), and, having substituted in it (10), we will obtain
dr (t ) + 4 dr (t) +—dr (t ) = f (t) t t
Differential equation (10) has solution
(10)
dr (t ) = e
r V2 . Vî
C1 cos—t + C2 sin—t 1 t 2 t
+ drf (t ),
(11)
(12)
in which coefficients C\, C2 are determined from original conditions dr (0), 8r(0), but the forced component
drf (t) = et sin V2 0
A t
f (t-e)-de,
(13)
corresponds to the convolution of pulse transition function of the inertial section of the second order, which is (11), and function f(t).
To obtain a rough estimation of the magnitude of component drf (t) at t ^^ without specifying functional dependence f(t), we will use ratios | f (t) |< f' and
. Vï
sin--
t
<1.
Thus, the inequality holds for any t > 0,
|drf (t 1 e Z*'f *-d e =
2VÎ"
242
f '.
(14)
Expression (14) shows that the magnitude of error (12) is limited by module. The problem on roughness degree of the obtained estimation for the component of error | drf (t) |, caused by disturbance, is interesting Consider a particular case, at f (t) = f' = const. Then the forced component takes the form
drf (t) = t- f' 6
and differs from the estimation from (14) less than by two times. From this we can conclude that the used simplifications in obtaining estimation |drf (t)| do not lead to a significant distortion of the result.
2
It may be concluded from (12) that at ideal measurements of the state vector at f ° 0 and non-zero initial conditions, control (9) provides asymptotic stability of motion relative to the reference trajectory. However, at non-zero, but limited disturbance | f (t)|< f \ the point surely does
not deviate from the reference trajectory more than by t2
magnitude .
5. Influence of measurement error on the accuracy of solution to a terminal problem
Now that we have equation (11) for the error of the state vector, we will analyze the influence of the error of measurement of the state vector of SINS on accuracy of operation of the closed control system. In particular, we are interested in terminal accuracy of solution of the control problem, i.e. there is estimation of the error of reaching the reference state vector at the final moment.
Let the state vector (r, r) be estimated by SINS with an error. Because dependences of errors of SINS of all components of navigation parameters are quite complicated, it does not seem possible to introduce them into consideration. We will use some generalized error of determining of navigation parameters.
We will consider a one-dimensional case. We will assume that initial measurement is value a(t) = r(t) + da, corresponding to the sum of true acceleration and constant of error da. Estimation of the state vector (r, r) in this case is performed in accordance with the principle of inertial notation by double integration of measurements.
The considered case is characteristic of SINS operation, in which error da consists of the error of accelerometer itself and the error of compensation of the vector of free fall acceleration, caused by the error of gyroscopes. In some modes of AV motion, such total error can be considered constant.
In view of the above, and considering that estimation of values r, r comes from solutions to equation r(t) = a(t), we obtain
r(t) = r(t) + da t,
r(t) = r (t) + da •-.
Subsequently, after substitution of expressions (10), (11) in the law of control (9), for the closed system, equation of implementation of reference motion takes the form
dr(t) + r(t) + £d r(t) =
4da t 3da
■t2 +f (t).
(15)
Provided that dr (0) = 0, dr(0) = 0, at the initial moment this equation (15) has the solution
? /a dam2 da 2
dr (t) =--1 +—t2
2 6
1 V2+ 1 . V2 , _ f..
—cos—t + —;= sin—t + dr (t).
2 t V2 t 1 w
Thus, we obtained the expression for determining a control error in the closed circuit, which includes a constant measurement error and a random component.
Analysis of expression (16) makes it possible to argue that unlike (12), where dr(t) is a magnitude, limited by module, in (16) absolute value of error dr(t) infinitely increases at an increase in time t, which is explained by accumulation of the error of identification of the state vector. It can be seen from (16) that the higher t, the less the influence da on dr ceteris paribus. It is possible to show that lim(drT (t)) = 0 holds for any final T
A similar conclusion is also true for the error of terminal speed.
drd
da
(t) =V2 ™
■2t . S t
■sin--T-da T.
t
The limit value of the error of terminal speed at t ^ ^ is equal to zero. Thus, at negligibly small dynamic disturbances f \ it is advisable to choose large parameter t, included in the law of control (9). This choice makes it possible to decrease the influence of the error of acceleration measurement on the accuracy of terminal control.
6. Discussion of results: experimental research into the influence of magnitude of parameter t on improvement of control accuracy
To verify the assumptions about a possibility for improvement of accuracy of AV control by SINS data, a computation experiment with the use of the developed software package, simulating operation of AFCS [12, 13], was conducted on some time interval. Here, we will present the results of modeling of steady rectilinear motion, at which only one coordinate changes.
Simulation of motion was conducted on a small time interval equal to 100 seconds, assuming that within this period, connection with a satellite is sure to be resumed. In the course of the experiment, accelerometer measurement errors da remained at the assigned level, for example, during plotting the curves, value da = 0.02 m/s2 in Fig. 2. The magnitude of prediction depth t changes, to demonstrate this, values of t were selected: 1 sec., 200 sec. and 500 sec.
100
\br(t,l)\ 80
§
i-T \br(t,200)\ 60
o t
o> 40
.1 •3 \br(t,500)\
o o 20
CJ
0
-10
20
40
60
t, sec
100
(16)
Fig. 2. Coordinate error when using the law of control with different values of t
The curves in Fig. 2 show the change in control error dr(t) during the modeled motion of an object. The figure shows a monotonous rise of the curves, i. e., error dr(t) in-
2
2 —t
2
creases, but with different dynamics. The rate of an increase depends on the magnitude of prediction depth t - an increase in t slows down the increase rate 5r(t). Indeed, at the final point of section of motion at t=1 sec at the 100th second of the flight, 5r(t)=100 m, at t=200 sec at the 100th second of the flight, the error decreased to 47.5 m, while at t= =500 sec, it decreased to 23 m.
The form of the curves in Fig. 2 proves the theoretical findings that at an increase in the value of interval t, the rate of an increase in the coordinate error of AV to the final point decreases.
In the second experiment, the possibility of weakening of the influence of errors of inertial measurement on control influence due to a change in the parameter of synthesis of t was explored. But, unlike the previous experiment, in which the value of t was fixed in each particular session of flight simulation, this parameter changed.
Fig. 3 shows the results of three series of experiments under the following conditions: t ranges from 0.1 to 10,000 sec; the time of AFCS control only by data of the inertial unit is fixed: T={100, 200, 300} sec; there is no dynamic disturbance (f(t)° 0); the error of the inertial unit measurement is non-zero and corresponds to error of MEMS-based accelerometers.
-50
5(100,t) 5(200,t) "10°
b(300,T)
-150
_ . - - -
y / /
/ /
experiments take time, might be difficult to carry out, and are expensive.
However, each sensor or, in the case of MEMS-based sensors, a particular model, have a technical passport with accuracy characteristics, given by manufacturer. It would be possible to use the data from the passport for the algorithmic compensation of measurement errors, but actual values of measurement errors do not comply with the passport data for a variety of reasons. In this case, algorithmic compensation will not improve accuracy of control. In comparison with the above, the proposed method makes it possible to improve control accuracy without conducting complicated experiments on calibration, but rather by appropriate choice of the parameter. It is enough to consider the generalized error of inertial measurements, which can be estimated based on the data of the technical passport. Roughness of this assessment will be compensated by selection of parameter t in the law of control, because the influence of the parameter was proved analytically and experimentally. Subsequently, it is possible to make a table of correspondence of accuracy of the inertial unit to the value of the parameter of synthesis of control actions.
Another method of application of the established dependence is to change parameter t during the control session only by SINS data. This will be done in the subsequent modeling experiments.
Analysis of the results of experimental studies makes it possible to argue that the data, coming from the MEMS sensors of medium accuracy, can be used for control and control accuracy will be high enough.
7. Conclusions
2X10"' 4X10-' 6X10"' prediction depth t, sec
8 X 10 IX 10"4
Fig. 3. Dependence of coordinate error at the terminal point on prediction
depth
Fig. 3 shows that at an increase in parameter t, there is a decrease in terminal error. Such an effect of an increase in terminal accuracy and, consequently, an increase in AFCS accuracy are caused by slowing down of an increase in terminal error with an increase in parameter t. We will note that in the motion interval of 100 seconds (the upper curve), which is permissible at a temporary loss of connection with satellites, the error decreases very fast and tends to zero. The results of the second experiment indicate a monotonous decrease in terminal error at an increase in t.
A method of synthesis of control functions, which unlike traditional methods contains the parameter, which makes it possible to regulate AFCS accuracy ceteris paribus, was proposed. This enables us to argue that data, coming from the MEMS sensors with high accuracy, can be used for control and control accuracy will be quite high.
To improve the accuracy of operation of inertial units, algorithmic compensation the primary measurements is widely used. For this, one conducts experiments on calibration of particular sensors and the inertial unit in the assembly, during which components of measurement errors are assessed. If each sensor of the unit is assessed and characterized in this way, results are sure to improve, but such
1. We present the procedure of deriving an analytical expression for calculation of values of the control function, which takes into account the value of the generalized error of measurements of the inertial module and the level of external dynamic disturbance. This approach made it possible to extend and modify the traditional method for solving the inverse problem of dynamics in the case of the terminal control problem. Thus, we received a modified law of control that can be applied in AFCS, which has a new component, such as the parameter, leveling the factors of decrease in accuracy control. The value of the so-called prediction was selected as the parameter.
2. Some aspects of the problem of using output information of the inertial navigation system as input data for synthesis of control influence on AV were explored. Analytical expressions for estimation of accuracy of object control at the point of destination of the motion were obtained. It was shown that, despite distorted navigation parameters, it is possible to maintain the accuracy of solution of the terminal problems within some time due to increasing errors of the inertial module. This possibility was obtained due to param-etrization of the law of control.
We provided a set of recommendations on selection of a changing parameter of control synthesis depending on the precision level of the measuring sensor, used in the inertial module. While selecting a variable value of prediction depth in the law of synthesis by the method of pursuing of leading point, it is recommended to consider the following:
- if the information-measuring inertial system provides obtaining of coordinates and speed of an object with a high degree of accuracy, i. e. inertial sensors have a small measurement errors, a selected prediction parameter should be minimal;
- when using low-precision sensors, the prediction parameter should be increased.
For determining a value of the parameter and final solution to the problems of adjustment of the algorithm of synthesis while working with certain inertial module, it is
necessary to conduct AFCS modeling using data on sensors' accuracy.
3. By using computer simulation, and analytically previously, it was shown that it is possible to decrease the influence of SINS measurement errors on the synthesis of values of control functions by selecting the value of prediction depth.
Subsequently, it is planned to use the presented algorithm for conducting of a full-scale experiment using a small unmanned aerial vehicle.
References
1. Revnivyh S. G. Tendencii razvitiya global'noy sputnikovoy navigacii // Mezhdunarodnaya konferenciya po integrirovannym nav-igacionnym sistemam: XIX mezhdunarodnaya konferenciya. Sankt-Peterburg: GNC RF CNII "Elektropribor", 2012. P. 266-276.
2. INS/CNS/GNSS Integrated Navigation Technology / Quan W., Li J., Gong X., Fang J. Springer, 2015. 372 p. doi: 10.1007/978-3662-45159-5
3. Chiang K.-W., Duong T., Liao J.-K. The Performance Analysis of a Real-Time Integrated INS/GPS Vehicle Navigation System with Abnormal GPS Measurement Elimination // Sensors. 2013. Vol. 13, Issue 12. P. 10599-10622. doi: 10.3390/s130810599
4. Krut'ko P. D. Obratnye zadachi dinamiki upravlyaemyh sistem. Lineynye modeli. Moscow: Nauka, 1987. 304 p.
5. Nebylov A., Tikhomirov M., Benzerrouk H. Development Stages and Prospects of Flight Control and Navigation Complexes for Aircraft // IFAC-PapersOnLine. 2017. Vol. 50, Issue 1. P. 6434-6439. doi: 10.1016/j.ifacol.2017.08.1032
6. Sabatini R., Moore T., Ramasamy S. Global navigation satellite systems performance analysis and augmentation strategies in aviation // Progress in Aerospace Sciences. 2017. Vol. 95. P. 45-98. doi: 10.1016/j.paerosci.2017.10.002
7. Miao L., Shi J. Model-based robust estimation and fault detection for MEMS-INS/GPS integrated navigation systems // Chinese Journal of Aeronautics. 2014. Vol. 27, Issue 4. P. 947-954. doi: 10.1016/j.cja.2014.06.012
8. Integrirovannaya s magnitnym kompasom BINS s vrashchayushchimsya izmeritel'nym modulem na MEMS / Ci U., Czichzhou L., Ya C., Czyan'e L. // Giroskopiya i navigaciya. 2016. Issue 3 (94). P. 3-13.
9. Vas'kov A. G., Klychnikov V. V. Sposoby kompleksirovaniya navigacionnyh sistem letatel'nyh apparatov // Nauchno-prakticheskie issledovaniya. 2017. Issue 2. P. 49-52.
10. Barbour N. M. Inertial Navigation Sensors // NATO RTO Lecture Series, RTO-EN-SET-116. 2010. URL: http://www.dtic.mil/ get-tr-doc/pdf?AD=ADA581016
11. Golovan A. A., Parusnikov N. A. Matematicheskie osnovy navigacionnyh sistem. Chast' I. Matematicheskie modeli inercial'noy navigacii. Moscow: Izd-vo MGU, 2011. 132 p.
12. Khatsko N. Ye. Prohramne zabezpechennia dlia modeliuvannia roboty systemy avtomatychnoho keruvannia polotom po vymiram inertsialnykh datchykiv // Materialy mizhnarodnoi naukovo-praktychnoi konferentsiyi «Suchasni metodyky, innovatsiyi ta dosvid praktychnoho zastosuvannia u sferi tekhnichnykh nauk». Liublin, 2017. P. 30-33.
13. Hac'ko N. E. Issledovanie vozmozhnosti upravleniya poletom po dannym inercial'nyh datchikov nizkogo klassa tochnosti // Prob-lemy mashynobuduvannia. 2013. Vol. 16, Issue 5. P. 32-35.