Материал: Лекция№10

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

ЛЕКЦИЯ 2

СПЕКТРАЛЬНЫЙ АНАЛИЗ НЕПЕРИОДИЧЕСКИХ СИГНАЛОВ.

Основная литература.

Волощук Ю.І. Сигнали та процеси у радіотехніці: Підручник для вищих навчальних закладів, том 1 – Харків: «Компанія СМІТ», 2003. с. 125 – 138.

¾ Пробразование Фурье

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

ки, приведенные на рис. 3 11. Здесь s(t)– сигнал, заданный на бесконечном интервале времени, но отличный от нуля на конечном интервале (t1, t2 ) . Выберем произвольный

интервал

Т, который больше интервала

 

(t1, t2 ) , и

получим периодический сигнал

Рисунок 3.11

sR (t) , который повторяет s(t) с периодом Т.

Это периодическое колебание можно записать в виде ряда Фурье:

 

sR (t) =

 

 

Cn e jnω1t ,

(3.37)

где ω1 = 2π/ T ;

 

 

n=−∞

 

 

 

 

 

t

 

 

 

 

1

 

Cn =

2 s(t)ejnω1t dt.

(3.38)

 

 

 

 

 

T t

 

 

 

 

 

 

 

 

 

 

1

 

 

 

Подставим (3.38) в (3.37), для удобств заменив обозначение переменной интегрирования. Тогда,

 

1 t2

 

sR (t) =

 

 

s( x)ejnω1x dx e jnω1t ,

 

 

T

t1

 

n=−∞

 

 

Учтем, что T = 2π/ ω1 :

2

 

 

 

1

t2

 

 

 

s

R

(t) =

 

s( x)ejnω1x dx e jnω1t ω .

(3.39)

 

 

 

 

 

 

1

 

 

 

 

2π n=−∞ t1

 

 

 

Из анализа графиков, приведенных на рис. 3. 11, следует, что для того, чтобы вернуться от сигнала sR (t) к сигналу s(t), необходимо осуществить

предельный переход вида

s(t) = lim sR (t) .

(3.40)

T →∞

 

Из (3.38) получаем, что при T → ∞ коэффициенты ряда Фурье Cn 0 ,

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

стремится к нулю, то есть спектр становится сплошным. Далее, при T → ∞ частота ω1 dω, nω1 →ω, а операция суммирования в (3.39) трансформи-

руется в операцию интегрирования. Таким образом, формула (3.39) принима-

ет вид

 

 

 

 

 

 

1

t2

 

 

s(t) =

e jωt

s(x)ejωx dx dω.

(3.41)

 

 

2π −∞

t

 

 

 

 

 

1

 

 

Интеграл в квадратных скобках является функцией только частоты ω,

то есть, можно ввести обозначение

 

 

 

 

t2

 

 

 

 

S(ω) = s(x)ejωxdx.

 

(3.42)

t1

Поскольку сигнал s(t) по условию отличается от нуля только в интервале (t1, t2 ) , интеграл (3.42) не изменится, если пределы интегрирования за-

менить на бесконечные. Кроме того, вернемся к переменной t. Получим:

 

 

 

 

 

 

 

 

 

 

S(ω) =

s(t)ejωt dt.

(3.43)

 

 

 

 

 

 

−∞

 

Теперь выражение (3.41) можно записать так:

 

 

 

 

 

1

 

 

 

s(t) =

S(ω)e jωt dω.

(3.44)

 

 

 

2π

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

U

Формула

S(ω) =

 

 

s(t)ejωt dt называется

формулой прямого

преобразования Фурье (ПФ),

−∞

 

 

 

 

 

 

 

 

 

1

 

 

 

 

U

Формула

s(t) =

 

 

S(ω)e jωt dω называется формулой обратно-

 

 

 

 

 

 

2π −∞

 

го преобразования Фурье (ПФО).

Функция S(ω) называется спектральной плотностью сигнала s(t).

(3.47)
являет-

3

Если в качестве переменной используется частота f [Гц], пара преобразований Фурье имеет вид

 

 

=

j2πft dt

 

S( f )

s(t)e

 

 

 

−∞

 

 

(3.45)

 

 

 

 

 

j2

πft

 

 

 

S( f )e

df

 

s(t) =

 

 

 

 

 

−∞

 

 

 

Пару преобразований Фурье (3.45) будем обозначать так: s(t) S( f ) или s(t) S(ω) , где двустороння стрелка указывает, на то, что S( f ) или S(ω) получены из s(t) с помощью ПФ, а s(t) получено из S( f ) или S(ω) с помощью ПФО.

Непериодический сигнал s(t) является Фурье преобразуемым, если он удовлетворяет некоторым требованиям – достаточным условиям Дирихле:

Функция s(t) однозначно определена для всех значений аргумен-

та t, имеет конечное число экстремумов и разрывов на любом конечном интервале времени.

Функция s(t) должна быть абсолютно интегрируемой, то есть:

 

| s(t) | dt < ∞.

(3.46)

−∞

Подчеркнем две важных особенности этих условий:

RУсловия Дирихле являются достаточными, но не являются необхо-

димыми, то есть существуют такие функции, для которых эти условия не выполняются, но они являются Фурье-преобразуемыми;

RУсловиям Дирихле удовлетворяют все сигналы с ограниченной энергией, то есть для которых

