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

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

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

§ 3.. Задача Коши для уравнения в частных производных второго порядка

21

Отсюда

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3.8)

 

 

φ(x)

 

ψ(3x) = C.

 

 

 

3

 

 

Из (3.6) и (3.8) находим

 

 

 

 

 

 

 

 

 

 

3

x2 +

3

 

 

9x2

3

 

φ(x) =

 

 

C,

ψ(3x) =

 

 

C.

 

4

4

4

4

 

Наконец, положив 3x = z, получаем

ψ(z) = z2 3C. 4 4

Требуемые функции φ и ψ найдены. Подставляя их в (3.5), мы получим искомое решение задачи:

u(x, y) =

3

(x + y)2

+

3

C +

(y − 3x)2

3

C = 3x2

+ y2.

4

 

4

 

 

 

 

4

4

 

Замечание 3. Осложнения, связанные с неодинаковостью аргументов функций φ и ψ в (3.6) и (3.7), устранятся, и решение поставленной задачи упростится, если произвольную функцию ψ в (3.5) мы представим в виде сложной функции ψ(ω(y − 3x)), где ψ — произвольная, а ω — определен-

ная функция,

причем последняя выбирается так, чтобы ψ(ω(y

3x)) при

e

 

ψ(x)

 

e

 

 

 

 

 

 

 

y = 0

принимала вид

 

 

 

 

 

 

следует положить

 

 

 

 

. В нашем случае, очевидно,

e

 

 

 

 

 

ω(z) =

z/3

 

общее решение уравнения (3.1) будет иметь вид

 

 

. Тогда

 

e

 

 

e (

 

)

 

 

 

e

(

3)

 

 

 

 

 

 

 

 

 

3

 

 

 

 

u(x, y) = φ(x + y) + ψ

y −

3x

=

φ(x + y) + ψ

x

 

 

y

,

 

 

 

 

 

 

а поэтому вместо (3.6) и (3.7) мы получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

 

 

 

2

 

 

 

 

 

 

 

и

 

 

 

 

 

 

 

 

φ(x) + e(x) = 3x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

e

 

 

 

1

e

 

 

 

 

 

 

 

 

 

 

φ(x)

3

ψ(x) = (φ(x)

3

ψ(x)) = 0.

 

 

 

 

 

Той же цели мы могли бы достичь, если бы в (3.3) вместо η = y −3x взяли η = x − y/3, воспользовавшись замечанием 1.

22

 

 

 

 

 

 

 

 

Глава I. Метод характеристик

Пример 3.2. Найти решение уравнения

 

 

 

 

x2

2u

2xy

2u

 

3y2

∂u

(3.9)

 

 

 

= 0,

∂x2

∂x∂y

∂y2

удовлетворяющее начальным условиям

 

 

 

 

 

 

 

 

 

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

y=1

= f0(x),

∂y

y=1

= f1(x),

(3.10)

 

 

 

 

 

 

 

 

 

 

 

 

где f0(x), f1(x) — заданные функции.

Решение. В области x ≠ 0, y ≠ 0 справедливо ∆ = 4x2y2 > 0, следовательно, в указанной области уравнение (3.9) является гиперболическим. Его характеристическое уравнение

x2dy2 + 2xy dx dy − 3y2dx2 = 0

распадается на два уравнения:

x dy − y dx = 0 и x dy + 3y dx = 0.

Функции y/x и x3y являются интегралами этих уравнений соответственно. Поэтому для приведения уравнения (3.9) к каноническому виду произведем замену переменных:

ξ =

y

 

и η = x3y.

 

x

 

 

 

 

 

 

 

 

 

 

 

 

В результате этой замены мы получим уравнение

 

 

2u

 

1 ∂u

(3.11)

 

 

 

 

 

 

 

 

= 0.

 

ξη

ξ

При ∂u/∂ξ = v оно принимает вид

 

 

 

 

 

 

 

 

 

∂v

 

v

 

 

 

 

 

 

= 0.

 

 

 

η

 

Это есть уравнение типа (2.3). В соответствии с (2.5) его общее решение

имеет вид w w

dvv 14 dηη = C(ξ),

§ 3.. Задача Коши для уравнения в частных производных второго порядка

23

откуда v = η1/4C1(ξ). Общее решение уравнения (3.11) получается из урав-

нения

∂uξ = η1/4C1(ξ)

интегрированием по ξ:

w

u = η1/4 C1(ξ)dξ + C2(η) = η1/4C3(ξ) + C2(η)

(функции C(ξ), C1(ξ), C2(η), C3(ξ) произвольные). Возвращаясь к пере-

менным x, y, получаем общее решение уравнения (3.9):

 

 

 

y

 

