ISSN 0868-5886
НАУЧНОЕ ПРИБОРОСТРОЕНИЕ, 2006, том 16, № 4, c. 106-111
= ОБРАБОТКА И АНАЛИЗ СИГНАЛОВ
УДК 621.391 © С. В. Дворников
ТЕОРЕТИЧЕСКИЕ ОСНОВЫ ПРЕДСТАВЛЕНИЯ СИГНАЛА В АНАЛИТИЧЕСКОМ ВИДЕ ФУНКЦИЯМИ ЕГО ОГИБАЮЩЕЙ
И ПОЛНОЙ ФАЗЫ
Излагаются теоретические аспекты представления сигнала в аналитическом виде функциями: огибающей сигнала и его полной фазой, представленными в частотно-временном и временном пространствах. Приводятся результаты теоретических исследований. Материалы статьи будут интересны специалистам, занимающимся синтезом частотно-временных распределений в приложении к задачам измерения параметров узкополосных нестационарных процессов в шумах высокой интенсивности.
ВВЕДЕНИЕ
занные между собой преобразованием Гильберта:
Методы частотно-временной обработки находят все более широкое практическое применение в различных областях радиотехники. Большинство работ, посвященных указанной проблематике, как правило, начинаются с анализа свойств функции Вигнера, которая базируется на понятии аналитического представления комплексного сигнала (аналитического сигнала). Практическое удобство его использования в алгоритмах синтеза частотно-временных распределений (ЧВР) несколько затенило теоретический аспект использования аналитического сигнала в качестве инструмента анализа сложных процессов. В связи с этим цель данной работы — в строгом теоретическом обосновании подхода, позволяющего обобщить традиционное представление аналитического сигнала функциями его огибающей и полной фазы (здесь и ниже выражение "функции огибающей и полной фазы" следует понимать как "функции, представляющие огибающую сигнала и полную фазу сигнала") в частотно-временном пространстве.
АНАЛИТИЧЕСКОЕ ПРЕДСТАВЛЕНИЕ СИГНАЛА ФУНКЦИЯМИ ЕГО ОГИБАЮЩЕЙ И ПОЛНОЙ ФАЗЫ
Совершенствование технологий цифровой обработки сигналов обусловило эффективность применения многих алгоритмов частотно-временного анализа, связанных с синтезом совместных представлений плотности распределения энергии в пространстве частота—время, которые в свою очередь базируются на понятии аналитического сигнала za ^) [1]:
Za ^) = Z ^) + ^), (1)
где г(^ и 7 (0 — в общем случае функции, свя-
z(t) = -! ] Mdr; z(t) = ^J ПJ t-T n j
1 } z(T)
dT.
n j t-t
(2)
Согласно (2), основное различие сигналов г) и ) состоит в том, что их спектральные составляющие сдвинуты относительно друг друга на п /2. В [1-3] подробно изучены свойства аналитического сигнала и рассмотрены особенности формирования на его основе ЧВР класса Коэна.
Представление сигнала в аналитической форме в первую очередь связано с удобством его описания. Так, если рассматривать, опираясь на понятие (1), в общем случае произвольный сигнал в виде квазигармонического колебания
za (t) = A(t )exp{ jQ(t)},
(3)
то его огибающую (рис. 1) можно представить как функцию
Л^) = 7^(07740 ,
а полную фазу рассчитать в соответствии с
)
6(t) = arctg
z (t)
При таком определении Л^) и в^) следует выполнение условий:
дг^) ал(о
dt dt d0(t)_ z(t)z(t) - z'(t)z(t)
dt
z 2(t) + z 2(t)
(4)
(5)
т. е. сигнал будет совпадать с функцией огибающей в точках, где в(^.) = 0, 2п, 4п,...
z(t)
--Д......A(() л--.. t
f (t) =
1 d0(t)
(6)
точек ti будет выполняться равенство dz(t) _ d[A(t)]
dt
dt
t = t,.
Рис. 1. Сигнал z (t) и функция его огибающей A(t)
Физический смысл указанных понятий заключается в том, что огибающая представляет функцию, связывающую совокупность точек сигнала, в которых фаза меняет знак. В свою очередь, фаза сигнала в^) представляет функцию, характеризующую закон изменения сигнала в точках его экстремума. Следовательно, мгновенную частоту сигнала (3) можно рассматривать как величину, пропорциональную скорости изменения фазы
2п dt
Из предложенного подхода следует, что
z (t) = A(t)cos(0(t)), а z(t) = A(t )sin(0(t)). (7)
Докажем выполнение условия (7) для равенств (4) и (5).
Утверждение 1
Для колебания вида z(t) = A(t)cos(0(t)), где A(t) — его огибающая, а Q(t) — фаза, справедливо выполнение условия: производная от функции сигнала в общем случае равна значению производной от функции его огибающей.
Доказательство
Распишем значение производной от функции исходного сигнала:
dz (t) = d[ A(t )cos(9(t))] = dt dt
d[A(t)] m, ... „ 4d[0(t)] . m/vi
= v n cos[0(t)] - A(t) —-———isin[0(t)]. dt dt
Учитывая, что функции cos и sin ортогональны, следовательно, всегда существуют такие точки tf, в которых cos[0(t,.)] = 1, а соответственно sin[0(t,)] = 0 . Отсюда следует, что для указанных
С учетом периодичности функции cos полученное равенство будет выполняться при в(и) = 0,2п,4п,...
Утверждение 2
Для колебания вида z(t) = A(t)cos(0(t)) (где A(t) — функция огибающей z(t), а 0(t) — полная фаза z(t)) и его комплексно-сопряженной формы Z(t) = A(t)sin(0(t)) справедливо выполнение следующего условия: производная от функции полной фазы равна отношению
z(t) z(t) - z (t) Z(t) z 2(t) + z 2(t) .
Доказательство
Распишем указанное условие с учетом аналитического описания исходного сигнала и его комплексно-сопряженной формы:
z (t) z(t) - z' (t) z(t) = z 2(t) + z 2(t) = = {A' (t)sin[0(t)] + A(t)0'(t)cos[0(t)]}A(t)cos[0(t)] -= A2(t )cos2[0(t)] + A2(t )sin2[0(t)]
_ {A' (t )cos[0(t)] - A(t )0 '(t )sin[0(t )]}A(t )sin[0(t)] = A2(t)cos2(0(t)) + A2(t )sin2[0(t)] =
{cos2[0(t)] + sin2[0(t)]}A2(t)0 '(t) .'.
A2(t ){cos2[0(t)] + sin2[0(t)]}
= (t).
Физическая сущность мгновенной частоты и огибающей принципиально отличается от частоты и амплитуды, используемых для представления сигнала с помощью ряда обобщенного описания [2]
z(t) = 2£AkcosI n-t + vk
k=0
или его интегральной формы z (t) = 2 J {A( f )cos[2nft + Ф( f)]}df.
(8)
(9)
Так, например, согласно ряду (8), определяемое колебание представляет совокупность гармоник, частота и амплитуда которых неизменны на всем протяжении оси ^ от , до ^, даже если само рассматриваемое колебание существует только в ограниченном временном интервале. Отсутствие же энергии вне сигнального интервала в этом случае
108
С. В. ДВОРНИКОВ
объясняется тем, что спектральные компоненты базисных гармоник в совокупности равны нулю.
С позиций описания сигнала на основе его мгновенной частоты и функции огибающей в каждый момент времени он будет иметь только одну частотную составляющую, определяемую выражением (6), в то время как его спектр может иметь множество различных спектральных компонент. Но в отличие от неизменных величин, определяющих сигнал — частот и амплитуд базисных функций, значение мгновенной частоты может непрерывно меняться во времени.
ДВУМЕРНЫЕ ПРЕДСТАВЛЕНИЯ ФУНКЦИИ ОГИБАЮЩЕЙ АНАЛИТИЧЕСКОГО СИГНАЛА
Развитие теории синтеза ЧВР однозначно связано с понятием комплексной плотности распре -
деления энергии, получаемой в результате комплексного произведения аналитического сигнала га ^) и его спектра ^ (/)
рк (/, 0 = г: (/у2п/. (10)
Распределение вида (10) называют ЧВР Риха-чека [4]. Так, на рис. 2 представлено идеальное ЧВР тестового сигнала, а на рис. 3 — распределение тестового сигнала, полученное на основе выражения (10).
Анализ результатов, представленных на рис. 3 показывает, что помимо истинных компонентов энергии тестового сигнала, расположенных вдоль линии его мгновенной частоты, на плоскости ЧВР присутствуют и так называемые интерференционные компоненты, представляющие результат совместного взаимодействия сигнальных отсчетов, составляющих анализируемый процесс в полосе анализа.
С целью подавления "паразитных" выбросов сигнальной энергии в [5] предложено использовать так называемые псевдораспределения, в которых для расчета "единичной" частотно-временной площадки в каждый момент времени задействуются только сигнальные компоненты, попадающие во временной интервал, ограниченный функцией окна, например Хемминга. Общность такого подхода базируется на рассмотрении псевдораспределения как некоторой сглаженной версии исходного комплексного ЧВР. На практике большее применение нашло распределение функции активной энергии, представляющее реальную часть ЧВР Рихачека [4]
Рл (/, t) = 0.5Яе[рд (/, t)]. (11)
0.8
0.6
0.4
0.2
1 / [кГц]
1
1
\
V
^ [мс]
0 1 2 3 4 5 6 7
Рис. 2. Идеальное частотно-временное распределение тестового сигнала
t, [мс]
1---------Д
0 1 2 3 4 5 6 7
Рис. 3. Частотно-временное распределение Риха-чека тестового сигнала
Так, на рис. 4 показан результат представления тестового сигнала на основе его псевдоформы.
Поскольку распределение энергии (11) в общем случае характеризуется мгновенным значением, то, очевидно, к нему применимо понятие огибающей. Однако в этом случае она будет представлена функцией двух переменных — частоты и времени. Тогда в соответствии с утверждением 1 аналитическое описание огибающей функции распределения энергии будет иметь следующий вид:
Ом(/, X) = ^/{Яе[(/, X)] +{1т[(/, X)] . (12)
Следует заметить, что полученные аналитические выражения для ЧВР имеют место только в отношении узкополосных сигналов, значение нижней границы спектра /н которых значительно превышает занимаемую ими полосу частот:
/h >AF .
(13)
В [3] уточнено само понятие узкополосного сигнала. Вместо традиционно принятого неравенства АР << /0 допускается
/0 > 1.5AF:
(14)
dt
dt
dt
Доказательство
Представим производную функции плотности распределения мгновенной энергии в следующем виде:
Р (/,х) = 2^ ОО8[0(Х)] X dt dt х[а,(/)со8(2л/) + Ъ(/)8ш(2п/Х)] -
-2 A(t )
d0(t )
dt
[at ( f ) cos(2nf ) - bt ( f ) sin(2n ft )] -
-4n ft • A(t )cos0(t ) x
x[at ( f ) sm(2nf ) - bt ( f ) cos(2nft )].
Функции базиса разложения cos и sin ортогональны, следовательно, всегда существуют такие точки tf , в которых функция cos[0(t¿)] = 1, а соответственно sin[0(t¿ )] = 0 . Тогда
dpM ( f, t ) =
= 2
dt
dA(t ) dt
[at ( f )cos(2nf ) + bt ( f)sm(2n ft )] -
-4л/X • А(Х)[а(/)8Ш(2Л/) - Ъ(/)С08(2я/X)].
При отсутствии частотных сдвигов вторая половина правой части равенства будет равняться нулю, в результате чего окончательно получим:
dPM ( f, t ) „dA(t )
dt
= 2-
dt
-X( f, t ).
где /0 — средняя частота в спектре сигнала. Корректность использования функции огибающей распределения в качестве адекватного представления плотности энергии в частотно-временном пространстве основывается на доказательстве следующего утверждения.
Утверждение 3
Для колебания вида г(X) = А(Х)со8[0(Х)] , где А(Х) и 0(Х) соответственно функции огибающей и фазы, значение производной его функции плотности распределения мгновенной энергии вида рм (/, X) = 2г(Х)Х(/, X) равно представлению производной огибающей колебания в двумерном тригонометрическом частотно-временном базисе разложения Х(/, X) = а (/) С08(2п/X) + Ъ (/) 8т(2тс/X)]:
dPм (/, X) = d[2г(X)Х(/, X)] = 2^4(0 Х(/ t)
Таким образом, в соответствии с утверждением (3) производная от функции распределения мгновенной энергии будет равна с точностью до постоянного множителя представлению производной от огибающей сигнала в двумерном тригонометрическом частотно-временном базисе разложения.
Особая привлекательность функции плотности распределения мгновенной энергии как инструмента анализа заключается в несмещенности оценок, получаемых на основе ее производной, в соответствии с утверждением 4.
Утверждение 4
Мгновенные оценки, получаемые на основе функции плотности распределения мгновенной энергии сигнала, являются несмещенными.
Доказательство
Допустим, что Рм (/, X) = 2гЦ)Х(/, X) — функция плотности распределения мгновенной энергии некоторого произвольного сигнала г (X), где
Х(/, X) = [а (/)С08(2я/) + Ъ (/)8Ш(2я/)] — двумерный тригонометрический частотно-временной базис разложения. Тогда
dРм (/, X)_^2г(t) Х(/, X)]_
dt
dt
= 2 •
*(t) X(f, t) + z(t) dX(f,t)
dt
dt
Решая данное уравнение, получим:
dPM ( f, t ) = dt
110
С. В. ДВОРНИКОВ
= 2
д£(о
[а,(/)ОС8(2п/0 + Ъ,(/)Яп(2п/0] -
-4п/• )[а,(/)*ш(2пА)-Ъ(/)ос8(2п/)].
Для мгновенных оценок вторая половина правой части равенства может быть приравнена нулю, причем без всякого смещения по частоте. В результате чего будем иметь
Фм (/,,) = Х(/,,).
ди
ди
Таким образом, оценка, полученная на основе производной от функции распределения мгновенной энергии, является несмещенной и численно равна произведению производной сигнала и функ-
ции тригонометрического частотно-временного базиса.
Проведенный в [3] анализ свойств распределения огибающей (12) показал, что приемлемые результаты его применения в качестве инструмента анализа удается получить только для узкополосных сигналов с низкой степенью нестационарности поведения их параметров. Указанное обстоятельство обусловлено достаточно слабо выраженной степенью концентрации энергии в пределах обрабатываемой частотно-временной области. Так, на рис. 5 представлено распределение Ом (/,,) тестового сигнала.
Анализ полученных результатов позволяет сделать вывод, что в пределах от 1 мс до 7 мс (область относительно высокой стационарности) распределение (12) локализует энергию сигнала вдоль линии его мгновенной частоты с высокой степенью концентрации. В то же время на интервале от нуля до 1 мс (область относительно низкой стационарности) степень концентрации сигнальной энергии такова, что сложно говорить о характере изменения мгновенной частоты.
Справедливости ради необходимо отметить, что оценивание энергетических параметров нестационарных процессов непосредственно на основе распределения мгновенной энергии
Рм (/,,) = 27(0К (/)ОС8(2П/0 + Ъ (/)81и(2п/0]
не лишено сложностей. Прежде всего — это наличие "паразитных" интерференционных компонент. Так, на рис. 6 представлено распределение рм (/,,) тестового сигнала. Даже на основе визуального анализа можно заключить, что в районе частоты / = 200 Гц на плоскости ЧВР отчетливо просматривается ложный пик энергии, обусловленный тем, что в пределах "единичной" площадки взаимодействовали не только близлежащие сигнальные компоненты, но и находящиеся на более отдаленных частотных позициях.
ЗАКЛЮЧЕНИЕ
Несмотря на указанные недостатки рассмотренных распределений, последние наряду с другими совместными частотно-временными описаниями класса Коэна находят широкое применение в практических приложениях радиотехники [6-7]. Рассмотренные же в работе ЧВР были выбраны с целью наглядного представления полученных теоретических результатов, поскольку они достаточно полно характеризуют роль и значение функции огибающей и полной фазы аналитической формы представления сигналов для измерения и оценивания их параметров.
СПИСОК ЛИТЕРАТУРЫ
1. Левин Б.Р. Теоретические основы статистической радиотехники. Кн. первая. Изд. 2-е, пере-раб. и доп., М.: Сов. радио, 1974. 552 с.
2. Гоноровский И. С. Радиотехнические цепи и сигналы: Учебник для вузов. 4-е изд., перераб. и доп. М.: Радио и связь, 1986. 512 с.
3. Алексеев А.А., Кириллов А.Б. Технический анализ сигналов и распознавание радиоизлучений. СПб.: ВАС, 1998. 368 с.
4. Коэн Л. Время-частотные распределения: Обзор // ТИИЭР. 1989. Т. 77, № 10. С. 72-121.
5. Claasen T.A.C.M., Meclenbrauker W.F.G. The Wigner distribution a tool for time-friquency signal analysis. Part 1, 2, 3 // Philips J. Res. 1980. V. 35. P.217-250, 276-300, 372-389.
6. Hlawatsch F., Papandreou-Suppappola A., Boun-dreaux-Bartels G.F. The power classes — quadratic time-frequency representations with scale covariance and dispersive time-shift covariance // IEEE Transactions on Signal Processing. 1999. V. 47, N 11. P. 3067-3083.
7. Дворников С.В., Железняк В.К., Комарович В.Ф. и др. Метод обнаружения радиосигналов на основе обработки их частотно-временных распределений плотности энергии // Информация и космос. 2005. № 4. С. 13-16.
Военная академия связи, Санкт-Петербург Материал поступил в редакцию 13.03.2006.
THEORY OF ANALYTIC SIGNAL PRESENTATION BY FUNCTIONS OF SIGNAL ENVELOPE AND TOTAL PHASE
S. V. Dvornikov
Military Communication University, Saint-Petersburg
Theoretical aspects of analytic signal presentation by functions of its envelope and total phase on time-frequency and time scales are offered. The results of theoretical researche are presented. This article may be interesting for specialists involved in time-frequency distribution synthesis applied to estimation of narrowband non-stationary process parameters in the presence of high intensity noise.