Научная статья на тему 'ИССЛЕДОВАНИЕ ОБЛАСТИ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ БАЗИСНОЙ ФУНКЦИИ ДВУХ АРГУМЕНТОВ ПРИ ПОСТРОЕНИИ Λ-ОРТОГОНАЛЬНОЙ БАЗИСНОЙ ФУНКЦИИ'

ИССЛЕДОВАНИЕ ОБЛАСТИ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ БАЗИСНОЙ ФУНКЦИИ ДВУХ АРГУМЕНТОВ ПРИ ПОСТРОЕНИИ Λ-ОРТОГОНАЛЬНОЙ БАЗИСНОЙ ФУНКЦИИ Текст научной статьи по специальности «Математика»

CC BY
18
3
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ЛИНЕЙНО НЕЗАВИСИМАЯ БАЗИСНАЯ ФУНКЦИЯ ДВУХ АРГУМЕНТОВ / Λ-ОРТОГОНАЛЬНАЯ БАЗИСНАЯ ФУНКЦИЯ / ОБЛАСТЬ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ / ПОКАЗАТЕЛЬ ТОЧНОСТИ

Аннотация научной статьи по математике, автор научной работы — Щербов Игорь Леонидович

Цель исследований - определение области задания аргумента τ и интервалов его дискретизации при построении алгоритма адаптивного нелинейного оптимального сглаживания многопараметрических данных траекторных измерений, позволяющего совместно реализовать пространственную и временную избыточность получаемых данных. Исследования проводились путем построения Λ-ортогональной базисной функции с целью получения независимых оценок коэффициентов сглаживающего полинома. Показано, что решение задачи по определению максимально правдоподобной оценки вектора коэффициентов сглаживающего полинома целесообразно осуществлять методом последовательных приближений. При построении Λ-ортогональной базисной функции максимально правдоподобная оценка вектора коэффициентов сглаживающего полинома достигается в 2-3 итерации. Из результатов исследований, приведенных в работе, следует, что показатель точности Qт как функция двух аргументов (τ0 - наименьшее значение аргумента и Δτ - интервал дискретизации аргумента τ) в широком диапазоне значений этих аргументов изменяется слабо, но резко увеличивается при Δτ ≪τ0 . При этом величины значений аргументов не должны превосходить соответственно максимально и минимально возможных чисел, которые можно без потери точности записать в разрядную сетку применяемой ЭВМ. При равномерном шаге дискретизации аргумента τ, аргумент τ0 целесообразно выбирать в средней части интервала, где qmin и qmax соответственно минимальные и максимальные числа, которые без потери точности можно записать в разрядную сетку ЭВМ. Приближение к краям интервала может при неблагоприятных условиях привести к возрастанию ошибок вычислений при определении вторичных параметров положения летательного аппарата из-за того, что основная матрица системы уравнений становится плохо обусловленной.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

STUDY OF THE AREA OF DETERMINATION OF PARAMETERS OF THE BASIS FUNCTION OF TWO ARGUMENTS IN CONSTRUCTING THE Λ-ORTHOGONAL BASIS FUNCTION

The purpose of the research is to determine the area of setting the argument τ and its discretization intervals when constructing an algorithm for adaptive nonlinear optimal smoothing of multi-parameter data of trajectory measurements, which makes it possible to jointly implement the spatial and temporal redundancy of the data obtained. The research was carried out by constructing a Λ-orthogonal basis function in order to obtain independent estimates for the coefficients of the smoothing polynomial. It is shown that it is advisable to solve the problem of determining the maximum likelihood estimate of the coefficient vector of the smoothing polynomial by the method of successive approximations. When constructing a Λ-orthogonal basis function, the maximum likelihood estimate of the coefficient vector of the smoothing polynomial is achieved in 2-3 iterations. It follows from the research results presented in the paper that the accuracy index Qт as a function of two arguments (τ0 is the smallest value of the argument and Δτ is the discretization interval of the argument τ) in a wide range of values of these arguments changes slightly, but increases sharply at Δτ≪τ0 . In this case, the values of these arguments should not exceed, respectively, the maximum and minimum possible numbers that can be written without loss of accuracy in the re grid of the computer used. With a uniform discretization step of the argument τ it is advisable to choose the argument τ0 in the middle part of the interval¸where qmin and qmax respectively, are the minimum and maximum numbers that can be written into the bit grid of a computer without loss of accuracy. In case of adverse conditions Approaching the edges of the interval can lead to an increase in calculation errors in determining the secondary parameters of the position of the aircraft due to the fact that the main matrix of the system of equations becomes ill-conditioned.

Текст научной работы на тему «ИССЛЕДОВАНИЕ ОБЛАСТИ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ БАЗИСНОЙ ФУНКЦИИ ДВУХ АРГУМЕНТОВ ПРИ ПОСТРОЕНИИ Λ-ОРТОГОНАЛЬНОЙ БАЗИСНОЙ ФУНКЦИИ»

