Навигация
Поиск
Разработки
LabVIEW в расчетах радиолиний систем передачи данных

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

В настоящее время все более широкое применение в России получают высокоскоростные беспроводные (радио) сети передачи данных. Зоны, обслуживаемые такими сетями, зачастую выходят за пределы одного населенного пункта, охватывая значительные территории пересеченной местности. При этом расположение пунктов приема - передачи далеко не всегда соответствует условиям, благоприятным с точки зрения распространения радиосигналов между этими пунктами. Кроме того, во многих случаях требуется организовать прием - передачу данных между стационарными (с одной стороны) и передвижными пунктами (с другой), что в свою очередь осложняет организацию обмена данными между ними. [1] В этих условиях для выбора технических решений (обеспечивающих требуемую скорость обмена данными в сети и разумных с точки зрения финансовых затрат на их реализацию) актуальной становится проблема оценки уровня затухания радиосигналов с учетом как реального рельефа местности (с застройкой), характеристик предлагаемого (или имеющегося) оборудования, так и погрешностей и нестабильности позиционирования антенных устройств на пунктах приема - передачи.

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

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

2 - технологических погрешностей, эксплуатационных и временных нестабильностей позиционирования антенных устройств, задаваемых углами отклонения главного лепестка диаграммы направленности от требуемого направления в горизонтальной и вертикальной плоскостях и углом поворота антенного устройства относительно этого направления;

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

Следует отметить, что в распространенных методиках проектирования радиолиний расчет затухания проводится лишь для свободного пространства (в лучшем случае с учетом поправок, определяемых в результате дополнительных измерений на местности) и рекомендуется таким образом размещать антенные устройства на пунктах приема - передачи, чтобы, по -возможности, большая часть первой зоны Френеля оставалась свободной от препятствий. При этом оговаривается, что реальное затухание сигнала, а, следовательно, и дальность радиопередачи (или, соответственно, скорость обмена данными) будет во многом определяться конкретными условиями на местности.

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

2. Предприятие, на котором внедрено решение

Данная программа разработана и используется в Институте проблем управления РАН.

3. Используемое оборудование и программное обеспечение National Instruments

Программа выполнена в среде программирования LabVIEW 7.0

4. Описание решения

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

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

1) Модуль синтеза объемных (3-х мерных) диаграмм направленности используемых антенн по заданным значениям ширины основного лепестка, уровней бокового излучения в горизонтальной и вертикальной плоскостях и обратного излучения, коэффициента перекрытия по частоте. Синтез осуществляется с использованием семейства 2-х параметрических Лу-функций первого рода с параметром (v), равным, соответственно, -1/2 и 1. [6]. Выбор функций, используемых в дальнейших расчетах, осуществляется по заданным значениям коэффициентов усиления антенн. В результате синтеза формируются массивы значений диаграмм направленности в любых (с заданным шагом) направлениях. Эти массивы используются для контрольного визуального отображения диаграмм направленности. Эти же массивы поступают и в управляющую программу, где используются для расчета затухания радиолинии с данными антеннами.

На рис. 1 показана лицевая панель модуля синтеза диаграмм направленности антенн.

2) Модуль формирования матрицы высот рельефа местности. В предлагаемой версии программы матрица высот рельефа формируется путем линейной интерполяции значений высот точек отдельных профилей рельефа, расположенных на выбранных расстояниях от вертикальной осевой плоскости. Сформированная таким образом матрица высот используется для контрольного визуального отображения подстилающей поверхности установленных размеров относительно эллипсоидов первой и любой другой (выбранной произвольно) зон Френеля (отображаемых относительно осевой линии, соединяющей антенные устройства, расположенные в заданных точках, на заданном расстоянии друг от друга и размещенные на выбранных высотах) Эта матрица поступает также в управляющую программу, где используется для расчета затухания радиолинии.

На рис. 2 показана лицевая панель модуля формирования матрицы высот.

3) Модуль управляющей программы.