u(x, y) = x3/4y1/4C3

(

 

) + C2(x3y).

(3.12)

x

Теперь наша задача состоит в том, чтобы из бесконечного множества функций C2 и C3 выбрать такие, при которых (3.12) будет удовлетворять условиям (3.10).

Когда в соответствии с (3.10) мы положим в (3.12) y = 1, аргументы функций C2 и C3 окажутся равными 1/x и x3. А это, как видно из предыдущего примера, вызовет определенные осложнения. Поэтому произвольные функции C2 и C3 мы, имея в виду замечание 3, заменим другими произ-

вольными функциями, положив

и C2(x3y) = ψ (ω2(x3y)),

C3

(x)

= φ (ω1

(x))

 

 

y

 

 

y

 

где функции φ и ψ произвольны, а ω1 и ω2 выберем так, чтобы аргументы функций φ и ψ при y = 1 оба оказались равными x. Для этого следует

положить

ω1(z) = z1 и ω2(z) = z1/3.

Действительно, тогда

 

 

y

 

x

 

C3 (

 

)

= φ (

 

)

и C2(x3y) = ψ ((x3y)1/3) = ψ (xy1/3).

x

y

Теперь общее решение уравнения (3.9) получает более подходящий вид

 

x

(xy1/3).

 

u(x, y) = x3/4y1/4

φ (y ) + ψ

(3.13)

24

 

 

 

 

 

 

 

 

 

 

 

 

Глава I. Метод характеристик

Отсюда находим ∂u/∂y:

 

 

 

 

 

 

 

 

 

 

∂u

 

 

 

1

 

 

 

 

x

 

x

 

 

 

=

 

x3/4y3/4φ

(

 

) − x7/4y7/4φ(

 

) +

 

 

∂y

4

y

y

 

 

 

 

 

 

1

xy2/3ψ

(xy1/3).

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

Полагая здесь и в (3.13) y = 1, из (3.10) получаем

 

 

 

 

 

 

 

 

 

 

 

x3/4φ(x) + ψ(x) = f0(x),

 

 

 

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

 

x3/4

φ(x) − x7/4φ(x) +

 

xψ(x) = f1

(x).

(3.14)

 

 

 

4

3

Исключая из этой системы уравнений функцию ψ(x), мы приходим к следующему уравнению относительно функции φ(x):

 

 

φ

 

x

) =

1

x3/4f

 

x

 

 

3

x7/4f

