Modeling of radio channels with Rayleigh fading and non-isotropic spectrum of Doppler spread

全文:

Введение Разработка и оптимизация алгоритмов современных и перспективных устройств передачи данных, ввиду их сложности, крайне затруднительна и малоэффективна без использования компьютерного моделирования широкого диапазона условий работы радиоканалов. Средства моделирования - имитаторы радиоканалов обеспечивают повторяемость и возможность управления воспроизводимыми условиями и характеристиками радиоканалов, отражающих существенные особенности поведения тестируемых устройств в реальных условиях. Кроме того, имитируя процессы прохождения сигналов по радиоканалам модель должна быть способна генерировать процессы со статистическими характеристиками аккуратно согласованными с набором некоторых эталонных характеристик: реализовывать референтную модель канала. Тем самым открывается возможность сопоставлять и исследовать характеристики устойчивости и надёжности передачи данных, которые обеспечивают разрабатываемые алгоритмы в сравнении с уже известными. Применение для целей моделирования средств цифровой обработки сигналов (ЦОС) в полной мере обеспечивает выполнение вышеперечисленных требований. Методы и алгоритмы передачи сигналов, успешно используемые при сравнительно низких скоростях движения мобильных объектов, нуждаются в модификации для работы с высокоскоростными объектами. Для разработки и тестирования новых средств обеспечения передачи информации в условиях высокой скорости движения необходимо иметь модели канала, учитывающие проявления разного вида замираний, возрастающего влияния эффекта допплеровского рассеяния [1; 8]. Процесс релеевских замираний с допплеровским рассеянием адекватно моделируется ARMA (англ. Autoregressive Moving-average Model - модель авторегрессии скользящего среднего) моделями гауссовских процессов, лежащих в основе представления релеевских замираний [2]. В данной статье рассматривается методика построения ARMA модели, имитирующей условия работы и эталонные характеристики каналов с релеевскими замираниями и допплеровским рассеянием в широком диапазоне значений параметров этих каналов. Основные эталонные характеристики референтной модели В теории и приложениях [2-3; 7] рассматриваются два типа допплеровского рассеяния: - изотропного рассеяния, предполагающего равномерную плотность распределения вероятности значений углов прихода сигнала, принимаемого на ненаправленную антенну движущегося объекта; - неизотропного рассеяния, проявляющегося при приеме на антенну с направленными свойствами сигнала с неравномерным распределением значений углов его прихода на движущийся объект. Референтная модель, имеющая аналитическое описание в замкнутом виде своих эталонных спектрально-корреляционных характеристик, описывающая оба типа рассеяния базируется на плотности распределения фон Мизеса-Тихонова углов прихода [4]: , (1) где - модифицированная функция Бесселя 0-го порядка; - среднее значение направления углов прихода; - параметр концентрации, задающий ширину рассеяния . Распределение (1) называют так же «круговым нормальным» распределением. Оно описывает непрерывное распределение углов прихода на окружности с периодом . Параметр соответствует математическому ожиданию, а параметр - величине, обратной дисперсии в обычном нормальном распределении. При распределение (1) превращается в равномерное распределение, порождающее изотропное допплеровское рассеяние: , (2) при , , детерминированное значение угла прихода. Все другие значения параметра порождают неизотропное допплеровское рассеивание. Графическое изображение (1) дано на рис. 1. Рис. 1. Распределение фон Мизеса-Тихонова Релеевские замирания это общепринятое описания для огибающей процесса: , (3) где - синфазная компонента, а - квадратурная компоненты, являющиеся гауссовскими процессами с нулевым математическим ожиданием. Спектральная плотность мощности (СПМ) огибающей процесса (3) соответствующая распределению (1), как установлено в [4] имеет вид: (4) где - максимальная частота допплеровского рассеивания. При СПМ (4) соответствует закону Джекса, то есть изотропной модели Кларка [2; 9]. В общем случае (4) несимметрично относительно , что означает проявление в корреляционной функции взаимной корреляции между квадратурными компонентами процесса в (1). Реализация модели Методами ЦОС этот ARMA процесс может быть реализован в виде рекурсивного цифрового фильтра (ЦФ), на входе которого действует - комплекснозначный белый гауссовский шум, обладающей равномерной СПМ . Главная особенность такого ЦФ - его коэффициенты должны быть комплексными числами. Реализация РЦФ с комплексными коэффициентами достаточно высокого порядка M удобно выполнять в виде каскадного соединения звеньев первого порядка, дискретизированного с частотой : , (5) где . При этом каждое звено реализуется в виде следующей структуры: , (6) где переменные с индексом r - действительные числа, с индексом i - мнимые числа. Видно, что для реализации одного отсчета требуется выполнить восемь умножений действительных чисел. Поскольку при этом выражение (4) в рамках ARMA моделей может быть получено только аппроксимацией дробно-рациональной функцией где - частотная характеристика ЦФ моделирующая (4), то синтез ARMA РЦФ можно осуществить только численными итерационными методами, в которых важную роль играет начальное приближение для параметров функции (5). Если начальное приближение выбрано «хорошим» - то есть достаточно близким к искомому конечному результату, то удаётся избежать многих проблем, присущих итерационным численным алгоритмам [6]. Хотя объекты аппроксимации это квадратичные, энергетические функции, критерий оптимальности, заключающийся в точности воспроизведения референтной модели, более логично полагать чебышевским, минимизирующим: , (7) где ; , (8) где A - масштабирующий множитель. Обозначив отклонение передаточной функции числителя в (8) как , а отклонение знаменателя как , получим, что изменение в (8) имеет вид: , (9) где и - соответственно, выражения для числителя и знаменателя в (8). В линейном приближении можно показать, что отклонения, вызванные изменением модуля и фазы коэффициентов в (9), разделяются и могут быть рассмотрены по отдельности, что значительно упростит процедуру синтеза ЦФ аппроксимирующего СПМ вида (4). Методика синтеза СФФ При создании модели средствами ЦОС для сокращения сложности реализации за счет уменьшения числа вычислительных операций применяют многоскоростную цифровую фильтрацию [5]. Требуемая СПМ получается на выходе каскадного соединения фильтра формирующего СПМ (СФФ), которая аппроксимирует (4) на частоте дискретизации , и последующего интерполирующего фильтра с коэффициентом повышения частоты дискретизации . Наиболее сложная проблема, возникающая при построении такого вида модели - это решение задачи синтеза СФФ: рекурсивного ЦФ с комплексно-значными коэффициентами. Как уже отмечалось выше особое внимание должно быть уделено определению начального приближения. Рассмотрим этот этап решения задачи синтеза СФФ подробнее. Первый шаг в процедуре синтеза состоит в оценке порядка начального приближения ЧХ в (5). Эту роль может выполнить оценка порядка эллиптического ЦФ, практически оптимального ЦФ в классе чебышевской аппроксимации модуля ЧХ приведённая в [5]: , (10) где - допустимая величина сигнала в полосе затухания; - допустимая величина пульсаций сигнала в полосе пропускания; - ширина нормированной переходной полосы. Второй шаг это определение начальных значений комплексных коэффициентов в (5). Для чего может быть использован алгоритм пошагового определения значений искомых коэффициентов из [2]. В состав квадратичной целевой функции алгоритма введена «барьерная» функция для соблюдения условий устойчивости синтезируемого СФФ. Рис. 2. Диаграмма нулей и полюсов передаточной функции РЦФ 12-порядка (пример из [2]) Далее полученное решение модифицируется: анализ структуры выражения (4) позволяет выделить знаменатель, который фактически соответствует СПМ изотропной модели и определяет симметричные граничные частоты в спектре . Такая симметрия может быть мотивом замены одиночных комплексных нулей и полюсов в окрестностях частот на комплексно-сопряженные пары этих точек, координаты которых могут быть выбраны усреднением координат, соответствующих комплексных нулей-полюсов решения полученного на втором шаге, описываемой процедуры, ближайших к точкам на окружности частот , (см. рис. 2) [2]. Произведенная замена имеет одним своим следствием сокращение числа выполняемых операций умножения из расчета два действительных умножения на один комплексный коэффициент. Другое следствие - появление отклонений в СПМ на выходе СФФ. На завершающем этапе синтеза СФФ выполняется поиск решения, в зависимости от условий синтеза, приближающегося или полностью удовлетворяющего чебышевскому критерию аппроксимации (условию альтернанса) в (7) [6]. Рассмотрим в качестве иллюстрации использования предлагаемой методики синтеза СФФ, порядок которого, в соответствии с формулой (10) равен 12 [2]. Для этого примера параметры СПМ из (4) есть , , . Вид СПМ и ее аппроксимация приведены на рис. 3. Рис. 3. Вид СПМ для фильтра 12-го порядка (пример из [2]) Характеристики приближения, полученные на втором шаге приведены на графиках рис. 4. Рис. 4. Вид теоретической и расчетной СПМ, и функция ошибок, соответствующих примеру [2]. Результаты последующей модификации, включающей создание вместо двух звеньев с комплексными коэффициентами вида (5)-(6) одного биквадратного звена РЦФ с действительными коэффициентами вида , (11) Использование двух звеньев вида (11), по одному для синфазной и квадратурной компонент (3) обеспечивает уменьшение числа операций умножения на 8 умножений на каждый отсчёт в (3). Модификация координат нулей и полюсов, обусловленная частичной заменой сомножителей в (5) на (11) выполняется использованием (9) для улучшения качества начального приближения. Рис. 5. Вид теоретической СПМ и модернизированной расчетной СПМ, и кривая ошибок. Сравнивая графики рис. 2-3 и рис. 4-5 легко видеть, что приближение рис. 4-5 намного ближе к оптимальному чебышевскому решению (по величине ошибки приближения и числу пульсаций). Приведенные графики получены с применением оригинальных программ в среде MATLAB, разработанных авторами статьи. Заключение Рассмотрено построение ARMA модели канала с релеевскими замираниями и неизоторопным допплеровским рассеянием. Референтная СПМ аппроксимируется по критерию чебышевского приближения квадратом передаточной функции рекурсмвного цифрового фильтра с комплексными коэффициентами. Синтез фильтра выполняется численными методами в виде многоэтапной и многоитерационной процедуры. При создании которой основное внимание уделено проблеме получения «хорошего начального» приближения, реализуемого с уменьшением числа арифметических операций.

作者简介

Sergey Eliseev

Povolzhskiy State University of Telecommunication and Informatics

Email: eliseev-sn@pguti.ru

Pavel Pesotsky

Povolzhskiy State University of Telecommunication and Informatics

Email: pvpesotsky@ya.ru

参考

补充文件