CC BY
49
^ Научно-технические ведомости СПбГПУ. Физико-математические науки № 1(189) 2014 ^ УДК 621.391
А.Г. Голубев, П.А. Молчанов
ОАО «Камчатский гидрофизический институт»
алгоритмы оценивания импульсной характеристики
МНОГОЛУЧЕВОГО КАНАЛА СВЯЗИ С ТРАНСФОРМАцИЕЙ РАБОЧЕГО
ДИАПАЗОНА ЧАСТОТ
При передаче данных в многолучевом канале связи, как правило, по реакции канала на испытательный импульс (ИИ) осуществляется предсказание формы каждого принимаемого информационного сигнала (ИС). Во избежание межсимвольной интерференции в последовательности ИИ целесообразна ее передача со скачкообразным изменением несущей частоты. При этом возникает необходимость оценивания импульсной реакции канала в полосе частот ИС по реакции канала на последовательность ИИ, передаваемых в разных диапазонах частот, в общем случае не совпадающих с полосой частот ИС. Настоящая статья посвящена алгоритмам такого оценивания.
МНОГОЛУЧЕВОЙ КАНАЛ, ИСПЫТАТЕЛЬНЫЙ ИМПУЛЬС, АЛГОРИТМ ОЦЕНИВАНИЯ, ПОГРЕШНОСТЬ ОЦЕНИВАНИЯ, ИМПУЛЬСНАЯ РЕАКЦИЯ, КОРРЕЛЯЦИЯ, АВТОКОРРЕЛЯЦИОННАЯ ФУНКЦИЯ.
Введение
Большинство каналов связи сравнительно большой протяженности характеризуется многолучевым распространением. Форма передаваемого сигнала при таком распространении искажается случайным образом, что приводит к существенному снижению эффективности когерентной обработки при приеме. В таких случаях применяется, в частности, так называемая «система с испытательным импульсом и предсказанием» (СИИП) [1, 2, 8]. Реализация СИИП предполагает несколько функций:
передачу информационных сигналов (ИС), соответствующих символам передаваемого алфавита;
передачу испытательных импульсов (ИИ), не содержащих информации;
оценивание по ИИ импульсной характеристики канала распространения (ИХК);
предсказание (прогноз) формы каждого из выходных ИС, т. е. ИС на выходе канала связи.
Указанный прогноз каждого из выходных информационных сигналов осуществляется посредством вычисления свертки соответствующего исходного ИС (т. е. сигнала на входе канала связи) с ИХК.
Чтобы обеспечить необходимую точ-
ность оценивания ИХК, испытательные импульсы должны передаваться многократно (непрерывно), однако при этом (особенно в случае их передачи при повторяющейся форме) происходят искажения формируемой оценки ИХК, вызванные эффектом, аналогичным межсимвольной интерференции (МСИ). Последняя имеет место при передаче через многолучевой канал информационных сигналов. Для исключения МСИ необходимо, например, передавать последовательность испытательных импульсов со скачкообразным изменением несущей частоты [3, 9]. При этом для накопления (осреднения) оценок ИХК, сформированных по разным испытательным импульсам (они передаются на разных частотах), необходимо «привести» эти оценки к одному частотному диапазону, а именно к диапазону (либо диапазонам) передачи информационных сигналов. Данная задача никем до нас не рассматривалась.
Настоящая статья посвящена решению указанной задачи, а именно разработке алгоритма оценивания ИХК с трансформацией диапазона частот, в котором по каждому испытательному импульсу формируется оценка ИХК (далее для краткости именуем их алгоритмами трансформации).
Модель многолучевого канала и используемая терминология
Моделями канала с рассеянием по времени могут служить соответствующие четырехполюсники с переменными параметрами, характеризующиеся импульсными характеристиками, зависящими от времени [4, 5, 10]. При этом принимаемый сигнал есть свертка переданного сигнала с ИХК. Импульсная характеристика четырехполюсника определяется как его реакция на входное воздействие в виде 5-функции; такое определение соответствует ИХК в бесконечной полосе частот. Применительно к многолучевому каналу, по которому передаются сигналы в ограниченной полосе частот, оправдано рассмотрение модели, содержащей в качестве входной цепи узкополосный фильтр, полоса пропускания которого совпадает с частотным диапазоном передаваемого сигнала [4]. Импульсная характеристика такого канала (его узкополосной модели) есть свертка ИХК с импульсной характеристикой указанного полосового фильтра; последняя же практически совпадает с автокорреляционной функцией (АКФ) передаваемого сигнала. Определенную таким образом ИХК будем называть УИХК (импульсная характеристика канала, определенная в узкой полосе частот).
Оценка УИХК в диапазоне частот испытательного импульса Б(п) (здесь и далее п — индекс аргумента дискретного времени) вычисляется как свертка принимаемого сигнала с этим импульсом (см., например, работу [3]). Нетрудно установить, что результат оценивания указанной УИХК (с точностью до реализации аддитивного шума) есть свертка ИХК с АКФ испытательного импульса Б(п).
Постановка задачи
Задача настоящей статьи — разработка и сравнение алгоритмов пересчета оценки УИХК, полученной в диапазоне частот испытательного импульса Б1(п), в оценку импульсной характеристики этого же канала в диапазоне частот ИС или испытательного импульса Б2(п). Последний передается в полосе частот, не совпадающей с таковой
для сигнала Б1(п).
Эта задача ставится следующим образом. Имеем следующие исходные соотношения:
У» = z(п) ® Б1(Ы - п) =
= АКФ» ® к(п) + N1);
у2(п) = *(п) ® Б2(М - п) = = АКФ2(п) ® Н(п) + N2),
где у1(п), у2(п) - оценки УИХК1 и УИХК2, которые были бы сформированы в случаях передачи испытательных импульсов Б1(п) и S2(n) при их взаимно-корреляционной обработке соответственно; *(п) — реализация сигнала на входе приемника; ® — операция свертки; М = т (/ , т — частота дискретизации каждого из результатов оценивания УИХК и длительность каждого из сигналов Б1(п) и Б2(п)); АКФ12 (п) — АКФ сигналов ^(п); к(п) — ИХК; ' Щп), Щп) — реализации шума или ошибки оценивания соответствующих УИХК.
Вместо оценки УИХК1 у1(п) решение задачи может быть также основано непосредственно на анализе наблюдаемой приемником реализации сигнала *(п).
Задана форма каждого из испытательных импульсов Б12(п), причем вместо одного из них можно рассмотреть и информационный сигнал. Фактически передан испытательный импульс Б1(п), и на основе взаимно-корреляционной обработки этого импульса сформирована оценка УИХК1 у 1(п); при этом форма ИХК к(п) неизвестна. Требуется вычислить оценку УИХК2, близкую по форме к той, которая была бы сформирована в случае передачи испытательного импульса Б2(п) на основе его взаимно-корреляционной обработки (это есть оценка у2(п)). Другими словами, требуется пересчитать оценку УИХК из диапазона частот испытательного импульса Б1(п) в аналогичную оценку УИХК (того же канала) в диапазоне частот испытательного импульса Б2(п); такую пересчитанную оценку обозначаем как у 1^2(п).
Варианты алгоритма трансформации (пересчета) рабочего диапазона частот
В рамках настоящей работы предлагаются и анализируются три варианта алго-
ритма трансформации рабочего диапазона частот (алгоритмы 1 — 3).
Алгоритм 1. Это алгоритм с промежуточным оцениванием непосредственно ИХК. Согласно этому варианту алгоритма, на основе заданной оценки УИХК1 у1(п) реализуется прямое вычисление оценки ИХК к(п) и далее вычисляется искомая оценка УИХК2 как свертка:
у^2(й) = АКФ2(й) ® к(п). (2)
При оценивании ИХК последняя рассматривается как решетчатая функция, представляющая собой последовательность бесконечно коротких импульсов, амплитуды которых и расположение на временной оси соответствуют амплитудам и относительным задержкам лучей, приходящих в точку приема. В этом случае формирование оценки ИХК сводится к следующим этапам:
нахождение оценки стандарта шума ст ш в реализации у 1(п);
вычисление величины порога обнаружения сигнала в луче П = к ст ш (к — регулируемая константа);
фиксация всех отсчетов массива у 1(п), превышающих порог П (т. е. принятие решений об обнаружении сигналов в лучах);
оценивание амплитуд этих отсчетов и их положения на оси времени (т. е. амплитуд и положения на оси времени сигналов, приходящих в точку приема в отдельных лучах).
Эмпирически установлена целесообразная величина константы: к = 4,5. Оцениванию стандарта шума предшествует отбор 10 % наибольших по уровню отсчетов из массива у 1 (п) и дальнейшее их удаление из этого массива. Затем искомая оценка стандарта шума Стш классически вычисляется по массиву оставшихся отсчетов. Оценивание амплитуд к. и временного положения лучей (отсчет времени выполняется относительно середины ИХК) производится на основе соответствующей параболической интерполяции следующим образом [6]:
к . = у1(п) - 0,25[у 1(п1 + 1) - у 1(п1 -1)]Ы;
(п + Ы).
Ы = о, 5[у1(п -1) - уДп + 1)] / [уДп + 1) -- 2у 1(п.,) + у1(п., -1)],
где п. — номер временного отсчета, соответствующего моменту обнаружения сигнала в 1-м луче.
Алгоритм 2. Это алгоритм, следующий из соотношений (1) с учетом теоремы о свертке [7]. При сформированной оценке УИХК1 у1(п) вычисление искомой оценки УИХК2 у 2(п) реализуется следующим образом:
у 1^2(п) = ОДПФЩк)};
У2(к) = Б,(к)у 1(к)/ Б2(к),
(4)
Ь =
Л
где ОДПФ{72(к)} — результат выполнения операции обратного дискретного преобразования Фурье над массивом спектральных отсчетов У2(к); Б1(к), £2(к), у1(к) — результаты вычисления операции дискретного преобразования Фурье от реализаций сигналов 5» и Б(п), а также оценки УИХК1 У 1(п).
Все спектры, входящие в соотношения (4), — это спектры соответствующих комплексных огибающих.
Следует отметить, что этот вариант алгоритма (в отличие от двух остальных) не связан с ограничением модели многолучевого канала. Ограничение вызвано предположением, что ИХК представляет собой решетчатую функцию.
Алгоритм 3. Этот вариант представляет собой комбинацию двух остальных. Он предусматривает вычисление оценки у^2(п)пр по формуле (3), но эта оценка рассматривается только как промежуточная (в отличие от алгоритма 2). Далее над ней выполняются все операции, которые в алгоритме 1 выполняются над оценкой у1(п) с получением оценки ИХК к (п). После этого искомая оценка у1^2(п) формируется по формуле (2).
Методика эмпирического анализа качества алгоритмов трансформации рабочего диапазона частот
Критерием качества алгоритма трансформации служит величина отношения
К = к1 / к0. Величины к0 и к1 определяются как
Е У У2 оИ к1 = , ^2 ^ ;
ЛЕ у 2^2(п)Е у 2 о(п)
, Г (5)
Е у '2 (п)у "2И к 0 = я . ДййЁйй
V п п
Здесь к1 — значение нормированной корреляции между оценкой УИХК2 у1^2(я), пересчитанной из диапазона частот сигнала S1(n) в диапазон частот сигнала S2(n), и ее условно эталонной формой у 2 0(п), которая бы имела место при оценивании непосредственно в диапазоне частот сигнала S2(n). (Условно эталонная величина оценки УИХК2 у 2 0(п) формируется путем вычисления свертки
у2 0(п) = АКФ2(п) ® к(п)
и добавления к ней реализации аддитивного шума); к0 — значение нормированной корреляции, вычисляемой (в отличие от значения к1) между двумя ее эталонными формами у '2(п) и у "2(п), совпадающими между собой с точностью до реализаций аддитивного шума.
Все упомянутые реализации аддитивного шума не корреллированы между собой.
В процессе анализа качества задается исходная ИХК к(п), нормированная так, что ее максимум равен единице. Затем вычисляется свертка этой ИХК с АКФ1 сигнала S1(n). Результат вычисления указанной свертки у1(п) суммируется с реализацией шума, и в итоге получается оценка у1(п) УИХК1 в диапазоне частот сигнала S1(n). Уровень (стандарт) шума ош варьируется. Далее моделируется собственно алгоритм пересчета, в результате чего формируется искомая оценка УИХК2, приведенная к диапазону частот сигнала S2(n) (т. е. оценка У1^(п)).
Представленное выше определение критерия качества алгоритмов трансформации рабочего диапазона частот обосновывается однозначной связью указанной величины К применительно к СИИП (в широком круге
вариаций этого метода передачи данных) с потерями в помехоустойчивости. Последние определяются выражением 10 ^К. Под потерями в данном случае подразумевается разница в пороговом отношении сигнал/ шум на входе приемника (декодера); она имеет место при оценивании УИХК непосредственно в требуемом диапазоне частот сигнала S2(n) (в случае передачи испытательного импульса или сигнала S2(n)), а также при использовании одного из рассматриваемых алгоритмов пересчета этой оценки из диапазона частот сигнала (испытательного импульса) S1(n) в требуемый диапазон частот испытательного импульса S2(n). Актуальность такого пересчета как альтернативы оценивания УИХК непосредственно в требуемом диапазоне частот оговорена во введении.
Результаты анализа качества алгоритмов трансформации
Анализ качества алгоритмов трансформации ИХК произведен в соответствии с изложенной выше методикой для двух вариантов исходной ИХК.
Рассматриваемые варианты исходной ИХК служат примерами, характеризующими сравнительно простую (рис. 1, а) и сложную с точки зрения решаемой задачи (рис. 1, б) ситуации. При этом вариант импульсной характеристики канала, типичной для сравнительно простой ситуации и содержащей малое количество лучей, выбран произвольно; при этом он вполне возможен на практике. Вариант же указанной характеристики, типичной для сравнительно сложной ситуации и содержащей большое количество лучей, был сформирован на основе реальных данных при временном разрешении примерно 3 мс, ввиду разнообразия возможного сочетания амплитуд и задержек сигналов в лучах.
Результаты вычисления критерия качества К применительно к исходным ИХК (см. рис. 1) приведены на рис. 2. Здесь и далее отношение сигнал/шум определяется как -20-1§ош (дБ).
На основании анализа графиков на рис. 2 можно сделать следующие выводы.
Алгоритмы 1 и 3 дают идентичные ре-
а)
0,0 — -1.0
0,0 Время, с
б)
0,6
0,0 -1,0
и
д
I
Время, с
Рис. 1. Варианты исходной импульсной характеристики канала при малом (а) и большом (б) количестве лучей, приходящих в точку приема
зультаты, с точки зрения указанного выше критерия качества.
Алгоритм 2 при малом количестве лучей дает результат несколько хуже, чем алгоритмы 1 и 3, за исключением области очень низких (ниже —15 дБ) отношений сигнал/ шум.
Алгоритм 2 при малых отношениях сигнал/шум значительно более эффективен, чем алгоритмы 1 и 3.
В случае большого количества лучей в области больших отношений сигнал/шум соотношение эффективностей рассматриваемых алгоритмов то же самое, что и при малом количестве лучей.
Для выбора предпочтительного алгорит-
а)
к
0,6
0,5
■20 -15 -10 -5 0 5 10 15 20
Сигнал/шум, дБ
ма для конкретных «промежуточных» случаев (типов ИХК) целесообразно провести соответствующие расчеты по предложенной методике.
Заключение
Разработаны варианты алгоритма оценивания импульсной характеристики многолучевого канала связи с трансформацией рабочего диапазона частот, т. е. пересчета оценки указанной характеристики в ее узкополосном представлении из одного диапазона частот в другой. Разработка позволяет реализовать передачу данных в соответствии с принципом СИИП при чередовании (скачкообразном изменении)
б)
1.0 I-1-1-1-1-1-1-1-
0,9
■20 -15 -10 -5 0 5 10 15 20
Сигнал/шум, дБ
Рис. 2. Результаты вычисления критерия качества К применительно к исходным ИХК, показанным на рис. 1, а (а) и б (б). Номера кривых соответствуют номерам алгоритмов
диапазонов частот передачи испытательных импульсов. Применение такого преобразования, как известно, позволяет исключать искажения формы этих импульсов из-за межсимвольной интерференции.
Приведены методика анализа качества алгоритмов трансформации и примеры расчетов по предложенной методике. Показано, что соотношение эффективностей рассмотренных алгоритмов определяется
как формой импульсной характеристики канала, так и отношением сигнал/шум, при котором осуществляется указанная трансформация. Установлено, что в зависимости от конкретных условий (количество принимаемых лучей и отношение сигнал/шум) предпочтительными могут быть как алгоритм 2, так и эквивалентные между собой (в свете рассмотренного критерия качества) алгоритмы 1 и 3.
СПИСОК ЛИТЕРАТУРЫ
1. Гольденберг Л.М., Кловский Д.Д. Метод приема импульсных сигналов, основанный на использовании вычислительных машин // Труды ЛЭИС. 1959. Вып.УП. С. 17-26.
2. Кловский Д.Д. Передача дискретных сообщений по радиоканалам; 2-е изд. М.: Радио и связь, 1982. 304 с.
3. Скляр Б. Цифровая связь. Теоретические основы и практическое применение; 2-е изд. М.: ИД «Вильямс». 2003. 1105 с.
4. Кейлос Т. Каналы с параметрами, изменяющимися во времени. Лекции по теории связи. Под ред. Е.Дж. Багдади. М.: Мир, 1964. С. 68-122.
5. Ван-Трисс Г. Теория обнаружения, оценок и модуляции. Т. 3. М.: Сов. радио, 1975. 664 с.
6. Оппенгейм Э. Применение цифровой обработки сигналов. М.: Мир, 1980. 552 с.
7. Голд Б., Рейдер Ч. Цифровая обработка сигналов. М.: Сов радио, 1973. 368 с.
8. Кловский Д.Д. Системы оптимального приема в каналах с эхо-сигналами // Труды учебных институтов связи. 1963. Вып.18.
9. Прокис Д.Д. Цифровая связь. М.: Радио и связь, 2000. 800 с.
10. Зюко А.Г., Кловский Д.Д., Назаров М.В., Финк Л.М. Теория передачи сигналов. М.: Связь, 1980. 288 с.
СВЕДЕНИЯ ОБ АВТОРАХ
ГОЛУБЕВ Анатолий Геннадиевич — доктор технических наук, заместитель по научной работе генерального директора ОАО «Камчатский гидрофизический институт». 197183, Россия, г. Санкт-Петербург, Липовая аллея, 9 agg300@mai1.ru
МОЛЧАНОВ Павел Александрович — кандидат физико-математических наук, начальник отдела комплексного проектирования ОАО «Камчатский гидрофизический институт». 684090, Россия, Камчатский край, г. Вилючинск, Академическая ул., 3 pasha@kgfi.ru
GolubevA.C., MolchanovP.A. ESTIMATION ALGORITHMS FOR IMPULSE RESPONSE OF MULTIBEAM COMMUNICATION CHANNEL WITH TRANSFORMATION OF THE WORKING FREQUENCY BAND.
In transmission messages, the prediction about the form of every information signal, as a rule, is carried out using response of communication channel on test impulses. To avoid errors in predictions of intersymbol interference it is useful to apply the jump of carrier frequencies from one test impulse to another one. In this case it is necessary to make an estimate of channel impulse response in frequency band of information signals from channel response on series of test impulses, which transmitted in different frequency bands. Our article focuses on the algorithms for this estimation.
MULTIBEAM COMMUNICATION CHANNEL, TEST PULSE, ESTIMATION ALGORITHM, ESTIMATION ERROR, PULSE RESPONSE, CORRELATION, AUTOCORRELATION FUNCTION.
REFERENCES
1. Goldenberg L.M., Klovsky D.D. Method computers. Proceeding of LEIC, 1959. No.VII. of receiving of impulse signals based on using of pp. 17-26 (rus).
2. Klovsky D.D. Transmission of discrete messages in radio channels. 2 edition, M.: Mir Publishers: Radio and communication. 1982. 304 p. (rus)
3. Sclyar B. Digital communication. Theoretical basis and practice. 2 edition., Moscow: «Viliams» Publishers. 2003. 1105 p. (rus)
4. Keilos T. Channels with time-varying parameters. Lectures on communication theory. Ed. E.D. Bagdady. Moscow: Mir Publishers, 1964. pp. 68-122. (rus)
5. Van-Triss G. Theory of detection, estimation and modulation. Vol. 3. Moscow: Soviet Radio Publishers, 1975. 664 p. (rus)
6. Oppengeim E. Implementation of digital signal processing. Moscow: Mir Publishers, 1980. 552 p. (rus)
7. Gold B., Reider C. Digital signal processing. Moscow: Soviet Radio Publishers, 1973, 368 p. (rus)
8. Klovsky D.D. System of optimal receiving in eho-channels. Proceeding of educational communication institutes. 1963. Vol.18. (rus)
9. Proakis J.G. Digital communications. Moscow: Radio and communication, 2000, 800 p. (rus)
10. Zyko A.G., Klovsky D.D., Nazarov M.V., Fink L.M. The theory of signal transmission. Moscow: Communication, 1980. 288 p. (rus)
THE AUTHORS
GOLUBEV Anatoly G.
Kamchatka Hydrophysical Institute JSС 197183, Lipovaya Ave. 9, St. Petersburg, Russia agg300@mail.ru
MOLCHANOV Pavel A.
Kamchatka Hydrophysical Institute JSС
684090, Akademicheskaya Str. 3, Vilyuchinsk, Kamchatka Krai, Russia pasha@kgfi.ru
© Санкт-Петербургский государственный политехнический университет, 2014
