Материал: конспект-лекций

Внимание! Если размещение файла нарушает Ваши авторские права, то обязательно сообщите нам

x=at, получаем:

 

 

 

s(at)

s(at)exp(-j

 

 

 

 

101

 

 

 

t) dt = (1/a)

s(x)exp(-jx /a) dx

 

 

 

 

s(at) (1/a) S( /a). (4.3.2')

Выражение (4.3.2') действительно при а>0. При а<0 происходит зеркальный поворот сигна-

ла относительно вертикальной оси, а замена переменной t=x/a вызывает перестановку пределов

интегрирования и, соответственно, изменение знака спектра:

 

s(at)

-(1/a) S(

/a).

(4.3.2'')

Обобщенная формула изменения аргумента:

 

 

s(at)

(1/|a|) S(

/a), a ≠ 0

(4.3.2)

Если под аргументом функции и ее спектра понимать определенные физические единицы,

например, время - частота, то отсюда следует: чем короче по своей длительности сигнал, тем шире по частоте его спектр, и наоборот. Это можно наглядно видеть на рис. 4.3.1. для сигналов s1(k) и s2(k) и их спектров S1( ) и S2( ).

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

ния сигналов и их спектров, но не изменяет самих сигналов и спектров. Так, при масштабе оси времен t=1 секунда, масштаб оси частот f=1/t=1 герц, а при t=1 мксек f=1/t=1 МГц (t=at, f=1/at, a=10-6).

4. Теорема запаздывания. Запаздывание (сдвиг, смещение) сигнала по аргументу функции на интервал to приводит к изменению фазочастотной функции спектра (фазового угла всех гар-

моник) на величину - to без изменения модуля (амплитудной функции) спектра. Применяя заме-

ну переменной t-to = x, получаем:

s(t-to)

s(t-to)exp(-j

t) dt=

s(x)exp(-j x)exp(-j to) dx= S(

)exp(-j to) (4.3.3)

Совершенно очевидно, что амплитуды гармоник сигнала при его сдвиге изменяться не

должны. С учетом того, что |exp(-j to)|=1, это следует и из (4.3.3):

 

 

 

|S(

) exp(-j to)| = |S( )|.

 

 

Фазовый спектр сдвигается на - to с линейной зависимостью от частоты:

 

 

 

S( ) exp(-j to) = R( ) exp[j(

( )] exp(-j to) = R( ) exp[j( ( )- to)].

(4.3.4)

 

 

 

 

 

 

 

 

 

 

Рис. 4.3.3. Изменение спектра сигнала при его сдвиге.

Пример двух одинаковых сигналов, сдвинутых относительно друг друга на to=1, и соответ-

102

ствующих данным сигналам спектров приведен на рис. 4.3.3.

Аналогично нетрудно показать, что сдвиг спектра в частотной области на 0 вызывает умножение сигнала на exp(j 0t):

S( - 0) s(t) exp(j 0t),

что эквивалентно модуляции сигналом комплексной экспоненты во временной области.

5. Преобразование производной (дифференцирование сигнала):

 

 

 

 

 

s(t) = d[y(t)]/dt = d[

Y( ) exp(j

t) d

 

 

 

 

 

 

 

 

 

 

 

 

 

j Y(

) exp(j

 

 

 

 

 

 

 

 

 

 

 

 

 

dt =

Y(

) [d(exp(j

t))/dt] d

 

 

 

 

 

 

t) d

 

j

Y( ).

(4.3.5)

Таким образом, дифференцирование сигнала отображается в спектральной области про-

стым умножением спектра сигнала на оператор дифференцирования сигнала в частотной обла-

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

Рис. 4.3.4. Спектры сигнала и его производной.

Пример сигнала, его производной и соответствующих им спектров приведен на рис. 4.3.4.

По изменению аргумента спектра (для четного исходного сигнала он был нулевым) можно видеть,

что для всех гармоник спектра появляется сдвиг фаз на

/2 (900) для положительных частот, и на -

/2 (-900) для отрицательных частот.

 

 

В общем случае, для кратных производных:

 

 

dn[y(t)]/dtn = (j

n Y( ).

(4.3.6)

При дифференцировании спектра функции соответственно получаем:

 

dn[S( )]/d n = (-jt) n

s(t).

 

6. Преобразование интеграла сигнала в частотной области при известном спектре сигнала может быть получено из следующих простых соображений. Если имеет место s(t) = d[y(t)]/dt j

Y( ) = S( ), то должна выполняться и обратная операция: y(t) = s(t) dt Y(

) = S( )/j . Отсю-

да следует:

 

s(t)dt (1/j )S( ).

(4.3.7)

103

Рис. 4.3.5. Сигналы и амплитудные спектры сигналов.

Оператор интегрирования в частотной области (1/j ) при >1 ослабляет в амплитудном спектре высокие частоты и при <1 усиливает низкие. Фазовый спектр сигнала смещается на -900

для положительных частот и на 900 для отрицательных. Пример модуля спектра сигнала и его ин-

тегральной функции приведены на рис. 4.3.5.

Формула (4.3.7) справедлива для сигналов с нулевой постоянной составляющей. При инте-

грировании сигналов с определенным значением постоянной составляющей С=const в правой ча-

сти выражения (4.3.7) появляется дополнительное слагаемое преобразования Фурье постоянной составляющей C, которое представляет собой дельта-функцию на нулевой частоте с весовым ко-

эффициентом, равным значению С:

 

 

 

 

 

 

 

 

 

(1/j

)S( ) + C·

(0).

 

 

7. Преобразование свертки сигналов y(t) = s(t) * h(t):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y(

) =

y(t) exp(-j t) dt =

 

 

s(

) h(t- ) exp(-j

t) d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Y(

 

) =

s(

) d

 

 

 

h(t-

) exp(-j

t) dt.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По теореме запаздывания (4.3.3):

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

h(t- ) exp(-j

t) dt = H(

) exp(-j

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда:

Y(

) =

 

H(

 

) s(

) exp(-j

 

) d

= H(

)·S(

).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

s(t) * h(t)

 

S(

) H( ).

 

 

dt

(4.3.8)

Рис. 4.3.6. Сигналы и амплитудные спектры сигналов.

Пример выполнения свертки в частотной области приведен на рис. 4.3.6. Отметим, что ча-

стотное представление H( ) импульсного отклика h(t) линейной системы (или соответствующей

104

линейной операции) имеет смысл частотной передаточной функции системы и позволяет опреде-

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

деление по частоте коэффициента пропускания частотных составляющих сигнала с входа на вы-

ход системы (операции).

Таким образом, свертка функций в координатной форме отображается в частотном

представлении произведением фурье-образов этих функций.

Это положение имеет фундаментальное значение в практике обработки данных.

Любая линейная система обработки данных (информационных сигналов) реализует опреде-

ленную операцию трансформации сигнала, т.е. выполняет операцию свертки входного сигнала s(t)

с оператором системы h( ). С использованием преобразования свертки эта операция может про-

изводиться как с динамической, так и с частотной формой представления сигналов. При этом об-

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

ласти, т.к. может быть на несколько порядков выше по производительности, чем во временной об-

ласти. Она представляет собой последовательность следующих операций.

1.

Перевод сигнала в частотную область: s(t) S( ).

 

2.

Умножение спектра сигнала на передаточную функцию системы: Y( ) = H( )·S(

).

Передаточная функция системы определяется аналогичным преобразованием h( ) H(

) или

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

менной области потребуются операторы h( ) с бесконечной импульсной характеристикой. 3. Перевод спектра обработанного сигнала во временную область: Y( ) y(t).

8. Преобразование произведения сигналов y(t) = s(t)·h(t):

Y(

) =

s(t) h(t) exp(-j

t) dt =

s(t) [(1/2

)

H( ') exp(j 't) d '] dt =

= (1/2 )

s(t)H( ') exp(-j(

- ')t) d 'dt = (1/2

)

H(

') d '

 

 

s(t) exp(-j( - ')t) dt =

= (1/2 )

H(

') S( - ') d ' = (1/2 ) H( ) * S( ).

(4.3.9)

Таким образом, произведение функций в координатной форме отображается в частотном представлении сверткой фурье-образов этих функций, с нормировочным множителем (1/2 ),

учитывающем несимметричность прямого и обратного преобразования Фурье функций s(t) и h(t)

при использовании угловых частот.

9. Производная свертки двух функций s'(t) = d[x(t) * y(t)]/dt. С использованием выражений

(4.3.6) и (4.3.8), получаем:

s'(t) = j [X( ) Y( )] = (j X( )) Y( ) = = X( ) (j Y( ). s'(t) = x'(t) * y(t) = x(t) * y'(t).

105

Это выражение позволяет выполнять вычисление производной сигнала с одновременным сглаживанием весовой функцией, которая является производной сглаживающей функции (напри-

мер, гауссиана).

10. Спектры мощности. Временная функция мощности сигнала в общей форме определя-

ется выражением:

w(t) = s(t) s*(t) = |s(t)|2.

Спектральная плотность мощности, соответственно, равна преобразованию Фурье произве-

дения s(t)·s*(t), которое отобразится в спектральном представлении сверткой Фурье-образов этих функций:

W(f) = S(f) * S*(f) = S(f) S*(f-v) dv. (4.3.10)

Но для всех текущих значений частоты f интеграл в правой части этого выражения равен произведению S(f)·S*(f), так как для всех значений сдвига v ≠ 0 в силу ортогональности гармоник

S(f) и S*(f-v) значения их произведения равны нулю. Отсюда:

W(f) = S(f) * S*(f) = |S(f)|2. (4.3.11)

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

нала по спектру мощности невозможно. Это означает также, что сигналы с различными фазовыми характеристиками могут иметь одинаковые спектры мощности. В частности, сдвиг сигнала не от-

ражается на его спектре мощности.

Для функций мощности взаимодействия сигналов в частотной области соответственно име-

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

Wxy(f) = X(f) Y*(f),

Wyx(f) = Y(f) X*(f), Wxy(f) = W*yx(f).

Функции мощности взаимодействия сигналов комплексные, даже если обе функции x(t) и y(t) вещественны, при этом Re[Wxy(f)] - четная функция, а Im[Wxy(f)] - нечетная. Отсюда полная энергия взаимодействия сигналов при интегрировании функций мощности взаимодействия опре-

деляется только реальной частью спектра:

 

 

 

 

 

Exy = (1/2 )

Wxy( ) d

Re[Wxy] d

 

 

 

 

 

 

 

и всегда является вещественным числом.

11. Равенство Парсеваля. Полная энергия спектра сигнала:

Es = W(f) df = |S(f)|2 df.

,

(4.3.12)

Так как координатное и частотное представление по существу только разные математиче-

Источник: https://studfile.net/preview/16436030/