К настройке алгоритма вычисления барометрической высоты летательного аппарата системой воздушных сигналов Текст научной статьи по специальности «Строительство и архитектура»

УДК 629.7.054.44

DOI: 10.14529/ctcr200106

К НАСТРОЙКЕ АЛГОРИТМА ВЫЧИСЛЕНИЯ БАРОМЕТРИЧЕСКОЙ ВЫСОТЫ ЛЕТАТЕЛЬНОГО АППАРАТА СИСТЕМОЙ ВОЗДУШНЫХ СИГНАЛОВ

В.И. Панферов1'2, Н.А.Тренин1, С.В. Панферов2, А.М. Хаютин1, Д.В. Пауков1, К.С. Халдин3

1 Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия им. проф. Н.Е. Жуковского и Ю.А. Гагарина», филиал в г. Челябинске, г. Челябинск, Россия,

2 Южно-Уральский государственный университет, г. Челябинск, Россия,

3 Targem Games, центр разработки компьютерных игр, г. Екатеринбург, Россия

Введение. В условиях повышения требований к точности определения высотно-скоростных параметров летательных аппаратов задача настройки алгоритмического обеспечения их систем воздушных сигналов (СВС) по реальным полетным данным является вполне актуальной. Цель исследования. Рассмотреть задачу оценки реального значения градиента температуры по высоте атмосферы для целей оперативной подстройки штатного алгоритмического обеспечения СВС. Материалы и методы. Традиционный подход заключается в использовании численного значения градиента температуры, принятого для международной стандартной атмосферы (МСА). Особенность предлагаемого подхода заключается в том, что градиент температуры по высоте определяется по текущим измерениям атмосферного давления и температуры, а не высоты и температуры. Результаты. Разработаны адаптивный и неадаптивный алгоритмы параметрической идентификации. Предложен усовершенствованный вариант адаптивного алгоритма, обеспечивающий должную настройку его фильтрующих свойств на конкретную статистическую обстановку. Проведена апробация алгоритмов как по характеристикам МСА с применением компьютерного генерирования случайных помех, так и по известным в литературе результатам реальных измерений, полученным с помощью метеорологической ракеты. Выполнено сравнение свойств адаптивного и неадаптивного алгоритмов. Показано, что многократная «прогонка» одних и тех же экспериментальных данных через адаптивный алгоритм приводит практически к одинаковой точности оценки градиента, что и при использовании неадаптивного алгоритма. Заключение. Реализация алгоритмов способствует повышению точности определения барометрической высоты воздушного судна (ВС), не требует какой-либо модернизации аппаратной части СВС, необходимо только некоторое дополнение алгоритмического обеспечения системы. Адаптивный алгоритм рекомендуется к применению в случае существенного дрейфа характеристик атмосферы для оперативной подстройки системы.

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

Постановка задачи

Известно [1], что для достаточно точного определения высоты воздушного судна (ВС) барометрическим методом необходимо знать, как реально по высоте атмосферы изменяется ее температура. Известно также, что в реальных условиях эта температура обычно убывает с ростом высоты. Причем скорость этого убывания непредсказуемо меняется как в течение суток, так и в течение года, а также она различна в разных пунктах и на разных высотах. Кроме того, в некоторых случаях в определенном диапазоне высот температура может не убывать, а, наоборот, возрастать. Таким образом, зависимость температуры T от высоты h, то есть T = T (И) , почти всегда неизвестна, да и непрерывно меняется. Поэтому при построении теории барометрических высотомеров используют такое понятие, как международная стандартная атмосфера (МСА). Это условное понятие. В МСА считается, что температура по высоте атмосферы изменяется в соответствии с формулой

(1)

где T0 - температура воздуха на некоторой высоте h0, у= 0,0065 К/м - градиент температуры для стандартной атмосферы. Также в МСА принято, что на уровне моря атмосферное давление составляет всегда 760 мм рт. ст. (1013,2 гПа), а температура равна +15 °C (288 К), т. е. считается, что