Этот модуль обеспечивает ввод в программу основных параметров проектируемой радиолинии (таких, как рабочая длина волны, поляризация сигнала и расстояние между пунктами обмена данными), ввод возможных отклонений в позиционировании антенных устройств и модуля коэффициента отражения в «блестящей» точке. Этот модуль обеспечивает ввод в программу предлагаемого размера зоны обзора и возможность выбора алгоритма расчета, оптимального для данных условий, определяющих размеры как зоны обзора, так и области, свободной от препятствий. Управляющая программа осуществляет расчет затухания с использованием всех поступающих данных. Расчет осуществляется для сигнала с вертикальной, или с горизонтальной поляризацией. Результаты расчета отображаются на цифровом и графическом индикаторах. На цифровом индикаторе отображается значения затухания радиолинии; на графическом индикаторе отображается годограф уровней сигналов при сканировании по эллипсоидам зон Френеля в выбранном объеме. На рис. 3 показана лицевая панель модуля управляющей программы.

5. Преимущества технологий National Instruments

Использование пакета LabVIEW позволило легко и компактно написать программы всех модулей, объединить их в единый комплекс и быстро отладить программу.

6., Перспективы внедрения и развития решения

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

Литература

[1] К. Ли Уильям. Техника подвижных систем связи. Перевод с английского. Москва «Радио и связь», 1985. McGraw - Hill Company New York, 1982.

[2] P. Фейнман, А. Хибс. Квантовая механика и интегралы по траекториям. Перевод с английского. Издательство «Мир» Москва 1968. McGraw - Hill Company New York, 1965.

[3] P. Фейнман. КЭД странная теория света и вещества. Перевод с английского. Москва «Наука» Главная редакция физико - математической литературы, 1988. Princeton, New Jersey: Princeton University Press, 1985.

[4] Ю.С. Легович, Ю.В. Максимов, Д.Ю. Максимов. Проектирование беспроводных сетей WiFi повышенной дальности. Материалы международной конференции CAD/CAM/PDM -2006, стр.64

[5] Ю.С. Легович, Ю.В. Максимов, Д.Ю. Максимов. Метод расчета затухания сигнала в беспроводных системах передачи данных с учетом рельефа местности и параметров антенн, Автоматика и телемеханика, 2008.

[6] Справочник по радиолокации. Перевод с английского в четырех томах. Том 2 «Радиолокационные антенные устройства». Москва «Советское Радио», 1977. McGraw-Hill Book Company, 1970.

Продолжение...

Контакты
an image
НПП Центральная лаборатория автоматизации измерений
111250 Москва, Энергетическая улица, д.7, офис 311
(495) 134-03-49
E-mail: Данный адрес e-mail защищен от спам-ботов, Вам необходимо включить Javascript для его просмотра.

Стенд для исследований рабочих процессов ДВС в динамических режимах

1. Введение

В настоящее время двигатели внутреннего сгорания используются во всех областях народного хозяйства: промышленность и сельское хозяйство, гражданское и дорожное строительство, транспорт, энергетика, нефтяная промышленность и т.д. Двигатели используемых в этих областях мобильных и стационарных агрегатов работают при резко изменяющихся нагрузках [1]. Это происходит при трогании с места автомобилей и тракторов, переключении передач, преодолении препятствий, выполнении дорожно-строительных и сельскохозяйственных работ [1]. По литературным данным переменные режимы для автомобилей составляют в городе 97%, на грунтовых дорогах 92%, на загородных магистралях 34% процентов всего времени работы [1 имеющиеся в ней ссылки]. Переменные режимы тракторных двигателей возникают при транспортных работах, которые составляют 40-45% от общего времени работы, так и при выполнении основных сельскохозяйственных операций. Если при транспортных работах двигатели загружены резкопеременными нагрузками, обусловленными в основном изменением профиля дорог, то для сельскохозяйственных операций характерен режим при котором нагрузка колеблется относительно какого-то среднего значения.

