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

Основные методы решения практических задач в курсе «Уравнения математической физики». Учебное пособие

.pdf
Скачиваний:
0
Добавлен:
12.08.2026
Размер:
1 Мб
Скачать

§ 1.. Преобразование Фурье и его свойства

161

Рассмотрим еще интеграл

w

Φ(z) = f(τ) sin z− x) dτ.

−∞

Функция Φ(z) нечетная для любого x по z и хотя, как можно показать, непрерывна по z, но не обязательно интегрируема на (−∞, ∞), однако интеграл от Φ(z) в смысле главного значения по Коши сходится и равен нулю:

 

 

 

 

 

 

 

 

 

 

 

v.p. w

Φ(z) dz = 0,

 

или

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

∞ ∞

 

 

 

 

 

 

1

w w

f(τ) sin z

x) dτ dz = 0.

(1.2)

 

 

 

 

 

 

 

 

 

 

−∞ −∞

 

 

 

 

 

Напомним, что

 

 

 

 

 

 

 

 

 

 

 

N

 

 

 

v.p. w

f(x) dx = lim

w

f(x) dx.

 

 

 

−∞

 

N→∞

−N

 

 

 

 

 

 

 

 

 

Умножая (1.2) на i и складывая с (1.1), получим интегральную формулу Фурье в комплексном виде:

 

 

f(x) = 21π w

f(τ)eiz−x) dτ dz.

(1.3)

−∞

Формула (1.3) позволяет построить прямое и обратное интегральные преобразования Фурье. Прямое преобразование функции f(x) — обозначим его F (z) — записывается так:

 

 

F (z) = 1 w

f(τ)eiτz dτ.

(1.4)

−∞

Возможны обозначения F (z) = fb(z) = V f и др. Из формулы (1.3) следует, что выражение

 

 

 

 

V 1F =

1

w

F (z)e−izx dz = f(x)

 

 

−∞

 

162

Глава V. Метод интегральных преобразований

является обратным преобразование Фурье.

Исходную функцию f(x) здесь и далее называют оригиналом или прообразом преобразования, F (z) — изображением или образом.

Перечислим некоторые свойства преобразования Фурье:

1. Если f(x) абсолютно интегрируема на (−∞, ∞), то F (z) непрерывна в каждой точке оси и стремится к нулю при z → ±∞ [12, с. 637].

2.Преобразование Фурье линейно и однородно, т. е. для любых чисел α,

βвыполняется равенство

V f + βg) = αV f + βV g.

3. Преобразование производных. Заметим, что в теории дифференциальных уравнений применяются интегральные преобразования, которые операцию взятия производных или определенный дифференциальный оператор переводят в более простую операцию. Без доказательства покажем, как работает преобразование Фурье.

Пусть f(x) и f(x) непрерывны и абсолютно интегрируемы на (−∞, ∞),

тогда

V (f(x)) = (−iz)V (f).

Если еще и f′′(x) удовлетворяет перечисленным выше требованиям, то

V(f′′(x)) = (−iz)2V (f) = −z2V (f).

4.Преобразование интеграла. Пусть f(x) абсолютно интегрируемая на

( ,

 

) функция, g(x) = wx f(ξ) dξ, тогда V (g) =

i V (f).

−∞

0

z

 

 

 

 

 

5. Сдвиг в оригинале:

V(f(x − ξ)) = eizξV (f).

6.Сдвиг в изображении:

V (e−ixξf(x)) = F (z − ξ).

Как следствие, получаем формулы

V (cos xξf(x)) = 12(F (z + ξ) + F (z − ξ)),

V (sin xξf(x)) = 21i(F (z + ξ) − F (z − ξ)).

§ 2.. Задача о теплопроводности бесконечного стержня

163

7. Преобразование свертки. Теоретически каждое интегральное преобразование во множестве преобразующих функций (оригиналов) имеет внутреннюю операцию, называемую сверткой, не всегда простую. Для преобразования Фурье свертка определяется так.

Определение 1.1. Пусть f(x) и g(x) абсолютно интегрируемы на (−∞, ∞). Под сверткой этих функций (обозначается f g) понимается следующая операция:

 

 

 

 

 

 

 

f g =

1

w

f(ξ)g(x − ξ) dξ =

1

w

f(x − t)g(t) dt = h(x).

 

 

−∞

 

 

−∞

 

Заметим, что h(x) также абсолютно интегрируема по оси x и

f g = g f.

Имеет место следующее свойство. Преобразование Фурье свертки равняется произведению преобразований компонентов, т. е.

V (f g) = V (f) · V (g).

§2. Задача о теплопроводности бесконечного стержня