T = 288 - 0,0065 (h - h0 ),

(2)

а если при этом h0 = 0, то высота ВС в этом случае называется абсолютной.

Заметим, что именно формула (2) «заложена» в конструкцию механических барометрических высотомеров [1]. Кремальера прибора позволяет конкретизировать только величину Н0) с помощью указания того давления Ро , какое есть в МСА на уровне Н0). Таким образом, температурная погрешность механических высотомеров обуславливается как тем, что его конструкцией «считается», что на уровне моря температура всегда равна 288 К, так и тем, что эта температура уменьшается на 6,5 градуса на каждый километр высоты.

В современных системах воздушных сигналов (СВС) используется следующая зависимость [1-10]

T = T0 - 0,0065 (h - h0),

(3)

т. е. имеется возможность указать по факту температуру Т0 , которая имеет место на уровне Н0) (Н0), как и в механических высотомерах, задается давлением р0)). Таким образом, методическая погрешность определения высоты СВС будет заметно меньше той погрешности, какая есть у механических высотомеров. Тем не менее остается погрешность, обусловленная тем, что в СВС по-прежнему «считается», что у= 0,0065 К/м, что в реальных условиях может быть совсем не так. Поэтому возникает задача разработки алгоритма оценки реального значения параметра у по тем данным, которые достаточно надежно измеряются в процессе полета с целью оперативной подстройки алгоритмического обеспечения СВС, задача совершенствования которого является вполне актуальной [2-10].

Решение задачи

Известно [1], что барометрическая высота ВС определяется в соответствии с формулой

р=^ f1-1^ h

у R

(4)

Здесь знак Т означает возведение в степень, а ^ и R - соответственно ускорение свободного падения и газовая постоянная атмосферного воздуха. Учитывая, что имеет место зависимость (1), данную формулу перепишем в виде

у ln

fp} Р0

= g ln R

f T ^ V T0 J

(5)

Далее предположим, что есть экспериментальные точки (рЭ; Т^) - это определенные в процессе полета реальные значения давления и температуры на разных высотах. Эта задача весьма просто решается с помощью СВС. После этого решим следующую задачу оптимизации:

'=1

г=1

у ln

