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

УРАВНЕНИЯ С ЧАСТНЫМИ ПРОИЗВОДНЫМИ в примерах и задачах-объединены

.pdf
Скачиваний:
0
Добавлен:
06.02.2025
Размер:
427.79 Кб
Скачать
x 6=,0 y 6=.0

уравнения (1.9).

ПРИМЕР 2. Решить уравнение

x

2 2u

− y

2 2u

+ 2x

∂u

− 2y

∂u

= 0,

(1.13)

 

∂x2

 

∂y2

∂x

∂y

в области

В указанной области справедливо = x2y2 > 0. Следовательно, в этой области уравнение (1.13) гиперболическое. Для приведения его к каноническому виду составим характеристическое уравнение x2dy2 − y2dx2 = 0, которое распадается на пару различных уравнений xdy −ydx = 0 è xdy + ydx = 0. Общие

интегралы этих уравнений соответственно таковы: xy = C è

y

= C. Ôóíê-

 

 

 

 

x

 

öèè æå ϕ1(x, y) = xy è ϕ2(x, y) = xy

являются их интегралами. Поэтому в

результате замены переменных

 

 

 

 

 

ξ = xy,

η =

y

 

 

(1.14)

x

 

 

 

 

 

уравнение (1.9) будет приведено к каноническому виду. Вычисляем частные производные по формулам (1.6)

 

 

∂u

 

 

 

 

∂u ∂ξ

 

 

 

 

 

∂u ∂η

 

 

 

 

∂u

 

 

 

∂u y

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

y −

 

 

 

 

 

 

 

 

 

 

 

 

 

 

| × 2x,

 

 

 

∂x

 

 

∂ξ ∂x

∂η ∂x

∂ξ

∂η

x2

 

 

∂u

 

 

∂u ∂ξ

 

 

 

 

∂u ∂η

 

 

∂u

∂u 1

| × −2y,

 

 

 

 

 

 

=

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

=

 

 

 

x +

 

 

 

 

 

 

 

 

 

 

 

 

∂y

 

∂ξ ∂y

 

 

∂η ∂y

∂ξ

∂η x

 

2u ∂2u

2

 

 

 

 

 

y2

 

 

2u

 

 

 

 

 

2u y2

 

∂u 2y

 

2

 

 

=

 

y

 

 

− 2

 

 

 

 

+

 

 

 

 

 

+

 

 

 

 

 

| × x

,

∂x2

∂ξ2

 

 

x2

∂ξ∂η

∂η2

x4

 

 

∂η

x3

 

 

 

 

2u ∂2u

2

 

 

 

 

 

 

 

2u

 

 

 

 

 

2u 1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

=

 

x

 

 

 

+ 2

 

+

 

 

 

 

 

 

 

 

 

 

| × −y

.

 

 

 

 

 

∂y2

 

 

∂ξ2

 

 

 

∂ξ∂η

∂η2

x2

 

 

 

 

 

 

 

Смешанную производную, естественно, не вычисляем. Подставляем найден-

11

ные производные в уравнение (1.13) и приведем подобные члены

2u 2

 

2

2

2

 

 

2u

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

(x

y

 

− x

y

) +

 

 

 

(−2y

 

− 2y

)+

 

 

 

 

 

 

 

 

∂ξ2

 

∂ξ∂η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u y2

 

y2

 

∂u

 

∂u 2y 2y 2y

 

 

 

 

 

 

+

 

(

 

 

) +

 

(2xy − 2xy) +

 

(

 

 

 

) = 0

 

 

 

 

 

 

∂η2

x2

x2

∂ξ

∂η

x

x

x

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

 

2

2u

2y ∂u

 

−4y

 

 

 

 

 

= 0

 

∂ξ∂η

x

∂η

èëè

 

2u

1 ∂u

 

 

+

 

 

 

= 0.

 

 

 

∂ξ∂η

2xy ∂η

С учетом, что ξ = xy, окончательно имеем:

 

2u

1 ∂u

 

 

 

+

 

 

 

= 0.

(1.15)

 

 

 

∂ξ∂η