Вкачестве иллюстрации рассмотрим задачу об интегрировании уравнения теплопроводности для неограниченного в оба конца тонкого однородного стержня с теплоизолированной боковой поверхностью.

Пример 2.1. В области −∞ < x < ∞, t > 0 найти решение уравнения

∂u

2u

 

(2.1)

 

= a2

 

,

 

∂x2

∂t

 

 

удовлетворяющее условию

u(x, 0) = φ(x).

(2.2)

Решение. Преобразуем по Фурье обе части уравнения (2.1) по переменной x, полагая V (u) = ub(z, t). Искомая функция u(x, t) обладает некоторыми предварительными свойствами, не связанными условиями вывода уравнения и физических особенностей решения. Так, автоматически

= −a2z2ub(z, t).

164 Глава V. Метод интегральных преобразований

u(x, t),

∂u

,

2u

,

∂u

непрерывны по x и t для x (−∞, ∞) и t > 0.

∂x

∂x2

∂t

Приходится считать, что lim u =

lim

∂u

= 0 для t > 0. Используя

 

 

 

 

 

 

 

x→±∞

x→±∞ ∂x

 

возможность дифференцирования несобственных интегралов по t, можно доказать справедливость равенства

()

 

∂u

 

 

V

 

=

 

(V (u)) =

 

ub(z, t).

∂t

∂t

∂t

Преобразованное по Фурье уравнение (2.1) примет вид

V (∂u) = a2V (2u), ∂t ∂x2

или

∂ub(z, t)

∂t

(2.3)

Уравнение (2.3) по существу считаем обыкновенным дифференциальным уравнением с параметром z, а решение ищется, удовлетворяющее начальному условию, которое получается преобразованием Фурье равенства (2.2):

V (u(x, 0)) = V (φ), или

ub(z, 0) = φb(z).

(2.4)

Решая уравнение (2.3), находим, что

u z, t

) =

C z e−a2z2t,

(2.5)

b(

( )

где C(z) — неизвестная функция параметра z. Полагая в (2.5) t = 0, получаем, что C(z) = φb(z).

В терминах образов преобразования Фурье получаем

ub(z, t) = φb(z)e−a2z2t.

Искомую функцию находим обратным преобразованием Фурье. Это можно сделать непосредственными вычислениями, но мы воспользуемся свойством свертки в обратную сторону, записав

u(x, t) = V 1(ub) = φ(x) g(x),

§ 2.. Задача о теплопроводности бесконечного стержня

165

где g(x) = V 1(e−a2z2t). Вычисляем:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

g(x) =

V 1(e−a2z2t) =

1

w

e−a2z2te−izx dz =

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

1

2

2

 

 

 

2

2

 

=

 

 

( w

e−a

z

t cos zx dx − i w

e−a

z

t sin zx dz).

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

−∞

 

 

 

Второй интеграл от нечетной функции равен 0. Сходимость обоих инте-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

гралов очевидна. В первом интеграле сделаем замену: ξ = az

t, z =

a

 

.

t

При этом пределы интегрирования не меняются:

 

 

 

 

 

 

 

 

 

 

 

 

 

g(x) =

 

1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

x .

 

 

 

 

 

 

 

 

 

 

w

eξ2 cos µξ dξ, µ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

t −∞

 

 

 

 

 

 

 

 

 

 

a

 

t

 

 

 

 

 

 

Вычислим методом дифференцирования по µ интеграл

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J(µ) =

w

eξ2 cos µξ dξ,

 

 

 

 

 

 

 

 

(2.6)

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

учитывая, что J(µ) вместе с J(µ) часто встречаются в задачах с исполь-

зованием преобразования Фурье. Имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

2

 

 

 

 

 

J(µ) = w

ξeξ

 

sin µξ dξ = w

sin µξ (

 

deξ

) =

 

 

 

 

2

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

ξ2

 

 

 

 

 

 

1

 

ξ2

 

 

 

 

 

 

1

 

 

 

 

 

 

=

 

e

 

sin µξ −∞

 

 

w

µe

 

cos µξ dξ

=

 

µJ(µ).

2

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Внеинтегральное выражение, очевидно, равно нулю. Таким образом, J(µ)

удовлетворяет уравнению J(µ) =

1

µJ(µ), интегрируя которое находим

2

J(µ) = Ceµ2/4. Постоянную C определяем из равенства

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

J(0) = w

 

eξ2 dξ =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

 

 

 

 

 

 

 

 

 

 

 

 

−∞

как интеграл Пуассона. В итоге J(µ) =

 

eµ2/4. А так как µ =

π

g x

V 1 e−a2z2t

1

1

 

 

πe4a2t

 

1 e4a2t .

( ) =

(

) =

 

·

a

 

 

 

 

 

x2

=

a

 

 

x2

 

 

 

 

 

 

 

 

 

t

2t

x

, то a t

166

 

 

 

 

 

Глава V. Метод интегральных преобразований

Искомая функция равна

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

x2

1

 

 

−x)2

 

