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

Размещено на http://www.allbest.ru//

28

ІSSN 0485-8972 Радиотехника. 2014. Вып. 177

Размещено на http://www.allbest.ru//

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

А.А. Костыря, канд. техн. наук

1. Введение и постановка задачи

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

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

В данной статье приводится описание и анализ модели системы, позволяющей производить локальное позиционирование движущихся объектов с точностью до 1 см [2, 3]. Исследуемая система использует ультразвук (УЗ) вместо электромагнитных волн, что и обеспечивает высокую точность определения местоположения объектов. Она состоит из УЗ меток - устройств, устанавливаемых на объектах позиционирования и имеющих в своем составе УЗ приемники, а также двух УЗ излучателей, устанавливаемых неподвижно в области позиционирования. Синхронизацией работы УЗ излучателей и УЗ меток управляет базовая станция, подключенная к персональному компьютеру. Ультразвуковая система локального позиционирования (УСЛП) использует ультразвуковой активный дальномерный метод определения координат. Обмен данными и синхронизация между базовой станцией и УЗ метками осуществляется по радиоканалу.

Целью проводимого математического моделирования работы УСЛП является анализ погрешностей измерения координат предлагаемой системы [4].

2. Создание математической модели

Для моделирования УСЛП используется система автоматизированного проектирования Mathcad 14.0 [5]. Исследуемое устройство состоит из следующих функциональных модулей: УЗ излучатель, звуковой канал, УЗ приемник, предварительный усилитель, детектор сигнала и пороговое устройство, определяющее временное положение УЗ сигнала [4]. В модели этим модулям соответствуют следующие блоки (рис. 1): генератор зондирующего сигнала (УЗ

излучатель), генератор шума (звуковой канал и внутренние шумы устройства), сумматор (звуковой канал, УЗ приемник), узкополосный фильтр (УЗ приемник), ограничитель амплитуды (предварительный усилитель УЗ метки), фиксатор уровня (аналогичный модуль) и пороговое устройство (компаратор УЗ метки).

Размещено на http://www.allbest.ru//

28

ІSSN 0485-8972 Радиотехника. 2014. Вып. 177

Размещено на http://www.allbest.ru//

Рис. 1. Схема математической модели УСЛП

2.1 Формирование сигнала и шума

В УСЛП применяется УЗ сигнал с несущей частотой 40 кГц и длительностью 1 мс.

УЗ зондирующий сигнал формируется в соответствии с выражением

(1)

где , - период сигнала; - амплитуда сигнала, - номер выборки; - количество выборок.

Таким образом, генерируется синусоидальный сигнал на интервале времени от до , в остальные моменты времени сигнал равен 0.

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

, (2)

где - СКО шума; - индекс реализации шума; - количество реализаций шума.

С использованием (1) и (2) формируется смесь сигнала и шума:

.

На рис. 2 представлены: графики зондирующего сигнала (рис. 2, а) и смеси сигнал-шум (рис. 2, б).



а

б

Рис. 2. Графики: а - сигнала; б - и смеси сигнал-шум

2.2 Моделирование полосового фильтра

В УСЛП применяются узкополосные УЗ приемники, имеющие резонанс на частоте 40 кГц. Для выполнения моделирования полосового фильтра, имеющего характеристику, аналогичную УЗ приемнику, сигналы и шум преобразовываются в частотную область с использованием функции Mathcad :

;

;

,

где - спектр зондирующего сигнала, - спектр шума, - спектр смеси сигнал + шум.

Функция выполняет дискретное преобразование Фурье двухмерного вектора :

,

где - число элементов в векторе , - мнимая единица.

Далее определяется частотная характеристика полосового фильтра:

,

где - частота несущей сигнала, - ширина полосы пропускания фильтра (относительно несущей частоты), - функция экспоненты Mathcad.

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

;

,

где - спектр принятого сигнала на выходе фильтра; - спектр шума, - спектр смеси сигнал + шум.

На рис. 3 представлены спектры сигнала и шума на входе фильтра (рис. 3, а) и соответствующие им спектры на выходе фильтра (рис. 3, б).

При помощи обратного преобразования Фурье определяются сигналы и шум на выходе фильтра во временной области:

,

для чего используется встроенная процедура:

;

;

,

где - принятый сигнал на выходе фильтра, - шум на выходе фильтра, - смесь сигнал + шум на выходе фильтра.

На рис. 4 представлены графики принятого сигнала (рис. 4, а), шума (рис. 4, б) и смеси сигнал + шум (рис. 4, в) на выходе фильтра.

2.3 Ограничение и фиксация сигналов

Вследствие того, что усилитель УЗ метки имеет ограниченный динамический диапазон, ограничение амплитуды сигнала устанавливается на уровне 0,5 В:

а

б

Рис. 3. Амплитудные спектры

а

б

в

Рис. 4. Сигналы на выходе фильтра

;

;

;

,

где - принятый сигнал после ограничения, - смесь сигнал + шум после ограничения, - уровень ограничения.

Функция в Mathcad имеет следующий синтаксис:

.

