Степанов В.Н. Ряды Фурье. Преобразования Фурье. Интеграл

advertisement
Министерство образования Российской Федерации
Омский государственный технический университет
РЯДЫ ФУРЬЕ
ПРЕОБРАЗОВАНИЯ ФУРЬЕ
ИНТЕГРАЛ ФУРЬЕ
Методические указания к практическим занятиям
для студентов технических вузов
Омск-2003
Составители: Геннадий Николаевич Бояркин, канд.физ.-мат.наук, профессор
Владимир Николаевич Степанов, канд.физ.-мат.наук, доцент
Любовь Степановна Рыженко, ст. преподаватель
2
1. РЯДЫ ФУРЬЕ ДЛЯ ФУНКЦИЙ С ПЕРИОДОМ 2 и 2
Рядом Фурье периодической функции f ( x ) с периодом 2 , определенной на
сегменте  , , называется ряд
a0 
(1)
  a n cosnx  b n sin nx  ,
2 n 1
где
1
a n   f ( x ) cos nx dx (n  0,1,2...),
(2)
 
bn 
1
 f ( x ) sin nx dx (n  1,2...).
 
(3)
Если ряд (1) сходится, то его сумма S( x ) есть периодическая функция с периодом 2 , т.е. S( x  2)  S( x ) .
Теорема Дирихле. Пусть функция f ( x ) на сегменте  ,  имеет конечное
число экстремумов и является непрерывной за исключением конечного числа точек разрыва 1-го рода (т.е. удовлетворяет так называемым условиям Дирихле). Тогда ряд Фурье этой функции сходится в каждой точке сегмента  ,  и сумма
этого ряда S( x ) вычисляется:
1) S( x )  f ( x ) во всех точках неразрывности f ( x ) , лежащих внутри сегмента
 , ;
2) Sx 0  
3) Sx  
1
f x 0  0  f x 0  0, где x 0 - точка разрыва 1-го рода функции f (x) ;
2
1
f    0  f   0 на концах промежутка, т.е. при x   .
2
В случае, когда f ( x ) - четная функция, ее ряд Фурье содержит только свободный член и косинусы, т.е.

a
f ( x )  0   a n cos nx ,
(4)
2 n 1
где
2
a n   f ( x ) cos nx dx (n  0,1,2...).
(5)
0
3
В случае, когда f ( x ) - нечетная функция, ее ряд содержит только синусы, т.е.

f ( x )   b n sin nx dx,
(6)
n 1
где
2
b n   f ( x ) sin nx dx.
(7)
0
Часто приходится разлагать в тригонометрический ряд функции периода, отличного от 2 . В этом случае, если f ( x ) - периодическая функция с периодом 2 ,
для которой выполняются на сегменте  ,  условия Дирихле, то указанная
функция может быть представлена в виде суммы ряда Фурье:
a0  
nx
nx 
   a n cos
 b n sin
,
2 n 1

 
где
an 
1
nx
f ( x ) cos
dx,

 

(8)
(9)
1
nx
b n   f ( x ) sin
dx.
(10)
 

В случае, когда f ( x ) - четная функция, как (4) – (5), ее ряд Фурье содержит
только свободный член и косинусы, т.е.
f (x) 
где
a0 
nx
  a n cos
dx,
2 n 1

(11)
2 
nx
a n    f ( x ) cos
dx .
 0

(12)
В случае, когда f ( x ) - нечетная функция, ее ряд Фурье содержит только синусы,
т.е.

f ( x )   b n sin
n 1
nx
dx,

где
(13)
2 
nx
b n    f ( x ) sin
dx.
(14)
 0

При разложении в ряд Фурье целесообразно придерживаться следующей схемы. Вначале проверяем, что данная функция удовлетворяет условиям Дирихле; затем вычисляем коэффициенты a n и b n по соответствующим формулам; подставляя их в ряд, получаем искомое разложение; наконец, основываясь на теореме Ди-
4
рихле, определяем, при каких x полученный ряд сходится к данной функции. Рассмотрим примеры разложения в ряд Фурье периодических функций.
1. Разложить в ряд Фурье функцию периода 2 , заданную на интервале
   x   формулой: f (x)  x (рис. 1).
y
-4π
-3
-2
-

0
2
3 4π
x
Рис. 1
Решение. Эта функция удовлетворяет условиям Дирихле и, следовательно, может быть разложена в ряд Фурье. Применяя формулу (7), найдем коэффициенты
x  u , sin nx dx  dv 



Фурье  интегрируем по частям
1
:
dx  du, v   cos nx 

n


2
2   x cos nx  1 
2
2 1

,
b n   x sin nxdx  
  cos nxdx    cosn  
0

n
n 

n
0

0


1
т.к.   cos nx dx   sin nx  0 .
n


n
Следовательно, ряд Фурье функции f ( x ) будет иметь вид
n
 sin x sin 2x sin 3x sin 4x


 1
2



 ... 
sin nx  ... .
2
3
4
n
 1

Так как функция f (x)  x удовлетворяет условиям Дирихле, то в любой точке
непрерывности f ( x ) сумма ряда равна значению функции. В точках  n и n
сумма ряда равна нулю. На рис. 2 показаны графики: функции f ( x ) и частичных
сумм ряда, содержащие 1, 2 и 3 члена. Из рисунка видно, как график частичных
сумм ряда приближается к графику функции f ( x ) при увеличении членов суммы.
5
sin 2x sin 3x 