УДК 004.942:519.254:519.876.5:621.396.96:629.7.018 DOI 10.18522/2311-3103-2022-6-106-116

И.Л. Щербов

ИССЛЕДОВАНИЕ ОБЛАСТИ ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ БАЗИСНОЙ ФУНКЦИИ ДВУХ АРГУМЕНТОВ ПРИ ПОСТРОЕНИИ Л-ОРТОГОНАЛЬНОЙ БАЗИСНОЙ ФУНКЦИИ

Цель исследований - определение области задания аргумента т и интервалов его дискретизации Дгу Дzz при построении алгоритма адаптивного нелинейного оптимального сглаживания многопараметрических данных траекторных измерений, позволяющего совместно реализовать пространственную и временную избыточность получаемых данных. Исследования проводились путем построения Л-ортогональной базисной функции с целью получения независимых оценок коэффициентов сглаживающего полинома. Показано, что решение задачи по определению максимально правдоподобной оценки вектора коэффициентов сглаживающего полинома целесообразно осуществлять методом последовательных приближений. При построении Л-ортогональной базисной функции максимально правдоподобная оценка вектора коэффициентов сглаживающего полинома достигается в 2-3 итерации. Из результатов исследований, приведенных в работе, следует, что показатель точности Qnt как функция двух аргументов (т0 - наименьшее значение аргумента и Ат - интервал дискретизации аргумента т) в широком диапазоне значений этих аргументов изменяется слабо, но резко увеличивается при Лт << т0. При этом величины значений аргументов не должны превосходить соответственно максимально и минимально возможных чисел, которые можно без потери точности записать в разрядную сетку применяемой ЭВМ. При равномерном шаге дискретизации аргумента т, аргумент т0 целесообразно выбирать в средней части интервала, где qmin и qmax соответственно минимальные и максимальные числа, которые без потери точности можно записать в разрядную сетку ЭВМ. Приближение к краям интервала может при неблагоприятных условиях привести к возрастанию ошибок вычислений при определении вторичных параметров положения летательного аппарата из-за того, что основная матрица системы уравнений становится плохо обусловленной.

Линейно независимая базисная функция двух аргументов; Л-ортогональная базисная функция; область определения параметров; показатель точности.

I.L. Shcherbov

STUDY OF THE AREA OF DETERMINATION OF PARAMETERS OF THE BASIS FUNCTION OF TWO ARGUMENTS IN CONSTRUCTING THE Л-ORTHOGONAL BASIS FUNCTION

The purpose of the research is to determine the area of setting the argument т and its discretization intervals Дтх, Д^у, Д^z when constructing an algorithm for adaptive nonlinear

optimal smoothing of multi-parameter data of trajectory measurements, which makes it possible to jointly implement the spatial and temporal redundancy of the data obtained. The research was carried out by constructing a Л-orthogonal basis function in order to obtain independent estimates for the coefficients of the smoothing polynomial. It is shown that it is advisable to solve the problem of determining the maximum likelihood estimate of the coefficient vector of the smoothing polynomial by the method of successive approximations. When constructing a Л-orthogonal basis function, the maximum likelihood estimate of the coefficient vector of the smoothing polynomial is achieved in 2-3 iterations. It follows from the research results presented in the paper that the accuracy index Qnt as a function of two arguments (т0 is the smallest value of the argument and Ат is the discretization interval of the argument т) in a wide range of values of these arguments changes slightly, but increases sharply at Лт << т0. In this case, the values of these arguments should not exceed, respectively, the maximum and minimum possible numbers that can be written without loss of accuracy in the re grid of the computer used. With a uniform discretization step of the argument

t it is advisable to choose the argument T0 in the middle part of the interval3where qmln and qmax respectively, are the minimum and maximum numbers that can be written into the bit grid of a computer without loss of accuracy. In case of adverse conditions Approaching the edges of the interval can lead to an increase in calculation errors in determining the secondary parameters of the position of the aircraft due to the fact that the main matrix of the system of equations becomes ill-conditioned.

Linearly independent basis function of two arguments; A-orthogonal basis function; area of determination of parameters; accuracy index.

Введение. Измерения, проводимые в процессе лётных испытаний, предназначены для определения и анализа реального положения летательных аппаратов. По результатам измерений и последующей их обработки должно быть принято решение о степени соответствия реальных характеристик летательных аппаратов, заданным на них требованиям.

Как правило, оценка лётных характеристик летательных аппаратов осуществляется не только по результатам летных испытаний, но и по результатам математического моделирования различных этапов и траекторий полета, когда речь идет о возможностях летательных аппаратов, проверка которых в процессе летных испытаний может быть связана с риском потери летательного аппарата.

Основная часть. В статьях [1, 2] приводятся результаты исследования линейно-независимой базисной функции (ЛНБФ) двух аргументов:

(t,r) =

