Дифференциальные уравнения. Учебное пособие-1
.pdf
4) u(0,t) A,
u(x,0) 0
Ответ:
U (x, p)
lim u(x,t) 0 ,
x
x 0 .
A |
e |
p x |
|
x |
||
|
|
|
; u(x,t) AErf |
|
. |
|
p |
|
|
|
2 t |
||
5) |
u(0,t) 0, |
u(l,t) A , |
|
|
|
|
|
|
|
|||||||
|
x |
|
0 x l . |
|
|
|
|
|
|
|
||||||
|
u(x,0) 0 |
|
|
|
|
|
|
|
|
|||||||
Ответ: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
U (x, p) |
Ach p x |
; |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
p ch |
pl |
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
( 1)n |
|
2n 1 2 |
|
|
|
|
|||
|
|
|
|
4 |
|
|
2 |
|
t |
|
(2n 1) |
|
||||
|
|
|
|
2l 2 |
|
|||||||||||
|
u(x,t) A 1 |
|
|
|
|
|
|
|
e |
|
|
|
cos |
|
x . |
|
|
|
|
2n 1 |
|
|
|
2l |
|||||||||
|
|
|
|
n 1 |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
101
СПРАВОЧНЫЙ МАТЕРИАЛ
1. Задача Штурма-Лиувилля.
X 2 X 0 – дифференциальное уравнение; X (0) X (l) 0 – граничные условия.
Разыскиваются числа n (собственные числа), при которых
существуют ненулевые решения дифференциального уравнения, удовлетворяющие граничным условиям, а также и сами ненулевые решения (собственные функции). Рассматриваются также задачи Штурма-
Лиувилля |
с граничными условиями |
|
|
|
|
|
X |
|
|
|
0 , |
|
X (l) 0 , |
X (0) |
(l) |
X (0) |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X (0) X (l) 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
Собственные функции, соответствующие граничным |
условиям |
||||||||||||
X (0) X (l) 0 , имеют вид X n (x) sin n x, |
|
|
n |
|
n |
, |
n 1,2,.... |
|
|||||
|
|
|
|
|
|
|
|
l |
|
|
|
|
|
Для остальных граничных условий собственные функции и |
|||||||||||||
собственные значения имеют вид ( n 0,1,2,... ): |
|
|
|
|
|
|
|
||||||
|
|
n |
|
n |
; |
|
|
|
|
|
|
||
X (0) |
X (l) 0 : X n (x) cos n x, |
l |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X (l) 0 : X n (x) cos n x, |
n |
|
|
(2n |
1) |
; |
|
|
|
|||
X (0) |
|
|
2l |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
n |
|
(2n |
1) |
|
. |
|
|
|
|||
X (0) X (l) 0 : X n (x) sin n x, |
|
|
2l |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2. Задачи для волнового уравнения.
а) Задача для волнового уравнения на отрезке 0,l с однородными граничными условиями.
utt а2uxx – дифференциальное уравнение;
u(x,0) (x), ut (x,0) (x) – начальные условия; u(0,t) u(l,t) 0 – граничные условия.
Рассматриваются также однородные граничные условия следующих видов:
ux (0,t) ux (l,t) 0 , u(0,t) ux (l, t) 0 , ux (0, t) u(l,t) 0 .
Решение этой задачи по методу Фурье получается в виде
u(x,t) X n (x) Tn (t) ,
n 1
где X n (x) – собственные функции задачи Штурма-Лиувилля с условиями,
соответствующими |
рассматриваемым |
граничным |
условиям; |
Tn (t) An cos a nt Bn sin a nt ; |
|
|
|
102
n – собственные числа задачи Штурма-Лиувилля;
An , Bn – коэффициенты, определяемые по начальным условиям.
б) Задачи для неоднородного волнового уравнения
utt а2uxx f (x,t) .
Решение задачи можно получить в виде разложения по собственным функциям задачи Штурма-Лиувилля
|
|
|
|
|
|
|
u(x,t) X n |
(x) un (t) , |
|||
|
n 1 |
|
|
|
|
|
|
|
2 |
2 |
|
|
a |
|
nun |
||
где un (t) – решения задач Коши un |
|
fn (t) , un (0) n , un (0) n ; |
|||
|
|
|
|
|
|
fn (t), n , n – |
коэффициенты |
разложений f (x,t) X n (x) fn (t) , |
|||
|
|
|
|
|
n 1 |
|
|
|
|
||
(x) X n (x) n |
, (x) X n (x) n . |
|
|
|
|
n 1 |
n 1 |
|
|
|
|
3. Задачи для уравнения теплопроводности
а) Задача для уравнения теплопроводности на отрезке 0,l с однородными граничными условиями.
ut а2uxx – дифференциальное уравнение; u(x,0) (x) – начальное условие;
u(0,t) u(l,t) 0 – граничные условия.
Рассматриваются также однородные граничные условия следующих видов:
ux (0,t) ux (l,t) 0 , u(0,t) ux (l, t) 0 , ux (0, t) u(l,t) 0 .
Решение этой задачи по методу Фурье получается в виде
u(x,t) X n (x) Tn (t) ,
n 1
где X n (x) – собственные функции задачи Штурма-Лиувилля с условиями,
соответствующими |
рассматриваемым |
граничным |
условиям; |
|||
T (t) С |
n |
e a2 2nt ; |
|
|
|
|
n |
|
|
|
|
|
|
n |
– собственные числа задачи Штурма-Лиувилля; |
|
||||
Cn |
– коэффициенты, определяемые по начальным условиям. |
|
||||
б) Задача для неоднородного уравнения теплопроводности ut а2uxx f (x, t) .
103
Решение задачи можно получить в виде разложения по собственным функциям задачи Штурма-Лиувилля
|
|
|
u(x,t) X n |
(x) un (t) , |
|
n 1 |
|
|
где un (t) – решения задач Коши un a2 2nun fn (t) , un (0) |
n ; |
|
|
|
|
fn (t), n – коэффициенты разложений f (x,t) Xn(x) fn(t), (x) Xn(x) n . |
||
|
n 1 |
n 1 |
4. Задачи для волнового уравнения и уравнения теплопроводности с неоднородными граничными условиями
u(0,t) A(t), u(l,t) B(t) .
Каждая из этих задач сводится к задаче с однородными граничными
условиями |
для |
функции |
v(x,t) u(x,t) w(x,t) , |
где |
w(x,t) A(t) |
B(t) A(t) x . |
|
|
|
|
l |
|
|
|
Решение получается в виде u(x,t) v(x,t) w(x,t) .
5. Задачи для уравнения Лапласа а) Краевая задача для уравнения Лапласа в круговом секторе
0 r R, 0 ( r, – полярные координаты, 2 ).
u urr 1r ur r12 u 0 – дифференциальное уравнение;
1)u (R, ) f ( ) – граничное условие (на окружности r R ),
2)u (r,0) u (r, ) 0 – граничные условия (на лучах 0, ).
Вместо условий 2) на лучах рассматриваются также условия
3) u (r,0) u (r, ) 0 , u (r,0) u (r, ) 0 , u (r,0) u (r, ) 0 .
Решение задачи по методу Фурье получается в виде
u(r, ) Rn (r) n ( ) ,
n 1
где n ( ) – собственные функции задачи Штурма-Лиувилля для
дифференциального уравнения 2 0 с условиями, соответствующими рассматриваемым граничным условиям вида 2) или 3);
Rn (r) Cnrn ; Cn – коэффициенты, определяемые по граничным условиям 1).
104
б) Задача Дирихле для уравнения Лапласа в круге 0 r R, |
0 2 |
||
( r, – полярные координаты). |
|
||
u urr 1 ur |
1 |
u 0 – дифференциальное уравнение; |
|
|
|||
r |
r 2 |
|
|
u (R, ) f ( ) – граничное условие.
Решение этой задачи по методу Фурье получается в виде
u(r, ) A0 An cos n Bn sin n r n ,
n 1
где An , Bn – коэффициенты, определяемые по граничному условию.
105
УКАЗАНИЯ К РЕШЕНИЮ ЗАДАЧ
Указания к задачам 2, 3 из сборника [29]
Указание к задаче 2. Найти общее решение уравнения, приведя его к каноническому виду 49uxx 14uxy uyy 14ux 2uy 0 .
Решение. Характеристическое уравнение имеет вид
|
|
49dy2 14dyx dx 2 |
0 |
(a |
|
49, a |
|
7, a |
1) |
|
|
||||||
|
|
|
|
|
|
|
11 |
|
12 |
|
|
22 |
|
|
|
||
Имеем: |
a2 |
a |
a 0, |
значит, |
уравнение |
имеет |
параболический |
тип. |
|||||||||
|
12 |
11 |
22 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Характеристическое уравнение равносильно уравнению 7dy dx 0 , |
его |
||||||||||||||||
решение есть |
7y x C . |
Значит, |
следует |
ввести новую переменную |
|||||||||||||
7y x . |
Переменную |
|
вводят |
произвольно, |
проверяя, |
однако, |
|||||||||||
выполнение условия независимости |
|
и : x y |
y x |
0 . Полагаем, |
|||||||||||||
например, |
x . |
Получаем: |
x 1, |
|
xx |
xy |
yy 0, |
x |
1, |
||||||||
y yy |
xx xy 0 ; далее, используя формулы (2.24), имеем: |
|
|
||||||||||||||
ux u |
u , uy 7u , uxy 7u 7u , uxx 2u u |
u , uyy 49u . |
|||||||||||||||
Подставляя последнее соотношение в исходное уравнение, получаем уравнение в канонической форме:
7u 2u 0.
При всяком фиксированном это есть обыкновенное линейное дифференциальное уравнение с постоянными коэффициентами; корни его
характеристического |
уравнения равны |
k1 0, k2 |
|
2 |
, поэтому общее |
||||
|
|
|
|
2 |
|
|
7 |
|
|
решение имеет вид |
u( , ) C1 ( ) C2 |
|
, |
где C1 ( ),C 2 ( ) – |
|||||
( ) exp |
7 |
|
|||||||
|
|
|
|
|
|
|
|
||
произвольные функции переменной . Возвращаясь к прежним
переменным, получаем общее решение уравнения:
u(x, y) C1 (7x y) C2 (7x y) ,
где C1 ( ) и C2 ( ) – произвольные дважды дифференцируемые функции одной переменной x 7y .
Указание к задаче 3. Найти общее решение уравнения, приведя его к каноническому виду
2uxx 5uxy 3uyy 0.
Решение. Имеем: a11 2, a12 5
2, a22 3, a122 a11a22 49
4 0 . Значит,
уравнение имеет гиперболический тип. Запишем характеристическое
106
уравнение |
dy |
|
1 |
|
5 |
|
25 |
( 3) |
|
|
1 5 |
|
7 |
. Интегрируя, находим два |
dx |
2 |
|
2 |
4 |
|
|
|
|||||||
|
|
|
|
|
|
|
2 2 |
|
2 |
|
семейства характеристик |
y 3x C1, y |
1 x C2 , где |
C1,C 2 – |
|
|
2 |
|
произвольные постоянные. Значит, следует ввести новые независимые
переменные |
y 3x, y |
1 x . Пересчитаем производные: |
||||||||||||||
|
|
3u 1 u , u |
|
|
2 |
|
|
|
|
|
|
5 u |
1 u , |
|||
u |
|
|
u u , u |
3u |
|
|
||||||||||
|
x |
|
2 |
|
|
y |
|
|
|
|
xy |
|
|
2 |
2 |
|
u |
xx |
9u |
3u |
|
|
1 u |
, u |
yy |
u |
u |
u . |
|
||||
|
|
|
|
4 |
|
|
|
|
|
|
|
|||||
Подставляя |
uxx , uxy , uyy |
|
в исходное уравнение, приведем его к |
|||||||||||||
каноническому виду u |
0, или |
|
|
|
|
|
|
|
|
|||||||
u 0.
Интегрируя по переменной , имеем u ( ) , где ( ) – произвольная
функция. Интегрируя еще раз, но уже по переменной , находим u( , ) ( )d ( ) или u( , ) ( ) ( ) ,
где ( ) ( )d , ( ) – произвольные функции. Итак, общее решение исходного уравнения имеет вид
u(x, y) ( y 12 x) ( y 3x) ,
где ( y |
1 x), ( y 3x) – |
произвольные |
дважды дифференцируемые |
|
2 |
1 x, y 3x |
|
функции переменных y |
соответственно. |
||
|
|
2 |
|
Указания к задачам 1-15 из сборника [13]
Указание к задаче 1. |
|
|
|
|
Решить |
начально-краевую |
задачу |
(нестационарную |
задачу |
теплопроводности для однородного тонкого стержня) |
|
|||
|
ut 9uxx ; |
|
(2.82) |
|
|
u(x,0) 31sin 3 x, |
x (0,9); |
(2.83) |
|
|
u(0,t) u(9,t) 0, |
t 0. |
(2.84) |
|
|
|
107 |
|
|
Условие (2.83) задает распределение температуры в начальный момент времени, условия (2.84) задают температуру (нулевую) на концах стержня. Требуется найти распределение температуры по длине стержня в произвольный момент времени.
Решение.
Согласно методу разделения переменных ищем частное решение уравнения (2.71) в виде u(x,t) X (x) T (t) . Найдем частные производные:
|
|
|
|
|
|
|
|
uxx X (x) T (t), ut X (x) T (t) . |
X (x) T (t) 9 X |
(x) T (t) . |
Поделив обе |
||||
Подставим в |
уравнение (2.71): |
||||||
|
|
|
|
|
|
|
|
части уравнения на 9X (x) T (t) , получим |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
X (x) |
|
T (t) |
const . |
(2.85) |
|
|
|
|
|
||||
|
|
X (x) 9T (t) |
|
|
|
||
Уравнения (2.85) равносильны системе: |
|
|
(2.86) |
||||
|
|
X (x) X (x) 0 , |
|
||||
|
|
|
|
|
|
|
(2.87) |
|
|
T (t) 9 T (t) 0. |
|
||||
|
|
|
|
|
|
|
|
Согласно граничным условиям X (0) X (9) 0.
Таким образом, требуется найти ненулевые решения уравнения (2.86), удовлетворяющие условиям X (0) X (9) 0. Уравнение (2.86) есть
линейное однородное уравнение второго порядка. Его характеристическое уравнение k 2 0 имеет корни k
. Рассмотрим три случая:
1) |
0. Общее решение уравнения (2.86) имеет вид |
X (x) Ax B , |
где |
|||||||||||
A, B const . Удовлетворяя условиям X (0) X (9) 0, |
получим |
A B 0 . |
||||||||||||
Следовательно, получили тривиальное решение X (x) 0 . |
|
|
|
|
||||||||||
2) |
0. |
Общее |
решение |
|
уравнения |
(2.86) |
имеет |
|
вид |
|||||
X (x) D e |
x D |
e x |
, |
D , |
D |
2 |
const . Удовлетворяя |
условиям |
||||||
|
1 |
2 |
|
|
1 |
|
|
|
|
|
|
|
|
|
X (0) X (9) 0, получаем |
D1 0 , |
D2 0 , X (x) 0 – |
опять тривиальное |
|||||||||||
решение. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3) |
0. |
Общее |
решение |
|
уравнения |
(2.75) |
имеет |
|
вид |
|||||
X (x) D3 cos x D4 sin |
|
x , |
D3 , D4 const . |
Удовлетворяя |
условиям |
|||||||||
X (0) X (9) |
0, получаем |
D3 0 , |
|
|
|
|
|
n |
2 |
где |
||||
D4 sin 9 0 . Отсюда |
|
, |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
n |
– произвольное натуральное |
число. Значит, |
X (x) D4 sin n x |
– |
это |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
9 |
|
|
|
решение уравнения (2.86), удовлетворяющее условиям X (0) X (9) 0. |
|
|||||||||||||
|
Общее решение уравнения (2.87) при 0 |
имеет вид T (t) Me t , |
||||||||||||
где M – произвольная константа. |
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
108 |
|
|
|
|
|
|
|
Таким образом, функция u(x,t) имеет вид |
|
|||
|
|
u(x,t) C e t sin |
x , |
(2.88) |
|
n |
2 |
|
|
где |
|
, C MD4 – произвольная постоянная. . |
|
|
|
9 |
|
|
|
В силу линейности уравнения (2.82), решением его является линейная комбинация частных решений вида
|
nx |
|
n |
2 |
|
|
u(x,t) Cn sin |
e |
|
|
t . |
(2.89) |
|
9 |
|
|||||
n 1 |
9 |
|
|
|
|
|
Это решение удовлетворяет граничным условиям (2.84).
Для нахождения коэффициентов Cn воспользуемся начальным условием
(2.83).
u(x,0) Cn sin nx |
31sin 3 x . |
|
|
|
|
n 1 |
9 |
|
Полученное соотношение есть разложение в ряд Фурье по синусам |
|||||||||||||||||
функции |
|
u(x,0) 31sin 3 x |
на |
отрезке |
0,9 . |
Найдем коэффициенты |
|||||||||||
разложения |
|
nx |
|
|
|
|
|
nx |
|
|
|
||||||
Сn |
2 |
9 |
|
|
|
dx |
62 9 |
sin 3 x sin |
dx |
|
|||||||
9 |
u(x,0)sin |
9 |
|
9 0 |
9 |
|
|||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|||||
|
31 9 |
|
|
3 x |
nx |
|
|
|
nx |
|
|
0, |
n 27; |
||||
9 |
|
cos |
|
|
cos 3 x |
9 |
|
dx |
31, |
n 27. |
|||||||
|
|
|
|
|
|
9 |
|
|
|
|
|
|
|||||
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Подставляя полученные значения в (2.89), находим решение задачи
(2.82)-(2.84)
u(x,t) 31sin 3 x e 81 2t . |
|
Указание к задаче 2. |
|
Решить начально-краевую задачу (теплопроводности) |
|
ut 2uxx ; |
(2.90) |
u(x,0) 31cos3 x cos 4 x; |
(2.91) |
ux (0,t) ux (6,t) 0. |
(2.92) |
Условия (2.92) означают, что концы стержня теплоизолированы. |
|
Решение.
Аналогично задаче 1 ищем частное решение уравнения (2.90) в виде u(x,t) X (x) T (t) . Получим систему:
X (x) X (x) 0 , |
(2.93) |
||
|
|
|
(2.94) |
T (t) 2 T (t) 0. |
|
||
|
(6) |
0. |
|
Согласно граничным условиям X (0) X |
|
||
|
|
|
|
109 |
|
|
|
Характеристическое уравнение для (2.93) имеет вид k 2 |
0 , его корни |
|||
k . Рассмотрим три случая: |
|
|
X (x) Ax B , где |
|
1) 0. Общее решение уравнения (2.93) |
имеет вид |
|||
A, B const . Удовлетворяя условиям |
X (0) X (6) 0, получим |
A 0 , |
||
|
|
|
|
|
X (x) B . При этом общее решение уравнения (2.94) при 0 имеет вид
T (t) E const . Тогда u(x,t) B E C0 . |
|
|
|
|
|
|
||||||||||||
2) |
0. |
Общее |
|
решение |
|
|
уравнения |
(2.93) |
имеет |
вид |
||||||||
X (x) D e x D |
e |
x |
, |
D , D |
2 |
const . |
Удовлетворяя |
условиям |
||||||||||
|
|
1 |
2 |
|
|
|
|
1 |
|
|
|
|
|
|
|
|||
|
|
|
0, получаем D1 |
D2 |
|
0 , X (x) 0 – тривиальное решение. |
||||||||||||
X (0) |
X (6) |
|
||||||||||||||||
3) |
0. |
Общее |
|
решение |
|
|
уравнения |
(2.93) |
имеет |
вид |
||||||||
X (x) D3 cos |
x D4 sin |
x , |
D3 , D4 const . |
Удовлетворяя |
условиям |
|||||||||||||
|
|
|
0, получаем D4 |
|
|
|
D3 sin 6 |
|
|
|
n 2 |
|||||||
|
0 , |
|
0 . Отсюда |
, |
||||||||||||||
X (0) |
X (6) |
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
nx |
|
6 |
|
где n – произвольное натуральное число. Значит, X (x) D |
cos |
– это |
||||||||||||||||
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0, |
|
|
|
|
|
||
решение |
уравнения |
(2.93) |
при |
|
|
удовлетворяющее |
условиям |
|||||||||||
|
|
|
0. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
X (0) |
X (6) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Общее решение уравнения (2.94) при 0 |
имеет вид T (t) Me 2 t , |
|||||||||||||||||
где M |
– произвольная |
константа. Тогда |
u(x,t) C e 2 t |
cos |
x , |
где |
||||||||||||
n |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
, C MD3 – произвольная постоянная. |
|
|
|
|
|
|||||||||||
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Решение уравнения (2.90), удовлетворяющее граничным условиям (2.92), можно записать в виде линейной комбинации полученных частных решений
|
|
|
nx |
|
n 2 |
|
|
|
2 |
t |
|||
u(x,t) C0 |
Cn cos |
e |
|
6 . |
||
|
n 1 |
|
|
6 |
воспользуемся начальным |
|
Для нахождения коэффициентов |
C0 , |
Cn |
||||
условием (2.91) |
nx |
|
|
|
|
|
|
31cos3 x cos 4 x . |
|||||
u(x,0) C0 Cn cos |
|
|||||
n 1 |
6 |
|
|
|
|
|
Полученное соотношение есть разложение в ряд Фурье по косинусам |
||||
функции u(x,0) 31cos3 x cos 4 x |
на отрезке |
0,6 . Найдем |
||
коэффициенты разложения C0 , Cn : |
|
|
||
C0 2 |
6 |
(31cos3 x cos 4 x)dx 0; |
|
|
6 |
0 |
|
|
|
110