2 sin x 


2
3 

sin 2x 

2 sin x 

2 

2 sin x
y
Рис. 2
2. Разложить в ряд Фурье функцию с периодом 2 , заданную на интервале
 ,  формулой

0,    x  0; a  x  

.
f ( x )  1,
0xa
1
 , x  0, x  a
2
Построим график функции (рис. 3).
y
1
-3
-2

-
2
3
x
Рис. 3
Решение. Функция удовлетворяет условиям Дирихле. Применяя формулы (2) и
(3), находим коэффициенты Фурье
a0 
1
1a
x a a
f
(
x
)
dx

dx

 ,


 
0
0 
an 
a
1
1a
1
sin na
f
(
x
)
cos
nx
dx

1

cos
nx
dx

sin
nx

,


 
0
n
n
0
6
bn 
a
1a
1a
11
cos nx


f
(
x
)
sin
nx
dx

1

sin
nx
dx


cos
nx



0
0
n
n
0
a
0

1  cos na
.
n
Разложение в ряд Фурье f ( x ) имеет вид
f (x) 
 sin na
a
1  cos na

 cos nx 
 sin na .
2 n 1 n
n
3. Разложить в ряд Фурье функцию с периодом T  4 , заданную на интервале
6, 0  x  2,
(0; 4) формулой f ( x )  
3x, 2  x  4.
Построим график функции (рис. 4).
-8 -6 -4 -2
0 2
4
6
8 10 12
Рис. 4
Решение. Пользуясь формулами (9) и (10), полагая   2 и разбивая интервал
интегрирования точкой x  2 на две части, поскольку в каждой из них функция задана различными формулами, получим
4
14
nx
12
nx
nx 
a n   f ( x ) cos
dx    6 cos
dx   3x cos
dx  
20
2
20
2
2
2

1  12
nx 2  2
nx
4
nx  4 
6
  sin
 3 sin
 2 2 cos
   2 2 1  cos n, n  0 .
2  n
2 0  n
2
n 
2  2 n 
При а – четном cos n  1 и x n  0 , при n – нечетном cos n  1 и a n 
7
12
.
n 2 2
При n=0
2
4
14
12
 1  2 3x
a 0   f ( x ) dx    6 dx   3x dx    6x 
20
20
2
 2 0 2

  15 ,

2
4
4
14
nx
12
nx
nx 
b n   f ( x ) sin
dx    6 sin
dx   3x  sin
dx  
20
2
20
2
2
2

1  12
nx 2  4
nx 2x
nx  4 
  cos
 3 2 2 sin

cos
 
2  n
2 0 n 
2
n
2  2 
1
nx 0 3  2x 
nx  4 4 2
nx 

 12 cos
    cos
dx 
   cos
2n
2 2 2  n 
2  2 2 n
2

0

1
121  cos n  32  4 cos2n  4 cos n  
2n

1
12  12 cos n  24  12 cosn   12   6 .
2n
2n
n
Искомое разложение данной функции имеет вид
f (x) 
15 12 
x 1
3x 1
5x

 2  cos  cos
 cos
 ... 
2  
2 9
2
25
2

6  x 1
2x 1 3x

  sin
 sin
 sin
 ... .

2 2
2
3
2

Оно справедливо во всей области определения данной функции: в интервале
(0, 2) сумма ряда S(x)  6 , в интервале (2, 4) - S( x )  3x . В точке разрыва x  2 ,
S( x) 
1
f (2  0)  f (2  0)  6 .
2
8
4. Разложить в ряд Фурье функцию с периодом T  2 , заданную на интервале
(1,1) формулой f ( x )  x (рис. 5).
у
1
-4
-3
-2
-1
0
1
2
3
4
х
Рис. 5
Решение. Функция удовлетворяет условиям Дирихле. Применяя формулу (12),
1
x2 1
(функция f ( x ) - четная), полагая   1 , получим a 0  2 x dx  2
 1,
2 0
0
 x
1
1
a n  2  x cos nx dx  2 sin nx  2 2 cos nx  
n 
0
 n
0
1
0, n  четные,
2

 2 2 cos n  1  
4

n 
2
 n  2 , n  нечетное.

f (x) 
1 4  cos2n  1 x
.
 
2  2 n 1 2n  12
ПРИМЕРЫ ДЛЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ
1. Разложить в ряд Фурье функцию с периодом 2 , заданную на интервале
 2x при    x  0,
 ,  уравнением f (x)  
.
3
x
при
0

x


.

5 10  cos x cos 3x cos 5x

 


 ... 
2
2
4
  12
3
5

 sin x sin 2x sin 3x




 ... .
2
3
 1

Ответ: f ( x ) 
2. Разложить в ряд Фурье функцию с периодом 2 , заданную на интервале
0, 2 формулой f (x)  x .
2


sin nx
Ответ: f ( x )   
.
2 n 1 n
9
3. Разложить в ряд Фурье функцию с периодом 2 , заданную на интервале
 ,  формулой y  sin x .
Ответ: f ( x ) 

2 4  cos 2x cos 4x
cos 2kx
 

 ... 
 ... .
   1 3
35
(2k  1)  (2k  1)

4. Разложить в ряд Фурье периодическую функцию f ( x ) с периодом 2 ,
заданную на интервале  ,  уравнением f ( x )    x .