1(

x .

 

 

 

4

 

 

 

 

 

 

(

 

 

0

( )

4

)

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

φ(x) =

1

x3/4f0

(x) +

3

w x7/4

(

1

f0(x) − f1(x))dx + C

(3.15)

 

 

 

 

4

4

4

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

 

1

 

 

3 x

 

(

1

(z) − f1(z))dz + C,

 

φ(x) =

 

x3/4f0

(x) +

 

xw0

z7/4

 

f0

(3.16)

4

4

4

x0 — произвольное фиксированное число. Как известно, интегралы в формулах (3.15), (3.16) могут отличаться друг от друга только на постоянное число, что, ввиду наличия в этих формулах произвольного постоянного C, ничего не меняет.

Из (3.14) и (3.16) находим

 

3

 

 

3

x

z7/4 (

1

(z))dz − Cx3/4.

 

 

 

x3/4 xw0

ψ(x) =

 

f0

(x)

 

 

f(z0) − f1

4

4

4

§ 3.. Задача Коши для уравнения в частных производных второго порядка

25

Подставляя найденные функции φ и ψ в (3.13), получаем решение поставленной задачи:

 

 

 

 

 

 

 

 

 

(

 

 

 

 

 

 

 

( )

 

 

 

 

 

 

u(x, y) = x3/4y1/4

1

 

x3/4y3/4f0

 

 

x

 

+

 

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(

)

 

 

 

 

 

 

x0

 

 

 

 

(

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

3 x/y

 

 

 

 

1

 

 

 

− f1(z) dz + C) +

3

 

 

 

+

 

 

 

 

 

w z7/4

 

 

 

 

f(z0)

 

f0

 

xy1/3

4

 

 

4

4

 

 

3

 

 

 

 

 

xy1/3

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x3/4y1/4

xw0

z7/4

(

 

 

 

 

f(z0) − f1

(z))dz − Cx3/4y1/4 =

4

4

 

1

 

 

x

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

yf0

(

 

) +

 

f0 (xy1/3) +

 

 

 

 

 

 

 

 

 

4

y

4

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

x/y

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

x3/4y1/4 xyw1/3

z7/4

(

 

 

f(z0) − f1

(z))dz.

 

 

 

 

4

4

 

 

 

 

Рассмотрим еще один пример.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пример 3.3. Найти решение уравнения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u

 

 

 

2u

 

 

 

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

2x

 

 

3y

 

+ 3

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x∂y

∂y2

∂y

 

 

 

 

удовлетворяющее начальным условиям

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

y=1 = x2 2,

 

∂y

y=1

= x3.

 

 

 

(3.17)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение. Составим характеристическое уравнение

2x dx dy − 3y dx2 = 0, (3.18)

(2x dy + 3y dx) dx = 0.

Приравниваем по очереди нулю сомножители. Равенство

2x dy + 3y dx = 0

26

Глава I. Метод характеристик

является дифференциальным уравнением с разделяющимися переменными. Разделяя их, получим

2dyy + 3dxx = 0.

Интегрируем каждое слагаемое

2 ln |y| + 3 ln |x| = ln C1,

отсюда

x3y2 = C1.

Далее, dx = 0 и x = C2 Первый вариант замены:

ξ = x3y2, η = x.

Функция ξ, как и в предыдущем примере, на прямой y = 1 принимает значение ξ|y=1 = x3, что мало согласуется с данными (3.17), но мы не меняем ее, а, как и в примере 3.1, запишем удобное общее решение. Это преобразование вполне будет оправдано простой заменой переменных. Вычисляем производные при ξ = x3y2, η = x:

 

∂u

 

 

∂u

∂u

∂u ∂u

 

 

 

 

 

=

 

 

3x2y2 +

 

 

,

 

 

 

=

 

 

2x3y,

 

×3

 

∂x

ξ

η

∂y

ξ

2u

 

2u

2u

 

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂x∂y

=

ξ2

6x5y3 +

ξη

2x3y +

ξ

6x2y,

 

×2x

 

 

 

 

 

 

 

 

2u

 

2u

 

 

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂y2

=

ξ2

4x6y2 +

ξ

2x3.

 

× − 3y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Подставим найденные производные в данное уравнение:

2u

(12x6y3 12x6y3) +

2u

4x4y +

ξ2

ξη

+

∂u

(12x3y − 6x3y + 6x3y) = 0.

 

ξ

После приведения подобных членов получаем

4x4y

2u

+ 12x3y

∂u

 

 

= 0,

ξη

 

 

 

ξ

§ 3.. Задача Коши для уравнения в частных производных второго порядка

27

сокращаем на 4x4y и меняем x на η:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u

 

 

3 ∂u

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξη

η ξ

 

Пусть v =

∂u

, тогда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ξ

 

 

 

∂v

 

 

 

 

3v

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

.

 

 

 

 

 

 

 

 

η

 

η

 

Воспользуемся равенством (2.5), поменяв в нем местами ξ и η:

 

 

 

w

dv

 

=

3 w

 

dη

+ ln C(ξ),

 

 

 

v

 

 

 

 

 

 

 

 

 

 

 

η

 

ln |v| = 3 ln |η| + ln C(ξ).

Теперь

 

C(ξ)

 

∂u

 

C(ξ)

 

 

 

 

 

 

 

v =

 

;

 

=

 

 

 

.

 

 

η3

ξ

η3

Далее

w C(ξ) dξ + φ(η) =

 

 

 

u, η) =

1

1

ψ(ξ) + φ(η).

η3

η3

Здесь C(ξ), φ(η) и ψ(ξ) — произвольные функции, w

ψ(ξ) = C(ξ) dξ.

Переходим к переменным x и y

u(x, y) = x13 ψ(x3y2) + φ(x).

Преобразуем ψ(x3y2) к произвольной функции от аргумента ξ1 = xy2/3:

ψ(x3y2) = ψ ((xy2/3)3) = g ((xy2/3)) = g1).

 

Итак,

(xy2/3) + φ(x).

 

1

(3.19)

u(x, y) = x3 g

Подберем функции g и φ, удовлетворяющие начальным условиям (3.17).

Найдем предварительно производную ∂u/∂y из равенства (3.19)

 

 

∂u

1

 

 

(xy2/3)x

2

 

 

 

=

 

gξ

 

 

y1/3.

(3.20)

 

∂y

x3

1

3

28

Глава I. Метод характеристик

Подставим y = 1 в равенства (3.19) и (3.20) и удовлетворим условиям (3.17):

1

g(x) + φ(x) = x2 2,

 

 

 

x3

 

1

g(x)x

2

= x3

, g(x) =

3

x5.

 

 

 

 

 

x3

3

 

 

 

 

 

 

2

Интегрируем последнее равенство:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x6

 

 

 

 

 

 

 

 

 

g(x) =

 

 

+ C.

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Далее,

 

 

 

 

 

 

x6

 

 

 

 

 

 

 

 

 

1

 

 

 

φ(x) = x2 2

 

