Preview

Journal of «Almaz – Antey» Air and Space Defence Corporation

Advanced search

Development of a methodological approach and software for information processing and experimental assessment of the accuracy of the “Sazhen-TA” measuring system

https://doi.org/10.38013/2542-0542-2020-1-96-106

Abstract

This paper presents algorithms for determining trajectory parameters using the information provided by the “Sazhen-TA” system. Issues associated with experimental assessment of the accuracy of the “Sazhen-TA” measuring system on the basis of reference data from satellite navigation equipment are considered. The results of mathematical modelling are presented.

For citation:


Simonov S.M., Kisin Yu.K., Nekrasov A.V., Shamshin M.V. Development of a methodological approach and software for information processing and experimental assessment of the accuracy of the “Sazhen-TA” measuring system. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2020;(1):96-106. https://doi.org/10.38013/2542-0542-2020-1-96-106

Introduction

With introduction of the Sazhen-TA trajectory measuring system into the measuring complex, a vital task is seen in the development of methodological approach and mathematical software support (MSS) of processing and experimental assessment of the measurement accuracy of this system.

The Sazhen-TA system, consisting of two sets, is deployed in the launch site area for taking measurements at the initial segment of aircraft (A/C) flight. The azimuth and elevation angle are measured in the visible-spectrum and infrared ranges. The tasks of determining A/C trajectory coordinates by the measured angular parameters and assessing the accuracy of the obtained trajectory are considered. The software is developed on the algorithmic programming languages PASCAL, C#, which allows to adapt its application to the operating system MSVS OS accepted and implemented in the Armed Forces of the Russian Federation.

The content of the paper is determined by the topicality of the task consisting in estimation of the Sazhen-TA system accuracy during measuring complex operation in the sea range conditions. Taking into account the nature of flight test objects, reference facilities are selected, based on which an experimental assessment of the Sazhen-TA system accuracy is made. At the first stage, this task of experimental assessment of the accuracy required development of a methodological approach and software for processing the Sazhen-TA system information on determining experimental trajectory, and on the second stage - creation of respective algorithms and programs for determining experimental accuracy of the Sazhen-TA measuring system.

Algorithms for determining aircraft trajectory by the end method and least-square method based on the Sazhen-TA system information

The task of determining trajectory coordinates is implemented as per two measuring facilities by the least-square method (LSM) [1] and the end method, using intersection of two beams determined by the azimuths and elevation angles measured by two Sazhen-ТА systems in accordance with [1–3].

Consideration is given to the task (α1, γ1, α2, γ2) → (x, y, z) of determining aircraft (A/C) coordinates by azimuths and elevation angles measured by two Sazhen-ТА systems. As pointed out in [2, 3], this task is equivalent to the task of determining A/C coordinates for two phase direction finders which measure direction cosines, as considered in the known monograph on the methods of statistical processing of trajectory measurements [1].

The association of direction cosines (cos θx, cos θz) with azimuth and elevation angle (α, γ) [2, 3] is determined by the formula:

cos θχ = sin γ · cos α,                                             (1)

cos θz = sin γ · sin α,

In [1], there are formulas for determining object coordinates for two phase direction finders which measure direction cosines. A/C coordinates in the local coordinate system (LCS) of the 1st measuring facility (MF) can be determined by the formula:

Thus, it is necessary to determine parameter D1 – distance between the A/C and standing point of the 1st MF.

The method of solving the task is illustrated in Fig. 1 [1], showing angles δ, φ, ψ used in subsequent calculations.

Fig. 1. Determining A/C coordinates as per two measuring facilities

The following designations are used in Fig. 1: O1X1Y1Z1 and O2X2Y2Z2 of the LCS of two measuring facilities numbered 1 and 2. Points O1 and O2 – location points of the measuring facilities. Point S – A/C current location point. Parameter D1 is determined as follows:

where b – distance between the two MF, b10 – unit vector (shown in Fig. 1 and calculated in accordance with [1]). Triangle angles O1O2S are determined from relations:

where Φ1, Φ2 – matrices of direction cosines for transition from the measuring facilities’ LCS to the Greenwich coordinate system (GCS) [1];
D01 = (cos θχ1, cos θy1, cos θz1), D2 = (cos θχ2, cos θy2, cos θz2).

The GCS and LCS coordinate systems are defined in accordance with [1, p. 194, 198].

Geocentric reference (Greenwich) coordinate system (GCS) is a 3D rectangular coordinate system OXгYгZг, with its datum point lying in the centre of Earth’s ellipsoid, axis OXr of which is located in the equatorial plane and directed towards zero (Greenwich) meridian, axis OZr is directed along the Earth’s rotation axis towards the North pole, and axis OYr complements the coordinate system to a right-handed one.

Local coordinate system (LCS) is a 3D rectangular coordinate system OXмYмZм associated with the Earth’s surface, with its datum point having geodetic coordinates: Вм (geodetic latitude), LM (geodetic longitude), hM (elevation above the surface of Earth’s ellipsoid); axis OYм is directed along the outward normal towards the Earth’s ellipsoid surface; axis OXм lies in the plane tangential to the Earth’s ellipsoid in point O and makes with the datum point meridian plane angle Ам (geodetic azimuth), which is reckoned clockwise from a direction to the North; axis OZм complements the coordinate system to a right-handed one.

For LCS of the Sazhen-TA system, Aм = 0.

The formulas of direction cosines association with azimuth and elevation angle (α1, γ1, α2, γ2) → (cos  θx1, cos θz1, cos θx2, cos θz2) are assigned by relation (1).

In this way, the task (α1, γ1, α2, γ2) → (x, y, z) of determining object coordinates by azimuths and elevation angles measured by two measuring facilities is solved as well. The output results are the trajectory parameters in the GCS.

Determination of trajectory parameters by azimuths and elevation angles measured by two measuring facilities is also performed by the LSM in the GCS.

The LSM task is solved iteratively, with a system of normal equations being solved at each step [1, 4, 5].

For each instant of the measurement period, trajectory coordinates to be assessed are determined iteratively by the formula:

where Xˆ(k+1), Xˆ(k) – assessments of trajectory coordinates in the GCS at iterations k +1 and k;

ΔΧ(k+1) – corrections to trajectory parameters at iteration k + 1;

A(k) – matrix of partial derivatives from the measured parameters of two MF as per assessed trajectory parameters in the GCS from the k-th iteration;