Ответ: f ( x )    2
 1n 1 sin nx .
n 1
n
5. Разложить в ряд Фурье периодическую функцию
при 0  x  1,
x
f(x+2)=f(x).
f (x)  
 2  x при 1  x  2,
Ответ: f ( x ) 
1 4  cos(2n  1) x
.
 
2  2 n 1 2n  12
6. Разложить в ряд Фурье периодическую функцию
при  3  x  0,
0
f(x+6)=f(x).
f (x)  
при 0  x  3,
x
2n  1 x  3   1 sin nx .
3 6 
1
Ответ: f ( x )   2 
cos

4  n 1 (2n  1) 2
3
 n 1 n
3
n
7. Разложить в ряд Фурье периодическую функцию
f (x)  10  x при 5  x  15, f (x  10)  f (x).
Ответ: f ( x ) 

10
n
  1
 n 1
sin
nx
5 .
n
10
2. РЯД ФУРЬЕ ДЛЯ НЕПЕРИОДИЧЕСКОЙ ФУНКЦИИ.
РЯДЫ ФУРЬЕ В КОМПЛЕКСНОЙ ФОРМЕ
Пусть y  f (x) – непериодическая функция, заданная на всей числовой оси. Так
как сумма тригонометрического ряда является периодической функцией, то
очевидно, что данная непериодическая функция не может быть разложена в ряд
Фурье. Но если функция задана на конечном интервале    x   , то для нее
можно построить ряд Фурье, который имел бы ее своей суммой на этом интервале.
Для этого рассматривают вспомогательную функцию ( x ) периода 2 ,
значения которой на интервале    x   совпадают со значениями функции f ( x )
(рис. 6).
четное
y
y
f (x)
y  f (x)
( x )
производная
-3  -2  - 

2 3
х
x
нечетное
Рис. 6
Рис. 7
Если для функции ( x ) выполняется условие теоремы Дирихле, то ее можно
представить соответствующим рядом Фурье. Этот ряд на интервале    x   во
всех точках непрерывности функции будет иметь своей суммой (x)  f (x) .
Иногда приходится иметь дело с функциями, заданными только в интервале
0  x   . В этом случае мы можем сначала продолжить по какому-либо закону
фукнцию на интервал    x  0 , а затем продолжить на всю числовую прямую
периодически с периодом 2 . Удобнее всего продолжить функцию на интервал
   x  0 четным или нечетным образом (рис. 7). В первом случае ряд Фурье
будет содержать только косинусы и свободный член. Во втором случае ряд Фурье
будет содержать только синусы.
Рассмотрим примеры разложения в ряд Фурье непериодической функции.
11
1. Разложить в ряд Фурье функцию f ( x ) , заданную на отрезке 0, 2 уравнением
x2
.
2
Решение. Функция может быть разложена в ряд Фурье бесчисленным количеством способов. Рассмотрим два наиболее важных варианта разложения.
f (x)  x 
y
1/2
-3 -2 -1 0
1
2
3
4
5
6
7
8
x
Рис. 8
А. Доопределим функцию f ( x ) на отрезке  2, 0 четным образом (рис. 8).
Имеем   2 .
 x2 1 3  2 2
1 2

a 0    x  x  dx    x   ;
2 
0
 2 6 0 3
2
0
1 
nx
2 
1  nx 2 2 2

1  x  sin nx dx 
a n    x  x 2  cos
dx   x  x 2  sin


2 
2
n 
2 
2 0 n 0
2
0
2

2 2
1  x sin nx dx .

n 0
2
Еще раз интегрируем по частям:
0
4
nx
4
nx
a n  2 2 1  x cos
 2 2  cos
dx 
n 
2 0 n  0
2
2

2


4
4
4
n
cos
n




1   1 .
2 2
2 2
2 2
n 
n 
n 
bn  0 .
Итак,
1 4  1   1
nx 1 8  1
1
1

f (x)   2 
cos
  2  2 cosx  2 cos2x  2 cos3x  ... .
2
3  n 1 n
2
3  2
4
6

n
12
Б. Доопределим функцию f ( x ) на отрезке 2, 0 нечетным образом (рис. 9).
y
1/2
-3 -2 -1
0 1
2
3
4
5
6
x
Рис. 9
0
1  nx
2 
1 
nx 2 2 2

1  x  cos nx dx 
b n    x  x 2  sin
dx    x  x 2  cos


2 
2
n 
2 
2 0 n 0
2
0
2
2 2
nx
2 
2
nx 2 2 2
nx
 dx 
  1  x cos
dx  1  x   sin
  sin
n 0
2
n 
n
2 0 0 n
2

4
nx 2
4 2 nx
 2 2 1  x sin

 sin 2 dx 
n 
2 0 n 2 2 0
0

8
nx
cos
3 3
n 
2
2
0



8
cos n  1  38 3 1   1n .
3
n 
n 
3
Итак, a n  0, a 0  0 .
8  1   10 
nx 16  x 1
3x 1
5x

f (x)  3 
sin
 3  sin
 3 sin
 3 sin
 ... .
3
 n 1
n
2
 
2 3
2
5
2