2ξ ∂η

 

Обозначим ∂u∂η = v(ξ, η). Тогда (1.15) принимает вид ∂v∂ξ + 2vξ = 0.

Это есть уравнение типа (1.9). На основании (1.12) его общее решение опре-

деляется равенством

R dvv +

21 R ξ = C(η). Отсюда

 

 

 

1

 

C(η)

 

ω(η)

 

 

 

v =

 

e

 

=

 

.

(1.16)

ξξ

(вследствие произвольности C(η) функция ω(η) также является произвольной). Решение уравнения (1.15) получается отсюда интегрированием по η:

R

ω(η)

 

 

 

 

ψ(ξ)

u

u =

ξ dη + ψ(ξ). Вместо произвольной постоянной мы написали здесь

произвольную функцию

 

, так как указанная функция удовлетворяет

уравнению (1.16) при произвольной ψ(ξ), в чем легко убедиться. Обозначив

ω(η)dη = ϕ(η), где функция ϕ(η) - снова произвольна, мы получаем общее

решение уравнения (1.15) в виде

1

 

. Возвращаясь к перемен-

 

 

 

u = ξ ϕ(η) + ψ(ξ)

R

 

12

∂u
∂ξ

ным x, y по формулам (1.14), выводим общее решение исходного уравнения (1.13): u = 1xy ϕ(xy ) + ψ(xy). Напоминаем, что функции ϕ è ψ здесь произвольны.

ПРИМЕР 3. Найти решение уравнения

2u

+ 2

2u

− 3

2u

= 0,

(1.17)

∂x2

∂x∂y

∂y2

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

u|y=0 = 3x2,

∂u

|y=0

= 0.

(1.18)

 

∂y

На всей плоскости ху справедливо

= 1 + 3 > 0. Следовательно, на всей

плоскости уравнение (1.17) является гиперболическим, а его характеристиче- ское уравнение dy2 −2dxdy−3dx2 = 0 распадается на два уравнения dy = −dx

è dy = 3dx. x + y = C è y −3x = C их общие интегралы, а функции x + y è y − 3x - интегралы этих уравнений. Поэтому к каноническому виду приводит следующая замена переменных:

ξ = x + y, η = y − 3x.

(1.19)

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

2u

 

 

= 0.

(1.20)

 

∂ξ∂η

 

Обозначив ∂u∂ξ = v(ξ, η), мы приводим его к виду ∂η∂v = 0, откуда заключаем,

что функция v îò η не зависит, то есть имеет вид: v = ω(ξ), ãäå ω(ξ) произвольная функция. Таким образом, для функции u мы имеем уравнение

= ω(ξ). Интегрируя по ξ, находим общее решение уравнения (1.20)

13

Z

u = ω(ξ)dξ + ψ(η) = ϕ(ξ) + ψ(ξ)

. Здесь обе функции ϕ(ξ), ψ(ξ) произвольные. Возвращаясь к переменным x, y, получаем общее решение уравнения (1.17):

u(x, y) = ϕ(x + y) + ψ(y − 3x).

(1.21)

Наша задача состоит в том, чтобы из бесконечного множества решений (1.21), выделить то, которое удовлетворяет условиям (1.18). Иными словами, мы должны найти вполне определенные функции ϕ è ψ, при которых функция u удовлетворяет условиям (1.18). Из (1.21) имеем ∂u∂y = ϕ0 (x + y) + ψ0 (y − 3x).

Полагая здесь и в (1.21) y = 0 , из (1.18) получаем:

ϕ(x) + ψ(−3x) = 3x2,

(1.22)

ϕ0 (x) + ψ0 (−3x) = 0.

(1.23)

Из последнего равенства желательно получить уравнение, в которое вместо производных ϕ0 (x), ψ0 (−3x) входили бы функции ϕ(x), ψ(−3x). Казалось

бы, этой цели можно достичь, изменив в (1.23) сумму производных на производную суммы функций. Но этому мешает то обстоятельство, что штрихи в (1.23) означают производные указанных функций по их аргументам ( ϕ0 (x) есть производная по х, ψ0 (−3x) - производная по -3х). Поэтому вначале в (1.23) перейдем к производным по одному и тому же переменному (например, по х).

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

