даишев лекция по дифурам
.pdfПравые части этих уравнений с точностью до первого порядка представим по формуле Маклорена:
x˙ = 2x + 8y + R1(x, y),
y˙ = −x − 3y + R2(x, y).
Следовательно, система первого приближения имеет вид:
x˙ = 2x + 8y,
y˙ = −x − 3y.
Характеристическое уравнение этой системы det (A − λE) = λ2 + λ+ 2 = 0, его корни λ1,2 = −1/2 ± i · 7/2. Легко видеть, что Reλ1 = Reλ2 = − 1/2 < 0. Решение x = 0, y = 0 является устойчивой точкой покоя.
ЛЕКЦИЯ 17
17.1. Уравнения в частных производных первого порядка
|
|
Определение. Выражение вида |
|
|
|
|
|
|
||||||||||||
Φ |
u(x |
, x |
, ..., x |
|
), |
|
∂u |
, |
∂u |
, ..., |
∂u |
, |
∂2u |
, |
∂2u |
, ..., |
∂mu |
|
= 0 |
|
n |
|
|
|
|
∂x12 |
∂x1∂x2 |
∂xnm |
|||||||||||||
|
|
1 |
2 |
|
|
∂x1 |
∂x2 |
∂xn |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(17.1)
называется дифференциальным уравнением в частных производных m-го порядка относительно неизвестной функции u(x) ≡ u(x1, x2, ..., xn). Как и в обыкновенных дифференциальных уравнениях, порядок старшей производной, входящей в уравнение, называется порядком этого уравнения. Рассмотрим некоторые примеры таких уравнений.
Пример 1. Рассмотрим дифференциальное уравнение
131
∂z(x, y) ∂x
Оно легко интегрируется, его решение имеет вид z = 12x2 + xy + ϕ(y), где ϕ(y) — произвольная функция y.
Пример 2.
∂2z = 1. ∂x∂y
Решение этого дифференциального уравнения также не вызы-
вает затруднений: |
∂ |
|
∂z |
= 1, откуда |
∂z |
= x + A(y) |
и, |
|
∂y |
||||||
|
∂x |
∂y |
|
|
|
||
следовательно, z = xy |
+ |
A(y) dy + B(x). |
Здесь A(y) |
и |
|||
B(x) — произвольные функции указанных в скобках аргументов.
Примеры показывают, что общее решение дифференциального уравнения в частных производных первого порядка содержит одну произвольную функцию, а уравнение второго порядка
— две произвольные функции, и так далее...
Предположения нуждаются в строгом обосновании. С. Ковалевской была доказана следующая теорема существования и единственности решения уравнения в частных производных.
Прежде чем сформулировать эту терему, напомним опре-
деление аналитической функции. |
|
|
|
|
|
Определение. Функция f (x1, x2, ..., xn) |
называется ана- |
||||
литической в окрестности точки |
◦ |
◦ |
◦ |
), |
если |
M0( x1, |
x2, |
... , xn |
|||
данная функция может быть разложена в степенной ряд, сходящийся в некоторой окрестности этой точки к самой этой функции.
132
17.2. Теорема Коши — Ковалевской существования и единственности решения уравнения в частных производных
Теорема Коши — Ковалевской. Существует единствен-
ное аналитическое в окрестности точки |
◦ |
◦ |
|
◦ |
) |
||||||||||||||
M0( x1, x2, ... , xn |
|||||||||||||||||||
решение уравнения |
m-го порядка, разрешённого относительно |
||||||||||||||||||
старшей производной по одной из переменных |
|
|
|
|
|||||||||||||||
∂mu |
= f x |
, ..., x |
|
, u, |
∂u |
, ..., |
∂u |
, |
∂2u |
, ..., ∂m−1 |
, ..., |
∂m |
, |
|
|||||
|
∂x1m |
n |
|
|
|
|
|
||||||||||||
|
1 |
|
|
|
∂x1 |
|
∂xn |
∂x12 |
∂x1m−1 |
|
∂xnm |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
удовлетворяющее условиям |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
u|x1=x◦1 = ϕ0(x2, ..., xn), |
|
|
|
|
|
|
|
|||||||||
|
|
|
∂u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|x1=x◦1 = ϕ1(x2, ..., xn), |
|
|
|
|
|
||||||||||
|
|
|
∂x1 |
|
|
|
|
|
|||||||||||
|
|
|
|
..................................., |
|
|
|
|
|
|
|||||||||
|
|
∂m−1u |
|x1=x◦1 = ϕm−1(x2, ..., xn), |
|
|
|
|
||||||||||||
|
|
|
∂x1m−1 |
|
|
|
|
||||||||||||
если функции ϕ0(x2, ..., xn), |
ϕ1(x2, ..., xn), |
... , ϕm−1(x2, ..., xn) |
|||||||||||||||||
являются аналитическими в окрестности точки M0 , а функция f в правой части уравнения является аналитической функцией в окрестности начальных значений своих аргументов.
Таким образом, решение определяется заданием начальных функций ϕ0(x2, ..., xn), ϕ1(x2, ..., xn), ... , ϕm−1(x2, ..., xn). Если мы их будем произвольно менять в классе аналитических функций, то получим совокупность аналитических решений исходного уравнения, зависящую от m произвольных функций. Иначе говоря, общее решение содержит m произвольных функций.
133
17.3. Линейные и квазилинейные уравнения в частных производных первого порядка
Определение. Квазилинейным неоднородным уравнением в частных производных первого порядка называется уравнение вида
X (x |
, x |
, ..., x |
|
, z) |
∂z |
+ X (x |
, x |
, ..., x |
|
, z) |
∂z |
+ ... |
||||||||||||
n |
|
|
n |
|
||||||||||||||||||||
1 |
1 |
|
2 |
|
|
|
|
∂x1 |
2 1 |
|
2 |
|
|
|
|
|
∂x2 |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
+X |
|
(x |
, x |
, ..., x |
|
, z) |
∂z |
= Y (x |
, x |
, ..., x |
|
, z). |
(17.2) |
|||||||||||
n |
n |
|
n |
|||||||||||||||||||||
|
|
1 |
|
2 |
|
|
|
|
|
∂xn |
|
1 |
|
2 |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Это уравнение, линейное относительно производных, может не быть, однако, линейным относительно неизвестной функции z(x1, x2, ..., xn). Если Y ≡ 0, а Xi = Xi(x1, x2, ..., xn), то
есть Xi не зависят от z, то уравнение |
n |
Xi(x1, x2, ..., xn) |
∂z |
= |
|
|
|
||||
|
∂x |
i |
|||
|
i=1 |
|
|
||
0 называется линейным однородным уравнением в частных производных первого порядка.
Для большей наглядности сначала целесообразно рассмот-
реть квазилинейное уравнение вида |
|
|
||||||
P (x, y, z) |
∂z |
+ Q(x, y, z) |
∂z |
= R(x, y, z), |
(17.3) |
|||
|
|
|||||||
|
|
∂x |
|
∂y |
|
|
||
где функции |
P (x, y, z), Q(x, y, z), R(x, y, z) задают непрерыв- |
|||||||
ное векторное поле |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F (x, y, z) = i P (x, y, z) + j Q(x, y, z) + k R(x, y, z). |
|
|||||||
Векторными линиями |
|
|
|
назы- |
||||
векторного поля F (x, y, z) |
||||||||
ваются линии, r = r(t), |
у которых касательный вектор в каж- |
|||||||
|
|
|
|
|
|
|
|
|
дой точке совпадает по направлению с вектором F (x, y, z) в |
||||||||
данной точке: |
dr |
|
|
k — const, t — параметр данной ли- |
||||
dt = kF , |
||||||||
нии. Векторные линии поля |
|
|
как известно, находятся |
|||||
F (x, y, z), |
||||||||
134
путём интегрирования системы обыкновенных дифференциальных уравнений
dx |
|
dy |
|
dz |
|
|
|
= |
|
= |
|
. |
(17.4) |
P (x, y, z) |
Q(x, y, z) |
R(x, y, z) |
||||
Поверхности, составленные из векторных линий, называются
F
N 
N F
Рис. 26. Векторная поверхность Рис. 27. Векторная трубка
векторными поверхностями (рис. 26) данного векторного поля. Из векторных линий построим векторную трубку (рис. 27). Боковая поверхность векторной трубки, как и любой векторной поверхности, характеризуется тем, что в любой ее точке выпол-
няется условие ортогональности |
N · F = 0, где N — вектор |
||||||||||||
нормали к боковой поверхности. |
|
|
|
|
|
|
|||||||
Если боковая поверхность векторной трубки определяется |
|||||||||||||
уравнением u(x, y, z) = 0, |
то |
|
|
j + |
∂z |
k, |
|
||||||
N = −−→ |
|
∂x i + ∂y |
|
||||||||||
|
|
∂u |
|
∂u |
|
∂u |
|
|
|||||
|
|
grad u = |
|
|
|
|
|
|
|
|
|
||
и условие ортогональности |
N · F = 0 даёт |
|
|
||||||||||
P (x, y, z) |
∂u |
+ Q(x, y, z) |
∂u |
|
+ R(x, y, z)∂u = 0 |
(17.5) |
|||||||
∂x |
∂y |
||||||||||||
|
|
|
|
|
|
|
∂z |
|
|||||
— линейное однородное уравнение в частных производных относительно функции трех переменных u(x, y, z).
135
Если боковая поверхность векторной трубки определяется
уравнением |
z = z(x, y), |
|
∂z |
|
∂z |
|
|
то N = −∂x |
i − |
∂y |
j + k и усло- |
||||
вие N · F |
= 0 после очевидных преобразований примет вид: |
||||||
P (x, y, z) ∂z |
+ Q(x, y, z) ∂z |
= R(x, y, z), что совпадает с уравне- |
|||||
∂x |
∂y |
|
|
|
|
|
|
нием (17.3).
Итак, чтобы найти уравнение боковой поверхности векторной трубки, необходимо проинтегрировать квазилинейное уравнение (17.3), когда мы ищем решение в явном виде z = z(x, y), или линейное уравнение (17.5), когда мы ищем решение в неявном виде u(x, y, z) = 0.
Для того чтобы решить уравнение (17.3) или (17.5), мы должны сначала найти векторные линии (так как векторные поверхности состоят из векторных линий), то есть решить вспомогательную систему обыкновенных дифференциальных уравнений (17.4):
dx |
= |
dy |
|
|
|
P (x, y, z) |
Q(x, y, z) |
dz
= R(x, y, z).
|
ψ1(x, y, z) = C1, — два независимых первых ин- |
|
|
Пусть (α) : |
|
|
ψ2(x, y, z) = C2, |
|
|
|
|
теграла системы (17.4), которые задают нам двухпараметрическое семейство векторных линий, называемых характеристиками уравнения (17.3) или (17.4). Но боковая поверхность векторной трубки состоит из однопараметрического семейства векторных линий. Для того чтобы из двухпараметрического семейства (α) с параметрами C1 и C2 выделить однопараметрическое семейство, мы должны задать зависимость (в общем случае произвольную) между этими параметрами: Φ(C1, C2) = 0.
136
Это означает, что искомое уравнение поверхности определяется соотношением
Φ (ψ1(x, y, z), ψ2(x, y, z)) = 0, |
(17.6) |
где Φ — произвольная функция своих аргументов, а (17.6) — общее решение уравнения (17.3).
Если требуется найти не произвольную векторную трубку
векторного поля |
|
а векторную поверхность, прохо- |
||||
F (x, y, z), |
||||||
|
|
|
|
ϕ1(x, y, z) = 0, |
|
|
|
|
|
|
|
|
|
дящую через заданную кривую |
|
|
(такая |
|||
(β) : |
ϕ2(x, y, z) = 0, |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
задача называется задачей Коши для дифференциального |
урав- |
|||||
нения в частных производных), то функция Φ в соотношении (17.6) уже не может быть произвольной: переменные x, y, z, входящие в это выражение, должны одновременно удовлетворять и условиям (α), и уравнению кривой (β). Иначе говоря, они одновременно должны удовлетворять системе уравнений
|
ϕ1(x, y, z) = 0, |
|
|
|
|
|
|
|
|
|
|
|
ϕ |
(x, y, z) = 0, |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ψ1(x, y, z) = C1, |
|
|
||
|
|
|
|
|
|
|
ψ2(x, y, z) = C2. |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Исключая x, y и z из этой системы, найдём конкретную
˜ |
Тогда |
связь между параметрами C1 и C2 : Φ( C1, C2) = 0. |
|
решением уравнения (15.3) будет конкретная функция |
|
˜ |
(17.7) |
Φ (ψ1(x, y, z), ψ2(x, y, z)) = 0, |
описывающая уравнение векторной поверхности, проходящей через заданную кривую (β).
137
|
Пример. Найти интегральную поверхность уравнения |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
∂z |
|
∂z |
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
x |
|
− y |
|
= 0, |
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
∂y |
∂x |
|
|
|
|
||||||
проходящую через кривую x = 0, |
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z = y |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
dx = |
dy |
= dz |
, первые интегралы этой |
|||||
Составим вспомогательную систему |
x |
|||||||||||||||||||
системы z = C1, x |
2 |
+ y |
2 |
= C2. |
|
|
−y |
|
|
|
0 |
|
||||||||
|
|
Следовательно, общим решением исходного |
||||||||||||||||||
уравнения является функция |
Φ *z, x2 + y2+ = 0, |
что эквивалентно z = |
||||||||||||||||||
f (x2 + y2) |
|
— поверхности вращения. |
x = 0, |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z = y |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Исключим x, y, z |
из системы |
|
|
|
|
|
|
получим C1 = C2 |
|||||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z = C1, |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
+ y |
2 |
= C2, |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
или |
z = x |
|
+ y |
|
— параболоид. |
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
ЛЕКЦИЯ 18
18.1. Общее решение линейного уравнения в частных производных первого порядка
Пусть задано линейное однородное уравнение первого порядка
X (x |
, ..., x |
|
) |
∂z |
+ X (x |
, ..., x |
|
) ∂z |
+ ... + X (x |
, ..., x |
|
|
) |
∂z |
|
= 0. |
||||||||||||||
|
|
|
|
|
∂xn |
|||||||||||||||||||||||||
1 1 |
|
|
|
n |
|
∂x1 |
2 1 |
|
|
n |
∂x2 |
|
n 1 |
|
|
n |
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(18.1) |
||||
Составим вспомогательную систему |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
dx1 |
|
|
|
|
|
= |
|
|
dx2 |
|
|
|
= ... = |
|
|
dxn |
|
|
|
|
. |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
X (x |
, x |
, ..., x |
n |
) |
X (x |
, x |
, ..., x |
n |
) |
X (x |
, x |
|
, ..., x |
n |
) |
|||||||||||||||
1 1 |
2 |
|
|
|
|
2 1 |
2 |
|
|
|
|
n 1 |
2 |
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(18.2) |
||||
138
Найдём n−1 функционально независимых первых интегралов:
|
ψ1(x1, x2, ..., xn) = |
C1, |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ψ |
(x |
, x |
, ..., x |
n |
) = |
C , |
|
|
2 |
|
1 |
2 |
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
............................... |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
ψn 1(x1, x2, ..., xn) = Cn 1. |
|||||||
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
− |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В n-мерном евклидовом пространстве с декартовыми координатами (x1, x2, ... , xn) эта система определяет (n − 1)-пара- метрическое семейство кривых, называемых характеристиками уравнения (18.1).
Теорема 1. Левая часть любого первого интеграла
ψ(x1, x2, ..., xn) = C
системы вспомогательных уравнений (18.2) является решением исходного линейного уравнения (18.1).
Доказательство. Вдоль любой интегральной кривой си-
стемы (18.2) выполняется тождество |
ψ(x1, x2, ..., xn) ≡ C. Это |
||||||||||||||||
значит, что dψ = |
n ∂ψ |
|
|
|
0. Из (18.2) следует, что |
|
|
||||||||||
|
|
|
|
dxi |
≡ |
dxi |
|||||||||||
|
|
|
|||||||||||||||
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|||
|
i=1 ∂x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
пропорциональны |
Xi(x1, ..., xn), и, следовательно, вдоль любой |
||||||||||||||||
интегральной кривой выполнено |
n |
∂ψ Xi(x1, x2, ..., xn) |
|
|
|||||||||||||
|
≡ |
0. |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
∂x |
i |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
i=1 |
|
|
|
|
|
|
|
Но интегральные кривые вспомогательной системы проходят че- |
|||||||||||||||||
рез каждую точку изменения переменных x1, |
x2, |
... , xn, |
а ле- |
||||||||||||||
вая часть тождества |
|
n |
∂ψ |
|
|
|
|
|
|
|
не зависит от |
||||||
|
|
|
|
Xi(x1, x2, ..., xn) |
≡ |
0 |
|||||||||||
∂x |
|
||||||||||||||||
|
|
i |
|
|
|
|
|
|
|
|
|
||||||
C1, C2, ... , Cn−1, |
|
i=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|||
и, следовательно, она не меняется при пере- |
|||||||||||||||||
ходе от одной интегральной кривой к другой. Это означает, что наше тождество справедливо не только вдоль интегральной кривой, но и во всей области изменения переменных x1, x2, ... , xn.
139
А это, в свою очередь, означает, что ψ( x1, x2, ... , xn ) есть решение уравнения в частных производных (18.1).
Очевидно, что
Φ (ψ1(x1, x2, ..., xn), ψ2(x1, x2, ..., xn), ..., ψn−1(x1, x2, ..., xn)) = C,
где Φ – произвольная дифференцируемая функция, является первым интегралом вспомогательной системы (18.2), так как вдоль интегральной кривой все функции ψi(x1, x2, ..., xn) ≡ Ci. Значит, и функция Φ обращается в постоянную вдоль интегральной кривой. Но тогда
z = Φ (ψ1(x1, x2, ..., xn), ψ2(x1, x2, ..., xn), ... , ψn−1(x1, x2, ..., xn))
является решением линейного уравнения (18.1).
Теорема. 2. Функция
z = Φ (ψ1(x1, x2, ..., xn), ψ2(x1, x2, ..., xn), ... , ψn−1(x1, x2, ..., xn)) ,
где Φ – произвольная дифференцируемая функция, является общим решением уравнения (18.1).
Доказательство. Мы должны доказать, что любое решение ϕ(x1, x2, ..., xn) уравнения (18.1) при конкретном выборе функции Φ может быть записано в виде ϕ = Φ (ψ1, ψ2, ... , ψn−1) . Пусть z = ϕ(x1, x2, ..., xn) – любое решение уравнения (18.1). Тогда выполнено
n |
∂ϕ X (x , x , ..., x ) = 0. |
(18.3) |
|||
|
|
i |
|
|
|
i=1 |
∂x |
|
i 1 2 |
n |
|
|
|
|
|
||
Докажем, что существует функция Φ такая, что
ϕ(x1, x2, ..., xn) = Φ (ψ1, ψ2, ..., ψn−1) ,
140