H – measurement vector (measured parameters (α1, γ1, α2, γ2);

Y(Xˆ(k)) – vector of measured parameters calculated as per assessment of coordinates at the k-th iteration;

W – inverse matrix to the covariance matrix of measurement errors Ky;

Ky – diagonal square matrix of the 4th order with dispersion of rated errors of the Sazhen-TA system on the diagonal.

The initial approximation Xˆ(0) is assigned as per the data of a priori estimated trajectory.

The matrix of partial derivatives from the measured parameters of two MF as per assessed trajectory parameters in the GCS has a block view:

where A1, A2 – respective matrices of partial derivatives for the first and second MF.

At that, each one of the said matrices has the following view:

where B – matrix of partial derivatives from the measured parameters as per trajectory parameters in the LCS (with MF index number omitted);

Ф – matrix of direction cosines of MF’s LCS association with the Greenwich coordinate system [1];

x, y, z – A/C coordinates in the LCS.

The elements of B matrix are calculated by the formulas:

When calculating estimated values of the vector elements of measured parameters Y(Xˆ(k)) as per coordinates at the k-th iteration, azimuth and elevation angle for each MF are determined by the formulas, with the iteration index omitted:

where x, y, z – A/C coordinates in the LCS (for a given MF).

When calculating arctangent, respective procedure of the programming language is applied, determining angle within an interval of (0,2).

During calculations, one should take into account possible transition for the estimated azimuth value from the 1st to the 4th quarter, and vice versa.

Transition in the measured parameters from azimuth and elevation angle to parameters in the form of direction cosines by formula (1) makes it possible to operate with continuous functions and ease the requirements to initial approximation in the iteration process. The programs implement processing of both types of measured parameters.

The accuracy of A/C trajectory parameters is determined by the covariance matrix:

where k corresponds to the number of the last iteration in formula (6).

The considered methods of determining trajectory parameters by the angular measurements are applicable both for the active and the passive segments of A/C flight.

It is also possible to determine A/C trajectory parameters by the results of measurements taken by one Sazhen-TA system – as the coordinates of intersection point of a beam, determined by azimuth and elevation angle, with the A/C flight plane. This method is applicable at a segment of measurements where trajectory parameters lie ‘almost’ in the same plane. To find trajectory segments lying in ‘almost’ the same plane, a special procedure is applied. In the said procedure, plane is determined by three selected trajectory points, and distances from other trajectory points to the selected plane are calculated. If for all the points of selected trajectory segment the calculated distances do not exceed small number ε, then the trajectory points lie ‘almost’ in the same plane. Number ε is selected at the level of errors of trajectory parameters determination by the LSM as per two measuring facilities in accordance with matrix Kx.

At the modelling stage, the calculations for selection of trajectory segment, the points of which lie ‘almost’ in the same plane, are made as per estimated trajectory. The estimated trajectory is provided by the A/C designer. During modelling, the parameters of this trajectory are taken as the true values.

Let us suppose that there is a trajectory segment with points lying in ‘almost’ the same plane, and there are measurements for this segment taken by one Sazhen-TA system. The end formulas for solving the task of determining trajectory by one Sazhen-TA system will be obtained by finding the coordinates of intersection point of beam (l1), determined by angles (α1, γ1), and the plane.

The plane can be defined by three points in the GCS: M1 = (X1, y1, Z1), M2 = (X2, y2, Z2), M3 = (X3, y3, Z3).

In a general plane equation A · x + B · y + C · z + D = 0, plane coefficients are determined in accordance with equation (8) [6, 7]:

In the LCS of the first measuring facility, beam (l1) lies on a straight line determined by equation in the parametric form:

In the GCS, straight-line equation of beam (l1) has the view:

Intersection of beam (l1), determined by equation (10), with plane (П) corresponds to parameter t1:

A/C coordinates in the GCS are determined by the formula:

The method of solving the task of determining A/C trajectory coordinates at a flight segment, the points of which lie ‘almost’ in the same plane, as per the data from one (first or second) Sazhen-TA system set, is illustrated in Fig. 2.

Fig. 2. Determining A/C trajectory coordinates at ‘almost’ flat segment of flight as per data from one (first or second) Sazhen-TA system set

The following abbreviations and designations are used in Fig. 2: OMF-1 – optical measuring facility No. 1 with LCS coordinates O1X1Y1Z1; EFLSC – Earth-fixed launch site coordinates; OcXcYcZc – EFLSC datum point and axes; M1, M2, M3, M4 – points which angular measurements of azimuth and elevation angle were taken in: (α1, γ1, α2, γ2, α3, γ3, α4, γ4).

Figs. 3–5 show deviations of A/C coordinates determined by three methods from the reference. The numbers of measurements (as well as in all other figures) are given on the horizontal axis, and the vertical axis shows deviations of A/C coordinates in km.

The plots of Figs. 3–5 demonstrate experiential coincidence of the results of trajectory determination by the LSM and beams intersection as per two Sazhen-TA systems, as well as acceptable accuracy of trajectory determination by one Sazhen-TA system at ‘almost’ flat segment of flight.

Fig. 6 shows the distances in metres to the plane from the current trajectory points when determining A/C trajectory coordinates at an ‘almost’ flat segment of flight.

Fig. 7 shows the angle between beam and plane and the angle of beams intersection in degrees. If the angle between beam and plane is close to zero, in this case, an error of trajectory determination by one Sazhen-TA system will be great. In this case, the task of trajectory determination by the measurements taken by one Sazhen-TA system is poorly visible [1, 4, 5]. It is recommended in the operation manual to ensure that the intersection angle of beams (lines of sight) from two Sazhen-ТА systems exceeds 10 degrees.

Fig. 6. Distances from plane to current trajectory points

Development of a methodological approach and software for experimental assessment of the accuracy of the Sazhen-TA measuring system

The efficiency of developed methodological approach and software is assessed on the basis of mathematical modelling. Experimental assessment of the accuracy is done at the level of measured parameters (azimuth and elevation angle) of the Sazhen-TA system and experimental trajectory parameters determined as per information of the Sazhen-TA system in accordance with the algorithms of the first section.

During mathematical modelling, estimated trajectory is taken as the true one. Based on this trajectory, angular measurements by two SazhenТА systems are modelled, with random errors of ‘measurements’ input by a pseudorandom number sensor with normal law of distribution and RMSD corresponding to the rated values, i. e. 5 arc seconds for the optical channel and 10 arc seconds for the infrared channel (all modelling calculations were performed with this RMSD level).

The methodology for assessing measurement accuracy with the use of a reference measuring facility is implemented in accordance with [1, Chapter 6]. Modelling of a reference trajectory is done by entering random errors into the estimated trajectory by means of pseudorandom number sensor with normal law of distribution and RMSD corresponding to the covariance matrix data obtained from the results of satellite navigation equipment information processing. Based on the reference trajectory data, the measured angular parameters of the Sazhen-ТА systems are calculated.

Random process vector is calculated:

where  – measured and reference vector, respectively, ti – time instant of taking the measurements, i – measurement number, j – measurement session number.

Vector components in (13) correspond to angular measurements (α1, γ1, α2, γ2) or experimental trajectory coordinates (x, y, z).

Fig. 8 illustrates implementation of the random process of deviations of experimental trajectory parameters by the LSM as per information of two measuring facilities from the reference data, in kilometres.

Fig. 9 illustrates implementation of the random process of deviations of two systems’ measurements from the reference data, in arc seconds.

The plots of Fig. 9 show that ‘random errors’ of the modelled angular measurements in a given implementation for two measuring systems do not exceed three RMSD of the rated values for the infrared channel, i. e. 10 arc seconds, even taking into account the reference trajectory errors.

Investigation into random processes  is carried out, same as in [1, Chapter 6], under assumption that they are stationary and ergodic, and assessment of statistical characteristics is performed within one implementation, therefore, further on index j (measurement session number) is omitted.

Assessment of statistical characteristics is performed at the level of mathematical expectation and RMSD.

Random processes  are represented in the form of three components of slowly varying (systematic) errors (SVE) and rapidly varying (random) errors (RVE), as well as the errors of reference data:

It is further assumed that, by the norm, condition  is fulfilled for error vectors and the reference data error can be disregarded.

Assessment   is done by approximation of experimental data from the LSM using algebraic polynomials or orthogonal Chebyshev polynomials [1, 8].

Smoothing out is performed independently by each vector coordinate .

Modelling of the introduction of slowly varying errors   is performed as per angular measurements using the following formulas:

where v – SVE linear variation coefficient assigned in the modelling program.

The influence of systematic errors of angular measurements on determination of trajectory parameters is effected after processing as per algorithms of the fi rst section through comparison with the trajectory parameters without systematic errors in the angular measurements.

Singling out of slowly varying errors in implementation of a random process is done through approximation with the 3rd grade orthogonal Chebyshev polynomials with the use of WinLTX software package developed by REC ETU, St. Petersburg.

The modelling results are given in Fig. 10, in kilometres.

The plots of Fig. 10 show considerable influence of systematic measurement errors on the experimental trajectory and efficient singling out of the systematic errors through approximation of deviations between the measured and the reference data by means of orthogonal Chebyshev polynomials.

Fig. 11 illustrates in kilometres root-meansquare deviations (RMSD) of experimental trajectory parameters by the LSM as per variance of covariance matrix Kx in accordance with (7).

Plots with designations RMSD X 2MF, RMSD Y 2MF, RMSD Z 2MF correspond to the RMSD obtained during processing of two measuring facilities.

Comparison of the plots in Figs. 11 and 8 shows that ‘random errors’ of the experimental trajectory in this implementation do not exceed three RMSD of experimental trajectory parameters by the LSM, even taking into account the reference trajectory errors.

Table 1 contains assessment of mathematical expectation and RMSD, in metres, of rapidly varying (random) errors of determining experimental trajectory as per the reference data.

Table 1

Assessment of mathematical expectation and RMSD of RVE as per reference data

 

X

Y

Z

Math expectation

-0.79

-0.51

-1.18

RMSD

7

118

59

The data in Table 1 correlate with the results of LSM measurements processing, taking into account the ‘reference’ trajectory errors.

When processing real information, experimental accuracy of the Sazhen-TA system is assessed applying satellite navigation equipment as data reference.

Conclusions

  1. Algorithms and programs for determining experimental parameters of the aircraft trajectory by azimuth and elevation angle, as measured by the Sazhen-TA system, are developed.
  2. Methodological approach and software for experimental assessment of the Sazhen-TA system accuracy with application of the reference data of satellite navigation equipment are developed.
  3. Performance capability of the proposed algorithms is confirmed through mathematical modelling.

References

1. Жданюк Б. Ф. Основы статистической обработки траекторных измерений. М.: Советское радио, 1978. 384 с.

2. Кисин Ю. К. О применении алгоритмов на основе метода наименьших квадратов и конечных формул в задачах обработки траекторных измерений // Вестник Концерна ВКО «Алмаз-Антей». 2016. № 3. С. 74–79.

3. Землянов А. Б, Кисин Ю. К. Определение параметров траектории летательного аппарата по минимальному набору угловых измерений // Морская радиоэлектроника. 2017. № 2. С. 56–58.

4. Мудров В. И., Кушко В. Л. Методы обработки измерений. 2-е изд. М.: Радио и связь, 1983. 304 с.

5. Барабанов О. О., Барабанова Л. П. Математические задачи дальномерной навигации. М.: Физматлит, 2007. 272 с.

6. Выгодский М. Я. Справочник по высшей математике. М.: Астрель – АСТ, 2006. 991 с.

7. Александров А. Д., Нецветаев Н. Ю. Геометрия: Учебное пособие. М.: Наука, Гл. ред. физ.-мат. лит., 1990. 672 с.

8. Хрущева И. В., Щербаков В. И., Леванова Д. С. Основы математической статистики и теории случайных процессов: Учебное пособие. СПб.: Лань, 2009. 336 с.


About the Authors

S. M. Simonov
Military Unit No. 09703

Sergey Mikhailovich Simonov – Cand. Sci. (Engineering), Deputy Head for R&D of the First Scientific Centre

Research interests: methods for polygon testing of aircrafts.



Yu. K. Kisin
Military Unit No. 09703

Yuri Konstantinovich Kisin – Cand. Sci. (Engineering), Senior Researcher, R&D Department

Research interests: methods for polygon testing of aircrafts.



A. V. Nekrasov
Military Unit No. 09703

Alexander Vladimirovich Nekrasov – Head of R&D Department

Research interests: methods for polygon testing of aircrafts.



M. V. Shamshin
Military Unit No. 09703

Mikhail Vladimirovich Shamshin – Junior Researcher, R&D Department

Research interests: methods for polygon testing of aircrafts.



Review

For citation:


Simonov S.M., Kisin Yu.K., Nekrasov A.V., Shamshin M.V. Development of a methodological approach and software for information processing and experimental assessment of the accuracy of the “Sazhen-TA” measuring system. Journal of «Almaz – Antey» Air and Space Defence Corporation. 2020;(1):96-106. https://doi.org/10.38013/2542-0542-2020-1-96-106

Views: 347


Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 License.


ISSN 2542-0542 (Print)