,

 

 

 

Имеем: dψ( 3x)

 

 

 

dψ( 3x)

·

 

 

d( 3x)

 

(−3x) · (−3)

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

= ψ

 

 

 

 

 

 

 

 

dx

d(−3x)

dx

 

 

 

 

 

откуда

0

(−3x) = −

1

dψ(−3x)

 

 

 

 

 

 

 

 

 

 

 

 

 

ψ

3

 

dx . Теперь (1.23) принимает вид

 

 

 

 

dϕ(x)

 

1

 

dψ(−3x)

=

d

(ϕ(x)

1

ψ(

3x)) = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx

 

 

3

 

 

3

 

 

 

 

 

 

 

 

 

dx

 

 

dx

 

 

 

14

Отсюда

 

 

ϕ(x) −

1

ψ(−3x) = C.

 

(1.24)

 

 

 

 

 

3

 

Из (1.22) и (1.24) находим

 

 

 

 

 

 

 

 

 

 

 

3

 

2

 

3

 

 

9x2

3

ϕ(x) =

 

x

 

+

 

C, ψ(−3x) =

 

 

 

C.

4

 

4

4

4

Наконец, положив −3x = z, получаем ψ(z) = z42 34C. Требуемые функции ϕ è ψ найдены. Подставляя их в (1.21), мы получим искомое решение задачи

 

3

 

2

 

3

(y−3x)2

3

2

 

2.

u(x, y) =

4

(x + y)

 

+

4C +

4

4C = 3x

 

+ y

 

ЗАМЕЧАНИЕ 3. Осложнения, связанные с неодинаковостью аргументов функций ϕ è ψ в (1.22) и (1.23), устранятся, и решение поставленной зада- чи упростится, если произвольную функцию ψ в(1.21) мы представим в виде сложной функции ψe(ω(y − 3x)), ãäå ψe - произвольная, а ω - определенная функция, причем последняя выбирается так, чтобы ψe(ω(y − 3x)) ïðè y = 0 принимала вид ψe(x). В нашем случае, очевидно, следует положить ω(z) = 3z . Тогда общее решение уравнения (1.17) будет иметь вид

u(x, y) = ϕ(x + y) + ψe(−y −33x) = ϕ(x + y) + ψe(x − y3),

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

ϕ(x) + ψe(x) = 3x2

è ϕ0 (x) − 13ψe0 (x) = (ϕ(x) − 13ψe(x))0 = 0.Той же цели мы могли бы достичь, если бы в(1.19) вместо η = y − 3x, мы взяли η = x − y3 , воспользовавшись ЗАМЕЧАНИЕМ 1.

15

ПРИМЕР 4. Найти решение уравнения

x

2 2u

− 2xy

2u

− 3y

2 ∂u

= 0,

(1.25)

 

∂x2

∂x∂y

 

∂y2

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

u|y=1 = f0(x),

∂u

|y=1

= f1(x),

(1.26)

 

 

∂y

ãäå f0(x) , f1(x) - заданные функции.

 

 

 

 

 

В области x 6= ,0 y 6= 0справедливо

 

= 4x2y2

> 0, следовательно, в

указанной области уравнение (1.25) является гиперболическим. Его харак-

теристическое уравнение x2dy2 + 2xydxdy − 3y2dx2 = 0 распадается на два уравнения

xdy − ydx = 0 è xdy + 3ydx = 0.

Функции y

x3y являются интегралами этих уравнений соответственно. По-

x è

 

этому для приведения уравнения (1.25) к каноническому виду произведем за-

мену переменных ξ = y

è

η = x3y. В результате этой замены мы получим

x

 

 

 

 

 

 

 

 

уравнение

 

 

2u

 

 

 

 

 

 

 

 

 

1 ∂u

= 0.

(1.27)

 

 

 

 

 

 

 

 

 

 

∂ξ∂η

4ξ ∂ξ

