Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Интеграл Фурье и его приложения. Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
09.08.2026
Размер:
331 Кб
Скачать

F(ΙΙ) ( ) =ei x0 f

′(x ) −f′(0) +i F(I) ( ) =

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

A

ei 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 (xb (x) и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

ei xdx = 2 sin a ;

 

Na ( ) =

(1.57)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η (x) =

1

sin a ei x dx.

(1.58)

a

π

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

С заменой здесь буквы а на букву b записываются такие же формулы для функции (1.56). С учетом того, что b a, имеем

bb

ei xdx = ηa (x)ei xdx =

b

 

 

 

 

 

 

b

 

 

 

 

 

 

 

 

1

b

 

sina

 

 

 

 

=

 

 

 

 

 

 

 

 

ei xdxei 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)eisxdx;

 

(1.65)

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f(x) =

 

F(s)eisxds.

 

(1.66)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

Вносим интеграл (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)eiπnxdx.

 

π

(1.67)

 

 

−1

 

24

Что мы видим? Разложение (1.63) будет однозначным, если функция F(s) есть интеграл Фурье от финитной функции с носителем x [–1, 1], т.е. если

1

 

F(s) = f(x)eisxdx.

(1.68)

−1

 

В этом случае коэффициенты в разложении (1.63) такие:

 

1

 

 

sin(πn s)

 

 

 

 

 

 

 

 

 

A(n) =

 

 

 

F(s)

 

 

 

ds = πFn).

(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)eisxdx,

(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

i x

Глава 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) +ei lu(l) +aU( ) =0.

Здесь

l

U( ) = u(x)ei xdx.

0

Разрешая (2.3) относительно U( ), записываем

= u(0) −ei lu(l) U( ) .

i +a

(2.3)

(2.4)

(2.5)

По спектральной плотности U( ) восстанавливаем решение исходного уравнения (2.1):

27

 

1

u(0) −ei lu(l)

 

u(x) =

 

 

ei xd =

i +a

 

 

−∞

 

 

 

1

ei x

 

 

 

1

=u(0)

d u(l)

i

ia

i

 

 

−∞

 

 

 

 

 

 

 

−∞

Обозначим

 

 

 

 

 

 

 

 

 

 

 

 

1

ei x

 

 

 

 

g(x) =

d .

 

 

i

ia

 

 

 

 

 

 

−∞

 

 

 

 

ei (xl)

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)

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;

xl −0, u(l −0) =u(l) =u(0)e−al;

x>l, u(x) =u(0)e−ax u(l)e−a(xl) =e−ax (u(0) −u(l)eal )=0.

Рассмотрим теперь неоднородное уравнение

u′+au =f(x), (2.10)

у которого правая часть f(x) – финитная функция с носителем x [0, l]. По той же методике представления решения инте-

28

гралом Фурье строим его как финитную функцию с тем же носителем.

Снова умножив обе части уравнения (2.10) на экспоненту ei x и проинтегрировав его по отрезку [0, l], получим

i U( ) −u(0) +ei lu(l) +aU( ) = F( ),

(2.11)

где

 

l

 

U( ) = u(x)ei xdx;

(2.12)

0

 

l

 

F( ) = f(x)ei xdx.

(2.13)

0

 

Из (2.11) находим

U( ) =

u(0) −ei lu(l) + F( )

.

i +a

 

 

Далее по формуле обращения:

 

 

1

 

 

u(0) −ei lu(l) + F( )

 

 

u(x) =

 

 

 

 

 

 

 

 

 

 

 

 

ei xd =

 

 

 

 

 

i +a

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

i x

 

 

1

i (xl)

=u(0)

 

 

e

d u(l)

e

d +

 

ia

i

ia

 

 

 

−∞

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

1

 

F( )

 

 

 

 

 

 

 

 

 

 

 

 

 

+

ei xd .

 

 

 

 

 

 

 

 

 

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′)eaxdx′.

(2.17)

0

 

Отметим, что выражение решения в виде (2.16) включает как бы две произвольные постоянные u(0) и u(l). На самом деле эти постоянные связаны условием

l

 

u(l) =u(0)e−al +e−al f(x′)eaxdx′.

(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)

−∞

30

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]