n
2. Разложить в ряд Фурье функцию f ( x ) , заданную на интервале 0,  уравнением f ( x )  x 2 .
Решение. Рассмотрим два возможных (из бесчисленных) способа разложения
этой функции в ряд Фурье на заданном интервале.
А. Будем полагать, что функция задана на отрезке длиной, равной периоду
T   , и периодически продолжить ее на всю числовую ось с этим периодом
(рис. 10).
13
y
-2 -
0

2
х
Рис. 10
Вычисляем коэффициенты Фурье полученной функции по общим формулам

(9), (10), полагая   .
2
3 

2
2x
2
a 0   x 2 dx 
 2 ;
0
3 0 3

 x2
2 2
2 x
1 
a n   x cos 2nx dx   2 cos 2nx  
 3  sin 2nx  
0
  2n
 2n 4n 
0

cos 2n 1
 2,
n2
n
n  0 ;

 x2
2 2
2 x
1 

b n   x sin 2nx dx   2 sin 2nx  
 3  cos 2nx    ;
0
  2n
n
 2n 4n 
0
2   1


f (x) 
   2 cos 2nx  sin 2nx  .
3 n 1 n
n

Б. Доопределим функцию f ( x ) на отрезке  , 0 четным образом и периодически на всю числовую ось. В данном случае T  2 . Вычисляем коэффициенты
Фурье полученной функции по формулам (5).
График функции представлен на рис. 11.
y
-

Рис. 11
14
x
2 2
2x 3
a 0   x dx 
0
3

2 2

;
3
0

x 2  u,
2
a n   x 2 cos nx dx 
0
du  2 x dx,
cos nx dx  dv

1
v  sin nx
n
u  x,


2  x2
1
4 
  sin nx   sin nx  2x dx     x sin nx dx 
 n
n0
n 0
0
du  dx,

0
4   1
1
 4 cos n

n 4




1
,
   x    cos nx    cos nx dx  
n2
n2
n   n
0 0n



dv  sin nx dx

1
v   cos nx
n
n  0.
3
 cos x cos 2 x cos 3x

n 1 cos nx
 4 2 

 ...   1
 ... .
Итак, f ( x ) 
2
2
2
3
2
3
n
 1

Заметим, что ряды Фурье, полученные в пп. А и Б, сходятся на отрезке 0,  к
одной и той же формуле f ( x )  x 2 , во втором случае вычислений нужно проводить
меньше, чем в первом.
Во многих случаях удобно использовать комплексную формулу ряда Фурье,
которую можно получить с помощью формул Эйлера:
e ix  cos x  i sin x .
eix  cos x  i sin x ;
Для функций с произвольным периодом T  2 ряд Фурье в комплексной форме имеет вид

f (x)   Cn e
in x

,
(15)
n  
где

1 
Cn 
f (x) e

2 
15
inx

dx .
(16)
3. Разложить в ряд Фурье функцию f ( x ) , заданную на интервале  ,  уравнением y  e  x .
Решение. В данном случае удобно использовать комплексную форму ряда
Фурье. По формуле (16)
1   x inx
1  1i n  x
e  1i n  x  e  1i n  x 
Cn 
dx  


 e e dx  2  e
2  
21  in   21  in  

e 1i n    e 1i n   e  ei n   e   e i n 
.


21  i n 
21  i n 
По формулам Эйлера
e  i n   cos n  i sin n   1 .
n
n

 1 e   e   
Следовательно, C n 
,
21  i n 
e
x

  Cn e
n  
e   e     1 e i n x
.


2 n  1  i n
n
inx
В интервале  ,  ряд представляет функцию e  x , а в точках x   его сум1
ма равна e   e  .
2
Заметим, что полученный ряд в комплексной форме можно преобразовать к
обычной тригонометрической форме ряда Фурье, для этого следует объединить
слагаемые с индексами n и  n и заменить в результате по формулам Эйлера показательные функции тригонометрическими:
n
n
inx

 1 e i n x  1 e i n x
 1  i n e i n x
n 1  i n e
u n  u n 

  1
2
1 in
 2   1
n
1 in
1 n
cos nx  n sin nx
, n  1,2,3,...
1  n2
при n  0
e  e
.
a0 
2
16

Следовательно,
e
x
e  e


 1   1n

cos nx  n sin nx  .
 
2
 2 n 11  n

ПРИМЕРЫ ДЛЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ
1. Разложить в интервале 0,   по синусам кратных дуг функцию


x
при
0

x

,

2
f (x)  
.

0 при
 x  .

2

Ответ: f ( x )   b n sin nx , где b 2 k   1
k 1
n 1
1
2
k
; b 2 k 1   1
2
2k
2k  1
.
2. Разложить в интервале 0,   по косинусам кратных дуг функцию


при 0  x  ,
cos x
2
f (x)  
 cos x при   x  .

2
Ответ: f ( x ) 
4 1 
n 1 cos 2nx 
   1
.

  2 n 1
4n 2  1 
3. Разложить в интервале (0,1) в ряд Фурье функцию f ( x )  2x .
Ответ: f ( x )  1 
2  sin 2nx
.

 n 1
n
4. Разложить в ряд Фурье функцию, заданную на интервале  ,  уравнением
 x при    x  0,
f (x)  
 0 при 0  x  .
Ответ:
 2  cos2n  1x 
n sin nx
.
 
   1
2
4  n 0 2n  1
n
n 1
17
5. Разложить в интервале (0, 1) в неполный ряд Фурье, содержащий только косинусы или только синусы, функцию
 0,3 при 0  x  0,5,