Исследования показывают, что учет динамических характеристик двигателя позволяет повысить эксплуатационную мощность, экономичность, надежность, долговечность и понизить токсичность и износ двигателя внутреннего сгорания [2]. Потенциальные возможности использования результатов исследования ДВС при неустановившихся режимах составляет повышение экономичности 15-20%, производительность машинно-тракторных агрегатов на 15-20%. Из всего вышесказанного очевидна актуальность разработки стендов для исследования двигателей внутреннего сгорания при неустановившейся нагрузке. Вариант такого стенда рассматривается в данной работе.

2. Схема стенда для исследований рабочих процессов двигателя внутреннего сгорания в динамических режимах

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

Схема стенда для исследований рабочих процессов двигателя внутреннего сгорания в динамических режимах приведена на рис. 1. Стенд работает следующим образом.

После запуска двигателя внутреннего сгорания 4 электротормоз 2 создает постоянный тормозящий момент, воздействующий на коленчатый вал двигателя внутреннего сгорания 4. После достижения двигателем внутреннего сгорания 4 установившегося режима имитатор типовых нагрузок 1 по сигналу автоматизированного блока управления 14 вырабатывает сигнал, приводящий к изменению момента сопротивления электротормоза 2. Изменение тормозящего момента приводит к появлению переходных процессов: изменяется крутящий момент на валу двигателя 4, изменяется частота вращения коленчатого вала двигателя 4, изменяется расход топлива и воздуха.

Схема стенда для исследований рабочих процессов двигателя внутреннего сгорания в динамических режимах

Рис. 1. Схема стенда для исследований рабочих процессов двигателя внутреннего сгорания в динамических режимах

Эти переходные процессы посредством датчиков преобразуются в электрические сигналы и поступают в блок сбора данных 10, который по команде автоматизированного блока управления 14 выдает на вход «сигнала отклика» блока идентификации параметров дифференциального уравнения один из четырех сигналов датчиков: датчика 5 частоты вращения коленчатого вала двигателя, датчика 6 момента сопротивления двигателя, датчика 7 расхода воздуха, датчика 8 расхода топлива, а на вход «сигнала» воздействия - сигнал с датчика 9 тока возбуждения электротормоза, при этом соответствующие сигналы преобразуются в цифровую форму. Сигналы, поступившие в блок идентификации параметров дифференциального уравнения 11, используются для определения параметров динамических моделей рабочих процессов двигателя внутреннего сгорания (коэффициентов дифференциального уравнения), при этом используются дифференциальных уравнений следующего вида: для регуляторной ветви регуляторной характеристики

где a1, a2, a3 - искомые коэффициенты дифференциального уравнения; ∆Ai - изменение наблюдаемого показателя двигателя внутреннего сгорания крутящего момента на валу двигателя, частоты вращения коленчатого вала двигателя, расхода топлива и воздуха при изменении нагрузки; AM - изменение момента сопротивления электротормоза; К - коэффициент пропорциональности.

Схема блока идентификации параметров дифференциального уравнения

Рис. 2. Схема блока идентификации параметров дифференциального уравнения

Схема блока идентификации параметров дифференциального уравнения 11 приведена на рис. 2, данный блок работает следующим образом, сигнал с входа сигнала воздействия: изменение момента сопротивление электротормоза 2, поступает на блок математической модели, выходной сигнал блока математической модели поступает на вход блока определения ошибки моделирования, где он сравнивается с сигналом, поступившим на вход сигнала отклика, который соответствует реакции двигателя на воздействие - изменение либо момента сопротивления двигателя, либо частоты вращения коленчатого вала, либо расхода топлива, либо расхода воздуха, на основании этого блок определения ошибки моделирования вырабатывает сигнал, соответствующий рассогласованию математической модели и исследуемого объекта, данный сигнал поступает на вход блока корректировки математической модели, изменяющий параметры блока математической модели с целью уменьшить ошибку моделирования, коэффициенты, соответствующие минимуму ошибки моделирования, принимаются в качестве искомых коэффициентов. Полученные коэффициенты дифференциального уравнения поступают в блок расшифровки результатов идентификации 12, который определяет порядок используемого для идентификации дифференциального уравнения и представляет коэффициенты дифференциальных уравнений в форме удобной для дальнейшего использования. После этого динамические характеристики либо выдаются на выход потребителю, либо запоминаются в запоминающем устройстве 13. Работа стенда автоматизирована и управляется автоматизированным устройством управления 14.

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