(00(tTx )(01(tTx )(02(tTx )-(Pm0(t'Tx )(m1(tTx )(m2(tTx ) (00(t,Ty )(01(tTy )(02(try )-Cpm0(tTy )(m1(try )(m2(try)

(00(tTz )(01(t'Tz )(02 (t' Tz )..(m0(t'Tz )(m1(tTz )(m2(tTz )

(1)

где

(t, r) = (t -10)0 r/0(t -10)0 rl (t -10)0 rf...(t - t0)mr0(t - t0)m t\ (t - t0)m rf;

тх, ту, т2 - независимая переменная по координатам X, Y, 2.

С учетом (1) описание вектор-функции г(^) с помощью полинома, определяющего положение ЛА и его координатных составляющих, имеет вид:

r(t,r, A) = (t,r)A -

(00(t, Tx )(01(t, Tx )(02(t, Tx )(m0(t, rx )(m1(t, Tx )(m2(t, Tx) (00(t,Ty)(01(t,Ty)(02(t,Ty) '(m0(t,Ty)(m1(t,Ty(m2(tTy) (00 (t,Tz )(01(t,T z )(02(t,T z )"'(m0(t,Tz (m1(tTz (m2(tT z)

A=

,2 2 ak(kl(trx) k = 01 = 0

my 2

k i 01 ?0 ak( (tTy)

,2 2 au(l(trz) k = 01 = 0

(t,rx ) A (t,ry ) A (f,rz ) A

(2)

ГДе AT \100а01а02 'ак0ак1ак2 -ат0ат1ат\

akl - коэффициенты сглаживающего полинома; I = 0, 1, 2;

m - степень сглаживающего полинома.

m о 2

mz 2

В представлении (2) координатные составляющие х(^тх,А), y(t,Ty,A), z(t,тz,A) будут формироваться с учетом всех составляющих вектора А.

Случайный характер траектории летательного аппарата вносит особенности по оптимальной оценки вектора А - коэффициентов сглаживающего полинома. При этом оптимальная оценка составляющих вектора коэффициентов сглаживающего полинома увязана с определением их состава.

Для этого необходимо получить статистически независимые оценки коэффициентов сглаживающего полинома. Так как случайный характер траекторий летательного аппарата трудно совместим с их высокой априорной определенностью, то прикладные методы информационной обработки данных траекторных измерений необходимо строить на основе статистических методов, которые не связанных с использованием априорной информации о распределении составляющих вектора А - коэффициентов сглаживающего полинома.

Так как векторы |-измерений первичных радиолокационных, кинотеодолитных данных измерений и коэффициентов сглаживающего полинома являются многомерными случайными величинами, то взаимное соответствие между ними будет определяться совместной плотностью вероятности [3].

Решение задачи по определению максимально правдоподобной оценки (МПО) вектора А - оценки вектора коэффициентов сглаживающего полинома целесообразно искать методом последовательных приближений, из-за нелинейной зависимости между вектором измерений первичных параметров случайной траектории от вектора коэффициентов сглаживающего полинома. Для совместной реализации пространственной и временной избыточности данных измерений, в работах [3-6] получен универсальный итеративный алгоритм:

Ау+1 = Ау + ДАУ = Ау+ОТлЬГ^ли - ^ Ау)]}, (3)

где J - матрица Якоби частных производных от измеряемых по вычисляемым параметрам; V - номер v-го приближения; JTдJv - основная матрица системы

уравнений на v-ом шаге приближения; Л - весовая матрица ошибок измерений.

Из выражения (3) видно, что СО достигается через ряд последовательных приближений.

Структура (1) может быть представлена таким образом:

Л0 X 1 X1 X1 . Ак X1 тЬк X1 Л* . X1 . А,™ X1 А,™ X1 А™ X1

А0 X 1 А0 X 1 А0 . X 1 . Ak X 1 тЧ X 1 Ак . X1 . А™ X1 А™ X1 т2..t.™ X1

Л0 X n А0 Xn А0 . Xn . лк Xn тЧ Xn т2> . Xn . Л™ Xn А™ Xn т2..t.™ Xn

А0 У 1 Л0 У1 А° . У1 . Af У1 Af У1 лк . У1 . А,™ У1 т11™ У1 т2™ У1

А0 У 1 А0 У 1 А0 . У 1 . Ak У 1 А4 У 1 Лк . У1 . Л™ У1 т11™ У1 т2..t.™ У1

А° У n А° У n Л0 . У n . А4 У n А4 У n Лк . У n . А™ У n т11™ У n т2..t.™ У n

z 1 д° z1 А° . z1 . Af z1 ?tk z1 Лк . z1 . Л™ z1 т11™ z1 т2™ z1

А0 z 1 А0 z 1 А0 . z 1 . Ak z 1 тЧ z 1 лк . z1 . А™ z1 тХ1™ z1 т2..t.™ z1

Л0 zn А0 zn А0 . zn . лк zn тЧ zn лк . zn . Л™ zn т1tm zn т2..t.™ zn

где - моменты времени на интервале сглаживания;

п - число точек на интервале сглаживания;

тх, ту, т2 - независимая переменная по координатам X, Y, Ъ, х=0, 1, 2 - индекс независимой переменной т по строке;

к=0, ..., т - степень сглаживающего полинома.

Известно [7-10], что для оценки вектора положения объекта в пространстве необходимо, чтобы определитель основной матрицы фТ ф системы уравнений не был равен нулю. Для этого необходимо, чтобы система базисных функций (4) являлась линейно независимой.

Из теории матричных вычислений [11-14] известно, что система векторов является линейно зависимой, если один из векторов линейно выражается через остальные векторы системы.

Из матрицы (4) видно, что она состоит из 3(ттах +1) столбцов и 3 • п строк

(где ттах - максимально возможная степень сглаживающего полинома, а

п - количество точек на интервале сглаживания), при этом система базисных функций будет линейно независимой, когда тж Ф ту Ф т2. Определим и исследуем область задания аргумента т, интервал его дискретизации дГх дГу дг_.

X•> у А

В работах [3, 4] рассмотрен способ построения Л-ортогональной базисной функции (Л-ОБФ) вида:

р00(^ V) Р02(^Рк1^..РтО^ Рт1(^ Рт2(<Г (5) с целью определения области задания интервала дискретизации аргумента т. Общий элемент р^ (^г) характерен тем, что недиагональные элементы основной

матрицы системы уравнений равны нулю.

В начале процесса построения Л-ОБФ примем следующие условия

р00(г,г) = ф00(Г /00 " ф00. В дальнейшем считаем:

Р01 (,г)= а00 01Р00 (*,Т) + Ф01 ^, г),

В данной формуле вспомогательный коэффициент определяется исходя из следующих условий:

1 шЛ]к1 = 0, (6)

где J - матрица Якоби, элементы которой = FPfcг; F - элемент матрицы проекций градиентов.

Для этого на базе ранее вычисленных матриц / и ф вычислим вектор-столбец /01 по формуле:

/01 = а00,01/00 + Ф01. Транспонируем матрицу /01 и умножим её справа на Л/00, а полученный

результат приравняем к нулю:

!Т Л

01Л 00 = а00,01"00Л 00 01Л 00

Л/™ ,+ ФТ = 0.

(7)

Из выражения (7) найдем значение неизвестного вспомогательного коэффициента а00,01:

a __Ф01Х700

00,01" JT / ■ / 00 00

Аналогично получим:

р02«,т) " «00,02Р00(^т) + «01,02 Р01(^т) + (2(',т)

и определим вспомогательные коэффициенты и .

Р «00,02 «01,02

I этого на базе известных матриц /

столбец:

Для этого на базе известных матриц JoQ, /01 и Ф^ вычислим вектор-

(10)

(11)

102 "«00,02/00 +«01,02/01 + Ф02. (9)

Транспонируем матрицу вектор-столбец , умножим её справа на Л/^ и,

приравняв полученный результат к нулю и согласно условия (6) получим:

Т Т Т _т

/02Л/00 " «00,02/00Л/00 + «01,02/01Л/00 + ф02Л/00 "0. (9)

Аналогично умножив / Т на Л/01 получим:

Т Т Т _т

702 Л/01 "«00,02,/00Л/01 +«01,02 /01Л/01 + Ф02 Л/01 " 0,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Согласно условия (6), выражения (9) и (10) упрощаются:

Т Т

«00,02/00Л/00 + Ф02Л/00 " 0

Т Т

«01,02 /01Л/01 + Ф02 Л/01 " а

Вычисляем коэффициенты из выражений (11)

« Ф02Л/00. « " Ф02Л/ 01 «00,02 " " т Л/ ' «01,02 " " /Т Л/ ■ /00 Л/00 01Л 01

Процесс повторяется до получения Р^От).

Если определена функция р и т) системы (5), то следующая функция

k (I—1)

Рк1 (^,т) будет найдена из предлагаемого нами рекуррентного соотношения:

k—1 2 к 1—1 Рк1 ^ "х^^Х^Х^ + £к А—«^1Р^,Т) + *к1 (Т

где

а фТ Л/х^ (12)

«х^М /Т. Л/ /

На базе ранее работ [15-17] можно показать, что между системами функций Р(^,т) и ((/, т) существует линейная связь:

Р(Г,т) "(Г,т)и,

где

и =

1 и00,01 и00,02

и.

00, к/

и,

00, т2

0

0

0

... и к/ ". и х^, т2

0 ... 1

верхняя треугольная матрица, диагональные элементы которой равны единице, а элементы, расположенные выше главной диагонали, вычисляются по формуле:

к /-1 к -1 2 и,, ,= ЁЕ и ^ а ,,+ ЕЕ и^ а ,, хКк/ р=кд±0 ХКРЧ РЧк р=0 Ч=0 ХКРЧ РЧ,к/

(13)

через вспомогательные коэффициенты а

Р4,к1

На основе ранее выполненных работ [12, 13] можно показать, что между (12) и (13) существует следующая связь:

1 и00,01 и 00,02

и

00, к/

и,

00, т2

и = 0

.. и

1 -а,

00,01

а

00,02

ХК, к/

00, к/

и

хК т2

00,т2

а ХК, к/ .. а хК, т2

0 ... 1

-1

Можно показать, что между векторами С и А существует связь , которая вытекает из тождественности следующих преобразований

Г( Ь, С) = (рС = (р1 С = (рии~ 1С = РА,

где I - единичная матрица; г - оценка вектора положения объекта в пространстве.

При решении задачи нелинейного сглаживания построение Л-ортогональной базисной функции необходимо производить на каждом шаге приближения к максимально правдоподобной оценке (МПО). Тогда на различных этапах последовательного приближения выражение (8) будет иметь вид:

Р1 = фи 1, Р2 = Р1и 2,

= ^у,

где V - номер последовательного приближения к МПО положения объекта в пространстве.

Из приведенных выше рассуждений следует общее выражение, позволяющее описывать линейную связь между Р и ф в виде:

Р^ = фП1и 2..и^=фиП,

0

0

1

0

1

0

0

0

0

1

0

0

1

0

0

0

где = ии2.. и у - матрица, полученная в результате перемножения и матриц

до v-го последовательного приближения включительно.

Для определения точности получения максимально правдоподобной оценки вектора оценки коэффициентов сглаживающего полинома, к которой приводит предложенный в работе [16] итеративный алгоритм через ряд последовательных приближений, предположим, что очередное вычисленное значение оценки вектора коэффициентов сглаживающего полинома совпало с истинным значением вектора

(Ау = А), и ошибки измерений относительно невелики. В практике траекторных

измерений данное допущение всегда выполняется.

Данные допущения позволяют определить отклонение вектора коэффициентов сглаживающего полинома на последнем шаге приближения формулой:

тТьГ )-1 тТ, у

AAV = Av- A = (JT AJ )-1 J T AAy

Так как систематические ошибки измерений нами исключены при обработке путем введения поправок, математическое ожидание ошибок будет иметь вид:

М (ДАу ) = (3Т Му ) -13Т ЛМ (Ду ),

а корреляционная матрица ошибок оценок коэффициентов сглаживающего полинома

К2 = М{[ДАу -М(ДАу)][ДАу -М(ДАу)]Т} = (3ТЛJу)-1.

А

Из выше приведенных рассуждений следует, что корреляционная матрица ошибок оценок коэффициентов сглаживающего полинома может быть получена обращением основной матрицы системы нормальных уравнений, вычисленной на последнем шаге приближения этих параметров [4, 18].

Для сохранения ранее введенных обозначений J для матрицы Якоби и <4 для вектора оценок коэффициентов сглаживающего полинома, полученных при независимом вычислении составляющих вектора оценок коэффициентов сглаживающего полинома при использовании ранее построенной Л-ОБФ, введем соответственно новые обозначения Ф и С для матрицы Якоби и МПО вектора коэффициентов сглаживающего полинома, получаемых при прямом решении системы уравнений с применением системы ЛНБФ [4].

Тогда вектор-поправка в новых обозначениях будет иметь вид

ДС у = (ФТ ЛФу )-1 ФТ ЛДу

и её можно рассматривать как результат решения матричной системы уравнений вида:

ФТ ЛФу ДСу = ФТ ЛД£, у , где Фу - матрица Якоби; Л - весовая матрица; Ду - отклонение вектора измерений от v-го приближения измеряемого вектора; ДСу - вектор-поправка коэффициентов полинома на v-м шаге приближения; фТЛФу - основная матрица системы уравнений на v-м шаге приближения.

Так как между вектором данных измерений £ и вектором вычисляемых параметров существует нелинейная зависимость, целесообразно искать Ф как матрицу последовательного преобразования [3]:

Ф = Fф,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где F - матрица пространственного преобразования (матрица проекций градиентов); ф - матрица временного преобразования.

При точном вычислении любым из методов недиагональные элементы равны нулю. На практике они отличаются от нуля из-за накопления ошибок вычислений, зависящих от способа построения ортогональной базисной функции, а это приводит к небольшой расстройке Л-ортогональности.

Приведем результаты исследования способа построения Л-ОБФ по показателю точности Qт.

Исследования производилось на различных интервалах сглаживания (15 и 25 точек) при различных степенях сглаживающего полинома (от 1 до 6) при удвоенной точности вычислений и обычной точности вычислений (для 25 точек на интервале сглаживания и 6-й степени сглаживающего полинома).

При обычной точности вычислений показатель точности Qт приблизительно в три раза хуже, чем в аналогичных условиях при удвоенной точности вычислений, время затраченное на Л-ортогонализацию системы базисной функции при выполнении одной итерации, одинаково, а объем занимаемой памяти на 1/3 меньше (табл. 1).

Таблица1

Показатель точности для структуры линейно независимых базисных функций двух переменных

Точек на интервале Степень ттах Показатель точности для структуры ЛНБФ двух переменных Точность вычисления на ЭВМ

15 6 9,31 • 10-14 Удвоенная

15 5 2,55^0-14 Удвоенная

15 4 1,1110-14 Удвоенная

15 3 6,25 10-15 Удвоенная

15 2 5,6910-15 Удвоенная

15 1 4,4910-15 Удвоенная

25 6 1, 0810-13 Удвоенная

25 5 5,3410-14 Удвоенная

25 4 2,1310-14 Удвоенная

25 3 8,9410-15 Удвоенная

25 2 9,3510-15 Удвоенная

25 1 7,1610-15 Удвоенная

25 6 4,77-10-4 Обычная

Метод диагонализации, улучшенный Якоби и переработанный Нейманом, позволяет вычислить собственные значения и собственные вектора действительной симметричной матрицы. В процессе вычислений заданная симметричная матрица не сохраняется, вычисленные собственные значения располагаются по главной диагонали заданной симметричной матрицы в убывающем порядке, а на недиагональных элементах располагаются промежуточные результаты вычислений. Вычисленные собственные вектора запоминаются по столбцам в том же порядке, что и собственные значения [12].

В связи с этим, непосредственно использовать результаты вычислений на ЭВМ, полученные для оценки точности вычислений по величине недиагональных элементов корреляционной матрицы ошибок оценок коэффициентов сглаживающего полинома, не представляется возможным [19].

Из работы [20] известно:

KÄ = Г - 1GtGTGt GT~ 1 = Г" 1 , (14)

где Г - диагональная матрица собственных значений основной матрицы системы T

Ф ЛФ, построенной с использованием структуры ЛНБФ [1];

G - матрица собственных векторов основной матрицы системы уравнений

ФТЛФ.

Однако в чистом виде матрица собственных значений не определяется. Воспользуемся для вычисления к развернутой частью выражения (14). Это позво-

A

лило учесть почти все ошибки вычислений для структуры ЛНБФ [3].

При построении Л-ОБФ максимально правдоподобная оценка вектора коэффициентов сглаживающего полинома достигается в 2-3 итерации.

Заключение. Экспериментальное исследование структуры ЛНБФ [3] на ЭВМ при нахождении области определения аргумента т0 и интервалов дискретизации Atj, ^Ty,^z полностью подтвердили теоретические выводы. Детально

структура ЛНБФ исследовалась в работах [15,-17, 21, 22].

Из результатов исследований, приведенных в этих работах, следует, что показатель точности Q,^ как функция двух аргументов (т0 - наименьшее значение аргумента и Дт - интервал дискретизации аргумента т) в широком диапазоне значений этих аргументов изменяется слабо, но резко увеличивается при .

При этом величины и не должны превосходить соответст-

венно максимально и минимально возможных чисел, которые можно без потери точности записать в разрядную сетку применяемой ЭВМ. При равномерном шаге дискретизации аргумента т аргумент целесообраз-

__ps

но выбирать в средней части интервала ( ^qmin,—^^-) (где qmin и qmах соответственно минимальные и максимальные числа, которые без потери точности можно записать в разрядную сетку ЭВМ). Приближение к краям интервала может при неблагоприятных условиях привести к возрастанию ошибок вычислений при определении вторичных параметров положения летательного аппарата из-за того, что основная матрица системы уравнений становится плохо обусловленной.

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Мильштейн А.В., Паслен В.В. Метод нелинейного сглаживания в обработке данных тра-екторных измерений // Сб. научных трудов. Вып. 28. - Донецк: Донецкий институт же-лезнождорожного транспорта, 2011. - С. 93-101.

2. Мильштейн А.В., Мотылев К.И., Паслен В.В. Исследование структур базисных функций // Сб. научных трудов. Вып. 29. - Донецк: Донецкий институт железнождорожного транспорта, 2012. - С. 23-30.

3. ОгоднийчукН.Д. Обработка траекторной информации. Ч. 1. - К.: КВВАИУ, 1981. - 141 с.

4. Огоднийчук Н.Д. О прикладных методах анализа траекторной информации // Сб. материалов НТК, посвященной 25-летию училища. Ч. 1. - К.: КВВАИУ, 1977. - С. 65-84.

5. Щербов И.Л. Исследование алгоритма адаптивного нелинейного оптимального сглаживания многопараметрических данных измерений // Информатика и кибернетика. - 2020.

- № 4 (22). - С. 5-12.

6. ЩербовИ.Л., Паслен В.В., Мотылев К.И., Михайлов М.В. Методы обработки данных измерений, которые владеют пространственной и временной избыточностью // Системные технологии: Региональный межвузовский сборник научных работ. Вып. 5 (46).

- Днепропетровск, 2006. - С. 95-100.

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

8. ОгоднийчукН.Д. Обработка траекторной информации. Ч. II. - К.: КВВАИУ, 1986. - 224 с.

9. Светозаров В.В. Основы статистической обработки результатов измерений. - М.: МИФИ, 2005. - 54 с.

10. СавчукВ.П. Обработка результатов измерений. - Одесса: ОНПУ, 2002. - 54 с.

11. Алгазин С.Д. Численные алгоритмы без насыщения в классических задачах математической физики. - М.: Научный Мир, 2002. - 155 с.

12. Грантмахер Ф.Р. Теория матриц. - М.: Физматлит, 2004. - 560 с.

13. Ланкастер П. Теория матриц. - М.: Наука. Глав. ред. физ.-мат. лит., 1962. - 280 с.

14. ВолковЕ.А. Численные методы. - М.: Наука. Глав. ред. физ.-мат. лит., 1987. - 248 с.

15. ЩербовИ.Л., Паслен В.В., Михайлов М.В., Мотътев К.И., Лебеденко Д.М., Антикуз А.Г. О построении ортогональных базисных функций // Туполевские чтения: Международная молодежная научная конференция, посвященная 1000-летию города Казани: Матер. конференции. Т. II. - Казань: Казан. гос. техн. ун-та, 2005. - С. 139-140.

16. Щербов И.Л., Паслен В.В. Обработка данных траекторного контроля с использованием ортогональных базисных функций // Вестник Академии гражданской защиты. - 2021.

- Вып. 1 (25).- С. 48-53.

17. Щербов И.Л. Информационная технология обработки данных траекторного контроля // Вестник Донецкого национального университета. Сер. Г: Технические науки. - 2021.

- № 1. - С. 71-77.

18. Огороднийчук Н.Д. Оценка параметров движения и точность их определения по данным коррелированных траекторных измерений. - К.: КВВАИУ, 1968. - 216 с.

19. Баклицкий, В.К., Юрьев А.Н. Корреляционно-экстремальные методы навигации. - М.: Радио и связь, 1982. - 256 с.

20. Дубровский С.А. Прикладной многомерный статистический анализ. - М.: Финансы и статистика, 1982. - 216 с.

21. Щербов И.Л., Паслен В.В., Фесенко Д.В., Михайлов М.В. Причины возникновения сбоев и методы их устранения // Международная молодежная научно-практическая конференция «Человек и космос». - Днепропетровск: НЦАОМУ, 2006. - С. 120.

22. Огоднийчук Н.Д., Паслен В.В., Велигдан С.В. Исследования на ЭВМ свойств систем ЛНБФ и Л-ОБФ как функции двух аргументов // Радиоэлектронное оборудование летательных аппаратов. Вып. 3. - К.: КВВАИУ, 1989. - С. 90-93.

REFERENCES

1. Mil'shteyn A.V., Paslen V.V.Metod nelineynogo sglazhivaniya v obrabotke dannykh traektornykh izmereniy [Method of non-linear smoothing in data processing of trajectory measurements], Sb. nauchnykh trudov [Collection of scientific papers]. Issue 28. Donetsk: Donetskiy institut zheleznozhdorozhnogo transporta, 2011, pp. 93-101.

2. Mil'shteyn A. V., Motylev K.I., Paslen V. V. Issledovanie struktur bazisnykh funktsiy [Investigation of the structures of basic functions], Sb. nauchnykh trudov [Collection of scientific papers]. Issue 29. Donetsk: Donetskiy institut zheleznozhdorozhnogo transporta, 2012, pp. 23-30.

3. Ogodniychuk N.D. Obrabotka traektornoy informatsii [Processing of trajectory information]. Part 1. Kiev: KVVAIU, 1981, 141 p.

4. Ogodniychuk N.D. O prikladnykh metodakh analiza traektornoy informatsii [On applied methods for analyzing trajectory information], Sb. materialov NTK, posvyashchennoy 25-letiyu uchilishcha [Collection of materials of the Scientific and Technical Committee dedicated to the 25th anniversary of the school]. Part 1. Kiev: KVVAIU, 1977. - S. 65-84.

5. Shcherbov I.L. Issledovanie algoritma adaptivnogo nelineynogo optimal'nogo sglazhivaniya mnogoparametricheskikh dannykh izmereniy [Study of the adaptive nonlinear optimal smoothing algorithm for multi-parameter measurement data], Informatika i kibernetika [Informatics and cybernetics], 2020, No. 4 (22), pp. 5-12.

6. Shcherbov I.L., Paslen V.V., Motylev K.I., Mikhaylov M.V. Metody obrabotki dannykh izmereniy, kotorye vladeyut prostranstvennoy i vremennoy izbytochnost'yu [Measurement data processing methods that have spatial and temporal redundancy], Sistemnye tekhnologii: Regional'nyy mezhvuzovskiy sbornik nauchnykh rabot [System Technologies. Regional interuniversity collection of scientific papers]. Issue 5 (46). Dnepropetrovsk, 2006, pp. 95-100.

7. Zhdanyuk B.F. Osnovy statisticheskoy obrabotki traektornoy informatsii [Fundamentals of statistical processing of trajectory measurements]. Moscow: Sov. radio, 1978, 384 p.

8. Ogodniychuk N.D. Obrabotka traektornoy informatsii [Processing of trajectory information]. Part II. Kiev: KVVAIU, 1986, 224 p.

9. Svetozarov V.V. Osnovy statisticheskoy obrabotki rezul'tatov izmereniy [Fundamentals of statistical processing of measurement results]. Moscow: MIFI, 2005, 54 p.

10. Savchuk V.P. Obrabotka rezul'tatov izmereniy [Processing of measurement results]. Odessa: ONPU, 2002, 54 p.

11. Algazin S.D. Chislennye algoritmy bez nasyshcheniya v klassicheskikh zadachakh matematicheskoy fiziki [Unsaturated numerical algorithms in classical problems of mathematical physics]. Moscow: Nauchnyy Mir, 2002, 155 p.

12. Grantmakher F.R. Teoriya matrits [Matrix Theory]. Moscow: Fizmatlit, 2004, 560 p.

13. LankasterP. Teoriya matrits [Matrix Theory]. Moscow: Nauka. Glav. red. fiz.-mat. lit., 1962, 280 p.

14. Volkov E.A. Chislennye metody [Numerical methods]. Moscow: Nauka. Glav. red. fiz.-mat. lit., 1987, 248 p.

15. Shcherbov I.L., Paslen V.V., Mikhaylov M.V., Motylev K.I., Lebedenko D.M., Antikuz A.G. O postroenii ortogonal'nykh bazisnykh funktsiy [About the construction of orthogonal basis functions], Tupolevskie chteniya: Mezhdunarodnaya molodezhnaya nauchnaya konferentsiya, posvyashchennaya 1000-letiyu goroda Kazani: Mater. Konferentsii [Tupolev readings: International youth scientific conference dedicated to the 1000th anniversary of the city of Kazan. Conference materials]. Vol. II. Kazan': Kazan. gos. tekhn. un-ta, 2005, pp. 139-140.

16. Shcherbov I.L., Paslen V.V. Obrabotka dannykh traektornogo kontrolya s ispol'zovaniem ortogonal'nykh bazisnykh funktsiy [Trajectory control data processing using orthogonal basis functions], VestnikAkademii grazhdanskoy zashchity [Vestnik of the Academy of Civil Protection], 2021, Issue 1 (25), pp. 48-53.

17. Shcherbov I.L. Informatsionnaya tekhnologiya obrabotki dannykh traektornogo kontrolya [Information technology of mathematical modeling of data processing of trajectory control], Vestnik Donetskogo natsional'nogo universiteta. Ser. G: Tekhnicheskie nauki [Vestnik of the Donetsk National University. Series G: Technical sciences], 2021, No. 1, pp. 71-77.

18. Ogorodniychuk N.D. Otsenka parametrov dvizheniya i tochnost' ikh opredeleniya po dannym korrelirovannykh traektornykh izmereniy [Evaluation of motion parameters and the accuracy of their determination according to data of correlated trajectory measurements]. Kiev: KVVAIU, 1968, 216 p.

19. Baklitskiy, V.K., Yur'ev A.N. Korrelyatsionno-ekstremal'nye metody navigatsii [Correlation-extremal methods of navigation]. Moscow: Radio i svyaz', 1982, 256 p.

20. Dubrovskiy S.A. Prikladnoy mnogomernyy statisticheskiy analiz [Applied multivariate statistical analysis]. Moscow: Finansy i statistika, 1982, 216 p.

21. Shcherbov I.L., Paslen V.V., Fesenko D.V., Mikhaylov M.V. Prichiny vozniknoveniya sboev i metody ikh ustraneniya [Causes of failures and methods of their elimination], Mezhdunarodnaya molodezhnaya nauchno-prakticheskaya konferentsiya «Chelovek i kosmos» [International Youth Scientific and Practical Conference "Man and Space"]. Dnepropetrovsk: NTSAOMU, 2006, pp. 120.

22. OgodniychukN.D., Paslen V.V., Veligdan S.V. Issledovaniya na EVM svoystv sistem LNBF i Л-OBF kak funktsii dvukh argumentov [Computer studies of the properties of the LNBF and Л-OBF systems as functions of two arguments], Radioelektronnoe oborudovanie leta-tel'nykh apparatov [Radio-electronic equipment of aircraft]. Issue 3. Kiev: KVVAIU, 1989, pp. 90-93.

Статью рекомендовал к опубликованию д.т.н., профессор А.В. Боженюк.

Щербов Игорь Леонидович - Донецкий национальный технический университет; e-mail:

scherbov@yandex.ru; r. Донецк, ДНР; проректор.

Shcherbov Igor Leonidovich - Donetsk National Technical University; e-mail: scherbov@yandex.ru;

Donetsk, DPR; vice rector.

i Надоели баннеры? Вы всегда можете отключить рекламу.