f (x)  
 0,3, при 0,5  x  1.
Ответ: а) f ( x ) 
б) f ( x ) 
1,2  cos x cos 3x cos 5x

n 1 cos 2n  1x





...


1
 ... ;

  1
3
5
2n  1

1,2  sin 2x sin 6x sin 10x
sin 2n  1x




...


...
 .
  1
3
5
2n  1
6. Разложить в ряд Фурье периодическую функцию f ( x )  e x при  2  x  2 ,
f ( x )  f ( x  4) .
n


e4  1  1
 1 
nx
nx 

2
2
cos

n

sin
Ответ: f ( x ) 

 .


2 2
2e 2  2
2
2 
n 1 4  n  
7. Разложить функцию f (x)  1 в интервале 0,   в неполный ряд Фурье, содержащий только синусы.
Ответ: f ( x ) 
2n  1x .
4  1
sin

 n 1 2n  1

3. ИНТЕГРАЛ ФУРЬЕ
Пусть функция (сигнал) y  f (t ) описывает некоторый периодический процесс.
С целью исследования этого процесса часто представляют функцию f ( t ) в виде
суммы постоянного члена и гармонических составляющих с частотами
0 , 1  20 ,... :
a0 

  a n cos n0 t  b n sin n0 t , 0  .
(17)
2 n 1

Совокупность коэффициентов Фурье периодической функции называется ее
спектром. Спектр периодической функции дискретный. С точки зрения физики
разложение в ряд Фурье можно трактовать как представление периодического сигнала в виде суммы гармонических колебаний.
f (x) 
18
Если на числовой оси  ,   задан непериодический сигнал f ( t ) , то для исследования такого процесса представим f ( t ) на промежутке  ,   в виде ряда
(17). За пределами рассматриваемого промежутка сумма тригонометрического ряда будет повторять значения функции f ( t ) в промежутке  ,  . После этого естественно сделать предельный переход при    . Оказывается, что в результате
такого предельного перехода произойдет качественный скачок. Непериодическая
функция, заданная на всей оси, представится в виде интеграла, который является
непрерывным аналогом ряда Фурье и представляет собой “сумму гармонических
составляющих”, частоты которых заполняют всю действительную полуось
0     . А именно, справедлива следующая теорема.
Теорема. (О представимости функции интегралом Фурье).
Если функция f ( t ) абсолютно интегрируема на всей числовой оси, т.е.

 f ( t ) dt   и удовлетворяет условиям Дирихле на любом конечном промежутке

этой оси, то при всех t имеет место равенство

1
f ( t )   d  f () cos t  d .
 0 
(18)
Если t 0 - точка разрыва первого рода функции f ( t ) , то левую часть формулы
f ( t 0  0)  f ( t 0  0)
(18) следует понимать как
, что мы всегда будем иметь в виду.
2
Формула (18) называется интегральной формулой Фурье, а ее правая часть –
двойным интегралом Фурье.
Обозначая
1
A()   f () cos   d и
 
1
B()   f () sin   d ,
 
(19)
запишем интегральную формулу Фурье (18) в виде

f ( t )   Acos t  Bsin t d .
(20)
0
Интегральная формула Фурье (20) аналогична разложению периодической
функции в ряд Фурье. Подынтегральная функция формулы (20) напоминает общий
член ряда Фурье, только здесь частота  , непрерывно изменяясь, пробегает все
значения от 0 до  , и потому суммирование заменяется интегрированием от 0 до
 . Функции A(), B , определенные формулами (19), аналогичными формулам
19
для коэффициентов ряда Фурье, дают закон распределения амплитуд и начальных
фаз в зависимости от частоты  .
Смысл интегральной формулы Фурье состоит в следующем: интегральная формула Фурье представляет непериодическую функцию как наложение гармоник с
непрерывной последовательностью частот.
С физической точки зрения это означает, что непериодический процесс уже
нельзя построить из гармонических колебаний только с определенными изолиро
ванными частотами n  n , теперь для его построения необходимы гармониче
ские колебания всех частот.
Рассмотрим частные случаи применения формулы (20).
1. Пусть f ( t ) - четная функция. Тогда A() 
2
 f () cos  d, B()  0 .
0
Формула (20) в этом случае принимает вид
f (t) 
2  

  f cos   d  cos t d .

 0 0

(21)
Запишем эту формулу в симметричном виде, положив f ( t ) 
2 
 FC  cos  t d ,
 0
2 
где FC () 
 f  cos   d .
 0
FC () называется косинус-преобразованием Фурье функции f ( t ) .
2. Пусть f ( t ) - нечетная функция. В этом случае A  0 ,
2
B()   f () sin   d ,
0
тогда
2  

f ( t )     f sin   d  sin t d .
 0 0

Введем синус-преобразование Фурье, положив FS () 
(22)
2 
 f  sin   d , тогда
 0
2 
 FS  sin  t d .
 0
3. Пусть f ( t ) определена только в интервале 0,   . Тогда ее можно представить при t  0 как формулой (21), так и формулой (22). Для этого следует функцию
формула (22) принимает вид f ( t ) 
20
f ( t ) продолжить на промежуток  , 0 так, чтобы она стала или четной или нечетной на всей действительной оси.
Теперь рассмотрим примеры представления непериодических функций интегралом Фурье.
1. Представить интегралом Фурье функцию
f (t )
1
1,  1  t  3,
 1