Лицевая панель программы

Рис. 3. Лицевая панель программы

3. Идентификация дифференциального уравнения второго порядка

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

Требуется, используя информацию о воздействие x(t) на исследуемый объект и об его отклике y(t), определить параметры 7V, Т2, к математической модели М описываемой дифференциальным уравнением третьего порядка (для регуляторной ветви регуляторной характеристики, мы будем рассматривать только его как наиболее общий):

Пронормируем дифференциальное уравнение

или

где

Рассмотрим способ идентификации по замкнутой схеме (рис.2). В этом случае задача нахождения коэффициентов дифференциального уравнения (5) сводится к задаче оптимального управления коэффициентами a2, a1, a0.

Пусть в результате измерения получена реализация выходного сигнала системы y(t)=f(t, a2, a1, a0) в ответ на входное воздействие x(t) (здесь a2, a1, a0 -неизвестные параметры системы). Зададимся параметрами модели a2м, a1м, a0м. Тогда входному сигналу x(t) будет соответствовать выходной сигнал модели yM(t)=fM(t, a2м, a1м, a0м). Следовательно, ошибку идентификации е можно определить согласно следующему выражению.

Очевидно, что ошибка ε будет минимальна при совпадении искомых коэффициентов ε с соответствующими параметрами исследуемого процесса a2м →a2 , a1м →a1, a0м →a0 => ε 20. Следовательно, задачу определения коэффициентов дифференциального уравнения, можно свести к задаче оптимального управления коэффициентами a2м, a1м, a0м. Оптимальное управление состоит в следующем: путем регулирования коэффициентами a2м, a1м, a0м найти такие коэффициенты дифференциального уравнения при которых ошибка ε 2=f(a2м, a1м, a0м,a2, a1, a0) будет минимальной.

Алгоритм идентификации (определения параметров) приведен на рис.4.

Входное воздействие поступает одновременно на вход исследуемого объекта и его модели. На основании их выходных сигналов y(t), yM(t) определяется ошибка моделирования ε. Ошибка ε используется для коррекции (регулировки) параметров модели a1м, a0м, с целью ее уменьшения.

В качестве искомых параметров модели примем параметры a1м, a0м для которых ошибка ε минимальна.

Таким образом, процесс определения параметров модели сводится к двум операциям:

• Моделирование;

• Определение ошибки моделирования и коррекция (регулировка) параметров модели.

Остановимся на этих этапах подробнее.

Алгоритм определения параметров модели

Рис.4. Алгоритм определения параметров модели

Моделирование (Численное решение дифференциального уравнения) Как отмечалось выше нам задано дифференциальное уравнение третьего порядка:

где x(t) - функция, определяющая воздействие на вход системы; y(t) - функция

(выходной отклик системы);

a2, a1, a0, К- постоянные коэффициенты.

При моделировании нам задаются следующие величины

3. начальные условия

1. воздействие x(t);

2. параметры a2, a1, a0 , К; ;

Требуется найти выходной сигнал y(t) системы при заданных выше условиях Разрешим уравнение (7) относительно старшей производной:

Такое дифференциальное уравнение (8) может быть смоделировано согласие следующей схеме рис. 5.

Структурная схема численного решения дифференциального уравнения

Рис. 5. Структурная схема численного решения дифференциального уравнения.

Входное воздействие x(t) умножается на весовой коэффициент К. После этого сигнал поступает на вход сумматора. На два других входа сумматора, подаются сигналы с выхода первого интегратора d2y/dt2, с выхода второго интегратора dy/dt выходной сигнал y(t) с соответствующими весовыми коэффициентами -a2, -a1, -a0. На выходе сумматора согласно выражению (8) получим третью производную выходного сигнала системы. Дважды интегрируя этот сигнал, получаем выходной сигнал системы - решение дифференциального уравнения.

Подпрограмма LabVIEW, выполняющая численное решение дифференциального уравнения согласно приведенному выше алгоритму изображена на рис.6.

