Толстых-Уравнматем-физики
.pdf2. Доказать, что общее решение уравнения теплопроводности
u a2 2ut x2
может быть представлено в виде
u f (a2t x) g(a2t x),
где f и g – произвольные функции. Привести примеры частных решений.
3. Проверить, что функции ln(x2 y2), x2 y2 , arctg y x
являются частными решениями уравнения Лапласа
2u 2u 0.x2 y2
4. Доказать, что выражение
u f (at x) g(at x), t
где f и g – произвольные функции, задает общее решение уравнения
2u |
a2 |
2u |
|
2 u |
|||
|
|
|
|
|
|
0. |
|
t2 |
x2 |
|
|
||||
|
|
t t |
Привести пример частного решения.
5. Найти общее решение u = u(x, y) уравнений:
|
u |
u |
|
|
2u |
|
2u |
|
|
а) |
|
cosx, б) |
|
x, в) |
|
|
0, г) |
|
1. |
|
y |
|
x2 |
x y |
|||||
|
x |
|
|
|
|
||||
Ответы: а) u sin x g y ; |
б) u xy f x ; |
|
|||||||
в) u xf y g y ; г) u xy f x g y . |
|
|
|||||||
Здесь f и g означают произвольные функции. |
|
11
2.Линейные дифференциальные уравнения
вчастных производных и их основные свойства
Линейное дифференциальное уравнение в частных производных содержит неизвестную функцию и ее частные производные по всем переменным только в первой степени, иначе говоря, частные производные входят в уравнение линейно (см. напр. (1.0)). Нетрудно видеть, что выше в п. 1 были рассмотрены только линейные уравнения.
Легко видеть также, что уравнения
3u |
u2, |
2u |
|
u |
|
u |
0, |
u |
u |
5 |
x3 |
x2 |
|
|
|
||||||
|
|
x y |
|
x |
являются нелинейными. Первое уравнение содержит искомую функцию во второй степени; второе – произведение частных производных; третье – произведение искомой функции на ее частную производную.
Общий вид линейного дифференциального уравнения в частных производных второго порядка относительно функции u = u(t, x) задает формула
|
2u |
|
2u |
|
2u |
|
u |
u |
|
|||
A |
|
|
B |
|
C |
|
D |
|
E |
|
Fu G, |
(2.1) |
t |
2 |
t x |
x2 |
|
|
|||||||
|
|
|
|
t |
x |
|
где A, B, C, D, E, F, G – функции аргументов t, x.
Если функция G тождественно равна нулю G 0 , то получаем уравнение
|
2u |
|
2u |
|
2u |
|
u |
u |
|
|||
A |
|
|
B |
|
C |
|
D |
|
E |
|
Fu 0, |
(2.2) |
t |
2 |
t x |
x2 |
|
|
|||||||
|
|
|
|
t |
x |
|
которое назовем однородным, а уравнение (2.1) – неоднородным.
Свойства однородного линейного дифференциального уравнения в частных производных
1 . Если u1 и u2 – частные решения уравнения (2.2), то их линейная комбинация c1u1+c2u2 (c1, c2 – произвольные постоянные) также решение этого уравнения.
Для доказательства достаточно применить свойства производных: производная суммы функций равна сумме производных, постоянный множитель можно вынести за знак производной (любого порядка), а затем определение решения уравнения.
12
Действительно,
u |
c |
u1 |
c |
|
u2 |
, |
u |
c |
u1 |
c |
|
u2 |
, |
|
|
|
|
|
|||||||||
t |
|
|
|
x |
|
|
|||||||
1 t |
2 t |
1 x |
2 x |
2u |
|
2u |
|
|
|
2u |
2 |
|
|
|
2u |
|
|
2u |
|
|
|
2u |
2 |
|
||||||||||||
|
|
|
c |
|
|
1 |
|
c |
2 |
|
|
|
|
|
, |
|
|
c |
|
|
1 |
c |
2 |
|
|
|
, |
|||||
|
2 |
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
2 |
|
|
2 |
||||||||||||||
t |
|
1 |
|
t |
|
|
|
t |
2 |
|
|
|
|
x |
1 |
x |
|
x |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
2u |
|
|
2u |
|
|
|
|
|
2u |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
c |
|
|
1 |
c |
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
t x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
1 t x |
|
|
2 t x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя эти значения частных производных в (2.2) и группируя слагаемые относительно u1 и u2,получаем:
с A |
2u1 |
B |
2u1 |
C |
2u1 |
D |
u1 |
E |
u1 |
Fu |
||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||
1 |
t2 |
|
|
t x |
|
|
x2 |
|
|
t |
|
|
x |
1 |
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
с2 A |
2u2 |
B |
2u2 |
C |
2u2 |
|
D |
u2 |
E |
u2 |
Fu2 0, |
|||||||||
|
|
x2 |
|
|
||||||||||||||||
|
|
|
t2 |
|
t x |
|
|
|
t |
|
x |
так как выражения в скобках равны нулю в силу того, что функции u1 и u2 являются решениями уравнения (2.2), то есть обращают его в тождество.
2 . Пользуясь методом математической индукции, можно показать, что любая конечная линейная комбинация (2.3) решений уравнения (2.2)
u c1u1 c2u2 ... cnun |
(2.3) |
также является решением уравнения (2.2). |
|
При n выражение (2.3) превращается в ряд: |
|
|
|
u cnun . |
(2.4) |
n 1 |
|
Вопрос о том, является ли ряд (2.4) решением данного уравнения (2.2), требует специального исследования. При этом часто используется признак Вейерштрасса (см., например, [3, 4, 8, 17]).
Признак Вейерштрасса. Если ряд (2.4) мажорируется сходящимся числовым рядом
13
un an, an ,
n 1
то ряд (2.4) сходится равномерно и абсолютно (иначе говоря, имеет место правильная сходимость ряда).
В этом случае можно показать, что (2.4) является решением уравнения (2.2), если ряды, полученные почленным дифференцированием, также обладают правильной сходимостью.
Например, функциональный ряд
u |
cosatsinx |
|
cos2atsin2x |
... |
cosnatsinnx |
... |
(2.5) |
|||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
2 |
|
|
|
|
|
|
|
|
|
|
|
|
22 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2n |
|
|
|
||||||
сходится равномерно и абсолютно, так как его члены ограничены |
|
|||||||||||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
(sin |
|
1, |
|
|
cos |
|
|
|
|
1): |
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
|
|
un |
|
|
|
cosnatsinnx |
|
|
|
|
|
|
|
1 |
; |
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2n |
|
|
|
|
|
|
|
|
|
|
2n |
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||
мажоранта – числовой |
|
|
|
|
|
|
|
ряд |
|
|
|
|
|
|
|
|
|
|
|
... |
|
|
|
|
... |
, |
представляющий |
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
|
4 |
|
2n |
||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
бесконечно убывающую геометрическую прогрессию: |
|
|
|
|||||||||||||||||||||||||||||||||||||||||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
b |
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|
|||||||||||
|
S |
|
|
|
S |
|
1 |
|
|
|
|
,b1 |
|
|
|
|
|
,q |
|
1; |
|
|
||||||||||||||||||||||||||||
|
|
n |
|
|
|
|
|
2 |
|
|
|
|||||||||||||||||||||||||||||||||||||||
|
n 1 2 |
|
|
|
|
|
|
|
|
|
|
|
|
1 q |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
cosnatsinnx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
u |
|
1. |
|
|
|
|
|
||||||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Заметим (учитывая упражнение 1.1), что функция u, определяемая рядом (2.5), является решением волнового уравнения (1.1).
По аналогии, с линейными неоднородными дифференциальными уравнениями (ЛНДУ) с постоянными коэффициентами, можно доказать следующий факт.
Теорема. Если функция uо – решение однородного уравнения (2.2), а функция uч частное решение неоднородного уравнения (2.1), то функция u uo uч решение уравнения (2.1).
Читателю предлагается доказать это утверждение самостоятельно.
14
Задачи для самостоятельного решения
1. Какие из ниже приведенных уравнений являются нелинейными?
а)
в)
2u |
|
sintx |
2u |
uln x, |
||||||
t2 |
|
x2 |
||||||||
|
|
|
|
|
|
|||||
|
u |
2 |
u u |
|
||||||
|
|
|
|
x |
|
|
|
|
|
0, г) |
t |
|
t |
|
x |
||||||
|
|
|
|
|
|
б) |
u |
tgu |
u |
tx, |
|
|
|||||
|
|
|
|
||||||||
|
t |
|
x |
|
|
|
|
|
|||
cost |
2u |
sinx |
|
2u |
|
2u |
0. |
||||
t2 |
t x |
x2 |
|||||||||
|
|
|
|
|
|
2. Является ли функция u, определяемая рядом
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
u |
|
|
|
|
|
2 |
u |
|
|
u |
cosant cosnx |
, решением уравнения |
|
|
a2 |
|
? |
|||||||||||||||||
4 |
|
t |
2 |
|
|
2 |
||||||||||||||||||
n 1 |
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
x |
||||||
3. Является ли функция u, определяемая рядом |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
e a2n2t sinnx |
|
e a2nt sinnx |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
a) u |
|
|
|
, б) u |
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
4 |
|
|
n |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
n 1 |
|
|
n 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
решением уравнения теплопроводности |
|
u |
a2 |
|
2u |
|
? |
|
|
|||||||||||||||
|
|
|
x2 |
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
|
|
||||
4. Является ли функция u ln(x2 |
y2) 3arctg |
y |
(x2 y2) |
|||||||||||||||||||||
|
||||||||||||||||||||||||
|
|
|
|
|
|
2u |
|
|
2u |
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
решением уравнения Лапласа |
|
|
0? |
|
|
|
|
|
|
|
|
|
|
|||||||||||
x2 |
|
y2 |
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ответы и советы:
1.Нелинейные уравнения а, г, в; линейное – б.
2.Да. Проверьте дифференцированием: выражение, стоящее под знаком суммы, является решением волнового уравнения; убедитесь в правильной сходимости данного ряда и рядов, полученных двукратным почленным дифференцированием.
3.а) Да; б) Нет.
4.Да.
15
3. Метод Эйлера решения линейных дифференциальных уравнений в частных производных
Математические модели различных прикладных задач часто сводятся к изучению специального вида уравнений в частных производных, для которых можно не только изучить свойства решений, но и найти их. К простейшему специальному виду относятся линейные однородные дифференциальные уравнения в частных производных с постоянными коэффициентами
A |
2u |
B |
2u |
C |
2u |
0, |
(3.1) |
|
t2 |
t x |
x2 |
||||||
|
|
|
|
|
где А, В, С – константы.
Методом Эйлера будем искать решения уравнения (3.1) в виде
u f (x t), |
(3.2) |
где f – произвольная функция (ср. метод Эйлера – метод характеристического многочлена для обыкновенных линейных однородных дифференциальных уравнений с постоянными коэффициентами).
Нужно подобрать неизвестную постоянную так, чтобы функция (3.2) была решением уравнения (3.1). Используя правило дифференцирования сложной функции, вычислим частные производные в
(3.1):
u |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
f |
|
|
|
|
|
|
|
||||||||||
|
|
(x t)(x t)t f |
(x t), |
|
||||||||||||||
t |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
u |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
f |
|
|
|
|
|
|
|
||||||||||
|
|
(x t)(x t)x f |
(x t), |
|
||||||||||||||
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
2u |
|
|
u |
|
|
|
2 |
|
||||||||||
|
|
|
|
|
|
|
|
|
f |
|
||||||||
t2 |
|
|
|
(x t)(x t)t |
f (x t), |
|||||||||||||
|
t |
t |
|
|
|
|
|
|||||||||||
2u |
|
|
|
u |
f |
(x t)(x t)x f (x t), |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
x2 |
|
|
|
|
|
|||||||||||||
|
x |
x |
|
|
|
|
|
|||||||||||
2u |
|
|
|
u |
f |
(x t)(x t)t f (x t). |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
t x |
|
|
|
|
|
|
|
|||||||||||
|
|
t |
x |
|
|
|
|
|
16
Все частные производные второго порядка имеют общий множитель f (x t), который, при подстановке частных производных в уравнение (3.1), можно вынести за скобки:
f |
|
2 |
(3.3) |
(x t)(A B C) 0. |
|||
Равенство (3.3) справедливо для любой функции f |
в том случае, если |
||
|
A 2 |
B C 0. |
(3.4) |
Уравнение (3.4), как и для обыкновенных линейных однородных дифференциальных уравнений, назовем характеристическим.
В зависимости от дискриминанта D B2 4AC квадратного уравнения, общее решение уравнения (3.1) может быть получено следующим образом.
1. D B2 |
4AC 0. Уравнение (3.4) |
имеет |
два |
различных |
|
вещественных |
корня |
1 2 . Им соответствуют |
частные |
решения |
|
линейного однородного дифференциального |
уравнения u1 f1(x 1t), |
||||
u2 f2(x 2t), |
где f1 |
и f2 – произвольные функции. Согласно свойству |
решений линейного дифференциального уравнения в частных производных (п. 2), функция u c1u1 c2u2 – также решение (3.1).
Итак, решение линейного однородного дифференциального уравнения
u f1(x 1t) f2(x 2t) |
(3.5) |
(произвольные постоянные c1 и c2, не ограничивая общности, для краткости записи можно опустить).
2. D B2 4AC 0.
В этом случае уравнение (3.4) имеет один вещественный корень
B (при A 0) кратности k = 2.
2A
По аналогии с обыкновенными линейными однородными дифференциальными уравнениями можно показать, что частными решениями уравнения (3.1) являются функции
u1 f (x t), u2 tf (x t).
17
Общее решение (3.1) |
|
u с1 f (x t) с2tf (x t) (с1 с2t) f (x t), |
(3.6) |
где f – произвольная функция.
Для доказательства справедливости формулы (3.6) достаточно
показать, что функция u2 tf (x t) |
является решением уравнения (3.1). |
|||||||||||||||
В самом деле, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
u2 |
|
f |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
t |
(x t) tf (x t), |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
u2 |
|
tf |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|||||||||
|
x |
(x t), |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
2u2 |
|
|
|
|
|
|
|
2 |
|
||||||
|
|
|
|
|
|
|
||||||||||
|
t2 |
|
|
f (x t) f |
|
(x t) tf |
(x t) |
|||||||||
2 f |
|
|
|
|
2 |
|
|
|
|
|
||||||
(x t) tf |
(x t), |
|
|
|||||||||||||
|
2u2 |
|
|
|
|
|
u2 |
|
|
|
|
|
|
|||
|
|
|
|
|
f |
|
t) tf |
|
||||||||
|
t x |
|
|
|
|
(x |
(x t). |
|||||||||
|
|
|
|
t x |
|
|
|
|
|
|
Подстановка найденных частных производных в уравнение (3.1) дает тождество
A(2 f |
|
2 |
|
|
|
|
|
) Ctf |
|
|
|
|
tf |
) B(f |
|
tf |
|
||||||
|
|
|
|
(2A B)f (A 2 B C)tf 0, |
|||||||
так как 2А В 0 |
при |
В |
, А 2 |
В С 0 при |
|||||||
|
|||||||||||
|
|
|
|
|
|
|
2А |
|
|
DB2 4AC 0.
3.D B2 4AC 0.
Вэтом случае корни характеристического уравнения (3.4)
комплексно-сопряженные |
i |
(i2 1, |
i – мнимая единица), а |
решения f (x ( i)t) |
– комплексные. Так |
как в физических задачах |
решение должно иметь реальный физический смысл, то комплексные решения не могут представлять для нас интереса. В качестве
18
вещественных решений рассматривают вещественную и мнимую часть комплексного решения
w u1 iu2 f1 x i t f2 x i t , f1, f2 – произвольные функции:
u c1u1 c2u2 ,
u1 Re f1 x i t f2 x i t , |
(3.7) |
||||
u |
Im f |
x i t f |
2 |
x i t . |
|
2 |
1 |
|
|
|
Физический смысл функций u1, u2 определяется содержанием реально решаемой задачи.
Задача 3.1. Найти общее решение уравнения Лапласа
|
2u |
|
|
2u |
0. |
(3.8) |
|
x2 |
y2 |
||||||
|
|
|
Решение. Функция, удовлетворяющая уравнению Лапласа, называется гармонической. Такое уравнение часто встречается при решении физических задач. Например, если рассмотреть стационарное распределение температуры по тонкой изотропной пластинке, изолированной от внешней среды (температура в точке М(x, y) пластинки зависит лишь от координат этой точки), то можно показать, что функция температуры u = f(x+ y) удовлетворяет уравнению Лапласа (3.8).
Очевидно, что уравнение Лапласа относится к типу линейных однородных дифференциальных уравнений в частных производных с постоянными коэффициентами. Его характеристическое уравнение2 1 0 имеет корни i.
Комплексное решение:
w f1(x iy) f2(x iy).
Вещественные решения:
u1 Re[ f1(x iy) f2(x iy)], u2 Im[ f1(x iy) f2(x iy)],
где f1, f2 – произвольные функции.
19
Например, если f |
(x iy) ex iy , f |
2 |
(x iy) (x iy)2 |
, то |
|
|
1 |
|
|
|
|
|
|
|
w ex iy (x iy)2. |
|
|
|||
Выделяем вещественную и мнимую часть функции w: |
|
|
||||
|
w exeiy x2 2xyi i2y2. |
|
|
|||
Применяя формулу Эйлера eiy cosy isin y и учитывая свойство i2 |
1, |
|||||
получаем, что w ex (cos y isin y) x2 |
2xyi y2 , откуда вещественная и |
|||||
мнимая части функции |
|
|
|
|
|
|
u ex |
cos y x2 y2, |
u |
2 |
ex sin y 2xy, |
|
|
1 |
|
|
|
|
|
и решение уравнения Лапласа может быть представлено в виде:
u c1 ex cos y x2 y2 c2 ex sin y 2xy .
Задача 3.2. Найти общее решение уравнения
2u |
|
2u |
12 |
|
2u |
0. |
|
t2 |
t x |
x2 |
|||||
|
|
|
Решение. Составляем характеристическое уравнение
2 12 0.
Его корни найдем, используя формулы Виета:
1 2 1, 1 2 12 1 4, 2 3.
Согласно (3.5), общее решение данного уравнения имеет вид:
u f1 x 4t f2 x 3t ,
где f1, f2 – произвольные функции.
Задача 3.3. Найти общее решение уравнения
2u |
8 |
2u |
16 |
2u |
0. |
|
t2 |
t x |
x2 |
||||
|
|
|
20