CC BY
50
ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И
УПРАВЛЕНИЕ
УДК621.396
АслановГ.К., КазибековР.Б., Набиулин Р.М., Тетакаев У.Р.
МОДЕЛИРОВАНИЕ ПРОЦЕССА ФОРМИРОВАНИЯ ДИАГРАММЫ НАПРАВЛЕНОСТИ ФАЗОВОГО УКВ АВТОМАТИЧЕСКОГО РАДИОПЕЛЕНГАТОРА В СРЕДЕ РАЗРАБОТКИ ПРИКЛАДНЫХ ПРОГРАММ LabVIEW
Aslanov G.K., KazibekovR.B., NabiulinR.M., Tetakaev U.R.
MODELING OF PATTERN FORMING PROCESS OF AUTOMATIC RADIO DIRECTION FINDER OF PHASE VHF IN THE DEVELOPMENT ENVIRONMENT OF LabVIEW APPLIED PROGRAMS
Разработана модель, демонстрирующая процесс формирования диаграммы направленности антенной системы аэродромного квазидоплеровско-го автоматического радиопеленгатора в среде разработки прикладных программ LabVIEW (Laboratory Virtual Instruments Engineering Workbench) компании National Instruments.
Ключевые слова: моделирование, радиопеленгатор, антенная система, диаграмма направленности.
In the article is developed the model demonstrating the forming process of pattern of antenna system of aerodrome quasidopler automatic radiodirection-finder station in the development environment of LabVIEW applied programs of National Instrument company.
Key words: modeling, radiodirection-finder, antenna system, pattern.
Введение. Разработка сложных радиотехнических систем требует сопоставления различных вариантов реализации для выбора оптимального из них. Эта задача может быть решена путем натурной реализации систем и сопоставления результатов испытаний, что является трудоемким, дорогостоящим и, что самое главное, длительным процессом.
Обойти указанные выше трудности можно моделированием процессов функционирования этих систем.
Постановка задачи. В части проведения испытаний автоматических радиопеленгаторов (АРП), моделирование обладает заведомо более широкими
49
возможностями. Так, замена натурного эксперимента для исследования процессов пеленгования, кроме уменьшения колоссальных затрат времени и средств, позволяет обеспечить повторяемость результатов эксперимента. При натурном эксперименте невозможно повторно обеспечить идентичность условий распространения радиоволн и характеристик подстилающей поверхности в районе размещения АРП, а также повторное нахождение воздушного судна в заданной точке пространства.
Методы исследования. Известно много методов обработки пеленгаци-онной информации. В связи с этим, возникает задача выбора оптимального из них по таким показателям, как быстродействие, точность определения пеленга и др. Эти задачи решаются путем моделирования процессов обработки пе-ленгационной информации.
Одним из факторов, влияющих на работу АРП являются когерентные источники излучения. В квазидоплеровских АРП можно искусственно формировать диаграмму направленности (ДН) в виде узкого луча [1].
Определение пеленга путем вращения диаграммы направленности в пространстве позволяет значительно уменьшить влияние когерентных источников излучения на точность пеленгования.
Как известно, сигнал на выходе радиоприемного устройства квазидопле-ровского АРП описывается выражением [2]:
е (у) = д со8
2к(! -1) 2жЯ
+ -
N Л
■С08Д 81И
2к(к -1)
N
-О
к = 1,Ы
(1)
где: А - амплитуда пеленгуемого сигнала,
к - номер временного интервала, в течение которого коммутируется /-й элемент антенной системы,
N - число элементов антенной решетки (АР),
01 и р1 - соответственно, пеленг и угол места на источник излучения, Л - длина волны пеленгуемого сигнала, Я - радиус АР.
В реальных условиях к основному сигналу на вибраторе могут добавляться отраженные сигналы, которые векторно складываются с ним [1,3]:
е(У) = ех(г) + е2(У) ,
где: е^) и е2(1:) - прямой и отраженный сигналы.
^ (^) = А1 008
2я-(/ -1) 2жЯ
+ -
N Л
-008Д 81П
'2ж(к -1) N
-о
e2 (t) = • A1 cos
2x(i -1) 2nR п . r л--cosp2 sin
N
Л
2ж(к -1) ^
• — <
N
К0 - коэффициент отражения от местного предмета;
в2, р2 - соответственно азимут и угол места отраженного сигнала;
Для корреляционной обработки формируется опорный сигнал вида:
e(m (t) = A cos
2n(i -1) 2nR( л.Г 2ж(к -1) .. ^ —~— л—— cos(A^ • j )sm| —^—~-A0^ m
N к \ N
где: Лви Лр - соответственно шаг сканирования по пеленгу и углу места. т - номер шага сканирования по пеленгу, т = 1, ; Nm - количество шагов сканирования по азимуту в, ^ = 360 / Лв;
j - номер шага сканирования по углу места, j = = 1, N;
Nj - количество шагов сканирования по углу места Д,
Nj = 60 / АД;
При значениях A6m и Afij равных 6i, Д\ и в2, Д2 будут иметь место максимумы ДН.
Для разработки программы, демонстрирующей работоспособность метода формирования диаграммы направленности, воспользуемся средой разработки прикладных программ LabVIEW (Laboratory Virtual Instruments Engineering Workbench) компании National Instruments. В ней используется интуитивно понятный язык графического программирования G.
Результаты эксперимента и их обсуждение. LabVIEW предоставляет усовершенствованные возможности для реализации оригинальных научно -исследовательских проектов за кратчайшие сроки, благодаря тесной интеграции с измерительными комплексами и большому числу встроенных функций для численной обработки и средств графического представления данных.
Кроме того, LabVIEW ускоряет процесс разработки за счет построения и отладки алгоритмов на основе данных реальных измерений, что очень важно для радиопеленгации.
В LabVIEW программа имеет два представления - передняя панель (Front Panel) и блок-диаграмма (Block Diagram) виртуального прибора (ВП, VI - Virtual Instruments). На передней панели размещаются объекты интерфейса пользователя, порты ввода, вывода.
В блок-диаграмме же сосредоточен исходный код программы, написанный на графическом языке программирования G.
На рис. 1. изображена передняя панель ВП.
Рисунок 1 - Передняя панель ВП «Метод формирования диаграммы
направленности»
В левом верхнем углу расположено поле ввода входных параметров для вычисления значений разностей фаз между центральным и кольцевыми вибраторами антенной системы, определяемых формулой:
v i = ( 2nR/A) cosfi cos(2n(i-1) /N - 0),
i = 1, N
(2)
Переключатель «Тум.» (тумблер) позволяет вносить в вычисленные разности фаз дополнительные ошибки, вызываемые наличием помех. На первой диаграмме «Phase on vibrators» показаны графики, отображающие выборки фаз вибраторов антенной системы. График, нарисованный штриховой линией получен по (2).
В эксплуатируемых АРП, фазовый детектор производит однозначное декодирование входного сигнала при изменении разностей фаз в пределах от 00 до 360°. Если девиация фазы входного сигнала выходит за пределы этого диапазона, то на выходе фазового детектора возникает фазовая неоднозначность, описываемая выражением:
V = Vu
Vu
360
■ 360
(3)
где: фд и фи - соответственно фаза сигнала на выходе фазового детектора и истинное значение фазы;
[фи/360.] - обозначают операцию антье (выделение целого числа из вещественного без операции округления).
На второй диаграмме «Phase on vibrators» сплошной линией изображены выборки фаз с вибраторов антенной системы с неоднозначностью - фд.
На второй диаграмме «Diagram», на графике «Etta», приведен результат корреляции опорного сигнала и пеленгуемого сигнала по углу в при известном Д На третьей диаграмме представлен результат корреляции по углу Д при известном в.
Значения найденных углов приведены в удобном для восприятия виде: кругового циферблата «Direction» для в и ползунка «Ugol Prihoda» для Д. На рис. 2. представлена блок - диаграмма основной программы ВП.
Рисунок 2 - Блок - диаграмма основной программы ВП
На блок диаграмме слева расположены порты ввода с названиями, в соответствии с названиями передней панели. Подпрограмма «Formula.vi» выдает массив вычисленных по (2) истинных значений фаз и по (3) значений фаз с неоднозначностью (direct and ambiguity phase).
В цикле «For» (на рисунке 2, помечен цифрой 1) производиться разрешение фазовой неоднозначности, возникающей после добавления фазовых составляющих, вызванных помехами и фазовой неидентичностью вибраторов АР.
Подпрограмма «Radiation pattern.vi» предназначена для формирования диаграммы направленности и нахождения азимута в и угла места ß.
На рис.3 представлена блок диаграмма подпрограммы «Radiation pat-tern.vi»,
Фд = Фи
Фи
360
360
На рис. 3а. расположены порты ввода входных параметров подпрограммы: шаг сканирования по углам в и Д входной массив фаз с неоднозначностью, частота пеленгуемого сигнала, радиус и число вибраторов АР.
Результат сканирования представляет собой 2-мерный массив, где строки - это номера шагов по Д, а столбцы - по в. Соответственно количество итераций главного цикла
б)
Рисунок 3 - Блок - диаграмма подпрограммы « Radiation pattern.vi» ВП
«For» вычисляется делением максимального возможного угла сканирования по ß (в настоящее время это 60°) на величину шага сканирования {2}. Здесь и далее в фигурных скобках приведены номера меток из рисунка 3.
В первом вложенном цикле « For», количество итераций для расчета в вычисляется делением максимального значения азимута в (360°) на величину
шага сканирования. У каждого цикла есть терминал счета количества итераций (символ «i» в синем квадрате) и, чтобы определить значение угла на каком-нибудь шаге сканирования, достаточно перемножить значение соответствующего шага на номер его итерации {3}.
Красным квадратом обведен второй вложенный цикл «For» в котором с помощью двух ВП «Sine Waveform.vi» {5} формируются сигналы eon(t) и e(t). Количество итераций в цикле «For» определяется размерностью входного массива фаз {4}. ВП «Sine Waveform.vi» является компонентом библиотеки «Waveform Generation», позволяющий получить цифровой синусоидальный сигнал LabVIEW формата - Waveform(DBL).
LabVIEW формата - Waveform(DBL) - это специальный кластер, который состоит из 3х типов данных: одномерного массива отсчетов сигнала (1 -d array of (DBL)), временного шага между отсчетами dt (DBL), начального значения времени первого отсчета t0 (Time Stamp). Значение амплитуды синусоидального цифрового сигнала в ВП « Sine Waveform.vi» по умолчанию равен 1.
В цикле «For» (очерчен красным квадратом), из всех входных параметров ВП «Sine Waveform.vi» на каждой итерации изменяются только значения фаз сигналов eon(t) и e(t){6}. В связи с чем, при каждой итерации на выходе ВП «Sine Waveform.vi» получаем цифровой синусоидальный сигнал частоты 5550 Гц, с амплитудой, равной 1 и фазой, соответствующей одному из элементов входного массива фаз, в соответствии с номером итерации.
На рис. 3 а) два рядом встречающихся значения 27750 и 40 - это параметры кластера «sampling info»: 27750 - частота дискретизации; 40 - количество генерируемых отсчетов.
Таким образом, на выходе цикла «For» накапливаются два двумерных массива синусоидальных сигналов e(t) и eon(t) с разными фазами, которые трансформируются в одномерные массивы {7}.
Верхний входной массив фаз поступает из вне в подпрограмму и соответствует фазам пеленгуемого сигнала e(t), а второй массив фаз формируется при помощи подпрограммы «Formula.vi» - соответствует массиву фаз опорного сигнала eon(t).
Сигналы e(t) и eon(t) перемножаются и суммируются, результат возводиться в квадрат. Это число, как результат корреляции фиксированного угла по Р и одного из углов 0 {8} поступает на выход первого вложенного цикла «For». Таким образом, на выходе этого цикла формируется массив значений ДН при сканировании по углу 0.
В полученном массиве (в диаграмме) производится определение максимума ДН по 0, значение которого поступает на выход основного цикла «For», которое в свою очередь, с целью увеличения амплитуды диаграммы направленности, возводится в квадрат.
Максимальные значения диаграммы направленности накапливаются на нижнем выходе основного цикла «Main for loop» и представляют собой диаграмму сканирования по углу Р {9}.
Из массива накопленных значений в для разных значений ß выбираем тот индекс, для которого в является максимальным.
На рис. 2. выходы подпрограммы «Radiation pattern.vi» подключены к соответствующим терминалам вывода найденных значений в и ß Весь код сосредоточен в цикле « While Loop», где условием остановки цикла является нажатие кнопки «stop».
Вывод. Исследования, проведенные на модели показали, что определение пеленга методом формирования диаграммы направленности [1] обладает более высокой помехозащищенностью, чем, например, находящийся в эксплуатации радиопеленгатор АРП-75 (обеспечивает работоспособность при уровне шума ниже 30 % от уровня основного сигнала, а рассматриваемый метод - при уровне ниже 90%).
Разработанная модель использовалась при выполнении опытно -конструкторских работ «Перспектива» в ОАО НПК «Русская радиоэлектроника».
Библиографический список:
1. Асланов Г.К., Магомедов К.Г., Дзюба А.П. Формирование диаграммы направленности АРП с помощью корреляционной обработки // Вестник Дагестанского государственного технического университета. Технические науки. - Махачкала: ДГТУ, 1997. - Вып. №1; С 99-101.
2. Асланов Г. К., Саидов А. С., Тагилаев А. Р. Проектирование фазовых автоматических радиопеленгаторов, М.: Радио и Связь, 1997.
3. Асланов Г.К., Гасанов О. И. Анализ причин возникновения аномальных ошибок в квазидоплеровских автоматических радиопеленгаторах // Научно-технические ведомости СПбГПУ №2 (76) 2009. Информатика. Телекоммуникации. Управление, С 87-93.
УДК 69.059.72
Байрамуков С.Х., Долаева З. Н., Омаров А. О.
О МЕТОДАХ ДИНАМИЧЕСКОГО ПРОГРАММИРОВАНИЯ ПРОЦЕССОВ КОМПЛЕКСНОЙ МОДЕРНИЗАЦИИ ЖИЛИЩНОГО ФОНДА
Bayramukov S.H., Dolaeva Z.N., Omarov A. O.
ABOUT THE METHOD OF DYNAMIC PROGRAMMING PROCESS OF INTEGRATED HOUSING MODERNIZATION
В данной статье проанализирована проблема обновления жилищного фонда. Введено более абстрактное и обширное понимание категории модернизации как комплексного развития жилой недвижимости. Построена ма-
56