Ïðè ∂u∂ξ = v оно принимает вид ∂η∂v 4vη = 0. Это есть уравнение типа (1.9). В

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dvv

41 dηη = C(ξ), îò-

соответствии с (1.12) его общее решение имеет вид

R

 

R

 

 

v = η

1

C1(ξ). Общее решение уравнения (1.27)

 

 

êóäà

 

 

4

 

 

 

 

 

 

 

 

 

 

получается из уравнения

 

1

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

1

∂u∂ξ = η4 C1

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

u = η4

C1(ξ)dξ+C2(η) = η4 C3(ξ)+C2(η)

 

 

C(ξ)

 

C1(ξ)

 

C2(η)

 

C3

(ξ) -

 

R

 

 

 

 

(функции

 

 

 

,

 

,

 

,

 

 

произвольные). Возвращаясь к пере-

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

3 1

 

y

(x3y). (1.28)

u(x, y) = x4 y4 C3

(

 

) + C2

 

 

 

x

 

16

Теперь наша задача состоит в том, чтобы из бесконечного множества функций

C2 è C3 выбрать такие, при которых (1.28) будет удовлетворять условиям (1.26).

Когда в соответствии с (1.26) мы положим в (1.28) y = 1, аргументы функ-

öèé C

è C

окажутся равными

1

 

x3. А это, как видно из предыдущего

2

3

 

x

è

 

 

 

 

 

примера, вызовет определенные осложнения. Поэтому произвольные функ- öèè C2 è C3 мы, имея в виду ЗАМЕЧАНИЕ 3, заменим другими произволь-

ными функциями, положив C3(xy ) = ϕ(ω1(xy )) è C2(x3y) = ψ(ω2(x3y)), ãäå

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

ω1(z) = z1

1

. Действительно, тогда мы будем иметь C3(xy ) = ϕ(xy ) è

è ω2(z) = z 3

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C2(x3y) = ψ(x3y3 ) = ψ(xy3 ). Теперь общее решение уравнения (1.25) получает

более подходящий вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

1

 

x

1

 

(1.281)

 

 

 

 

 

 

u(x, y) = x4 y4

ϕ(

 

) + ψ(xy3 ).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂y :

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂u 1

3

3

x

7

7 0 x

1

2

1

 

 

 

=

 

x4 y

4

ϕ(

 

) − x4 y4 ϕ (

 

) +

 

xy

3

ψ(xy3 ).

 

 

∂y

4

y

y

3

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

 

 

3

ϕ(x) + ψ(x) = f0(x),

 

 

 

 

x4

 

 

1

3

 

7

0

1

0

 

 

 

x4

ϕ(x) − x4

ϕ (x) +

 

 

xψ (x) = f1

(x).

(1.29)

4

3

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

ϕ0 (x) = 14x34 f00 (x) − 34x74 f1(x).

17

Отсюда

ϕ(x) = 4x

4 f0(x) +

4

Z

x4

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

(1.30)

1

 

3

3

 

7

1

 

(слагаемое с f00 (x) мы проинтегрировали по частям). Мы увидим далее, что

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

ϕ(x) = 4x

4 f0(x) +

4

x

z4

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

(1.31)

Zx0

1

 

3

3

 

7

1

 

x0 - произвольное фиксированное число. Как известно, интегралы в формулах (1.30), (1.31) могут отличаться друг от друга только на постоянное число, что ввиду наличия в этих формулах произвольного постоянного С ничего не меняет.

Из (1.29) и (1.31) находим

ψ(x) = 4f0(x) −

4x4

x

z4

(4f(z0) − f1(z))dz − Cx4 .

Zx0

3

 

3

3

 

7

1

3

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

u(x, y) = x4 y4

(4x4 y4 f0

 

 

 

 

 

 

x

z

4

(4f(z0) − f1(z))dz + C)+

(y ) + 4 Zx0

3

1

1

 

 

3

3

x

3

 

 

y

 

7

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

+4f0(xy3 ) −

 

 

 

 

xy 3

z4

(4f(z0) − f1(z))dz − Cx4 y4 =

4x4 y4 Zx0