f ( t )   , t  1, t  3,
2
0, t  3, t  1.
1/2
-1
0
3
t
Рис. 12
График функции представлен на рис. 12.
Решение. Данная функция на любом конечном промежутке числовой оси
 ,   удовлетворяет условиям Дирихле. Очевидно, что f (t ) является абсолютно
интегрируемой функцией на всей числовой оси, так как

3
3

1
1
 f ( t ) dt  1  dt  t  4   .
Следовательно, данная функция может быть представлена интегралом Фурье; по
формуле (17) имеем

1
1 3
f ( t )   d  f () cos  t  d   d 1  cos  t   d 
 0 
 0 1

1   sin t   3
1   sin t  3  sin t  1
d


d 


0

0

1

2  cost  1  sin 2
d .

0

В точках разрыва f ( t ) , т.е. при t  1 и t  3 , полученное представление сохра1
1
няется, т. к. в этих точках f ( t  0)  f ( t  0)  f ( t ) .
2
2
В частности, при
t  3 имеем f (3) 
sin x

dx  .
x
2
0

находим значение интеграла: 
21
2  cos 2  sin 2
1
d  , отсюда легко

0

2
Таким образом, в результате решения основной задачи – представления заданной функции интегралом Фурье – мы смогли вычислить интеграл от функции, первообразная которой через элементарные функции не выражается. И еще одна характерная особенность. Как видно из данного примера, интеграл Фурье может
представлять функцию, которая на разных промежутках числовой оси задается
разными аналитическими выражениями.
2. Представить интегралом Фурье заданную на всей оси функцию
 t  2,  2  t  0,

f ( t )   t  2, 0  t  2,

0,
t  2.

f (t )
2
-2
2
t
Рис. 13
График функции представлен на рис. 13.
Решение. Условия представимости данной функции интегралом Фурье
выполняются. Эта функция – четная, поэтому воспользуемся формулой (21):
 
 2



2 
2


2   cos   d   cos  t d 
f (t) 
f cos   d cos  t d 


 0  0

0 0




2   sin   2 cos   2  sin   2 
2  1  cos 2


 2 


cos  t d 
 cos t d   0
 0
 0
2 0

2
0

4  sin  


 cos  t d .
 0  

2


 sin  
В частности, полагая t  0 , найдем 
 d  .
2
  
0
2

3. Представить интегралом Фурье функцию
f (t )
sin t,
f (t)  
 0,
n - целое число.
Рис. 14
График функции представлен на рис. 14.
22
0  t  n,
t  0, t  n.
Решение. Функция f ( t ) непрерывна на всей числовой оси и абсолютно интегрируема, поскольку

n

0
 f ( t ) dt   sin t dt  n   .
Следовательно, возможно представление этой функции интегралом Фурье. Согласно формуле (18)
f (t) 
1  n
 d  sin   cos t  d .
0
0
Используя равенство
sin   cos t   
1
sin    t    sin    t  ,
2
интегрируя и производя несложные преобразования, получим
1  cos  t   1 cos t  n 
f (t)  
d.
0
1  2
n 1
4. Представить интегралом Фурье функцию, продолжив ее на всю числовую
ось: а) четным, б) нечетным образом.

cos t ,
f (t)  
 0,

f (x)
0
/2
0t

,
2

t .
2
t
Рис. 15
График функции представлен на рис. 15.
Решение. Функции, получаемые продолжением f ( t ) на  , 0  четным и нечетным образом, удовлетворяют условиям теоремы. Для четного продолжения
воспользуемся формулой (21)


2   2
f ( t )    cos   cos   d  cos  t d .

 0  0

Вычислим отдельно внутренний интеграл:
23



1 2
1 2


cos


cos


d


cos
1



d




 cos1   d 
2
2
0
0
0
2
  
sin
  
sin 1    2 sin 1    2
2 2 



2 1  
1   0
2 1   0


  
sin     cos   cos   cos  
2 2 
2 
2 
2 .
 

21  
21   21   1  2 
Следовательно, при четном продолжении

cos 
2
2 cos  t d .
f (t) 
 0 1  2


Если f ( t ) продолжена нечетным образом, то следует применить формулу (22).
Получим

 2

2 

f (t) 
cos


sin


d

 sin  t d .
 0  0




Производя аналогичные вычисления, придем к следующему представлению
функции f ( t ) :

f (t) 

2
0


2 sin  t d .
1  2
sin
Таким образом, для t  0

f (t) 

2
0


 cos


2 sin  t d  2
2 cos  t d .
1  2
 0 1  2
sin

24
5. Найти косинус- и синус-преобразования Фурье функции
f ( t )  exp  at , a  0, t  0 .
Решение.
2
FC () 
 f cos   d , следовательно, используя формулу для вычисления
0
циклического интервала, получим

2  a
2  e a
2
a


FC  
e
cos


d




a
cos




sin



.


0
  a 2  2 
 a 2  2
0
Аналогично

2  a
2  e a
2



FS  
e
sin


d




a
sin




cos



.



2
2
2
2
0
  a   

a


0

Если теперь к полученным функциям применить обратные косинус- и синуспреобразования Фурье, то найдем
2a  cos  t
d  e at , t  0,

2
2
 0a 
и
2   sin  t
d  e at , t  0.

