Интеграл Фурье и его приложения. Учебное пособие
.pdfF(ΙΙ) ( ) =e−i x0 f |
′(x ) −f′(0) +i F(I) ( ) = |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
= |
|
A |
e−i x0 − |
A |
+2i |
|
A |
|
|
|
x |
|
|
− |
i x0 |
= |
||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
sin |
|
|
0 |
e |
|
2 |
||||||||||||
x |
x |
|
x |
2 |
|
|
||||||||||||||||||
|
|
|
0 |
|
|
|
0 |
|
|
0 |
|
|
|
|
) |
|
|
|
|
|
|
|
||
|
A |
|
|
−i x0 |
|
|
( |
|
i x0 |
|
− |
i x0 |
|
− |
i x0 |
|
||||||||
= |
|
|
−1 |
+ |
|
2 |
|
−e |
|
2 |
|
|
2 |
|
=0. |
|||||||||
|
|
e |
|
|
e |
|
|
|
|
|
|
e |
|
|
|
|
|
|||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Как и должно быть!
Упражнения
1.Получить выражения автокорреляционной функции одиночного импульса прямоугольной формы (1.35) в явном виде и в виде разложения в интеграл Фурье. Записать равенство Парсеваля. Сделать чертежи.
2.Получить выражения автокорреляционной функции гауссова импульса (1.44) в явном виде и в виде разложения в интеграл Фурье. Записать равенство Парсеваля. Сделать чертежи.
3.Пусть f(x) – финитная функция, т.е. функция, отличная от нуля лишь в пределах некоторого отрезка [a, b]. Записать две формы явного выражения ее автокорреляционной функции. Показать четность автокорреляционной функции. Представить автокорреляционную функцию в виде интеграла Фурье. Записать равенство Парсеваля.
4.Пусть f(x) и g(x) – финитные функции, определенные на одном и том же отрезке [a, b]. Записать две формы явного выражения корреляции функций f(x) и g(x). Представить эту корреляцию в виде интеграла Фурье. Определить отрезок, на котором корреляция функций f(x) и g(x) отлична от нуля.
5.Пусть f(x) – финитная функция, f′(x) – ее производная,
тоже финитная. Записать различные формы явного выражения корреляции функции f(x) и ее производной f′(x). Представить эту корреляцию в виде интеграла Фурье и вычислить его
вточке x = 0. Записать различные формы явного выражения автокорреляционной функции производной f′(x). Получить
выражение спектра автокорреляционной функции производной f′(x) через спектр самой функции f(x) (с учетом четности). Записать равенство Парсеваля.
21
1.6. Еще о некоторых свойствах спектральных плотностей финитных функций
Пусть |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1, |
|
x |
|
|
|
<a |
|
|||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ηa (x) = |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=a; |
|
|||||
|
|
|
|
|
, |
|
|
|
|
|
|
|
x |
|
(1.55) |
|||||||||||||
2 |
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
0, |
|
|
|
x |
|
>a |
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1, |
|
|
x |
|
|
|
<b ≤a |
|
|||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ηb (x) = |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
=b . |
|
||||||||
|
|
|
, |
|
|
|
|
x |
|
(1.56) |
||||||||||||||||||
2 |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0, |
|
x |
|
>b |
|
||||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||||||||||
Тогда ηb (x) =ηa (x)ηb (x) и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
a |
e−i xdx = 2 sin a ; |
|
|||||||||||||||||||||||||
Na ( ) = |
∫ |
(1.57) |
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
η (x) = |
1 |
∞ sin a ei x dx. |
(1.58) |
|||||||||||||||||||||||||
a |
π |
∫ |
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−∞
С заменой здесь буквы а на букву b записываются такие же формулы для функции (1.56). С учетом того, что b ≤ a, имеем
bb
∫e−i xdx = ∫ηa (x)e−i xdx =
−b |
|
|
|
|
|
|
−b |
|
|
|
|
|
|
|
|
|
1 |
b |
∞ |
|
sina ′ |
|
|
|
|
||||||
= |
|
∫ |
|
∫ |
|
|
|
|
|
|
ei xdx′ e−i xdx = |
||||
π |
|
′ |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
−b |
|
−∞ |
|
|
|
|
|
|
|
|
|
||
|
1 |
∞ |
|
sina ′ |
b |
|
−i( − ′)x |
|
′ |
||||||
= |
|
∫ |
|
|
|
|
|
|
|
∫ |
e |
|
dx d = |
||
π |
|
|
|
′ |
|
||||||||||
|
|
|
|
|
|
|
|
|
|||||||
|
|
−∞ |
|
|
|
|
|
|
|
−b |
|
|
|
|
|
=2 ∞∫ sina ′ sinb(−′) ′ π −∞ ′ −′ d .
22
Таким образом:
1 |
∞ |
sina ′ sinb( − ′) |
|
sinb |
|
||
|
∫ |
|
|
|
d ′= |
|
. |
π |
′ |
− ′ |
|
||||
|
−∞ |
|
|
|
|
|
|
Положим
=l−l′,
тогда последнее равенство может быть записано так:
1 |
∞ |
sina(l−x) sinb(l′−x) |
|
sinb(l−l′) |
|
|
||
|
∫ |
|
|
|
dx = |
|
. |
(1.59) |
π |
l−x |
l′−x |
l−l′ |
|||||
|
−∞ |
|
|
|
|
|
|
|
Обратим внимание на то, что интегралы в левых частях последних равенств не зависят от величины параметра а, лишь бы выполнялось неравенство a ≥ b > 0. Более того, после замены а = b, u = al, u′ = al′, s = ax этот параметр вообще сокращается и равенство (1.59) приобретает вид
1 |
∞ |
sin(u −s) sin(u′−s) |
|
sin(u −u′) |
|
|
||
|
∫ |
|
|
|
ds = |
|
. |
(1.60) |
π |
u −s |
u′−s |
u −u′ |
|||||
|
−∞ |
|
|
|
|
|
|
|
Если здесь аргументам u и u′ присвоить дискретные значения:
u = πn; u′= πm,
где n, m – целые числа,
то получим равенство
1 |
∞ |
sin(πn −s) sin(πm −s) |
1, n =m |
|
||
|
∫ |
|
|
|
ds = 0, n ≠m, |
(1.61) |
π |
πn −s |
|
πm −s |
|||
|
−∞ |
|
|
|
|
|
говорящее о том, что функции
|
1 |
|
sin(πn −s) |
(1.62) |
||
|
|
|
|
|
|
|
|
|
|
πn −s |
|||
|
|
π |
|
|
|
|
образуют ортонормированную систему во всем бесконечном промежутке –∞ < s < ∞.
23
Найдем представления коэффициентов разложения некоторой функции F(s) по этой системе. Имеем
F(s) =
где
A(n) =
|
1 |
|
|
∞ |
sin(πn −s) |
|
||||
|
|
|
∑ An |
, |
||||||
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
||||||
|
|
|
π n=−∞ |
πn −s |
||||||
|
|
|
|
|
|
∞ |
|
|
|
|
|
1 |
|
|
∫ F(s) |
sin(πn −s) |
ds. |
||||
|
|
|
|
|
|
|
||||
|
|
π |
−∞ |
πn −s |
||||||
(1.63)
(1.64)
Пусть F(s) есть спектральная плотность некоторой функции f(x). Другими словами, пусть имеет место пара преобразований Фурье:
|
|
|
|
|
|
|
|
|
|
|
|
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F(s) = ∫ f(x)e−isxdx; |
|
(1.65) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f(x) = |
|
∫ F(s)eisxds. |
|
(1.66) |
|||||||||||
|
|
|
|
|
|
|
2π |
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
Вносим интеграл (1.65) в интеграл (1.64): |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
1 |
|
∞ |
∞ |
|
|
|
|
−isx |
sin (πn − s) |
|
||||||
A(n) = |
|
|
|
|
|
∫ |
∫ |
f |
(x)e |
|
dx |
|
|
|
πn − s |
ds = |
|||||||
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
π |
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
−∞ −∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
∞ |
|
|
∞ |
sin (πn − s) |
|
|
−isx |
|
|
||||||||
= |
|
|
|
|
|
|
∫ |
f (x) |
∫ |
|
|
πn − s |
|
e |
|
ds dx = |
|||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
π |
−∞ |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
||||||
|
|
|
/ = s − πn, s = πn + , ds = d / |
|
|||||||||||||||||||
|
|
1 |
|
|
|
∞ |
|
|
−iπnx |
∞ sin |
|
−i x |
|
|
|||||||||
= |
|
|
|
|
|
|
∫ f (x)e |
|
|
|
|
∫ |
|
e |
|
|
d dx. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
π |
−∞ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
||||
В силу (1.58) получаем
1 |
|
||
A(n) = |
|
∫f(x)e−iπnxdx. |
|
π |
(1.67) |
||
|
|
−1 |
|
24
Что мы видим? Разложение (1.63) будет однозначным, если функция F(s) есть интеграл Фурье от финитной функции с носителем x [–1, 1], т.е. если
1 |
|
F(s) = ∫f(x)e−isxdx. |
(1.68) |
−1 |
|
В этом случае коэффициенты в разложении (1.63) такие:
|
1 |
|
∞ |
|
sin(πn −s) |
|
|
|
||||
|
|
∫ |
|
|
|
|
||||||
A(n) = |
|
|
|
F(s) |
|
|
|
ds = πF(πn). |
(1.69) |
|||
|
|
|
πn |
−s |
|
|||||||
|
π |
|||||||||||
|
|
|
−∞ |
|
|
|
|
|
|
|||
Внося (1.69) в (1.63), запишем |
|
|
|
|
|
|||||||
|
|
|
|
|
∞ |
|
sin(πn −s) |
|
|
|
||
|
F(s) = ∑ F(n) |
. |
|
(1.70) |
||||||||
|
|
|
|
|
||||||||
|
|
πn −s |
|
|||||||||
|
|
|
|
|
n=−∞ |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
Преобразование (1.69) само по себе тоже примечательно. Заменяя в нем аргумент pn на непрерывную переменную s, перепишем его так:
|
1 |
∞ |
sin( −s) |
|
|
F(s) = |
|
∫ F( ) |
|
d . |
(1.71) |
π |
−s |
||||
|
|
−∞ |
|
|
|
Итак, для функции F(s), представленной интегралом (1.68), имеем два, по сути дела, тождественных преобразования (1.70) и (1.71).
Если рассмотреть финитную функцию fl(x) на отрезке [–l, l] произвольной длины, т.е. если вместо (1.68) взять преобразование
l |
|
Fl (s) = ∫fl (x)e−isxdx, |
(1.72) |
−l
то разложения (1.70), (1.71) надлежащими заменами приводятся к разложениям
∞ |
π |
sin(πn −ls) |
|
|
||
Fl (s) = ∑ Fl |
n |
|
; |
(1.73) |
||
πn −ls |
||||||
n=−∞ |
l |
|
|
|
||
|
|
|
|
|
||
25
F (s) = |
1 |
∞ |
F ( ) |
sinl( −s) |
d . |
(1.74) |
|
∫ |
|
||||
l |
π |
l |
−s |
|
||
|
|
−∞ |
|
|
|
|
Еще важно отметить, что если в преобразовании (1.72) перейти к пределу при l → ∞, то формула (1.74) сохранится. При этом имеют место следующие предельные переходы:
|
|
|
|
|
|
|
|
1, u >s |
|
|||
|
|
|
|
1 |
u |
sinl( −s) |
|
|
|
|||
|
|
|
|
|
1 |
|
|
|||||
lim |
|
|
|
|
d = |
|
, u = s; |
|
||||
π−∞∫ |
−s |
|
|
|||||||||
l→∞ |
2 |
|
|
|||||||||
|
|
|
|
|
|
|
|
0, u < s |
|
|||
|
|
|
|
|
|
|
|
|
|
|
||
|
1 sinl( −s) |
|
|
0, ≠ s |
|
|||||||
lim |
|
|
|
|
|
|
=d( −s) = |
. |
||||
π −s |
||||||||||||
l→∞ |
|
|
∞, =s |
|
||||||||
(1.75)
(1.76)
Этот предел следует понимать в обобщенном смысле, и он представляет собой обобщенную d-функцию, упоминавшуюся в предыдущем пункте.
Посредством d-функции равенство (1.74) при предельном переходе l → ∞ запишется так:
∞ |
|
∫ F( )d( −s)d = F(s). |
(1.77) |
−∞
Это равенство считается справедливым для любой непрерывной функции F(s).
Упражнения
Записать разложения (1.73), (1.74) для спектров: а) гауссова импульса (1.42);
б) двух одинаковых прямоугольных импульсов противоположного знака (1.48);
в) треугольного импульса (1.50).
26
Глава 2 Решение обыкновенных
линейных дифференциальных уравнений
2.1. Линейное уравнение первого порядка
Пусть дано линейное однородное уравнение первого порядка с постоянным положительным коэффициентом
u′+au =0. |
(2.1) |
Его аналитическое решение на всем бесконечном промежутке от –∞ до +∞ таково:
u =u(0)e−ax. |
(2.2) |
Применение преобразования Фурье позволяет записать решение этого уравнения в виде финитной функции, т.е. в виде отличного от нуля решения лишь на некотором конечном отрезке [0, l]. Наверное, можно применить и термин «кусочное решение».
Умножаем обе части уравнения (2.1) на экспоненту e
и интегрируем его по отрезку [0, l]. Эта процедура переводит уравнение (2.1) в его алгебраический Фурье-образ:
i U( ) −u(0) +e−i lu(l) +aU( ) =0.
Здесь
l
U( ) = ∫u(x)e−i xdx.
0
Разрешая (2.3) относительно U( ), записываем
= u(0) −e−i lu(l) U( ) .
i +a
(2.3)
(2.4)
(2.5)
По спектральной плотности U( ) восстанавливаем решение исходного уравнения (2.1):
27
|
1 ∞ |
u(0) −e−i lu(l) |
|
|
u(x) = |
|
∫ |
|
ei xd = |
2π |
i +a |
|||
|
|
−∞ |
|
|
|
1 |
∞ |
ei x |
|
|
|
1 |
∞ |
|||
=u(0) |
∫ |
d −u(l) |
∫ |
||||||||
2πi |
−ia |
2πi |
|||||||||
|
|
−∞ |
|
|
|
|
|
|
|
−∞ |
|
Обозначим |
|
|
|
|
|
|
|
|
|||
|
|
|
|
1 |
∞ |
ei x |
|
|
|||
|
|
g(x) = |
∫ |
d . |
|||||||
|
|
2πi |
−ia |
||||||||
|
|
|
|
|
|
−∞ |
|
|
|
|
|
ei (x−l)
−ia d .
(2.6)
(2.7)
Теперь функция (2.6) записывается так: |
|
u(x) =u(0)g(x) −u(l)g(x −l). |
(2.8) |
Интеграл (2.7) представляет собой разрывную функцию. Его значения известны и легко вычислимы, например, по вычету подынтегральной функции. Имеем
|
|
|
|
e−ax, x >0 |
|
||
|
1 |
∞ |
ei x |
|
1 |
|
|
|
∫ |
|
|
|
|||
g(x) = |
|
|
d = |
|
, x =0 . |
(2.9) |
|
2πi |
−ia |
|
|||||
|
−∞ |
2 |
|
|
|||
|
|
|
|
0, x <0 |
|
||
|
|
|
|
|
|
|
|
В соответствии с этими значениями кусочное решение (2.8) по отдельным промежуткам распределяется так:
x<0, u(x) =0;
x→0 +0, u(0 +0) =u(0); 0 <x <l, u(x) =u(0)e−ax;
x→l −0, u(l −0) =u(l) =u(0)e−al;
x>l, u(x) =u(0)e−ax −u(l)e−a(x−l) =e−ax (u(0) −u(l)eal )=0.
Рассмотрим теперь неоднородное уравнение
u′+au =f(x), (2.10)
у которого правая часть f(x) – финитная функция с носителем x [0, l]. По той же методике представления решения инте-
28
гралом Фурье строим его как финитную функцию с тем же носителем.
Снова умножив обе части уравнения (2.10) на экспоненту e–i x и проинтегрировав его по отрезку [0, l], получим
i U( ) −u(0) +e−i lu(l) +aU( ) = F( ), |
(2.11) |
где |
|
l |
|
U( ) = ∫u(x)e−i xdx; |
(2.12) |
0 |
|
l |
|
F( ) = ∫f(x)e−i xdx. |
(2.13) |
0 |
|
Из (2.11) находим
U( ) = |
u(0) −e−i lu(l) + F( ) |
. |
|
i +a |
|||
|
|
Далее по формуле обращения:
|
|
1 |
|
|
∞ |
u(0) −e−i lu(l) + F( ) |
|
|
||||||||||
u(x) = |
|
|
|
∫ |
|
|
|
|
|
|
|
|
|
ei xd = |
||||
2π |
|
|
|
|
|
i +a |
|
|||||||||||
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
|
|
||
|
1 |
∞ |
|
|
|
|
i x |
|
|
1 |
∞ |
i (x−l) |
||||||
=u(0) |
∫ |
|
|
e |
d −u(l) |
∫ |
e |
d + |
||||||||||
2π |
|
−ia |
2πi |
−ia |
||||||||||||||
|
|
|
−∞ |
|
|
|
|
|
|
|
|
|
−∞ |
|
|
|
||
|
|
|
|
|
|
|
|
|
1 |
|
∞ |
F( ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
∫ |
ei xd . |
|
|
|||||||
|
|
|
|
|
|
|
2πi |
−ia |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
−∞ |
|
|
|
|
|
|
|
|
(2.14)
(2.15)
Внося сюда (2.13) и используя уже введенное обозначение (2.7), записываем данное кусочное решение в виде
u(x) =u(0)g(x) −u(l)g(x −l) +
l |
|
+∫f(x′)g(x −x′)dx′=u(0)g(x) −u(l)g(x −l) + |
(2.16) |
0
x
+∫f(x′)g(x −x′)dx′.
0
29
Здесь тоже в соответствии с определением функции g(x) будет u(x) = 0, когда x [0, l], а на отрезке [0, l] имеем известную формулу
x |
|
u(x) =u(0)e−ax +e−ax ∫f(x′)eax′dx′. |
(2.17) |
0 |
|
Отметим, что выражение решения в виде (2.16) включает как бы две произвольные постоянные u(0) и u(l). На самом деле эти постоянные связаны условием
l |
|
u(l) =u(0)e−al +e−al ∫f(x′)eax′dx′. |
(2.18) |
0 |
|
То есть решение уравнения первого порядка имеет только одну произвольную постоянную, как и должно быть.
2.2. Линейное уравнение второго порядка
Однородное линейное уравнение второго порядка с постоянными коэффициентами записывается так:
u′′+2au′+b2u =0. |
(2.19) |
Его общее решение, аналитическое на всей оси, дается формулой
y =e−ax (C1 coskx +C2 sinkx), |
(2.20) |
||
где |
|
||
|
|
|
|
k = b2 −a2 , |
(2.21) |
||
считаем, что a2 < b2.
С помощью постоянных С1, С2 можно удовлетворить любым граничным и начальным условиям.
Покажем, как работает преобразование Фурье для решения начальных и граничных задач на конечном отрезке [0, l]. Ищем решение уравнения (2.19) в виде интеграла Фурье:
∞
u(x) = 1 ∫U(x)ei xd . (2.22)
2π −∞
30