3

 

1

 

 

3

3

1

 

7

 

1

 

 

 

 

 

3

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

4

y

 

4

 

 

 

4

 

 

 

 

 

 

 

4

 

 

 

 

1

 

 

 

1

Zxy 3

 

7

 

 

1

x

 

3

 

 

3 3

 

 

y

 

1

 

 

 

=

 

yf0(

 

) +

 

 

f0(xy3 ) +

x4 y4

1

z4 (

f(z0)

 

f1(z))dz.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В случае гиперболического уравнения его частное решение может быть найдено также по заданию искомой функции на паре независимых характеристик этого уравнения (характеристикой уравнения (1.1) называется

18

кривая ϕ(x, y)=C, ãäå ϕ(x, y) есть какой-нибудь интеграл характеристическо-

го уравнения (1.4), а С - любая постоянная. Характеристики ϕ1(x, y) = C1,

ϕ2(x, y) = C2 называются независимыми, если независимы интегралы ϕ1(x, y),

ϕ2(x, y), т.е. если для них выполняется неравенство (1.3)).

Рассмотрим еще один пример. ПРИМЕР 5. Найти решение уравнения

2x 2u − 3y2u + 3∂u = 0, ∂x∂y ∂y2 ∂y

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

u|y=1 = x2 − 2,

∂u

|y=1

= x3.

(1.32)

∂y

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

−2xdxdy − 3ydx2 = 0,

(1.33)

(2xdy + 3ydx)dx = 0,

Приравниваем по очереди нулю сомножители. Равенство 2xdy + 3ydx = 0

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

Разделяя их, получим

2dyy + 3dxx = 0.

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

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

отсюда

x3y2 = C1.

19

Далее dx = 0 è x = C2 Первый вариант замены ξ = x3y2, η = x.

Функция ξ ,как и в предыдущем примере, на прямой y = 1 принимает зна-

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

 

 

 

 

∂u

 

 

∂u

 

 

 

 

 

 

 

 

 

 

∂u

 

 

∂u

 

 

∂u

 

 

 

 

 

 

 

=

 

 

 

3x2y2 +

 

 

 

;

 

 

 

 

 

 

=

 

 

 

 

 

2x3y;

| × 3

 

 

 

 

∂x

∂ξ

∂η

∂y

∂ξ

 

 

 

2u

 

 

2u

 

 

 

 

 

 

 

 

 

 

 

2u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂u

 

 

 

 

 

 

=

 

 

 

6x5y3

+

 

 

 

 

 

2x3y +

 

 

 

6x2y;

 

| × 2x

 

 

∂x∂y

 

∂ξ2

∂ξ∂η

∂ξ

 

 

 

 

 

 

 

2u ∂2u

4x6y2 +

 

 

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

2x3.

 

 

 

 

| × −3y

 

 

 

 

 

 

 

 

∂y2

∂ξ2

 

 

∂ξ

 

 

 

 

 

 

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

 

2u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u

 

 

 

 

 

 

 

 

 

 

∂u

 

 

 

 

 

 

 

 

 

(12x6y3 − 12x6y3) +

 

 

 

4x4y +

 

 

 

 

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

 

∂ξ2

∂ξ∂η

 

∂ξ

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4

 

 

 

 

2u

3

 

∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

4x

y

 

 

 

 

 

 

+ 12x

y

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ξ∂η

 

 

 

 

 

 

 

 

 

 

 

∂ξ

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2u

3 ∂u

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+

 

 

 

 

 

 

 

 

 

= 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η ∂ξ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

∂ξ∂η

 

 

 

 

 

 

 

 

 

 

 

 

Пусть v = ∂u∂ξ тогда ∂η∂v = −

3v

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

η

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

ξ è η,

 

 

 

 

 

 

 

 

 

Z

 

 

 

v = −3 Z

 

 

 

 

 

η

 

 

+ lnC(ξ);

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dv

 

 

 

 

 

 

 

 

 

 

dη(ξ)

 

 

 

 

 

 

 

ln|v| = −3ln|η| + lnC(η).

20