зволяет существенно упростить задачу анализа устойчивости широтно-импульсной системы даже в случае высокого порядка ее линейной части.
Библиографические ссылки
1. Цыпкин Я. З. Теория линейных импульсных систем. М. : Физматгиз, 1963.
2. Цыпкин Я. З., Попков Ю. С. Теория нелинейных импульсных систем. М. : Наука, 1973.
A. M. Miheenko, S. S. Abramov, I. I. Pavlov TO THE ANALYSIS OF PULSE-DURATION SYSTEM WITH A FEEDBACK
Results of research ofpulse-duration system (PDS), covered by the chain of a negative feedback are shown. On the basis ofsymptom-free method of decrease of degree of order of a linear part of a system, the technique of converting PDS to equivalent nonlinear pulse-amplitude system for which the acquainted methods of research are suitable is offered.
Keywords: stability, pulse-duration systems, complex linear chain, nonlinear distortion.
©Михеенко А. М., Абрамов С. С., Павлов И. И., 2010
3. Полов К. П. Условия устойчивости усилителя в режиме Б с обратной связью // Радиотехника. 1971. N° 6.
4. Полов К. П. К исследованию устойчивости усилителя в режиме Б с обратной связью // Радиотехника. 1974. № 1.
5. Ицхоки Я. С. Приближенный метод анализа переходных процессов в сложных линейных цепях. М. : Сов. радио, 1969.
УДК 621.393.3
Е. В. Малинкин
ИНВАРИАНТНАЯ ПЕРЕДАЧА СИГНАЛОВ С ОБРАБОТКОЙ В ЧАСТОТНОЙ ОБЛАСТИ
Предложен способ передачи сигналов с последующей обработкой в частотной области. Для реализации такой обработки разработана структура, основанная на прямом и обратном преобразовании Фурье с элементом деления, позволяющим осуществить деление информационных и обучающих сигналов.
Ключевые слова: помехоустойчивость, инвариант, инвариантная относительная амплитудная модуляция, вероятность попарного перехода, отношение сигнал/шум.
В [1] была предложена относительная амплитудная модуляция (ОАМ). Ее суть заключается в том, что предыдущая посылка служит опорной для последующей, а их отношение на приеме, несмотря на воздействие мультипликативной помехи с коэффициентом передачи канала К, равно отношению этих посылок на передаче. ОАМ позволяет существенно снизить влияние мультипликативной помехи. Однако влияние аддитивной помехи при этом не устраняется.
Автор данной статьи воспользовался идеей ОАМ, но при этом опорные посылки одинаковой величины выделялись им в последовательность одинаковых обучающих сигналов с их дальнейшим усреднением для снижения воздействия аддитивной помехи, а информационные сигналы различной амплитуды выделялись в информационную последовательность. В результате отношение сигналов информационной последовательности к значению усредненных обучающих сигналов также снижало воздействие мультипликативной помехи. Такой способ формирования сигналов на передаче и их демодуляции на
приемной стороне позволяет улучшить помехоустойчивость предложенной структуры как минимум на два порядка по сравнению с классическими системами с ОАМ.
Данная статья посвящена синтезу алгоритма обработки информации инвариантной системой передачи информации в частотной области и анализу полученных технических характеристик.
Постановка задачи. Пусть имеется аналоговый канал связи с полосой пропускания отf до f, по которому передаются амплитудно-модулированные прямоугольной огибающей сигналы. Последовательности сигналов информационной и обучающей частей объединены в блоки.
Передаваемые сигналы подвергаются действию следующих помех: мультипликативной, описываемой изменением коэффициента передачи канала связи Щ) на определенной частоте, и аддитивной, представляющей собой белый шум с некоррелированными отсчетами и дисперсией б2 (при моделировании дисперсия принимает значение, равное 1). Помимо этого на принимаемый сиг-
нал воздействуют амплитудно-частотные и фазочастотные искажения канала связи.
Будем считать, что канал связи подвержен гладким общим замираниям и поэтому можно выделить интервалы стационарности свойств канала, в частности интервалы, на которых коэффициент передачи канала связи и амплитуда отсчетов амплитудно-частотной характеристики (АЧХ) постоянны, т. е. не меняются, а также что длительность передаваемого блока не превышает длительность интервала стационарности.
Для борьбы с аддитивной помехой введено накопление с усреднением обучающих сигналов. Для этого арифметически складываются одноименные отсчеты обучающих сигналов и затем эта сумма делится на количество слагаемых.
Необходимо синтезировать алгоритм обработки сигналов, устойчивый к воздействию комплекса помех (мультипликативной и аддитивной), а также структуру системы передачи информации, позволяющую демодулиро-вать сигналы передачи в соответствии с предложенным алгоритмом с последующим анализом количественных характеристик этой системы.
Синтез алгоритма обработки сигналов. Для решения поставленной задачи потребуется, во-первых, сформировать сигнал передачи (противоположная сторона) и, во-вторых, синтезировать алгоритм обработки на приеме.
Сигнал передачи сформируем блоками (рис. 1). Каждый блок будет состоять из обучающего сигнала (пилот-сигнала) и информационного сигнала. Сигнал обучения будет одинаковым на всех последующих блоках. Для простоты формирования обучающий сигнал Яоб(пТ) представим в виде отсчетов равной амплитуды. В соответствии с законами цифровой обработки сигналов обучающему сигналу Яоб(кТ) соответствует энергетический спектр Яоб( /к).
Энергетический спектр информационного сигнала определим как
5 тф (к ) = 5 об (/к) Я мод (/к), (1)
где Яоб(/к) - энергетический спектр сигнала обучения; Я (/А) - энергетический спектр сигнала передачи. Таким образом, энергетический спектр сигнала передачи на входе блока ОБПФ будет иметь вид
[ Яб (/к) •Я „ (/к) ],;
ІЛ № Я*. №>]_; к;
[5, (>к^ (>к)]; -.
В блоке ОБПФ производится преобразование отсчетов энергетических спектров в отсчеты сигнала передачи Яинф(пТ), а затем в модуляторе с помощью АМ-модуля-ции происходит модуляция информационного сигнала.
Первыми за время Тоб передаются радиоимпульсы одинаковой амплитуды, представляющие собой сигналы обучающей последовательности. Значения амплитуд отсчетов обучающих сигналов на выходе канала связи впоследствии запоминаются в элементе памяти и служат для демодуляции совместно со значениями амплитуд отсчетов информационных сигналов. Кроме того, значения амплитуд отсчетов сигналов обучающей последовательности на выходе канала связи используются и при определении коэффициента передачи канала связи. Затем за время Т передается последовательность радиоимпульсов различной амплитуды, представляющая собой последовательность информационных сигналов. Их количество равно количеству обучающих сигналов, а время передачи информационной последовательности равно времени передачи обучающей последовательности, и поэтому Тинф = Тоб. В дальнейшем вновь следует передача обучающей последовательности, а за ней информационной и так далее до конца передачи всего блока.
После приема обучающих сигналов значения амплитуд отсчетов обучающих сигналов в элементе памяти арифметически складываются и усредняются по количеству слагаемых для борьбы с аддитивной помехой, что повышает помехоустойчивость всей системы передачи информации.
Отсчеты информационной и обучающей последовательности сигналов и их энергетические спектры характеризуются следующими особенностями (рис. 2). Следует отметить, что амплитуды отсчетов обучающего сигнала во временной области на входе канала связи равны между собой, как и значения амплитуд отсчетов энергетических спектров обучающих сигналов, а амплитуды отсчетов сигналов информационной последовательности на входе канала связи отличаются своей величиной во временной области и имеют различные значения амплитуд отсчетов энергетических спектров сигналов. Совокупности N отсчетов информационной части во временной области соответствует N отсчетов его энергетического спектра. Аналогично совокупности N отсчетов обучающего сигнала во временной области соответствует N отсчетов его энергетического спектра.
Тогда энергетический спектр для сигнала информационной последовательности /-го блока может быть представлен в виде
^вых (Д) = ^ (Д)• к.(]к)• И1 (]к) + N(]к), (2)
где к - номер гармоники; З1. (/к) - энергетический спектр сигнала на выходе канала связи на /-м блоке; З (/к) -энергетический спектр сигнала на входе канала связи на /-м блоке; к(/к) - энергетический спектр мультипликативной помехи на /-м блоке; И(/к) - передаточная харак-
Рис. 1. Структура передающего устройства:
БПФ - быстрое преобразование Фурье; ОБПФ - обратное преобразование Фурье
теристика канала на і-м блоке обработки; N (/к) - спектральная плотность мощности аддитивного шума сигнала информационной последовательности на і-м блоке обработки.
Энергетический спектр для сигнала обучающей последовательности і-го блока находится аналогично:
х,..< /к)=«.б. (/к )х
хк( /к )• Н(/к) + ,У( /к ), (3)
где «б (/к) - энергетический спектр обучающего сигнала на выходе канала связи; «іоб вх(/к) - энергетический спектр обучающего сигнала на входе канала связи; N (/к) - спектральная плотность мощности аддитивного шума сигнала обучающей последовательности на і-м блоке обработки.
Следует заметить, что после усреднения сигналов обучения значения амплитуд отсчетов N (/к) будут много меньше значений амплитуд отсчетов N (/к) за счет эффекта накопления с усреднением. Значения амплитуд отсчетов N (/к) больше, чем значения амплитуд отсчетов Ж2і(/'к) в 1 ООО раз, т. е. N (к) = 1 ОООЫ2(/к), при количестве усреднений, равном 1 ООО [5].
Как было отмечено ранее, модулирующий параметр на передаче модулирует отношение энергетических спектров информационной и обучающей частей. Тогда для демодуляции сигналов приема необходимо поделить энергетический спектр информационной части на энергетический спектр обучающей части. В результате получим:
(/к) =
(/к) (/к)
= ^ (/к)•к, (/к)•И, (/к) + N (/к) (4)
= ( /к)•к ( )•И, ( ) + NЛ ).
В силу того что значения амплитуд отсчетов энергетического спектра аддитивного шума в обучающей последовательности существенно меньше значений амплитуд отсчетов энергетического спектра аддитивного шума информационной последовательности, то величиной N (/к) в выражении (4) можно пренебречь [3]. Тогда выражение (4) после элементарных преобразований сводится к виду
«мод.вых (/к ) = ЯШы* (/к) +
_________________N1i^ (/к)
«іоб.вх (/к )• к (/к)• Н< (/к) ’
(5)
где З/модвых(/к) - энергетический спектр модулирующей последовательности на выходе канала связи на /-м блоке; З.модвх(/'к) - энергетический спектр модулирующей последовательности на входе канала связи на /-м блоке.
По теореме Парсеваля, энергия сигнала, вычисленная во временной области, равна энергии сигнала, вычисленной в частотной области. Поэтому отношение энергии информационного сигнала к энергии обучающего сигнала сохраняется как во временной, так и в частотной области. Следовательно, результаты обработки во
ад
Информационная
последовательность
Обучающая по следовательно сть
О, 1, 2,
N-1 0, 1, 2,
N-1 0, 1, 2,
N-1 0, 1, 2,
N-1 п
і-й блок
г+1-й блок
б
Рис. 2. Информационные и обучающие последовательности сигналов (а) и их энергетические спектры (б)
а
временной и частотной области дают одинаковый результат, в связи с чем оценку помехоустойчивости произведем во временной области. Для этого обратимся к рис. 3, показывающему принцип такой обработки сигналов.
В точке а информационные и обучающие сигналы представлены в виде отсчетов во временной области. На выходе БПФ (точка б) сигнал представлен отсчетами энергетического спектра, в точке в - отсчетами энергетического спектра информационной части, в точке г - отсчетами энергетического спектра обучающей части, в точке д - отсчетами энергетического спектра модулирующей последовательности совместно с отсчетами энергетического спектра аддитивного шума, величина которого оценивается по формуле (5). В точке е сигнал представлен отсчетами модулирующей последовательности во временной области.
Для расчета вероятности ошибочного приема во временной области введем понятие инварианта.
С учетом используемой относительной амплитудной модуляции оценку значения инварианта можно найти следующим образом [2]:
• =
=1
об N • Жу. • 5
об
X Зоб " об
/=1
Будем оценивать вероятность перебора инвариантов, выбрав первым в паре самый наименьший по значению инвариант - ГЫУ1 = 1, а в качестве второго инварианта перебирая значения от 2 до 6. При переходе одного инварианта в другой происходит ошибочный прием, когда вместо первого инварианта в паре сравниваемых принимается второй с определенной вероятностью и, наоборот, вместо второго инварианта в паре сравниваемых инвариантов принимается первый с той же вероятностью.
Произведем расчет вероятности ошибочного приема. Для этого воспользуемся известным подходом [2]:
(6)
где Рпер - вероятность перехода первого инварианта во второй и наоборот; Р1 - вероятность появления первого инварианта; Р2 - вероятность появления второго инварианта; первый интеграл - это вероятность появления второго инварианта, когда послан первый; второй интеграл
- вероятность появления первого инварианта, когда послан второй; хр - пороговое значение, необходимое для вычисления Рпер, при известных Р1 и Р2.
Величина хр определяется с помощью наилучшей байесовской оценки путем минимизации Р по г . Полага-
пер р
ем, что появление первого и второго инвариантов равновероятно и поэтому выбираем Р1 = Р2 = 0,5.
Для дальнейших вычислений, исходя из выражения (6), необходимо знать аналитические выражения для Ж1(х) и Ж2(х).
Для когерентного приема с синусоидальной поднесу-щей расчет функции Ж1(х) и Ж2(х) известен и приведен в [2].
В нашем случае прямоугольная огибающая сигнала выделяется с помощью синхронного детектора, а значит, помеха имеет нормальное распределение [4]. Поэтому мы можем использовать тот же подход для нахождения аналитического выражения плотности вероятности оценки инварианта [2].
Итак, величина оценки инварианта в рассматриваемой нами системе определяется по выражению N
=
у (к • шу.+х(і)) і=1
І N
- 5о,
(7)
т УУ(к5об + П(т /'))
т т=1 / =1
' т=1 /=1
где в числителе ГЫУ7 - первый передаваемый инвариант, Х(/) - /-е значение гауссовой помехи; в знаменателе 5'об -значение обучающего сигнала, "п(т,/) -/-е значение гауссовой помехи в т-й реализации сигнала Зоб, к - коэффициент передачи канала связи, N- число отсчетов, взятых по огибающей ГЫУ, или З Ь - число обучающих сигналов.
Для расчета Рпер необходимо знать математическое ожидание и дисперсию числителя и знаменателя выражения (7). Для их расчета воспользуемся следующими выражениями:
Рис. 3. Структурная схема инвариантной системы передачи информации
- математическое ожидание числителя (7)
т = к • N • ШУ,;
числ ,
- дисперсия числителя (7)
Б = N • ст2,
(8)
(9)
где ст2 - дисперсия гауссовой помехи;
- математическое ожидание знаменателя (7) после преобразований
т = к • N (10)
знам 7 4 '
- дисперсия знаменателя (7) после преобразований
Б =-
(11)
Тогда выражение для плотности вероятности оценки инварианта с учетом (8)—(11) будет иметь вид [2]:
1
W ( 2 ) = |
х е
2лст ст
-¥ ! 2
хх26 (х-т )2 е 2"а2 |х|дх, (12)
; т - количество сигналов об-
к2 МУ,2
1 2 3 4 5 6 7
9 10 11 к
10'
10'
10'
10'
10'
10'
10'
10“:
■10
■20
■30
■40
■50
■60
■70
10
■90
-100
2
и
10'
1 пер
Рис. 4. Вероятность попарного перехода одного инварианта в другой при следующих заданных условиях: к = 1; ШУ1 = 1; МУ2 = 2, 3, ..., 10; кривая 1 - когерентный прием, инвариантная относительная амплитудная модуляция, прямоугольная огибающая; кривая 2 - классическая относительная амплитудная модуляция
Следует отметить, что особенностью любой инвариантной системы, основанной на принципе инвариантной относительности амплитудной модуляции, является то, что по каналу передаются амплитудно-модулированные сигналы, образованные и Передача этих сигна-
лов на основе классических алгоритмов обработки информации обеспечивает невысокую помехоустойчивость [4]. И только после обработки этих сигналов в соответствии с алгоритмом частного по выражению (7) будет получена оценка инварианта, являющегося по сути числом, а не сигналом. Вероятность попарного перехода одного инварианта в другой при больших отношениях сигнал/шум определяется величинами 10-30...10-80 (см. рис. 4, 5). При тех же значениях отношения сигнал/шум вероятность ошибочного приема единичного символа в классических системах лежит в пределах 10-6...10-10.
12345678 А
где Ст1 = ^/Бчисл работки.
Расчет Рпер проводился численно аппроксимацией формулы (12). Предложенная система сравнивалась с классической когерентной системой. Вероятность попарного перехода вычислялась в обоих случаях для одинаковых значений й-отношения сигнал/шум по формуле
У к 2 МУ2 А/
й2 = ^_____________=_________
N А/ст2 ст2 '
Пороговые значения рассчитывались путем минимизации Рпер в формуле (6). Для к = 1 и МУ1 = 1, ГЫУ2 = 2,
3, 4, 5, 6 вычисления дают результат = 1,5; 2; 2,5; 3; 3,5. Для к = О,7 и ГЫУ1 = 1, ГЫУ2 = 2, 3, 4, 5, 6 получены значения г = 1,5; 2; 2,5; 3; 3,5. р
Результаты моделирования приведены ниже (рис. 4, 5).
Рис. 5. Вероятность попарного перехода одного инварианта в другой при следующих заданных условиях: к = О,7; ШУ1 = 1;
ШУ2 = 2, 3, ..., 1О; кривая 1 - когерентный прием, инвариантная относительная амплитудная модуляция, прямоугольная огибающая; кривая 2 - классическая относительная амплитудная модуляция
По мнению автора, помехоустойчивость исследуемой инвариантной системы необходимо сравнивать с помехоустойчивостью аналогичных инвариантных систем, что и будет сделано в последующих работах.
Таким образом, предложена инвариантная когерентная система передачи информации и определены ее качественные характеристики. Проведенный анализ показывает, что эта система обладает высокой помехоустойчивостью: вероятность попарного перехода в инвариантной системе как минимум на два порядка меньше вероятности ошибки классического алгоритма с относительной амплитудной модуляцией. Поэтому предложенную систему следует использовать в телекоммуникационных системах, системах телеуправления и других системах, предъявляющих высокие требования к помехоустойчивости.
Библиографические ссылки
1. Петрович Н. Т. Передача дискретной информации в каналах с фазовой манипуляцией. М. : Сов. радио, 1965.
ст
2. Инвариантный метод анализа телекоммуникацион- 3. Левин Б. Р. Теоретические основы статистической
ных систем передачи информации : монография / радиотехники. 3-е изд. М. : Радио и связь, 1989.
В. Б. Малинкин, Е. И. Алгазин, Д. Н. Левин, В. Н. Попанто- 4. Теплов Н. Л. Помехоустойчивость систем передачи
нопуло. Красноярск, 2006. дискретной информации. М. : Связь, 1964.
E. V Malinkin
INVARIANT TRANSMISSION OF SIGNALS WITH PROCESSING IN FREQUENT FIELD
The work offers a method of transmission of signals with the subsequent processing in the frequency field. In order to realize such kind of processing, the structure based on direct and inverse Furier transformation with division element that makes it possible to realize the quotient of the current and subsequent processing block is offered.
Keywords: noise immunity, invariant, invariant relative amplitude modulation, probability of pairwise transition, signal/noise relation.
© Малинкин Е. В., 2010
УДК629.7.05
О. В. Зайцев
МЕТОДИКА ВОССТАНОВЛЕНИЯ НЕСУЩЕЙ ЧАСТОТЫ СИГНАЛА ТЕЛЕМЕТРИЧЕСКИХ СИСТЕМ ПРИ ЦИФРОВОЙ ОБРАБОТКЕ
Рассмотрены вопросы восстановления уровня несущей частоты телеметрического сигнала при цифровой обработке в тракте автоподстройки частоты и вычисления порога для принятия решения об истинности принятого символа информации, поступающей с космических аппаратов и ракетоносителей. Описаны методика и алгоритм вычисления уровня несущей частоты и значения порога для принятия решения на основе гистог-раммной обработки сигнала с выхода частотного детектора.
Ключевые слова: система управления, обработка сигналов, телеметрия.
Контроль правильности отработки полетного задания бортовыми системами ракетоносителей (РН) и космических аппаратов (КА) осуществляется на основе обработки телеметрической информации (ТМИ) о состоянии большинства агрегатов и элементов объекта [1]. Для оптимального приема и обработки ТМИ, поступающей на наземные средства приема, необходимо восстановить несущую частоту радиосигнала и вычислить порог принятия решения для достоверного определения принимаемого символа информации (нуля или единицы).
В настоящее время в цифровых трактах приемных устройств для восстановления несущей радиосигнала с частотной манипуляцией чаще всего используются традиционные методы: быстрое преобразование Фурье (БПФ) или линейная фильтрация информационного параметра [2]. Применение БПФ требует большого количества производимых операций для получения необходимого результата и, следовательно, значительного времени для проведения анализа сигнала. А при линейной фильтрации весьма сложно получить несмещенную оценку обрабатываемого параметра. Кроме того, традиционные методы при обработке сигнала на фоне сильных шумов не позволяют достоверно определять несущую радиосигнала [3], что значительно сни-
жает чувствительность цифровых радиоприемных трактов.
Методики и алгоритмы восстановления уровня несущей сигнала, альтернативные БПФ и методу линейной фильтрации, в настоящее время по-прежнему остаются слабо освещенными в отечественной и зарубежной литературе [4-6]. Поэтому разработка методики и алгоритма восстановления несущей радиосигнала на основе ги-стограммного метода, требующего проведения значительно меньших по сравнению с БПФ вычислительных затрат, является весьма актуальной.
Суть этого метода заключается в следующем. Диапазон возможных уровней сигнала на выходе частотного детектора разбивается на оптимальное число контрольных уровней, или интервалов. Сигнал после частотного детектора подается на блок построения гистограммы уровней. В этом блоке значения амплитуды сигнала сравниваются со значениями контрольных уровней, при этом фиксируется число попаданий значений амплитуды А - А
сигнала в интервал А к = — --—, где Ак и Ак-1 - значе-
ния соседних контрольных уровней. Затем формируется одномерный вектор, т. е. вектор числа попаданий уровня сигнала в интервал «А(к) - гистограмма Н(к). Параллель-