1490
.pdfДля функции ξ (t ), |
интегрируемой на отрезке [−T / 2; T / 2], |
можно за- |
||||||
писать |
|
|
|
|
|
|
|
|
T / 2 ξ |
(t )dt = |
0 |
ξ (t )dt +T / 2ξ (t )dt =T / 2 |
ξ (t )+ξ (−t ) dt , |
(1.24) |
|||
∫ |
|
∫ |
|
∫ |
∫ |
|
|
|
−T / 2 |
|
−T / 2 |
0 |
0 |
|
|
|
|
поскольку при замене переменной t на −t выполняются условия |
|
|||||||
|
|
∫0 |
ξ (t )dt = − ∫0 |
ξ (−t )dt = T∫/ 2ξ (t )dt . |
|
|||
|
−T / 2 |
+T / 2 |
|
0 |
|
|
||
Из выражения (1.24) следует, что |
|
|
|
|
||||
|
|
|
T / 2 |
ξ (t )dt, если функция |
ξ (t ) четная; |
|
||
T∫/ 2 |
ξ (t )dt = 2 ∫0 |
(1.25) |
||||||
−T / 2 |
|
|
|
если функция ξ (t ) нечетная. |
|
|||
|
|
|
0, |
|
В состав основной тригонометрической системы входят как четные функ-
ции (1, cos 2πT k t , …), так и нечетные функции (sin 2Tπ t , …, sin 2πT k t ).
Если функция ξ (t ), заданная на интервале [−T / 2; T / 2], четная, то ее
разложение в ряд Фурье имеет вид
∞ |
|
2πk t , |
|
|
|
ξ (t )= a0 + ∑ak |
cos |
|
|
(1.26) |
|
k =1 |
|
T |
|
|
|
т.к. если функция ξ (t ) четная, то произведение ξ (t )cos |
2πk |
t |
есть четная |
||
|
|||||
|
|
|
T |
|
функция, а произведение ξ (t )sin |
2πk |
t – нечетная функция. В этой связи |
|||||||||||||
|
|||||||||||||||
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
||
|
|
a0 = |
1 |
|
T∫/ 2 ξ (t )dt = |
2 |
T∫/ 2ξ (t )dt; |
|
|||||||
|
|
|
|
|
|
||||||||||
|
|
|
|
T −T / 2 |
T 0 |
|
|
||||||||
ak = |
2 |
T∫/ 2 |
ξ (t )cos |
2πk tdt = |
|
4 |
T∫/ 2 |
ξ (t )cos |
2πk tdt; |
||||||
|
|
||||||||||||||
|
T −T / 2 |
|
|
|
|
|
T |
|
T 0 |
|
T |
||||
|
|
bk |
= |
2 |
T∫/ 2 |
ξ (t )sin 2πk tdt = 0. |
|
||||||||
|
|
|
|
||||||||||||
|
|
|
|
T −T / 2 |
|
|
T |
|
|
22
Если функция ξ (t ) нечетная, то ξ (t )cos |
2πk |
t и ξ (t )sin |
2πk |
t есть не- |
|||||||||||||
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
T |
|
четная и четная функции соответственно. Поэтому |
|
|
|
||||||||||||||
|
|
|
|
a0 = |
|
1 |
T∫/ 2 ξ (t )dt = 0; |
|
|
|
|||||||
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
T −T / 2 |
|
|
|
|
|
|
|
|
|
||
|
|
ak |
= |
2 |
T∫/ 2 |
ξ (t )cos |
2πk tdt = 0; |
|
|
|
|||||||
|
|
|
|
||||||||||||||
|
|
|
|
T −T / 2 |
|
T |
|
|
|
|
|
||||||
bk = |
2 |
T∫/ 2 |
ξ (t )sin |
2πk tdt = |
4 |
T∫/ 2 |
ξ (t )sin |
2πk tdt |
|||||||||
|
|
||||||||||||||||
|
T −T / 2 |
|
|
|
|
T |
|
T 0 |
|
|
T |
|
|
||||
и ряд Фурье для нечетной функции имеет вид |
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
∞ |
2πk t . |
|
|
|
|
|
||||
|
|
|
ξ (t )= ∑bk sin |
|
|
|
(1.27) |
||||||||||
|
|
|
|
|
|
|
k =1 |
T |
|
|
|
|
|
1.4. Ряд Фурье в комплексной форме. Спектр периодической функции
Перепишем для удобства выражение (1.11) для тригонометрического ряда Фурье
|
|
|
|
∞ |
|
|
|
|
∞ |
|
|
|
2πk t |
|
||
|
|
ξ (t )= a0 + ∑ak |
cos 2πk t + ∑bk sin |
|
||||||||||||
|
|
|
|
k =1 |
|
T |
|
|
k =1 |
|
|
|
|
T |
|
|
и, подставив в него известные соотношения |
|
|
|
|
|
|
|
|
|
|||||||
cosα = |
1 |
exp( jα)+ exp(− jα) , |
|
sinα = |
1 |
exp( jα)−exp(− jα) , |
||||||||||
2 |
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
2 j |
|
|
|
|
|||
получим |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ a |
|
|
2πk |
|
|
|
|
|
2πk |
t |
|
− |
||
ξ (t )= a0 + ∑ k |
exp j |
T |
t |
+ exp − j |
|
T |
|
|
||||||||
|
|
k =1 2 |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
∞ b |
|
2πk |
|
|
|
− j |
2πk |
|
= |
|
|
|||
|
− j∑ k exp |
j |
T |
t −exp |
T |
|
t |
|
|
|||||||
|
|
k =1 2 |
|
|
|
|
|
|
|
|
|
|
|
23
∞ 1 |
|
2πk |
|
+ |
1 |
(ak + |
|
− j |
2πk |
|
. (1.28) |
|
= a0 + ∑ |
(ak − jbk )exp j |
T |
|
t |
2 |
jbk )exp |
T |
t |
||||
k =1 2 |
|
|
|
|
|
|
|
|
|
|||
Из выражений (1.21) и (1.22) |
следует, что |
|
|
|
|
|
ak − jbk = 1 2 T
= T1
ak + jbk = 1 2 T
= T1
T∫/ 2
−T / 2
T∫/ 2
−T / 2
T∫/ 2
−T / 2
T∫/ 2
−T / 2
|
2πk |
t − |
j sin |
2πk |
|
|||
ξ (t ) cos |
T |
|
T |
t dt = |
||||
|
|
|
|
|
|
|||
|
− j |
2πk |
|
|
(1.29а) |
|||
ξ (t )exp |
|
T |
|
t dt; |
||||
|
|
|
|
|
|
|
|
|
|
2πk |
t |
+ |
j sin |
2πk |
|
||
ξ (t ) cos |
T |
|
T |
t dt = |
||||
|
|
|
|
|
|
|||
|
j |
2πk |
t |
|
|
(1.29б) |
||
ξ (t )exp |
|
T |
|
dt. |
|
|||
|
|
|
|
|
|
|
|
Рассматривая соотношения (1.29а, б) можно сделать следующие выводы:
1. Выражение 0,5(ak + jbk ) может быть получено из выражения
0,5(ak − jbk ) при изменении знака индекса k : k → −k . Тогда, если ввести обозначение
|
|
|
0,5(a |
k |
− jb )=C& |
, |
|
|
|
||||||
|
|
|
|
|
|
|
|
|
k |
k |
|
|
|
|
|
то можно записать |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,5(a |
+ jb |
)=C& |
|
. |
|
|
||||||
|
|
|
|
|
|
|
k |
|
k |
|
(−k ) |
|
|
|
|
2. Коэффициент a0 |
можно записать в виде |
|
|
||||||||||||
|
a0 |
= |
|
1 |
T∫/ 2 |
ξ (t )exp( j0)dt . |
|
|
|||||||
|
|
|
|
||||||||||||
|
|
|
|
|
T −T / 2 |
|
|
|
|
|
|
|
|||
Этот коэффициент есть частный случай выражения для значений C& |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(±k ) |
C&(±k ) ≡C&k |
= |
|
1 |
T∫/ 2 |
ξ (t )exp |
− j 2πk t |
dt , |
(1.30) |
|||||||
|
|
||||||||||||||
|
|
|
T −T / 2 |
|
|
|
|
|
T |
|
|
в случае k = 0 . Следовательно, выражение (1.28) можно переписать в виде
24
∞ |
exp |
|
2πk t . |
|
ξ (t )= a ∑C&k |
j |
(1.31) |
||
k =−∞ |
|
|
T |
|
Соотношение (1.31) есть разложение периодической функции ξ (t ) в
ряд Фурье, обладающий комплексными коэффициентами C&k , которые оп-
ределяются выражением (1.30), или, иначе, рядом Фурье в комплексной форме. Этот ряд является более простым по структуре и имеет то преимущество, что все коэффициенты разложения определяются единственной формулой (1.30). Рассматривая физический смысл ряда (1.31), нетрудно видеть, что в разложении присутствуют как положительные гармоники с индексом k , так и отрицательные с индексом −k . Хотя отрицательные частоты физически не существуют, их использование создает определенные удобства для анализа, поскольку операции с комплексными числами значительно проще. Нетрудно видеть также, что ряд Фурье в действительной форме может быть получен путем отбрасывания гармоник с отрицательными частотами и удвоения амплитуд гармоник с положительными частотами (для значений k ≥1).
Относительно определения значений коэффициентов C&k для случая k < 0 (т.е. для отрицательных частот) нетрудно показать, что для действи-
тельных функций ξ (t ) справедливо следующее: поскольку для положи-
тельных значений k
C&k |
= |
1 |
T∫/ 2 |
ξ (t )exp |
− j |
2πk t dt , |
|
|
|||||||
|
|
T −T / 2 |
|
|
T |
|
то для отрицательных величин (−k )
|
1 |
T / 2 |
|
2πk |
|
|
1 |
T / 2 |
|
|
2πk |
|
|
|
|
C&(−k ) = |
|
∫ |
ξ (t )exp j |
|
t dt = |
|
∫ |
ξ (t )exp |
− j |
|
t dt |
= C&k |
. (1.32) |
||
T |
T |
|
T |
||||||||||||
|
−T / 2 |
|
|
T |
−T / 2 |
|
|
|
|
|
|
Таким образом, для действительных функций времени комплексные коэффициенты Фурье (комплексные амплитуды) гармоник с отрицательной
25
частотой определяются комплексно-сопряженными коэффициентами Фурье соответствующих гармоник с положительной частотой.
Сопоставляя амплитудно-фазовую форму ряда Фурье (1.10б), его тригонометрическую (1.11) и комплексную (1.31) формы, нетрудно видеть,
что амплитуды Ak и фазы ϕk формы (1.10б), коэффициенты ak |
и bk триго- |
|||||||||||||||||
нометрической формы (1.11) и комплексные коэффициенты |
C&k формы |
|||||||||||||||||
(1.31) связаны соотношениями |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
C&k |
= ReC&k |
+ j ImC&k |
= 0,5[ak − jbk ], |
|
|
||||||||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ReC& |
k |
= 0,5a , |
ImC& |
k |
= 0,5b ; |
|
(1.33а) |
|||||||||
|
|
|
|
|
|
k |
|
|
|
|
|
|
k |
|
|
|||
|
|
|
|
C& |
|
= 0,5 |
a2 |
+b2 |
= 0,5A ; |
|
(1.33б) |
|||||||
|
|
|
|
|
|
|||||||||||||
|
|
|
|
k |
|
k |
k |
|
|
|
|
k |
|
|
||||
|
|
argC&k |
= −arctg bk |
|
|
|
|
Ak sinϕk |
|
|
|
|||||||
|
|
= −arctg |
|
=ϕk . . |
(1.33в) |
|||||||||||||
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
ak |
|
|
|
|
Ak cosϕk |
|
|
|||||
Набор комплексных коэффициентов C&k |
= |
|
C&k |
|
exp(− jϕk ), характеризую- |
|||||||||||||
|
|
|||||||||||||||||
щих |
комплексные |
|
амплитуды |
гармонических |
составляющих |
|||||||||||||
|
2πk |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
называется ком- |
|
exp − j |
T |
t периодической функции (сигнала) ξ (t ), |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
плексным спектром этой функции [1, 7]. Набор модулей комплексных ам-
плитуд C&k , характеризующий частотную зависимость амплитуд гармоник называется амплитудным спектром функции ξ (t ). Набор величин ϕk , ха-
рактеризующий зависимость начальной фазы гармоник от частоты, назы-
вается фазовым спектром, а набор величин Pk =C&kC&k = C&k 2 , характери-
зующий частотное распределение мощности гармоник называется спектром мощности.
Комплексный спектр функции можно наглядно представить графиче-
ски. Если расположить векторные амплитуды C&k гармонических состав-
26
ляющих разложения в трехмерном пространстве, совместив действитель-
ную часть ReC&k с осью OX (при этом мнимая часть ImC&k будет автома-
тически совмещена с осью OY ), а номера гармоник расположить вдоль оси OZ , то будет получено трехмерное разложение действительной функции
ξ (t ) в комплексный ряд Фурье (см. рис. 1.1).
|
ImC&k |
|
|
|
|
|
ImC& |
& |
|
ImC& |
f |
|
1 |
ImC2 |
|
||
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
3 f1 |
ϕ3 |
|
|
|
|
|
|
|
2 f1 |
|
ϕ |
2 |
ReC&3 |
|
ϕ1 |
|
|
|
|
|
f1 |
|
|
& |
|
|
ReC& |
|
|
ReC2 |
|
0 |
1 |
|
|
|
|
|
|
|
|
|
ReC&k
Рис. 1.1. Графическое представление комплексного спектра
1.5. Примеры разложения периодических функций в ряд Фурье. Некоторые свойства спектра периодических функций
Рассмотрим ряд примеров разложения периодических функций в ряд Фурье и продемонстрируем на этих примерах некоторые свойства спектра.
1.5.1. Разложение в ряд Фурье периодической последовательности прямоугольных импульсов
Зададим периодическую последовательность прямоугольных импульсов, следующих с периодом T и имеющих длительность τ (рис. 1.2). Для упрощения расчетов амплитуду импульсов примем за единицу. Считая, что
27
один из импульсов последовательности расположен симметрично относительно начала координат (см. рис. 1.2.), представим данную последовательность комплексным рядом Фурье (1.31) и определим коэффициенты разложения
|
|
1 |
T / 2 |
|
|
|
|
2πk |
|
|
|
1 |
τ / 2 |
|
|
2πk |
|
||
C&k |
= |
|
∫ |
ξ (t )exp |
− j |
|
|
t dt = |
|
|
∫ |
exp |
− j |
|
t dt = |
||||
T |
|
T |
T |
T |
|||||||||||||||
|
|
−T / 2 |
|
|
|
|
|
|
|
|
−τ / 2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
2πk |
τ / 2 |
|
|
|
|
|
|
|
||
|
|
|
|
|
1 |
sin |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
= |
|
|
T |
|
|
. |
|
|
|
|
(1.34) |
||||
|
|
|
|
q |
|
|
2πk |
τ / 2. |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
ξ (t )
τ
−τ 2 0 τ 2 |
t |
T
Рис. 1.2. Периодическая последовательность прямоугольных импульсов
Величина q =T /τ называется скважностью. Набор коэффициентов C&k
представляет собой комплексный спектр импульсной последовательности,
изображенный на рис. 1.3. Опуская общий множитель 1/ q , |
видим, что |
||||||||
функция |
|
|
|
|
|
|
|
||
|
|
|
2πk |
τ / 2 |
|
|
|
|
|
|
sin |
T |
|
|
|
sin x |
|
||
|
|
|
|
|
|
= |
(1.35) |
||
|
|
2πk |
τ |
/ 2 |
|
x |
|||
|
|
|
|
|
|||||
|
|
|
T |
|
|
|
|
характеризует распределение амплитуд и фаз частотного (спектрального) представления анализируемого сигнала. Амплитуды спектральных (гармонических) составляющих могут быть изображены в виде так называемых «спектральных линий», расположенных вдоль оси частот, а изменение
28
огибающей комплексной амплитуды спектральных коэффициентов определяется функцией (1.35), известной в математике как первый замечательный предел. Значение этого предела в точке k = 0 есть единица. Таким образом, спектр периодической функции является дискретной функцией частоты.
C&K |
∆f |
|
sin x / x |
||
|
|
|
|
3 |
4 |
k |
1 |
1 |
1 |
2 |
|
|
|
|
||||
T |
T |
|
|
|
|
Рис. 1.3. Спектр периодической последовательности прямоугольных импульсов
Для определения спектра последовательности, изображенной на рис. 1.2. была использована комплексная форма ряда Фурье, но нетрудно показать, что, с учетом четности анализируемой функции, коэффициенты представления этой функции тригонометрическим рядом Фурье имеют вид
|
|
|
a0 |
= |
|
1 |
τ∫/ 2 |
dt =1/ q; |
|
bk |
= 0; |
|
|
|
(1.36а) |
|||
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
T −τ / 2 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
2 |
τ / 2 |
2πkt |
|
1 |
|
|
2πkτ / 2 |
|
2πkτ / 2 |
|
|
|||||
ak |
= |
|
∫ cos |
|
|
dt = |
|
sin |
|
|
|
/ |
|
. |
(1.36б) |
|||
T |
|
T |
q |
T |
|
T |
||||||||||||
|
|
−τ / 2 |
|
|
|
|
|
|
|
|
|
|
Таким образом, использование комплексной формы ряда Фурье предпочтительнее, поскольку сокращает объем вычислений.
Отметим, также, что функция (1.36) является знакопеременной и не может быть спектром амплитуд, который представляет собой набор модулей коэффициентов C&k ≥ 0 . Поскольку отрицательные лепестки спектра ха-
29
рактеризуются отрицательными величинами C&k , эти величины можно за-
писать в виде (−1) C&k = C&k exp{jπ}. Таким образом при переходе функции sin x / x в отрицательную область фазы спектральных линий, входящих в состав отрицательного лепестка спектра, претерпевают скачок на 1800 относительно фаз спектральных линий положительных лепестков. На рис. 1.4 изображены спектры амплитуд и фаз периодической последовательности прямоугольных импульсов.
C&k
1
|
|
|
|
|
k |
|
|
|
|
|
|
0 1 2 3 |
|
|
|
||
|
argC&k |
|
|
|
|
π
0
Рис.1.4. Амплитудный и фазовый спектры периодической последовательности прямоугольных импульсов
Продолжим анализ спектрального разложения (1.34) с целью установления связи между шириной основного лепестка спектральной функции и длительностью одиночного импульса последовательности. За основной лепесток примем часть графика функции шириной ∆f , заключенную меж-
ду первыми нулями функции sin x / x |
(рис. 1.3). |
Поскольку первый нуль функции |
sin x / x расположен в точке x =π , то |
полуширина основного лепестка (т.е. расстояние от начала координат до первого нуля) определяется из условия
30
2π |
k |
τ =π , |
|
||
|
|
|
|||
|
T 2 |
|
|||
откуда следует, что |
|
||||
|
|
k |
τ =1. |
(1.37) |
|
|
|
|
|||
|
|
T |
|
||
Величину k /T можно обозначить как ∆f , |
т.к. эта величина определяет |
||||
верхнюю частоту спектрального разложения, |
отвечающую первому нулю |
спектра (основному лепестку). Назовем величину ∆f шириной основного
лепестка спектра импульсной последовательности. |
|
Тогда выражение (1.37) может быть записано в форме |
|
∆fτ =1, |
(1.38) |
из которой следует, что произведение длительности импульса на ширину главного лепестка его спектра есть постоянная величина.
Выражение (1.38) есть одна из форм принципа неопределенности, из которого следует, что ширина спектра и длительность импульса обратно пропорциональны.
В дальнейшем будет показано, что в основном лепестке спектра сосредоточено более 90% энергии сигнала и поэтому в большинстве случаев ширину основного лепестка называют шириной спектра сигнала.
1.5.2. Теорема смещения для разложения в ряд Фурье
Рассмотрим разложение в ряд Фурье последовательности прямоугольных импульсов с единичной амплитудой, длительностью τ и периодом T (рис. 1.5). Начало первого импульса последовательности в данном случае совпадает с началом координат, т.е. данная последовательность смещена во времени на величину τ / 2 относительно последовательности ξ (t ), изо-
31