Минимизация функционала (коррекция коэффициентов дифференциального уравнения)

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

Ошибку моделирования ε можно представить, как функцию нескольких переменных

ε 2=f(a2м, a1м, a0м,a2, a1, a0)

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

ε 2=f(a2м, a1м, a0м) .

Подпрограмма численного решения дифференциального уравнения в LabVIEW

Рис. 6. Подпрограмма численного решения дифференциального уравнения в LabVIEW.

Следовательно, можно использовать широко известные методы поиска минимума функции нескольких переменных. В качестве такого метода выберем «метод покоординатного спуска» (метод Гаусса - Зейделя).

Цикл, осуществляющий минимизацию функционала согласно алгоритму Гаусса - Зейделя в LabVIEW

Рис. 7. Цикл, осуществляющий минимизацию функционала согласно алгоритму Гаусса - Зейделя в LabVIEW

Суть этого метода заключается в следующем:

1. Выбираем некоторые начальные коэффициенты a20, a10, a00

2. После этого мы фиксируем два коэффициента, a1, a0 и производим поиск минимума ошибки ε при изменении коэффициента a2. Коэффициент a2 соответствующий минимуму принимается, за новое приближенное значение коэффициента a2.

3. Далее повторяется процедура аналогичная описанной в п.2. Но при этом фиксируется параметры a2 и a0, и ищется минимум функции при изменении a1.

4. Далее повторяется процедура аналогичная описанной в п.2. Но при этом фиксируется параметры a2 и a1, и ищется минимум функции при изменении a0.

5. Повторяется п. 2

6. Процесс останавливается если при выполнении п.2, п.З и п.4 коэффициенты a2, a1, a0 изменились не более чем на определенную величину ∆. И эти коэффициенты принимаются в качестве искомых параметров модели

Цикл, осуществляющий минимизацию функционала согласно алгоритму Гаусса - Зейделя в LabVIEW, изображен на рис. 7.

Все описанные алгоритмы были реализованы в программной среде LabVIEW, что позволило производить автоматизированное определение искомых параметров.

4. Выводы

Описанный в данной работе стенд может быть использован при исследовании рабочих процессов двигателей внутреннего сгорания в динамических режимах (в условиях эксплуатации). На предлагаемом стенде обеспечивается возможность определения динамических характеристик рабочих процессов двигателя в виде дифференциальных уравнений, описывающих реакции изменения частоты вращения коленчатого вала двигателя, момента сопротивления двигателя, расхода воздуха и расхода топлива при изменении момента сопротивления на коленчатом валу двигателя, при этом процесс измерения и обработки информации автоматизирован. Программное обеспечение реализовано в программной среде LabVIEW. Выбор данной программной среды для реализации приведенного выше алгоритма обусловлен следующими важными факторами:

• На языке графического программирования LabVIEW, именуемого "G", программы создаются в виде блок-схемы, которые являются естественной формой проектирования.

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

• LabVIEW позволяет создавать удобный пользовательский интерфейс.

• Наглядность обработки сигналов.

Литература

1. Юлдашев А.К., Хайрутдинов И.Н. Стенды для исследования двигателей внутреннего сгорания при неустановившихся нагрузках. - Казань: Изд-во «Фэн», 2002.

2. Использование динамических характеристик рабочих процессов тепловых двигателей для проектирования, эксплуатации, диагностики и ремонта: Материлы международного научно-технического семинара, Казань, 22-25 сентября 2003 года.

3. Евдокимов Ю.К., Линдваль В.Р., Щербаков ПИ. LabVIEW для радиоинженеров: от виртуальной модели до реального прибора. Практическое руководство для работы в программной среде LabVIEW. - М: ДМК Пресс, 2007. - 400 с.

4. Патент на полезную модель (РФ) № 66526 Стенд для исследования рабочих процессов двигателя внутреннего сгорания в динамических режимах/ Юлдашев А.К., Евдокимов Ю.К., Синицкий С.А., Денисов ЕС. и д.р. 10.9.2007 Бюл. №25

..............................................................................................................................