2
2
0 a 
Применяя косинус- и синус-преобразования Фурье, можно получить таблицу
значений несобственных интегралов, зависящих от параметра. Однако основное
назначение косинус- и синус-преобразований Фурье состоит в применениях к решению задач математической физики.

6. Найти функцию ( t ) , если  () sin   d  e  ,   0 .
0
2 
e , как это видно из представленного уравнения,

является синус-преобразованием Фурье функции ( t ) . Поэтому и на основании
формулы обратного синус-преобразования Фурье находим
Решение. Функция F 
2
t  



0

2 
2 e 
 sin  t  t cos  t   2  2 t .
e sin  t d  2

 t 1
 t 1
0
25
sin  
cos  t d .

0

7. Вычислить интеграл   
Решение. Замечая, что
sin   
  cos   d , можем записать

0
1
1  

      cos   d   cos  t d .

 0 0

1, 0  t  ,
Если ввести функцию f ( t )  
то предыдущее равенство примет
0
,
t


,
t

0
,

вид
, 0  t  ,
1
1  

      f () cos   d   cos  t d  f t  , т.е.   

 0  

0, t  , t  0.
ПРИМЕРЫ ДЛЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ
Представить интегралом Фурье функции:
1, 0  t  1,
 1
1. f ( t )   , t  0, t  1,
2
0, t  0, t  1.

Ответ: f ( t ) 

2
0
sin

1  2 t 
 cos
2
2
d .

 1  t , 0  t  1,

2. f ( t )   1  t ,  1  t  0,

0, t  1, t  1.


2   sin 
Ответ: f ( t ) 
sin  t d .
 0 2

26
t  1,
 t ,
3. f ( t )  
 0,
t  1.

2  sin   cos  1
Ответ: f ( t ) 
cos  t d .
0
2

 et ,

4. f ( t )   0,
 e t ,

t  0,
t  0,
t  0.

2  sin  t
Ответ: f ( t ) 
d .
 0 1  2

3

2
t

3
,
0

t

,

2
5. Найти косинус-преобразование Фурье функции f ( t )  
3
0,
 t  .

2
3 

2 cos  1
2 
2
.
Ответ: FC  

2


6. Найти синус-преобразование Фурье функции f ( t )  a  t , 0  t  , a  0 .
Ответ: FC  
2
ln a
.

 ln a 2  2

7. Найти функцию t  , если    cos   d 
0
Ответ: t  
1
.
1  2
2
1

.
 1 t2

 1  cos h  
8. Вычислить интеграл   
 sin  t d .



0


 , 0  t  h,
Ответ:    2
0, t  h.
27
4. ПРЕОБРАЗОВАНИЕ ФУРЬЕ
Интегральную формулу Фурье можно записать в виде
f (t) 

1 
d

f  e j  t   d .


2 

(23)
Это есть комплексная форма интеграла Фурье. Введем функцию

F   f  e  j   d ,
(24)
1 
f (t) 
f e j  t d .

2 
(25)

то согласно (23) получим
Переход от F  к F по формуле (24) называется прямым преобразованием
Фурье. Восстановление f ( t ) по F с помощью формулы (25) называется обратным преобразованием Фурье. Функция F называется спектральной функцией
или спектральной плотностью сигнала f ( t ) . Функция F называется амплитудным спектром, функция     arg F называется фазовым спектром функции
f (t ) .
F и   - спектральные характеристики сигнала f ( t ) соответственно амплитудная и фазовая F - четная, а   - нечетная функция.
F  F  exp j  .
Рассматривая интеграл Фурье как предельный случай ряда Фурье при    ,
можно заметить, что спектральные линии в пределе сливаются. Поэтому амплитудный спектр непериодической функции будет сплошным и его изображают непрерывной линией.
Формулы (24) и (25) показывают, что если известна спектральная плотность
сигнала f ( t ) , то можно восстановить сигнал f ( t ) , и, наоборот, по известному сигналу f ( t ) можно определить его спектральные характеристики. Таким образом,
описания процессов временными функциями (сигналами) и спектральными функциями равноправны. При решении конкретных задач, связанных с распространением сигналов, используют ту или иную форму представления, исходя из простоты
математического анализа.
Рассмотрим важные для практики примеры нахождения спектральной плотности и спектральных характеристик непериодических сигналов.
28
1. Единичная функция t  изображается графиком, как показано на рис. 16.
Единичная функция t  определяется сле0, t  0,
дующим образом: t   
. Если по1, t  0.
пытаться вычислить спектральную плотность единичной функции t  “напрямую”,
возникает затруднение, связанной с тем, что
эта функция не является абсолютно интегрируемой.
( t )
1
0
t
Рис. 16
В этом случае умножают заданную функцию на затухающую экспоненту
e   t ,   0 . Вычислив спектральную плотность функции t e   t , искомую спектральную плотность находят предельным переходом при   0 .

F  lim   e
 0
 
e
 j 

d  lim  e   j   d 
 0

0

1
1
1
e   j    lim

.
 0   j 
 0   j 
j
0
 lim
Таким образом, амплитудная характеристика единичной функции F 
1

изображается графиком, как показано на рис. 17.
F()
F() 
1


0
Рис. 17
2. Прямоугольный импульс.

h ,
Сигнал, определяемый выражением: S( t )  
0,


t  ,
2