(

 

+ C).

x3

4

Подставляем найденные выражения для g(x) и φ(x) в (3.19), заменив аргумент в g(x) на xy2/3:

1

 

(

xy2/3

)

6

 

C

+ x2 2

x3

C

u(x, y) =

 

 

 

 

+

 

 

 

.

x3

4

 

x3

4

x3

После небольших упрощений получаем окончательный результат:

u(x, y) = x3(y4 1) + x2 2.

4

§ 4. Задачи с данными на характеристиках

Помимо условий задачи Коши, для уравнения (1.1) на практике возникают другие требования. В частности, для уравнения гиперболического типа ставится задача нахождения частного решения по значениям на двух характеристиках разных семейств. Приведем характерный пример.

Пример 4.1. Найти решение уравнения

 

 

2u

 

2u

 

1 ∂u

(4.1)

 

 

+ y

 

+

 

 

 

= 0 (y < 0),

 

∂x2

∂y2

 

 

 

 

 

2 ∂y

 

если u(x, y) = φ1(x) на характеристике x − 2

 

= 0

(L1)

−y

и u(x, y) = φ2(x) на характеристике x + 2

 

= 1,

(L2)

−y

§ 4.. Задачи с данными на характеристиках

29

где φ1(x) и φ2(x) — функции, заданные соответственно на

отрезках

0 6 x 6 1/2 и 1/2 6 x 6 1, причем φ1 (1/2) = φ2 (1/2) (рис. 4.1).

 

Решение. Из условия видно, что решение ищется в полуплоскости y < 0 в области, ограниченной линиями L1, L2 и y = 0.

Y

0

0,5

1

X

L2 L1

Рис. 4.1

√ √

Легко видеть, что подстановкой ξ = x−2 −y, η = x+2 −y уравнение

(4.1) приводится к каноническому виду

 

 

 

 

 

2u

 

 

 

 

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

ξη

 

 

 

 

Поэтому (см. пример 2.1) его общее решение таково:

 

 

 

 

 

 

(4.2)

u(x, y) = θ1(x − 2 −y) + θ2(x + 2

 

−y),

где θ1 и θ2 — произвольные функции. Наша задача состоит в том, чтобы

выбрать такие θ1 и θ2, при которых функция (4.2) удовлетворяет условиям

задачи.

 

На характеристике L1, т. е. при 2

 

= x, мы по условию имеем

 

−y

(4.3)

θ1(0) + θ2(2x) = φ1(x).

Аналогично на L2, т. е. при 2

 

= 1 − x, справедливо

 

−y

(4.4)

θ1(2x − 1) + θ2(1) = φ2(x).

Из (4.4), полагая 2x − 1 = z, получаем

()

θ1(z) = φ2

z + 1

θ2(1),

(4.5)

2

30

 

 

 

 

 

 

 

 

 

 

 

 

Глава I. Метод характеристик

а из (4.3) при z = 2x получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

 

 

 

 

 

 

 

 

θ2(z) = φ1

(

 

) θ1(0).

 

 

 

 

2

Требуемые функции θ1 и θ2 найдены. Подставляем их в (4.2):

u(x, y) = φ2 (

 

2

 

) + φ2 (

+

2

 

) 1(0) + θ2(1)).

2

y + 1

 

x

 

 

 

 

 

 

x

 

2 y

 

Подставляя в (4.3) x = 1/2, получим θ1(0) + θ2(1) = φ1 (1/2), так что искомое решение задачи приобретает вид

u(x, y) = φ1

(

 

+ 2

 

 

)

+ φ2 (

 

 

2

) φ1

(2).

x

 

x

2

 

 

2

 

y

 

 

 

y + 1

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Задачи

Определить тип заданного уравнения в заданной области.

I 1.1. (y + 1)2u 2 2u + x2u ∂u = 0 ∂x2 ∂x∂y ∂y2 ∂y

в прямоугольнике 1 < x < 3, 0 < y < 1.

I 1.2.

y

2u

+ x

2u

+ 2(x + y)

2u

= 0

∂y2

∂x2

∂x∂y

в круге x2 + (y − 6)2 < 1.

I 1.3. 2xy 2u + x2 2u + y2 2u − x∂u + y∂u = 0 ∂x∂y ∂y2 ∂x2 ∂y ∂x

в квадрате |x| < 1, |y| < 1.

I 1.4.

(x + y)

2u

+ (x − y)

2u

+ xu = 0

∂x2

∂y2

в круге (x − 5)2 + y2 < 1.

I 1.5. (x + 1)2u 2x2 ∂u + (y − 3)2u + u = 0

∂x2 ∂x ∂y2

в квадрате 0 < x < 1, 0 < y < 1.

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