Tlim→∞ T | x(t) |2 (t)dt K, где К – положительное вещественное число.

В общем случае спектральная плотность S( f ) (или S(ω) ) является

комплексной функцией частоты, то есть:

S( f ) =| S( f ) | e jθ( f ) .

Тогда, если s(t) S( f ) и s(t) является вещественной функцией времени, функция частоты S( f ) должна удовлетворять свойствам симметрии, аналогичным отмеченным ранее для коэффициентов ряда Фурье, а именно:

| S( f ) |=| S(f ) | и θ( f ) = −θ(f )

Модуль спектральной плотности вещественного сигнала s(t)

ся четной функцией частоты, а аргумент – нечетной функцией частоты

4

R В общем случае две функции частоты | S( f ) | и θ( f ) опреде-

ляют сплошной спектр непериодического сигнала s(t).

R Функция частоты | S( f ) | называется амплитудным спектром сиг-

нала s(t).

R Функция частоты θ( f ) называется фазовым спектром сигнала s(t).

Запишем формулы прямого и обратного преобразований Фурье в тригонометрической форме. Воспользуемся формулой Эйлера (А.1) – формула 1 в приложении к первой части учебника, стр. 521. (3.43) примет вид

 

S(ω) =

s(t) ej ωt dt = s(t) cos (ωt) dt j s(t) sin (ωt) dt =A(ω) j B(ω).

 

−∞

−∞

−∞

Здесь A(ω)

– ПФ четной составляющей сигнала, B(ω) – ПФ нечетной со-

ставляющей сигнала (Рис. 3.3). Теперь выражения для амплитудного и фазового спектров сигнала s(t) можно записать так:

| S(ω) |= A2

(ω) + B2

(ω),

 

 

(3.48)

θ(ω) = −arctg B(ω) .

 

A(ω)

 

 

 

Из (3.48) следует такое свойство спектральной плотности:

U Если сигнал s(t) является вещественной и четной функцией

времени, то его спектральная плотность будет вещественной функцией частоты, следовательно, фазовый спектр θ (ω ) ( или θ ( f )) может

принимать только два значения 0 или ± π.

U Если сигнал s(t) является вещественной и нечетной функцией времени, то его спектральная плотность будет мнимой функцией частоты, следовательно, фазовый спектр θ (ω ) ( или θ ( f )) может прини-

мать только два значения +π / 2 или π / 2 .

Приведем к тригонометрической форме ПФО (3.44).

 

1

 

 

1

 

 

s(t) =

 

S(ω)e jωt dω=

| S(ω) | e j[ωt (ω)]dω=

 

 

 

2π−∞

 

2π −∞

 

 

 

 

1

 

 

1

 

=

 

 

| S(ω) | cos[ωt (ω)] dω+ j

 

| S(ω) | sin[ωt (ω)] dω.

 

 

 

 

 

 

 

2π −∞

 

 

2π−∞

Второй интеграл правой части этого выражения равен нулю.

Упражнение 3.7

Докажите справедливость последнего утверждения.

Получим:

s(t) =

1

(3.49)

 

| S(ω) | cos[ωt (ω)] dω.

 

π

 

 

 

0

 

5

Отметим, что в (3.49) уже нет составляющих с отрицательными частотами, то есть мы вернулись к физической частоте.

Переход от комплексно-экспоненциальной формы ПФ к тригоно-

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

Отметим одно замечательное свойство ПФ – свойство площадей.

В (3.43) подставим ω=0. Получим S(0) = s(t)dt. В формулу ПФО под-

 

−∞

ставим t=0 и получим, что s(0) =

S( f ) df . Из геометрической интерпрета-

 

−∞

ции интегрирования и приведенных соотношений следует:

U Если s(t) S( f ) , то на частоте f (или w), равной нулю, значе-

ние спектральной плотности численно равно площади функции s(t).

U Если s(t) S( f ) , то площадь функции S( f ) численно равняется значению сигнала в момент времени t =0.

Найдем связь спектров импульсного сигнала и периодической последовательности этих импульсов.

Пусть имеется один импульс – сигнал s1(t) длительностью τ. Повторим этот импульс с периодом T > τ так, чтобы получить периодический сигнал sR(t). Этот периодический сигнал можно представить в виде ряда Фурье с коэффициентами

Cn =

1

 

 

sR (t)ej2πnf1t dt,

f1 =

1

.

(3.50)

 

 

 

 

 

 

T

τ

 

 

 

 

T

 

 

 

 

 

 

 

 

По формуле ПФ пары (3.45) найдем значение спектральной плотности

импульса на частоте nf1.

 

 

 

 

 

 

 

 

 

 

 

 

 

(t)ej2πnf1t dt =

sR (t) ej2πnf1t dt.

 

S(nf1) = s1

(3.51)

 

 

−∞

 

 

 

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Сравнивая формулы (3.50) и (3.51), получаем:

 

 

 

 

C

n

 

=

1

S(nf ) = f S(nf ).

(3.52)

 

 

 

 

 

 

 

 

 

 

 

 

T

1

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из (3.52) следует, что в случае, когда размерность сигнала, например [B], то размерность спектральной плотности будет [B/Гц] или [B c]. Кроме того, из анализа выражения (3.52) можно сделать такие выводы:

U Модуль спектральной плотности одиночного импульса и оги-

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