t  ,
2
находит широкое распространение как в технике, так и в теории сигналов и цепей.
Прямоугольный импульс высотой h , длительностью  изображен на рис. 18.
29
S( t )
h
-/2
0
/2
Рис. 18
Применяя формулу (24), находим спектральную плотность этого импульса.

2
F  h  e  j   dt  
2
h 
e
j  
j 

2
e
j 
2
  
 sin 2 
 2h

  sin
.
 h 
 
 
2



 2 
Далее на рис. 19 представлены графики спектральной плотности, амплитудного
и фазового спектров прямоугольного импульса. На графике фазового спектра каждая перемена знака F учитывается приращением фазы на  .
()
F()
F()
h
h
4

2

8

6

2


 -
2

6


Рис. 19
3. Треугольный импульс.
s( t )
h
-/2
0
/2
t
Рис. 20
График функции представлен на рис. 20.
30
 
h 1 
 
 
 

S( t )  h 1 

 

0,



t 

,   t  0,
 
2
2
t 

, 0t ,
 
2
2


t  ,t   .
2
2
Решение. Вычисляем спектральную плотность F .


2 
t   j  t
t   j  t


F   h 1 
e
dt   h 1 
e
dt 


 
 
0
2


2
2
0

0
e  j  t 0 2h e  j  t
e  j t 2


 h


 j t  1
h

 j    2   j 2
 j 0
2
 


sin
 j t


2
2h e
h
4

 .
 j t  1  ... 




  j 2
2 
0



4


График спектральной плотности F() изображен на рис. 21.
2
F()
h
2
4

0
8


Рис. 21
4. Колоколообразный импульс.

t2
S( t )  A  e 2a ,    t   . Этот импульс совпадает по форме с графиком нормального (гауссовского) закона распределения вероятностей и называется также
гауссовским импульсом. Колоколообразный импульс и его спектральная плотность
изображены на рис. 22.
2
F()
S( t )
0

t
Рис. 22
31
Будем находить спектральную плотность данного импульса. По формуле (24)
имеем

F   A  e

t2
2a 2
 e  j  t d .

Для вычисления интеграла удобно в подынтегральной функции дополнить показатель степени до квадрата суммы
2
 t

 t 2
 t2


2
2
  2  j  t     2  j  t  d   d    
 d  d2  ,

 2a


 2a

 2a

где величина d определяется из условия
j t  2
j a
t
.
 d , т.е. d 
2
2a
Таким образом, выражение для F приводится к виду
F  A e

d2
e
 t


 d 
 2a 
2
dt .

Перейдем к новой переменной x 
t
 d , получим
2a
F  A e d
2

2
2 a  e  x dx .


Так как  e  x dx   , то окончательно F  A 2 a e
2


a 2 2
2
 Be

2
2 b2
1
, где b  ,
a
B  2 a A .
Полученный результат имеет важное значение для теории сигналов. Оказывается, что гауссовский импульс и его спектр выражаются одинаковыми функциями и
обладают свойством симметрии: для получения одной из них по заданной другой
достаточно совершить замену t на  и наоборот.
32
5. Волновой цуг. Так называют функцию, определяемую равенством:

h  sin 0 t ,

S( t )  
0,

S( t )
 / 20
t
t 

,
0
t 

.
0
Рис. 23
График функции представлен на рис. 23.
Рассматриваемый сигнал играет в теории связи большую роль. Находим его спектральную плотность.

0
F  h 


sin 0 t e
0
 j t
0
e  j t


dt  h

j


sin

t


cos

t

0
0
0
 j2  0 2

0
j 
  j 
  2 h 0 j
h
  2 h 0 j

0
0
 2
e



e




sin


sin
.


0
0
2
2
0  2  0 2
0
0  2 
 0  
ПРИМЕРЫ ДЛЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ
Найти спектральные плотности и спектральные характеристики следующих непериодических сигналов.
1. Косинусоидальный импульс.
t


h cos  , t  2 ,
S( t )  

0,
t

.

2
Ответ: F 

2 .
2
cos
2 h

  
2
  

33
2. Экспоненциальный импульс.
e  a t ,
S( t )  
0,
t  0,
t  0.
1
Ответ: F 
a 2  2
.
3. Линейно-экспоненциальный импульс.
h t e  h t ,
S( t )  
0,
t  0,
t  0.
2h 
 j arctg
h
h 2  2
e
.
Ответ: F  2
2
h 
2

sin t, t  ,
4. f ( t )  
t  .

0,
Ответ: F 
4 sin  
.
 4  2 

sin 2 m t, t  ,
5. f ( t )  
m  натуральное число.
0
,
t


,


Ответ: F 
4m j
 sin  .
4m 2  2
34
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Гоноровский И.С. Радиотехнические цепи и сигналы. – М.: Радио и связь,
1986. – 512 с.
2. Колобов А.М. Избранные главы высшей математики. – Ч. 1. Ряд Фурье. Интеграл Фурье. Операционные исчисления. – Минск: Высшая школа, 1985. – 220 с.
3. Ефимов А.В. Математический анализ (специальные разделы) – Ч. 1. Общие
функциональные ряды и их приложение. – М.: Высшая школа, 1980. – 279 с.
4. Будак Б.М., Фомин С.В. Кратные интегралы и ряды. – М.: Наука, 1965. –
– 607 с.
35
36
Download