( Э ^ Рэ_

Р0

- g ln R

( T э ^

->min.

(6)

Здесь п - количество экспериментальных точек. Таким образом, требуется найти такое значение параметра у , при котором сумма квадратов невязок между правой и левой частями уравнения (5) с подставленными в него экспериментальными данными будет минимальной.

Решая задачу (6), определим, что оптимальное значение параметра у следует вычислять по соотношению

2

У=

I

g '=1

ln

э

т

ln

f Э"\ Ро

R

Ijln

i=1

рЭ Л

Ро

(7)

Рассмотренная задача - это задача параметрической идентификации (параметрической настройки) модели Земной атмосферы на экспериментальные данные по таким ее характеристикам, как температура, давление и высота (уравнения (1), (4) (5)). Решение (7) - это решение данной задачи неадаптивным методом наименьших квадратов (МНК) [11]. Однако представляет интерес решение задачи и адаптивным методом наименьших квадратов, который позволяет уточнять значение идентифицируемого параметра у по мере получения новой экспериментальной информации [11-14]. Существо адаптивного метода заключается в следующем.

Пусть известно значение идентифицируемого параметра на i - м шаге идентификации. Пусть

получена новая экспериментальная точка (рТч\) . Требуется определить «новое» откорректированное значение параметра уг-+1 через у 1 и (рЭ+1; ТЭ1) . Формально задача записывается так:

(ргэ+1; ТЭ1; у,) ^уг+1. (8)

Запишем выражение для квадрата локальной невязки между правой и левой частями уравнения (5) - локальное значение критерия I.

I,+i =

У,ln

Ро

- g ln R

f тЭ \ T+i

T

Вычислим производную

51

i+1

ду,

=2

У,ln

( РЭ+1 ^ Ро

-g ln R

( тЭ Л

т

i+1

т

ln

( Рйл Ро

(9)

(1Q)

и будем определять «новое» откорректированное значение параметра по следующему соотношению:

У,+1 = У,- S,+i

д I

,+1

дУ,

(11)

где 51+1 - некоторый коэффициент. Оптимальное значение этого коэффициента найдем из условия достижения минимума локальной невязкой после подстановки в ее выражение формулы (11), т. е. 1

S,

,+1ОПТ

2

2 '

(12)

ln

V v

Ро

/У

Если использовать полученное оптимальное значение коэффициента 5г+1ОПТ. что алгоритм коррекции будет иметь вид

то окажется,

У,+1 = У,

У,ln

э \ Р+1

Ро

- g ln R

f тЭ \ T+1

v то У

ln

( РЭ Л Р+1

Ро

(13)

Если упростить это соотношение, то получим

У ,+1 =

g ln ( т э T+1 Л

R V T У

f э Л

ln Р,+1

V ро У

(14)

2

2

Как видно из формулы (14), алгоритм коррекции не имеет «памяти» и с точки зрения помехоустойчивости такое свойство алгоритма (14) считается плохим [11-14]. Данную формулу можно было бы получить, и напрямую отыскивая минимум локальной невязки (9), однако проделанные выкладки позволяют сформировать структуру алгоритма, имеющую большую помехоустойчивость. Для этого несколько откорректируем формулу для 5г-+1 и будем полагать, что

°'5 (15)

+1 - "

(

а +

1п

V V

Ро

/V

где а - некоторый параметр, 0 <а«х> - используется для повышения помехоустойчивости алгоритма. В этом случае алгоритм коррекции будет иметь вид

У/+1 - Уг

Э ( т э V Т+1 ( Э V

у г 1п Рг+1 - £ 1п 1п Рг+1

V р0 V К V Т0 1 V Р0 V

а +

1п

( Э ^

Рг+1

(16)

V V

Ро

/V

Если сравнить соотношения (13) и (16), то вполне очевидно, что по алгоритму (16) параметр у будет исправляться от шага к шагу заметно осторожнее. В новом вычисленном значении этого

параметра уг+ будет содержаться информация как о прежнем его значении уг, так и «новая» информация об экспериментальных значениях давления и температуры (рЭ_1; Т+1) .

2

Особенности предлагаемого подхода

Если искать решение задачи определения градиента температуры по высоте у по экспериментальным данным, как говорится «в лоб», то в первую очередь возникает идея напрямую использовать для этого уравнение Т - Т0 -у(Ъ - Ъ0). Однако это уравнение в данном случае непригодно для решения задачи, так как для этого необходимо иметь достаточно точные экспериментальные значения высоты Ъ , что в принципе невозможно, поскольку вся теория измерения барометрической высоты основана на предположении о том, что у уже известно. Радиовысотомер для этой цели также непригоден, так как даже при горизонтальном полете он будет выдавать переменный сигнал из-за изменения рельефа пролетаемой местности. Особенность предлагаемого подхода к решению задачи заключается в том, что для идентификации параметра у , необходимого для подстройки алгоритмического обеспечения СВС, используются измеренные той же СВС значения давления и температуры на разных высотах. Кроме того, такой подход позволяет не нарушать свойство автономности СВС.

Апробация алгоритмов идентификации

Для апробации адаптивного и неадаптивного алгоритмов идентификации были разработаны специальные компьютерные программы.

Если с помощью алгоритма (7) обработать данные, содержащиеся в характеристиках СТА, то вычисленное значение у в точности равно его заявленному в литературе значению для СТА. Кроме того, если эти же данные использовать для апробации алгоритма (13), то сразу же на первой итерации ошибка в задании начального значения у исправляется и последующие вычисляемые его значения в точности равны 0,0065 К/м .

Далее при исследовании свойств предлагаемых алгоритмов идентификации на эти же данные накладывались помехи, генерируемые датчиком случайных чисел, качество работы алгоритма (7) в этом случае получается также практически удовлетворительным. Так, при максимальной помехе в определении давления в ±50 Па и помехе в определении температуры в ±5 °С погрешность в определении у не превысила по модулю 0,000022 К/м, что весьма удовлетворительно.

Алгоритм (13), или что то же самое (14), в данной ситуации настолько неприемлем, насколь-

ко велики помехи в сигналах давления и температуры. Причем погрешность в определении температуры оказывает большее влияние на погрешность оценки у, чем погрешность измерения давления. Так, при максимальной случайной погрешности в определении давления в ±20 Па и случайной помехе в определении температуры в ±1 °С максимальная погрешность оценки у по модулю составила 0,0023 К/м . При увеличении случайной помехи в измерении температуры до ±2 °С максимальная погрешность оценки у по модулю выросла до 0,0046 К/м . Если же случайную погрешность в измерении температуры оставить на уровне ±1 °С, а случайную погрешность в измерении давления уменьшить с ±20 до ±10 Па, то погрешность в определении у практически

не изменится - ее изменение составило всего лишь 9 • 10"7 К/м .

Работа алгоритма (16) иллюстрируется рис. 1.

8 7,5

О 2 4 6 8 10 12

Номер точки

Рис. 1. Отслеживание численного значения градиента температуры Fig. 1. Tracking the numerical value of the temperature gradient

На рис. 1 кривая, отмеченная «треугольниками», - это оценка градиента температуры по алгоритму (13), или что то же самое по алгоритму (14), или что то же самое по алгоритму (16), при параметре а=0 ; три другие кривые - это оценки градиента температуры по алгоритму (16), при этом для кривой, отмеченной «кружками», а = 0,25; для кривой, отмеченной «квадратиками», а = 0,5; для кривой без отметок (марок) а = 1,0. Здесь погрешность в измерениях давления и температуры генерировалась с помощью датчика случайных чисел в диапазонах [-50, + 50] Па и [—1, +1] К соответственно.

Как видно из рис. 1, подбирая численное значение параметра а для конкретной статистической обстановки, можно добиться удовлетворительного качества работы адаптивного алгоритма (16). Как видно из рис. 1, алгоритм (16) с а = 0,25 вполне приемлем для отслеживания параметра у в имеющейся статистической обстановке. Увеличение численного значения параметра а приводит к замедлению скорости отслеживания параметра у .

Для апробации алгоритмов идентификации были использованы также и экспериментальные данные работы [15]. По этим данным найденное неадаптивным МНК (алгоритм (7)) численное значение градиента температуры составило у = 0,00693 К/м, что подтверждает возможность отличия реальных значений этого параметра от значения МСА. Кроме того, экспериментальные данные этой же работы обрабатывались также и адаптивным алгоритмом (16), результаты представлены на рис. 2.

Рис. 2. Отслеживание численного значения градиента температуры по экспериментальным данным Fig. 2. Track the numerical value of the temperature gradient from experimental data

Здесь также для кривой, отмеченной «кружками», а = 0,25; для кривой, отмеченной «квадратиками», а = 0,5; для кривой без отметок (марок) а = 1,0 .

Как видно из рис. 2, здесь статистическая обстановка такова, что ошибка в задании начального значения параметра у практически удовлетворительно исправляется уже по 6 точкам, причем лучшая скорость отслеживания наблюдается при а = 0,25. Как известно [16], общих рекомендаций по выбору параметра а дать нельзя, необходима «...подгонка его к конкретной статистической обстановке»[16, с. 64].

В работе также исследовался тот факт, что многократная прогонка одних и тех же экспериментальных данных через адаптивный алгоритм позволяет добиться того же качества идентификации, что и при использовании неадаптивного алгоритма [11]. В данном случае применялись двукратная и трехкратная прогонка 6 экспериментальных точек. На рис. 3 приведены результаты двукратной прогонки экспериментальных данных.

7,5 — —— —— —— —— —— —

>

4,5

4 ---1-1--

О 2 4 6 8 10 12

Номер точки

Рис. 3. Двукратная прогонка экспериментальных данных через адаптивный алгоритм Fig. 3. Double run of experimental data via adaptive algorithm

Как видно из рис. 3, при тех значениях параметра а , что и для случаев на рис. 1 и 2, все три кривые к концу уже второй прогонки практически «сливаются», а отслеживаемое значение параметра у практически полностью совпадает со значением, найденным неадаптивным МНК. На рис. 4 приведены результаты трехкратной прогонки экспериментальных данных.

Рис. 4. Трехкратная прогонка экспериментальных данных через адаптивный алгоритм Fig. 4. Three-fold run of experimental data via adaptive algorithm

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

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

Выводы

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

Литература

1. Панферов, В.И. Авиационные приборы и пилотажно-навигационные комплексы: учеб. пособие: в 3 ч. /В.И. Панферов, Н.А. Тренин, А.М. Хаютин. - Челябинск: филиал ВУНЦВВС «ВВА», 2018. - Ч. I. - 145 с.

2. Козицин, В.К. Система воздушных сигналов вертолета на основе свободно ориентированного приемника давлений: автореф. ... дис. канд. техн. наук / В.К. Козицин. - Казань: Изд-во КГТУ им. А.Н. Туполева, 2006. - 23 с.

3. Василюк, Н.Н. Повышение точности ракетного барометрического высотомера с учётом аэродинамики и итогов радиолокационных наблюдений в лётных экспериментах: автореф. ... дис. канд. физ-мат. наук /Н.Н. Василюк. - М. : Физтех-полиграф, 2004. - 23 с.

4. Солдаткин, В.В. Алгоритмы формирования и обработки информации системы воздушных сигналов вертолета на основе неподвижного многофункционального аэрометрического приемника / В.В. Солдаткин // Фундаментальные и прикладные проблемы техники и технологии. -2013. - № 3 (299). - С. 110-120.

5. Крылов, Д.Л. Функциональная схема, алгоритмы обработки информации и погрешности системы воздушных сигналов самолета с неподвижными невыступающими приемниками потока / Д.Л. Крылов, В.М. Солдаткин // Новые технологии, материалы и оборудование российской авиакосмической отрасли: сб. докл. Всерос. науч.-практ. конф. - Казань: Изд-во Академии наук Республики Татарстан. - 2016. - С. 570-574.

6. Алгоритмическое обеспечение системы воздушных сигналов одновинтового вертолета с неподвижным аэрометрическим приемником на характерных режимах эксплуатации / В.В. Солдаткин, И.Ф. Мингазов, А.Р. Мустафин, Д.Н. Нурутдинова // Труды Международной научно-технической конференции «Перспективные информационные технологии». - Самара: Изд-во Самарского научного центра РАН. - 2018. - С. 881-885.

7. Ефремова, Е.С. Алгоритмическое обеспечение вихревой системы воздушных сигналов дозвукового самолета / Е.С. Ефремова // Труды Международной научно-техн. конф. «Перспективные информационные технологии». - Самара: Изд-во Самарского научного центра РАН. - 2018. -С. 812-816.

8. Солдаткин, В.В. Автоматическая подстройка измерительных каналов системы воздушных сигналов вертолета / В.В. Солдаткин // Вестник КГТУ им. А.Н. Туполева. - 2004. - № 2. -С. 26-29.

9. Kaletka, J. Evaluation of the HelicoPter Low AirsPeed System Lassie / J. Kaletka // Journal of American HelicoPter Society. - 1983. - No. 4. - Р. 35-43. DOI: 10.4050/JAHS.28.35

10. Leishman J.G. Airdynamic Interactions Bethcen a Rotor and a Fuselage in Forward Flight / J.G. Leishman, B.I. NaiPei // J.Amer. HelicoPter. Soc. - 1990. - Vol. 35, no. 3. - P. 22-31. DOI: 10.4050/JAHS.35.22

11. Растригин, Л.А. Современные принципы управления сложными объектами / Л.А. Рас-тригин. - М. : Сов. радио, 1980. - 232 с.

12. Панферов, В.И. Об адаптации системы управления температурным режимом протяжных термоагрегатов / В.И. Панферов, Г.П. Кулаченков //Изв. вузов. Черная металлургия. - 1990. -№ 10. - С. 72.

13. Identifying linear models from time domain simulations / J.J. Sanchez-Gasca, K. Clark, N.W. Miller et al. // IEEE ComPuter APPlications in Power. - APr. 1997. - Vol. 10, no. 2. - P. 26-30. DOI: 10.1109/67.582448

14. Hemerly, E.E. PC-based Packages for identification, oPtimization, and adaPtive control / E.E. Hemerly // IEEE Control Systems Magazine. - Feb. 1991. - Vol. 11, iss. 2. - P. 37-43. DOI: 10.1109/37.67674

15. Ракетные исследования атмосферы / П.П. Алексеев, Е.А. Бесядовский, Г.И. Голышев и др. //Метеорология и гидрология. - 1957. - № 8. - С. 3-13.

16. Салыга, В.И. Идентификация и управление процессами в черной металлургии / В.И. Са-лыга, Н.Н. Карабутов. - М. : Металлургия, 1986. - 192 с.

Панферов Владимир Иванович, д-р техн. наук, профессор, профессор кафедры авиационных комплексов и конструкций летательных аппаратов, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия им. проф. Н.Е. Жуковского и Ю.А. Гагарина», филиал в г. Челябинске; профессор кафедры информационно-аналитического обеспечения управления в социальных и экономических системах, Южно-Уральский государственный университет, г. Челябинск; tgsiv@mail.ru.

Тренин Николай Александрович, канд. воен. наук, начальник 13 кафедры авиационных комплексов и конструкций летательных аппаратов, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия им. проф. Н.Е. Жуковского и Ю.А. Гагарина», филиал в г. Челябинске, г. Челябинск.

Панферов Сергей Владимирович, канд. техн. наук, доцент кафедры градостроительства, инженерных сетей и систем, Южно-Уральский государственный университет, г. Челябинск; panferovsv@susu.ru.

Хаютин Андрей Михайлович, канд. техн. наук, заместитель начальника 13 кафедры авиационных комплексов и конструкций летательных аппаратов, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия им. проф. Н.Е. Жуковского и Ю.А. Гагарина», филиал в г. Челябинске, г. Челябинск.

Пауков Дмитрий Викторович, преподаватель 13 кафедры авиационных комплексов и конструкций летательных аппаратов, Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия им. проф. Н.Е. Жуковского и Ю.А. Гагарина», филиал в г. Челябинске, г. Челябинск.

Халдин Константин Сергеевич, канд. техн. наук, программист, Targem Games, центр разработки компьютерных игр, г. Екатеринбург; lar3811@yandex.ru.

Поступила в редакцию 4 ноября 2019 г

DOI: 10.14529/ctcr200106

TO THE SOLUTION OF THE PROBLEM

ON THE DETERMINATION OF BAROMETRIC HEIGHT

OF THE AIRCRAFT SYSTEM OF AIR SIGNALS

V.I. Panferov1'2, tgsiv@mail.ru, N.A. Trenin1,

S.V. Panferov2, panferovsv@susu.ru, A.M. Hayutin1, D.V. Paukov1,

K.S. Haldin3, lar3811@yandex.ru

1 Russian Air Force Military Educational and Scientific Center "Air Force Academy named after Professor N.E. Zhukovsky and Y.A. Gagarin", Chelyabinsk branch, Chelyabinsk, Russian Federation,

2 South Ural State University, Chelyabinsk, Russian Federation,

3 Targem Games, Computer Games Development Center, Ekaterinburg, Russian Federation

Introduction. In conditions of increasing requirements for the accuracy of determining the altitude and speed parameters of aircraft, the task of tuning up the algorithmic support of their air signal systems (SHS) based on real flight data is quite relevant. Aim. Consider the task of assessing the real value of the temperature gradient along the height of the atmosphere for the purpose of operational adjustment of the standard algorithmic support of SHS. Materials and methods. The traditional approach is to use the numerical value of the temperature gradient, adopted for the international standard atmosphere (ISA). A feature of the proposed approach is that the temperature gradient in height is determined by current measurements of atmospheric pressure and temperature, and not height and temperature. Results. Adaptive and non-adaptive parametric identification algorithms are developed. An improved version of the adaptive algorithm is proposed, which ensures proper tuning of its filtering properties to a specific statistical situation. Algorithms were tested both according to the characteristics of the ISA using computer-generated random interference, and according to the results of real measurements known in the literature obtained using a meteorological rocket. The properties of adaptive and non-adaptive algorithms are compared. It is shown that repeated "running" of the same experimental data through an adaptive algorithm leads to almost the same accuracy of the gradient

estimate, as when using a non-adaptive algorithm. Conclusion. The implementation of the algorithms improves the accuracy of determining the barometric altitude of the aircraft (AC), does not require any modernization of the hardware of the SHS, only some addition of the algorithmic support of the system is necessary. The adaptive algorithm is recommended for use in the case of a significant drift of atmospheric characteristics for the operational adjustment of the system.

Keywords: barometric altitude, air signal system, temperature gradient, identification algorithms, experimental data, filtering properties.

References

1. Panferov V.I., Trenin N.A., Hayutin A.M. Aviatsionnyye pribory i pilotazhno-navigatsionnyye kompleksy: uchebnoye posobiye: v 3 ch. Ch. I. [Aircraft Instruments and Flight-Navigation Complexes: A Manual: at 3 Parts. Part. I.]. Chelyabinsk, Branch of VUZS VVS Air Force Publ., 2018, 145 p.

2. Kozitsin V.K. Sistema vozdushnykh signalov vertoleta na osnove svobodno oriyentirovannogo priyemnika davleniy: avtoref. dis. kand. tekhn. nauk [Helicopter Air Signal System Based on Freely Oriented Pressure Receiver. Abstract of Cand. Sci. Diss.]. Kazan, KSTU Publ., 2006, 23 p.

3. Vasilyuk N.N. Povysheniye tochnosti raketnogo barometricheskogo vysotomera s uchetom aerodinamiki i itogov radiolokatsionnykh nablyudeniy v letnykh eksperimentakh: avtoref. dis. kand. fiz.-mat. nauk [Improving the Accuracy of the Rocket Aarometric Altimeter, Taking into Account Aerodynamics and the Results of Radar Observations in Flight Experiments. Abstract of Cand. Sci. Diss.]. Moscow, Fiztekh-Polygraph Publ., 2004, 23 p.

4. Soldatkin V.V. [Algorithms for the Formation and Processing of Information System Helicopter Air Signals Based on a Fixed Multi-functional Aerometric Receiver]. Fundamental and Applied Problems of Engineering and Technology, 2013, no. 3 (299), pp. 110-120. (in Russ.)

5. Krylov D.L. [Functional Diagram, Information Processing Algorithms and Errors of the Aircraft's Air Signal System with Fixed, Non-Protruding Flow Receivers]. Novyye tekhnologii, materialy i oborudovaniye rossiyskoy aviakosmicheskoy otrasli: sb. dokl. Vserossiyskoy nauch.-prakt. konf. [New Technologies, Materials and Equipment of the Russian Aerospace Industry: Sat. Report All-Russian Scientific and Practical. Conf.]. Kazan, 2016, pp. 570-574. (in Russ.)

6. Soldatkin V.V. [Algorithmic Support of the Airborne Signal System of a Single-Rotor Helicopter with a Fixed Aerometric Receiver at Characteristic Operating Conditions]. Trudy Mezhdunarodnoy nauchno-tekhn. konf. "Perspektivnyye informatsionnyye tekhnologii " [Proceedings of the International Scientific and Technical. Conf. "Perspective Information Technologies"]. Samara, 2018, pp. 881- 885. (in Russ.)

7. Efremova E.S. [Algorithmic Support of the Vortex System of Air Signals Subsonic Aircraft]. Trudy Mezhdunarodnoy nauch.-tekhn. konf. "Perspektivnyye informatsionnyye tekhnologii" [Proceedings of the International Scientific and Technical. Conf. "Perspective Information Technologies"]. Samara, 2018, pp. 812- 816. (in Russ.)

8. Soldatkin V.V. [Automatic Adjustment of the Measuring Channels of the Helicopter Air Signal System]. Bulletin of KSTU, 2004, no. 2, pp. 26-29. (in Russ.)

9. Kaletka J. Evaluation of the Helicopter Low Airspeed System Lassie. Journal of American Helicopter Society, 1983, no. 4, pp. 35-43. DOI: 10.4050/JAHS.28.35

10. Leishman J.G., Naipei B.I. Airdynamic Interactions Bethcen a Rotor and a Fuselage in Forward Flight. J. Amer. Helicopter. Soc, 1990, vol. 35, no. 3, pp. 22-31. DOI: 10.4050/JAHS.35.22

11. Rastrigin L.A. Sovremennyye printsipy upravleniya slozhnymi ob"yektami [Modern Principles of Control of Complex Objects]. Moscow, Sov. Radio Publ., 1980, 300 p.

12. Panferov V.I. [On the Adaptation of the Temperature Control System of the Continuous Thermoaggregates]. News of Higher Education Institutions. Ferrous Metallurgy, 1990, no. 10, p. 72.

13. Sanchez-Gasca J.J., Clark K., Miller N.W., Okamoto H., Kurita A., Chow J.H. Identifying Linear Models from Time Domain Simulations. IEEE Computer Applications in Power, Apr. 1997, vol. 10, no. 2, pp. 26-30. DOI: 10.1109/67.582448

14. Hemerly E.E. PC-Based Packages for Identification, Optimization and Adaptive Control. IEEE Control Systems Magazine, Feb. 1991, vol. 11, iss. 2, pp. 37-43. DOI: 10.1109/37.67674

15. Alekseev P.P., Besyadovsky E.A., Golyshev G.I., Izakov M.N., Kasatkin A.M., Kokin G.A., Livshits N.S., Masanova N.D., Shvidkovsky E.G. [Rocket Atmospheric Research]. Meteorology and Hydrology, 1957, no. 8, pp. 3-13. (in Russ.)

16. Salyga V.I., Karabutov N.N. Identifikatsiya i upravleniye protsessami v chernoy metallurgii [Identification and Control of Processes in the Steel Industry]. Moscow, Metallurgy Publ., 1986, 192 p.

Received 4 November 2019

ОБРАЗЕЦ ЦИТИРОВАНИЯ

К настройке алгоритма вычисления барометрической высоты летательного аппарата системой воздушных сигналов / В.И. Панферов, Н.А. Тренин, С.В. Панферов и др. // Вестник ЮУрГУ. Серия «Компьютерные технологии, управление, радиоэлектроника». - 2020. - Т. 20, № 1. - С. 60-70. DOI: 10.14529/ctcr200106

FOR CITATION

Panferov V.I., Trenin N.A., Panferov S.V., Hayu-tin A.M., Paukov D.V., Haldin K.S. To the Solution of the Problem on the Determination of Barometric Height of the Aircraft System of Air Signals. Bulletin of the South Ural State University. Ser. Computer Technologies, Automatic Control, Radio Electronics, 2020, vol. 20, no. 1, pp. 60-70. (in Russ.) DOI: 10.14529/ctcr200106