Аналитические сигналы

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

Аналитические сигналы

Введение

аналитический сигнал вещественный казуальный

Аналитический сигнал - это один из способов комплексного представления сигнала, который применяется при анализе сигналов и систем их обработки. Он позволяет ввести в анализ понятия огибающей и мгновенной частоты сигнала.

1. Понятие аналитического сигнала

Комплексное представление вещественных сигналов. При математическом анализе вместо вещественных сигналов с целью упрощения математического аппарата преобразования данных удобно использовать эквивалентное комплексное представление сигналов. Так, например, в теории электрических цепей вещественная запись синусоидального напряжения

u(t) = U_o cos (w_ot+j)

заменяется комплексной формой записи:

, , при этом .

Произвольный сигнал s(t), заданный на определенном участке временной оси (как конечном, так и бесконечном), имеет комплексную двустороннюю спектральную плотность S(w). При обратном преобразовании Фурье раздельно реальной и мнимой части спектра S(w) сигнал s(t) разделяется на четную и нечетную составляющие, которые являются двусторонними относительно t = 0, и суммирование которых полностью восстанавливает исходный сигнал. На рис. 1.1 приведен пример сигнала (А), его комплексного спектра (В) и получения четной и нечетной части сигнала из реальной и мнимой части спектра (С).

Рис. 1.1. Сигнал, спектральная плотность сигнала, четная и нечетная составляющие.

Аналитический сигнал. Можно выполнить обратное преобразование Фурье и в другой форме - раздельно для положительных и отрицательных частот спектра:

s(t) =S(w)·exp(jwt) dw + S(w)·exp(jwt) dw.----------(1.1)

Информация в комплексном спектре сигнала является избыточной. В силу комплексной сопряженности полную информацию о сигнале s(t) содержит как левая (отрицательные частоты), так и правая (положительные частоты) часть спектра S(w). Аналитическим сигналом, отображающим вещественный сигнал s(t), называют второй интеграл выражения (1.1), нормированный на p, т.е. обратное преобразование Фурье спектра сигнала s(t) только по положительным частотам:

z_s(t) = (1/p)S(w) exp(jwt). (1.2)

Дуальность свойств преобразования Фурье определяет, что аналитический сигнал z_s(t), полученный из односторонней спектральной функции, является комплексным и может быть представлен в виде:

z_s(t) = Re z(t) + j·Im z(t). (1.2')

Рис. 1.2. Сигналы z(t) и z*(t).

Аналогичное преобразование первого интеграла выражения (1.1) дает сигнал z_s*(t), комплексно сопряженный с сигналом z(t):

z_s*(t) = Re z(t) - j·Im z(t),

что наглядно видно на рис. 1.2 при восстановлении сигналов по односторонним частям спектра, приведенного на рис. 1.1-В.

По рисунку 1.2 можно видеть, что при сложении функций z_s(t) и z_s*(t) мнимые части функций взаимно компенсируются, а вещественные части, с учетом нормировки только на р, а не на 2р, как в (1.1), в сумме дают полный исходный сигнал s(t):

[z_s(t)+z_s*(t)]/2 = Re z(t) = = (1/2p)S(w) cos wt dt = s(t).

Отсюда, реальная часть аналитического сигнала z_s(t) равна самому сигналу s(t).

Для выявления характера мнимой части сигнала z_s(t) выполним перевод всех членов функции (1.2') в спектральную область с раздельным представлением по положительным и отрицательным частотам (индексами - и +) реальных и мнимых частей спектра:

Z_s(w) = A_-(w) + A₊(w) + jB_-(w) + jB₊(w) + j [A'_-(w) + A'₊(w) + jB'_-(w) + jB'₊(w)],

где индексами A' и B' обозначены функции преобразования Im(z(t)). В этом выражении функции в левой части спектра (по отрицательным частотам) должны взаимно компенсировать друг друга согласно определению аналитического сигнала (1.2), т.е.:

B'_-(w) = A_-(w), A'_-(w) = -B_-(w).

Отсюда, с учетом четности вещественных A'_-(w) и нечетности мнимых B'_-(w) функций спектра, следуют равенства:

B'₊(w) = - A₊(w), A'₊(w) = B₊(w).

Но эти четыре равенства есть не что иное, как преобразование Гильберта в частотной области спектра функции Re z(t) A(w)+jB(w) в спектр функции A'(w)+jB'(w) Im z(t) умножением на сигнатурную функцию -jsgn(w). Следовательно, мнимая часть аналитического сигнала z_s(t) является аналитически сопряженной с его действительной частью Re z(t) = s(t) через преобразование Гильберта. Эта часть аналитического сигнала получила название квадратурного дополнения сигнала s(t):

Im z(t) = = TH[s(t)] = s(t) _* hb(t), (1.3)

hb(t) = 1/(рt),

z_s(t) = s(t) + j. (1.4)

где индексом обозначен сигнал, аналитически сопряженный с сигналом s(t), hb(t) - оператор Гильберта.

Таким образом, квадратурное дополнение сигнала s(t) представляет собой свертку сигнала s(t) с оператором 1/(рt) и может быть выполнено линейной системой с постоянными параметрами:

= , (1.3')

Аналитический сигнал зависит от действительного аргумента, является однозначным и дифференцируемым. На комплексной плоскости он отображается вектором, модуль и фазовый угол которого изменяются от аргумента, а проекция сигнала на вещественную ось равна значению исходного сигнала s(t).

Какой-либо новой информации аналитический сигнал не несет, так как получен линейным преобразованием из исходного сигнала и представляет собой его новую математическую модель.

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

Z_s(w) = z_s(t) exp(-jwt) dt.

Эта функция, с учетом определения аналитического сигнала по выражению (1.2), должна быть отлична от нуля только в области положительных частот, где ее значения (в силу нормировки на p, а не на 2p) должны быть равны удвоенным значениям спектральной плотности сигнала s(t):

Z_s(w) = (1.5)

С другой стороны, при непосредственном преобразовании Фурье левой и правой части формулы (1.4) аналитического сигнала z_s(t), получаем:

Z_s(w) = S(w) + j. (1.6)

Данное выражение действительно для всей частотной оси (от - до +) и должно быть равно выражению (1.5). А это означает, что левая часть спектра (1.6) (отрицательные частоты w) должна быть обращена в ноль, аналогично формированию каузальной функции из ее четной и нечетной части. Это может быть выполнено следующим образом.

Рис. 1.3.

Если левые части спектра сигнала S(w) умножить на -1, обнулить реальную часть на частоте w=0 и оставить без изменения правые части спектра, то будут получены функции, показанные пунктиром на рис. 1.3), которые дают нули в левой части спектра при сложении с исходной функцией S(w) и увеличивают в 2 раза правые части спектра. Такая операция может быть выполнена умножением спектра S(w) на сигнатурную функцию sgn(w):

sgn(w) = (1.7)

Однако при этом реальная часть новой функции sgn(w)·S(w), как это можно видеть на рис. 1.3, становится нечетной, а мнимая часть четной, что не соответствует статусу спектральных функций. Для восстановления статуса полученный результат нужно дополнительно умножить на -j. Применяя для левой и правой части частотных аргументов индексирование соответственно w_l--и--w_r, можно записать подробные выражения для спектров:

S(w) = Re S(w_l) + j·Im(w_l) + Re S(w_r) + j·Im(w_r),

= j·Re S(w_l) - Im(w_l) - j·Re S(w_r) + Im(w_r).

При умножении квадратурной функции на j (для выражения в (1.6)):

j·= -Re S(w_l) - j·Im(w_l) + Re S(w_r) + j·Im(w_r).

Отсюда нетрудно видеть результат:

Z_s(w) = S(w) + j = = 2·Re S(w_r) + j·2·Im(w_r) = 2·S(w_r),

что полностью соответствует выражению (1.5). В краткой форме:

= = -jsgn(w)S(w), (1.8)

Hb(w) = -jsgn(w) = (1.9)

Таким образом, спектральная плотность аналитически сопряженного сигнала образуется из спектра S(w) исходного сигнала s(t) умножением на функцию -jsgn(w). Это обеспечивает при суммировании S(w) + j удвоение амплитуд частотных составляющих в области положительных частот и их взаимную компенсацию в области отрицательных частот.

Рис. 1.4.

Из выражения (1.8) в спектральной области непосредственно следует соответствующая связь функций s(t) и во временной области:

= s(t) _* hb(t), (1.10)

s(t) = -_* hb(t). (1.11)

где hb(t) = TF[-jsgn(w)] = 1/(pt) - обратное преобразование Фурье функции -jsgn(w).

Пример преобразования сигнала x(t) оператором Гильберта для формирования аналитического сигнала z_x(t) = x(t) + j· приведен на рис. 1.4.

Частотную характеристику оператора Гильберта (1.9) можно записать и в следующем виде:

Hb(w) = |Hb(w)|exp(jj_h(w)), где |Hb(w)| = 1.

Hb(w) = -jsgn(w) = , (1.12)

Если спектр функции x(t) также представить в форме

S(w) = |S(w)|exp(jj_s(w)),

то выражение (1.8) преобразуется к следующей форме:

= |S(w)|exp(jj_s(w))exp(jj_h(w)) = |S(w)|exp[j(j_s(w)+j_h(w))], (1.8')

т.е. модуль |S(w)| - амплитудный спектр сигнала как результат преобразования Гильберта сигнала s(t), не изменяется и остается равным амплитудному спектру сигнала s(t). Фазовый спектр сигнала (начальные фазовые углы всех гармонических составляющих сигнала) сдвигается на -90^о при w > 0 и на 90^о при w < 0 относительно фазового спектра сигнала s(t). Но такой фазовый сдвиг означает не что иное, как превращение косинусных гармоник в синусные, а синусных в косинусные. Это можно наглядно видеть на единичной гармонике. Так, если x(t) = cos(2pf_ot), то имеем следующее преобразование Гильберта через частотную область:

(t) = TH[x(t)] TF[TH[x(t)]] = -j sgn(f)[d(f+f_o)+d(f-f_o)]/2.

(f) = -j[-d(f+f_o)+d(f-f_o)]/2 = j·[d(f+f_o)-d(f-f_o)]/2.

Но последнее уравнение - спектр синусоиды. При обратном преобразовании Фурье:

(t) = TF^-1[(f)] = sin(2pf_ot).

При x(t) = sin(2pf_ot) аналогичная операция дает (t) = -cos(2pf_ot). Знак минус демонстрирует отставание (запаздывание) выходного сигнала преобразования, как операции свертки, от входного сигнала. Для гармонических сигналов любой частоты с любой начальной фазой это запаздывание составляет четверть периода колебаний. На рис. 1.5 этот сдвиг на четверть периода для единичной гармонической составляющей (несущей частоты радиоимпульса) виден достаточно наглядно.

Таким образом, аналитический сигнал, по существу, представляет собой сумму двух ортогональных сигналов, все гармонические составляющие которых сдвинуты по фазе на 90⁰ друг относительно друга.

2. Примеры применения аналитических сигналов [1,2].

Огибающая и мгновенная фаза сигналов. Допустим, что имеем зарегистрированный радиоимпульсный сигнал x(t) с несущей частотой ?_o, который содержит определенную информацию, заключенную в огибающей сигнала u(t) и его фазе j(t):

x(t) = u(t) cos (w_ot+j(t)). (2.1)

Требуется выделить информационные составляющие сигнала

Запишем выражение (2.1) в другой форме:

x(t) = a(t)cos(w_ot) + b(t)sin(w_ot), (2.2)

где функции a(t) и b(t) называются низкочастотными квадратурными составляющими сигнала x(t):

a(t) = u(t) cos jt, b(t) = u(t) sin jt.

u(t) =, tg j(t) = b(t)/a(t).

С использованием преобразования Гильберта из сигнала x(t) можно сформировать аналитически сопряженный сигнал (t). Математическую форму сигнала (t) получим из выражения (2.2) с учетом свойства модуляции преобразования Гильберта:

(t) = a(t)sin(w_оt) - b(t)cos(w_ot).

z(t) = x(t) + j(t).

Квадрат модуля сигнала z(t):

|z(t)|² = x²(t)+²(t) = a²(t)[cos²(w_ot)+sin²(w_ot)] + b²(t)[cos²(w_ot)+sin²(w_ot)] = u²(t).

Отсюда, огибающая u(t) и мгновенная фаза f(t) сигнала x(t):

u(t) =. (2.3)

f(t)--=--w_ot+j(t) = arctg[(t)/x(t)]. (2.4)

j(t)--=--f(t)-----w_ot.

Рис. 2.2.

Мгновенная частота сигнала определяется по скорости изменения мгновенной фазы:

df(t)/dt = . (2.5)

Для амплитудно-модулированных сигналов с одной несущей частотой эти результаты достаточно очевидны (см. рис. 1.5). Но выражения (2.3-2.5), полученные из общих соображений, остаются действительными и для любых произвольных сигналов.

На рис. 2.2. представлен сигнал, сложенный двумя гармониками:

x(t) = a(t)cos(w₁t) + b(t)cos(w₂t).

Квадратурное дополнение и аналитический сигнал:

(t) = a(t)sin(w₁t) + b(t)sin(w₁t).

z(t) = x(t) + j(t).

Огибающая такого сигнала, как это можно видеть на рисунке 2.2, должна вычисляться по формуле (2.3). При этом для данного сигнала получаем:

u(t) =,

что может существенно отличаться от функции .

Мгновенная фаза сигнала, график которой приведен на рис. 2.3, зависит от времени нелинейно:

f(t) = .

Рис. 2.3. Рис. 2.4.

Мгновенная частота сигнала (рис. 2.4) также имеет нелинейную зависимость от времени, причем ее значения могут существенно превышать даже суммарное значение частот, составляющих сигнал:

w(t) = .

Аналогичная методика определения огибающих, мгновенных значений фазы и частоты применяется и для анализа случайных процессов.

Огибающие модулированных сигналов. В качестве примера применения огибающих рассмотрим связь форм относительно узкополосных радиосигналов с формой модулирующих сообщений.

Амплитудная модуляция. Уравнение модулированного сигнала:

x(t) = U_o[1+ms(t)]cos w_ot, s(t) 1, m 1

Квадратурное дополнение и аналитический сигнал:

(t) = U_o[1+ms(t)]sin w_ot, z_x(t) = x(t) + j(t).

Огибающая сигнала x(t):

u(t) = |z_x(t)| = U_o[1+ms(t)],

т.е. точно повторяет форму модулирующего сообщения (см. рис. 2.5)

Рис. 2.5. Амплитудная модуляция.

Балансная модуляция. Уравнение модулированного сигнала, приведенного на рис. 2.6:

Рис. 2.6. Балансная модуляция.

x(t) = U_os(t)cos w_ot,

Квадратурное дополнение, аналитический сигнал, огибающая сигнала x(t):

(t) = U_os(t)sin w_ot, z_x(t) = x(t) + j(t), u(t) = |z_x(t)| = U_o|s(t)|.

Огибающая сигнала x(t) существенно отличается от модулирующего сообщения, но связана с ним простым соотношением.

Анализ каузальных систем. Каузальная (физически осуществимая) линейная система задается односторонним импульсным откликом h(t), t 0, и имеет частотную характеристику H(f):

H(f) = X(f) - jY(f),

Осуществим обратное преобразование Фурье для всех частей выражения раздельно:

h(t) = x(t) + y(t),

x(t) =X(f) cos(2pft) df,

y(t) =Y(f) sin(2pft) df,

где x(t) и y(t) - четная и нечетная части функции h(t). Нечетная функция y(t) в каузальной системе однозначно связана с четной функцией x(t):

y(t) = sgn(t)x(t). (2.6)

Осуществляя обратное преобразование Фурье обеих частей равенства (2.6) при известном преобразовании сигнатурной функции (sgn(t) -j/(pf)), получаем:

TF[y(t)] = (-j/pf) _* X(f) = (-j/p)[X(u)/(f-u)] du.

Отсюда:

Y(f) = (1/p)[X(u)/(f-u)] du = ТН[X(f)],

т.е. мнимая часть спектра импульсного отклика каузальной системы (и любой каузальной функции) является преобразованием Гильберта действительной части спектра. Соответственно, уравнение для определения действительной компоненты спектра по мнимой части:

X(f) = -ТН[Y(f)] = -(1/p)[Y(u)/(f-u)] dv.

ЛИТЕРАТУРА

1. Баскаков С.И. Радиотехнические цепи и сигналы Учебник для вузов. - М. Высшая школа, 1988.

2. Бендат Дж., Пирсол А. Прикладной анализ случайных данных. - М.: Мир, 1989. - 540 с.

25. Сергиенко А.Б. Цифровая обработка сигналов. - СПб.: Питер, 2003. - 608 с.

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

...