АНАЛИТИЧЕСКИЕ СИГНАЛЫ

Нуль и единица - от Бога, все остальное - дело рук человеческих.

Леопольд Кронекер. Немецкий математик, XIX в.

Содержание: 10.1. Понятие аналитического сигнала. Комплексное представление вещественных сигналов. Аналитический сигнал. Спектральная плотность аналитического сигнала. 10.2. Примеры применения аналитических сигналов. Огибающая и мгновенная фаза сигналов. Мгновенная частота. Огибающие модулированных сигналов. Анализ каузальных систем. Литература.

Введение

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

10.1. Понятие аналитического сигнала [1,25].

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

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

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

,, при этом .

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

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

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

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

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

z_s(t) = . (10.1.2)

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

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

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

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

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

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

При сложении функций z_s(t) и z_s*(t) с учетом нормировки в (10.1.2) только на 1/π, а не на 1/2π, как в (10.1.1), мы обязаны получить полный исходный сигнал s(t):

s(t) = [z_s(t)+z_s*(t)]/2 = Re z(t).

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

Реальная и мнимая части спектра произвольных каузальных сигналов связаны преобразованием Гильберта. Оно позволяет производить определение любой части частотной характеристики каузальной функции, действительной или мнимой, путем свертки другой ее части с оператором Гильберта 1/pf. Аналогично, мнимая часть аналитического сигнала z_s(t) является аналитически сопряженной с его действительной частью Re z(t) = s(t) через преобразование Гильберта, и называется квадратурным дополнением сигнала s(t):

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

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

z_s(t) = s(t) + jЧ . (10.1.4)

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

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

= , (10.1.3')

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

Почему именно оператор Гильберта применяется для получения квадратурного дополнения сигнала? Какую физическую операцию он выполняет? Ответ на этот вопрос может быть получен при рассмотрении спектра аналитического сигнала.

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

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

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

Z_s(w) = (10.1.5)

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

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

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

Рис. 10.1.3.

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

sgn(w) = (10.1.7)

Однако при этом реальная часть новой функции sgn(w)·S(w), как это можно видеть на рис. 10.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 (для выражения в (10.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),

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

= = -jЧ sgn(w)Ч S(w), (10.1.8)

Hb(w) = -jЧ sgn(w) = (10.1.9)

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

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

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

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

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

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

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

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

Hb(w) = -jЧ sgn(w) = , (10.1.12)

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

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

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

= |S(w)|Ч exp(jj_s(w))Ч exp(jj_h(w)) = |S(w)|Ч exp[j(j_s(w)+j_h(w))], (10.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). Знак минус демонстрирует отставание (запаздывание) выходного сигнала преобразования, как операции свертки, от входного сигнала. Для гармонических сигналов любой частоты с любой начальной фазой это запаздывание составляет четверть периода колебаний. На рис. 10.1.5 этот сдвиг на четверть периода для единичной гармонической составляющей (несущей частоты радиоимпульса) виден достаточно наглядно.

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

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

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

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

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

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

x(t) = a(t)Ч cos(w_ot) + b(t)Ч sin(w_ot), (10.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) получим из выражения (10.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) =. (10.2.3)

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

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

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

df(t)/dt = . (10.2.5)

Рис. 10.2.1.

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

На рис. 10.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).

Рис. 10.2.2.

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

u(t) =,

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

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

f(t) = .

Рис. 10.2.3. Рис. 10.2.4.

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

w(t) = .

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

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

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

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

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

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

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

u(t) = |z_x(t)| = U_oЧ[1+mЧ s(t)],

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

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

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

x(t) = U_oЧs(t)Ч cos w_ot,

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

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

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

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

Анализ каузальных систем. Каузальная (физически осуществимая) линейная система задается односторонним импульсным откликом 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). (10.2.6)

Осуществляя обратное преобразование Фурье обеих частей равенства (10.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 с.

к оглавлению

Знаете ли Вы, что "тёмная материя" - такая же фикция, как черная кошка в темной комнате. Это не физическая реальность, но фокус, подмена.
Реально идет речь о том, что релятивистские формулы не соответствуют астрономическим наблюдениям, давая на порядок и более меньшую массу и меньшую энергию. Отсюда сделан фокуснический вывод, что есть "темная материя" и "темная энергия", но не вывод, что релятивистские формулы не соответствуют реалиям. Подробнее читайте в FAQ по эфирной физике.

НОВОСТИ ФОРУМА

Рыцари теории эфира