Если условие выполняется, принимает значение , иначе .

Для того чтобы поднять крутизну фронта принимаемого сигнала и увеличить размах последнего, применяется фиксатор уровня [6]. Его схема приведена на рис. 5.

Рис. 5. Принципиальная схемафиксатора уровня

Фиксатор уровня смещает постоянную составляющую сигнала таким образом, что нижние пики импульсов на его выходе находятся на нулевом уровне (см. рис. 6, б).

Для того чтобы промоделировать работу фиксатора, определим сигнал на емкости C1, когда на входе присутствует принятый сигнал:

,

где - напряжение разряда конденсатора в течение одной выборки.

Сигнал на выходе фиксатора находится как разница между сигналом на его входе и напряжением на C1:

.

Аналогичным образом моделируется шум на выходе фиксатора уровня:

;

.

Смесь сигнал + шум после фиксации:

;

.

Графики смеси сигнал + шум до и после фиксации представлены на рис. 6, а и б соответственно). На рис. 6, а пунктирной линией показан сигнал на емкости фиксатора.

а

б

Рис. 6. Смесь сигнал + шум до и после фиксатора

2.4 Оценка временного положения сигнала

Для устранения влияния многолучевости на работу УСЛП временное положение принятого сигнала определяется по его переднему фронту.

Так как размах полезного сигнала на выходе фиксатора уровня составляет 1 В, порог обнаружения фронта сигнала устанавливаеется на уровне 0.5 В: . Временное положение фронта сигнала оценивается следующим образом. На первом этапе создается массив с индексами выборок, амплитуда которых превышает пороговое значение :

.

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

.

Аналогично оценивается временное положение смеси сигнал + шум:

;

.

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

,

где - момент излучения зондирующего сигнала.

Рис. 7. Истинное и оцененное временное положение сигналов

На рис. 8 представлены значения оценки временного положения принятого сигнала при различных шумовых реализациях. На этом же графике сплошной линией отмечена оценка временного положения зондирующего сигнала.

Рис. 8. Значения оценки временного положения принятого сигнала при наличии шумов

2.5 Погрешности оценки временного положения сигнала

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

.

Систематическая погрешность находится как разность между оценкой временного положения смеси сигнал + шум и истинным положением сигнала:

.

Численно абсолютная и систематическая погрешность находится как среднее арифметическое (функция ) векторов и соответственно:

;

.

Находим СКО абсолютной и систематической погрешностей, используя функцию :

;

.

В результате проведенного моделирования в пакете Mathcad полученные значения погрешностей составляют:

- абсолютная погрешность оценки временного положения сигнала мкс;

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

- СКО абсолютной погрешности мкс;

- СКО систематической погрешности мкс.

Выводы

Математическая модель УСЛП, созданная в пакете Mathcad, позволяет определить основные характеристики системы, в частности погрешности измерения временного положения сигнала, а также оперативно изменять параметры схемы УСЛП и отслеживать влияние данных изменений на конечный результат работы системы.

Результатом работы математической модели являются погрешности определения временного положения принятого сигнала. С учетом усредненного значения скорости звука в воздухе (331 м/с) значения погрешностей измерения дальности УСЛП в данном случае составляют:

абсолютная погрешность измерения дальности см;

систематическая погрешность измерения дальности  см;

СКО абсолютной погрешности измерения дальности см;

СКО систематической погрешности измерения дальности см.

Очевидно, что наибольшее влияние на точность работы УСЛП оказывает систематическая погрешность измерения дальности, но учитывая постоянство данной погрешности, ее можно устранить программными методами.

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

Список литературы

позиционирование сигнал фильтр

1. Новиков А.А., Амелина М.А. Конспект лекций по курсу Математическое моделирование в электронике - Смоленск : Моск. энерг. ин-т (Технический университет), 2006.

2. Пат. 71860 Україна, МПК (2012.01) G01S 11/00. Спосіб ультразвукового локального позиціонування рухомих об'єктів / Лівнов В.О., Коваль Ю.О., заявник та власник Харк. нац. ун-т радіоелектр. - № u2012 01318, подано 08.02.2012, опубл. 25.05.2012 бюл. №14.

3. Современные проблемы радиотехники и телекоммуникаций «РТ-2011»: Материалы 7-й междунар. молодежной науч.-техн. конф., 11 - 15 апреля 2011 г. - Севастополь : Изд-во СевНТУ, 2011. - 479 с.

4. Ливнов В. А. Анализ погрешностей и неопределенностей системы позиционирования с ультразвуковым измерительным и электромагнитным информационным каналами / В. А. Ливнов, Ю. А. Коваль // Радиотехника. - 2013. - Вып. 175. - С.96 - 100

5. Product & Service Advantage - Mathcad [Электронный ресурс] / PTC Inc - Режим доступа: \WWW/ URL: http://blogs.ptc.com/product/mathcad/ - 11.12.13 г. -- Загл. с экрана.

6. Ицхоки Я. С., Овчинников Н. И. Импульсные и цифровые устройства. - М. : Сов. радио, 1972. - С.243.

Размещено на Allbest.ru

...