u(x, t) = φ(x)

 

 

 

e

 

=

 

 

w

φ(τ)e

 

dτ.

(2.7)

 

 

 

4a2t

 

4a2t

 

a

 

2t

 

 

 

2a

 

πt −∞

 

 

 

 

§ 3. Синус-, косинус-преобразования Фурье

Названные преобразования используются при решении начальнокраевых задач для уравнения теплопроводности в случае полубесконечного стержня, для определенности x (0, ∞), при таких краевых усло-

виях u(0, t) = µ(t), либо ∂u∂x(0, t) = ν(t). Начальное условие вида (2.2): u(x, 0) = φ(x).

Если краевые условия однородны, то поставленная задача обычно сводится к случаю x (−∞, ∞), т. е. рассмотренному ранее. Если u(0, t) = 0, то искомую функцию u(x, t) продолжают для x < 0 нечетным образом, по-

лагая

 

 

 

u˜(x, t) = {

u(x, t), x > 0,

φ˜(x) = {

φ(x), x > 0,

−u(−x, t), x < 0,

φ(−x), x < 0.

u˜(x, t) находится по формуле (2.6), в которой u(x, t) = u˜(x, t) при x > 0. Если же ∂u∂x(0, t) = 0, то u(x, t) продолжается четным образом:

{

u˜(x, t) =

u(x, t), x > 0, u(−x, t), x < 0.

Пусть краевые условия неоднородные. Рассмотрим варианты

u(0, t) = µ(t) и ∂u(0, t) = ν(t), ∂x

а u(x, t) является решением уравнения (2.1) в области x > 0, t > 0. К решению подобных задач привлекают косинусили синус-преобразования Фурье.

§ 3.. Синус-, косинус-преобразования Фурье

167

3.1. Косинус-преобразование Фурье

Полагаем в интегральной формуле (1.1) f(x) четной функцией, после несложных рассуждений получаем

 

 

 

 

 

 

 

∞ ∞

 

 

 

 

 

f(x) =

2

w w f(τ) cos zτ cos zx dτ dz.

 

 

π

 

 

0

0

 

 

 

Из последнего равенства выводим прямое преобразование

 

 

 

 

 

 

 

π

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fc(z) =

2

w f(τ) cos zτ dτ = Vc(f)

(3.1)

и обратное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

w Fc(z) cos zx dz =

f(x), x > 0;

 

 

 

 

 

{ f(−x), x < 0.

 

 

π 0

 

 

 

 

Приведем без вывода формулу преобразования второй производной f′′(x) в предположении, что lim f(x) = lim f(x) = 0:

x→∞

x→∞

 

 

 

 

 

 

 

 

 

2

 

 

 

Vc(f′′) =

 

f(0) − z2Vc(f).

(3.2)

π

3.2. Синус-преобразование Фурье

Полагаем теперь в формуле (1.1) f(x) нечетной на (−∞, ∞) функцией. Имеем

 

 

 

 

 

 

∞ ∞

 

 

 

 

f(x) =

2

w w f(τ) sin zτ sin zx dτ dz.

π

 

 

 

 

 

0

0

 

 

Прямое преобразование введем равенством

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Fs(z) =

 

 

 

2

w f(τ) sin zτ dτ = Vs(f),

 

 

 

 

 

 

 

 

 

π

0

 

 

обратное запишется так:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

f(x),

x > 0,

 

w Fs(z) sin zx dz = {

−f(−x),

x < 0.

π

0

ν(τ)ea2z2τ dτ + C(z) .
U(z, t) = e−a2z2t
)

168

 

 

Глава V. Метод интегральных преобразований

 

Построим преобразование f′′(x):

 

 

 

 

 

 

 

 

 

2

 

 

 

 

Vs(f′′) =

 

· z · f(0) − z2Vs(f).

(3.3)

 

π

Замечание. Исходя из формул (3.2), (3.3), при краевом условии u(0, t) = = µ(t) следует применять синус-преобразование, а если ∂u∂x(0, t) = ν(t), — то косинус-преобразование. Первые слагаемые в (3.2), (3.3) должны быть известны.

Пример 3.1. В области x > 0, t > 0 найти решение уравнения (2.1)

∂u

2u

, удовлетворяющие условиям

 

 

 

= a2

 

 

 

 

∂x2

 

 

∂t

 

 

 

 

 

 

 

 

 

 

∂u

x=0= ν(t),

(3.4)

 

 

 

 

 

 

 

 

 

∂x

 

 

 

 

 

 

 

 

(3.5)

 

 

 

 

u(x, 0) = 0,

 

 

 

 

 

lim u = lim

∂u

= 0.

(3.6)

 

 

 

 

 

 

 

x→∞

x→∞ ∂x

 

 

Решение. Учитывая замечание, применяем косинус-преобразование Фурье. Преобразуем уравнение (2.1), обозначив Vc(u) = U(z, t):

 

 

 

∂u

 

 

 

 

 

 

 

2u

 

 

 

Vc (

 

)

= a2Vc (

 

).

 

∂t

∂x2

 

Далее, с учетом равенств (3.2) и (3.4) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂U

 

2 ∂u

 

 

 

 

 

 

 

 

 

 

 

 

(z, t) = a2 (

 

 

 

 

 

 

x=0−z2U(z, t))

,

 

∂t

 

 

π

 

∂x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂U

2

 

 

= −a2z2U − a2

 

· ν(t).

∂t

π

Полученное уравнение имеет общее решение вида

(−a22 wt

π

(3.7)

(3.8)

0

§ 4.. Преобразование Лапласа

169

Из начального условия (3.5) после его преобразования получаем начальное требование для (3.7) Vc(u(x, 0)) = U(z, 0) = 0, которое удовлетворяется при C(z) = 0. Таким образом,

 

 

 

 

 

 

 

 

 

wt ν(τ)e−a2z2(t−τ) dτ.

U(z, t) =

a2

 

2

 

 

 

 

 

 

 

 

π

0

 

Искомое решение находим обратным преобразованием:

 

a2

t

 

 

 

 

2 2

 

u(x, t) =

2

 

w (w

ν(τ)e−a z

(t−τ) dτ)cos zx dz =

π

00

=2πa2 wt ν(τ)(we−a2z2(t−τ) cos zx dz)dτ.

00

Внутренний интеграл заменой z =

 

ξ

 

 

 

сводится к J(µ). Равенство

(2.6) при µ =

x

примет вид

 

 

 

a t − τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a t − τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

w e−a2z2(τ) cos zx dz =

 

1

 

 

w eξ2 cos µξ dξ =

0

 

 

 

 

 

 

 

 

 

 

a

 

t − τ 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

1

 

w

eξ2 cos µξ dξ =

 

 

 

 

1

 

e

 

 

=

 

 

 

 

π

4a2(t−τ) .

 

2a t − τ −∞

 

 

 

 

 

 

 

 

 

 

2a t − τ

 

 

Окончательно

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

u(x, t) =

 

a w

 

ν(τ)

e

 

dτ.

 

 

 

 

 

 

 

 

4a2(t−τ)

 

 

 

 

 

 

 

 

 

 

−√

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

π

 

t − τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

§ 4. Преобразование Лапласа

4.1. Функция-оригинал. Функция-изображение

Определим класс функций, на котором можно построить интегральное преобразование Лапласа. В качестве функций-оригиналов определим множество комплекснозначных функций f(x) вещественного переменного t, −∞ < t < ∞, удовлетворяющих условиям:

170

Глава V. Метод интегральных преобразований

1.f(t) однозначна и кусочно непрерывна для любого t, −∞ < t < ∞.

2.f(t) 0 для t < 0.

3.При t → +функция f(t) растет не быстрее показательной функции. Это означает, что найдутся M > 0 и a > 0, при которых имеет место оценка

|f(t)| 6 Meat, t > 0.

(4.1)

В конкретных задачах условие 1 может быть усилено требованием существования у f(t) производных до определенного порядка.

Преобразованием Лапласа функции f(t) назовем несобственный интеграл с комплексным параметром p = x + iy

 

F (p) = w e−ptf(t) dt.

(4.2)

0

 

Функцию F (p), естественно, назовем функцией изображения или образом преобразования, тогда f(t) будем называть прообразом. Требования на f(t) обеспечивают существование F (p) в некоторой области комплексной плоскости p. Для обозначения операции преобразования воспользуемся вариантом

f(t) : F (p).

Для F (p) имеет место приведенная ниже теорема.

Теорема 4.1. Пусть f(t) удовлетворяет неравенству (4.1), тогда F (p) определена на множестве Re p > x0 > a, непрерывна и, более того, аналитическая по p в этой области.

Приведенная теорема подчеркивает важность постоянной a в оценке (4.1), и в ряде случаев удается указать точную границу полуплоскости определения F (p).

Например, для оригиналов со степенным ростом (tν, ν > 0) F (p) аналитическая в полуплоскости Re p > 0, если f(t) = e(α+iβ)t, — то в области

Re p > +α.

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