ЛЕКЦИЯ 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)e− jnω1t dt. |
(3.38) |
|||||||
|
|
||||||||
|
|
|
T t |
|
|||||
|
|
|
|
||||||
|
|
|
|
|
1 |
|
|
|
|
Подставим (3.38) в (3.37), для удобств заменив обозначение переменной интегрирования. Тогда,
∞ |
|
1 t2 |
|
|
sR (t) = ∑ |
|
|
∫s( x)e− jnω1x dx e jnω1t , |
|
|
||||
|
T |
t1 |
|
|
n=−∞ |
|
|
||
Учтем, что T = 2π/ ω1 :
2
|
|
|
1 |
∞ |
t2 |
|
|
|
|
s |
R |
(t) = |
∑ |
|
s( x)e− jnω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)e− jωx dx dω. |
(3.41) |
|||
|
||||||
|
2π −∞ |
t |
|
|
||
|
|
|
1 |
|
|
|
Интеграл в квадратных скобках является функцией только частоты ω, |
||||||
то есть, можно ввести обозначение |
|
|
||||
|
|
t2 |
|
|
|
|
|
S(ω) = ∫s(x)e− jωxdx. |
|
(3.42) |
|||
t1
Поскольку сигнал s(t) по условию отличается от нуля только в интервале (t1, t2 ) , интеграл (3.42) не изменится, если пределы интегрирования за-
менить на бесконечные. Кроме того, вернемся к переменной t. Получим:
|
|
|
|
|
|
|
∞ |
|
|
|
S(ω) = |
∫s(t)e− jωt dt. |
(3.43) |
||||
|
|
|
|
|
|
−∞ |
|
|
Теперь выражение (3.41) можно записать так: |
|
|||||||
|
|
|
|
1 |
∞ |
|
||
|
|
s(t) = |
∫S(ω)e jωt dω. |
(3.44) |
||||
|
|
|
2π |
|||||
|
|
|
|
|
−∞ |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
U |
Формула |
S(ω) = |
|
|
∫s(t)e− jωt dt называется |
формулой прямого |
||
преобразования Фурье (ПФ), |
−∞ |
|
|
|||||
|
|
∞ |
|
|
||||
|
|
|
1 |
|
|
|
|
|
U |
Формула |
s(t) = |
|
|
∫S(ω)e jωt dω называется формулой обратно- |
|||
|
|
|
||||||
|
|
|
2π −∞ |
|
||||
го преобразования Фурье (ПФО).
Функция S(ω) называется спектральной плотностью сигнала s(t).
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) e− j ω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)e− j2πnf1t dt, |
f1 = |
1 |
. |
(3.50) |
||||||
|
|
|
|
||||||||||||
|
|
T |
τ |
|
|
|
|
T |
|
||||||
|
|
|
|
|
|
|
|||||||||
По формуле ПФ пары (3.45) найдем значение спектральной плотности |
|||||||||||||||
импульса на частоте nf1. |
|
|
|
|
|
|
|
|
|
|
|||||
|
|
∞ |
|
(t)e− j2πnf1t dt = ∫ |
sR (t) e− j2π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 Модуль спектральной плотности одиночного импульса и оги-
бающей линейчатого спектра периодической последовательности, полученной повторением этого импульса, совпадают по форме и отличаются только масштабным множителем.