- •Введение
- •1.1 Задача, приводящая к уравнению с частными производными
- •1.2 Определение дифференциального уравнения с частными производными
- •1.3 Решение УЧП
- •2.2 Типы линейных (квазилинейных) УЧП второго порядка
- •2.4 Простейшие виды задач для УЧП
- •2.4.1 Задачи для гиперболических УЧП
- •2.4.2 Задачи для параболических УЧП
- •2.4.4 Требования к начальным и граничным условиям и правым частям УЧП
- •2.4.3 Задачи для эллиптических УЧП
- •3 Гиперболические УЧП
- •3.1 Задача Коши для бесконечной струны
- •3.1.1 Решение задачи Коши методом Даламбера
- •3.1.2 Физическая интерпретация формулы Даламбера
- •3.2 Смешанная задача для свободных колебаний струны с закрепленными концами
- •3.3 Смешанная задача для вынужденных колебаний струны с закрепленными концами
- •3.4 Смешанная задача для неоднородных УЧП, НУ и ГУ
- •4 Параболические УЧП
- •4.1 Смешанная задача для однородного уравнения теплопроводности и однородных ГУ (задача об остывании стержня)
- •4.2 Смешанная задача для неоднородного уравнения теплопроводности и однородных НУ и ГУ
- •4.4 Смешанная задача для неоднородных УЧП, НУ и ГУ
- •4.5 Задача Коши для бесконечного стержня
- •4.6 Физический смысл фундаментального решения
- •5 Эллиптические УЧП
- •5.1 Задача Дирихле для уравнения Лапласа в круге
- •5.2 Интегральная формула Пуассона
- •6 Применение преобразования Лапласа для решения задач математической физики
- •Список использованных источников
v(x,t)= v1(x,t)+ v2 (x,t).
Функция v1(x,t) - решение задачи (см. пункт 4.1)
УЧП: vtt = a2vxx (0 < x < L, t > 0);
НУ: v(x,0)=ϕ~(x) (0 ≤ x ≤ L);
ГУ: v(0,t)= 0 , v(L,t)= 0 (t ≥ 0).
Функция v2 (x,t) - решение задачи (см. пункт 4.2)
УЧП: vtt = a2vxx + ~f (x,t) (0 < x < L, t > 0);
НУ: v(x,0)= 0 (0 ≤ x ≤ L); vt (x,0)= 0 (0 ≤ x ≤ L);
ГУ: v(0,t)= 0 , v(L,t)= 0 (t ≥ 0).
4.5 Задача Коши для бесконечного стержня
Задача Коши для бесконечного стержня без источников тепла имеет
вид:
УЧП: ut = a2uxx (−∞ < x < +∞, t > 0);
НУ: u(x,0)=ϕ(x) (−∞ < x < +∞).
75
Применим метод разделения переменных. Для этого будем искать ее решение в виде
u(x,t)= X (x)T (t).
Подставив это решение в УЧП, получим
X (x)T ′(t)= a2 X ′′(x)T (t).
Разделим полученное уравнение на a2 X (x)T (t).
′ |
|
′′ |
|
|
|
T (t) |
= |
X (x) |
= −µ . |
(24) |
|
a2T (t) |
X (x) |
|
|||
Из соотношения (24) получаем два обыкновенных дифференциальных уравнения:
T ′+ µa2T = 0 , |
(25) |
X ′′+ µX = 0 . |
(26) |
Решив обыкновенное дифференциальное уравнение (25) методом разделения переменных, получим
T = Ce−µa2t . |
(27) |
Так как из физических соображений, очевидно, что температура точек стержня при t → +∞ не может неограниченно возрастать, то тогда µ ≥ 0 .
Для упрощения вида получаемых результатов, введем обозначение
µ = λ2 .
Рассмотрим теперь уравнение (26)
76
X ′′+λ2 X = 0. |
(28) |
При λ = 0 оно имеет общее решение (пункт 3.2)
X = C1x +C2 . |
(29) |
Из физических соображений, очевидно, что при x → ∞ температура точек стержня не может неограниченно возрастать. Отсюда C1 = 0 , и решение (29) примет вид
X = C2 . |
(30) |
При λ ≠ 0 уравнение (27) имеет общее решение (пункт 3.2)
X = C1 cosλx +C2 sinλx . |
(31) |
В это общее решение при λ = 0 входит и решение (30). Из (27) и (31) получаем решение УЧП
u(x,t)= X (x)T (t)= Ce−a2λ2t (C1 cosλx +C2 sinλx)=
= e−a2λ2t (C C1 cosλx +C C2 sinλx)=
►Введем обозначения A = C C1 , B = C C2 .◄
= e−a2λ2t (Acosλx + Bsinλx).
77
Полученное решение зависит от λ , тогда в дальнейшем будем использовать обозначения uλ (x,t). Кроме того, произвольные постоянные для каждого значения λ могут быть различны, поэтому выберем для них обозначения A(λ) и B(λ)
uλ (x,t)= e−a2λ2t (A(λ)cosλx + B(λ)sinλx).
Так как стержень бесконечен, то для уравнения (26) нет граничных условий как в пункте 3.2. Тогда на параметр λ не будут накладываться дополнительные ограничения, и он может принимать любые значения от −∞ до + ∞.
Так как при любом значении λ (−∞,+∞) функция uλ (x,t) является решением однородного УЧП, то естественно предположить, что его решением также будет являться и
|
|
|
+∞ |
+∞ |
|
|
|
|
|
|
u(x,t)= ∫uλ (x,t)dλ = ∫ e−a2λ2t (A(λ)cosλx + B(λ)sinλx)dλ. |
(32) |
|||||||
|
|
|
−∞ |
−∞ |
|
|
|
|
|
Найдем функции A(λ), |
B(λ), используя НУ u(x,0)=ϕ(x) . |
|
|||||||
|
|
|
+∞ |
|
|
|
|
|
|
|
|
u(x,0)= ∫(A(λ)cosλx + B(λ)sinλx)dλ =ϕ(x). |
|
(33) |
|||||
|
|
|
−∞ |
|
|
|
|
|
|
Представим функцию ϕ(x) интегралом Фурье [9, с. 583]: |
|
||||||||
1 |
+∞ |
+∞ |
|
|
+∞ |
|
|
|
|
|
|
∫ |
∫ϕ(ξ)cosλξdξ cosλx + |
∫ϕ(ξ)sinλξdξ |
sinλx dλ =ϕ(x). (34) |
||||
|
2π |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
−∞ −∞ |
|
−∞ |
|
|
|
||
Из (33) и (34) имеем
78
|
|
|
|
1 |
|
+∞ |
|
|
1 |
+∞ |
|
|
|
|
||
A(λ)= |
|
|
|
∫ϕ(ξ)cosλξdξ , A(λ)= |
∫ϕ |
(ξ)cosλξdξ . |
|
(35) |
||||||||
2π |
2π |
|
||||||||||||||
|
|
|
−∞ |
|
|
−∞ |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Подставим выражения (35) в (32) |
|
|
|
|
|
|
|
|||||||||
+∞ |
2 2 |
|
1 |
+∞ |
|
|
|
1 |
|
+∞ |
|
|
||||
u(x,t)= ∫ |
e−a λ t |
|
|
∫ϕ(ξ)cosλξdξ cosλx + |
|
|
∫ϕ(ξ)sinλξdξ sinλx dλ = |
|||||||||
|
2π |
2π |
||||||||||||||
−∞ |
|
|
|
|
|
−∞ |
|
|
|
−∞ |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
1 |
|
+∞ |
|
+∞ |
|
|
|
|
|
|
|
|
|
= |
|
|
|
|
∫ e−a2λ2t |
∫ϕ(ξ)(cosλξ cosλx +sin λξ sin λx)dξ dλ = |
|
||||||||
|
|
2π |
|
|
||||||||||||
|
|
|
−∞ |
|
−∞ |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
►Используем формулу косинуса разности
cos(α − β)= cosα cos β −sinαsin β .◄
|
1 |
|
+∞ |
|
+∞ |
|
||
= |
|
∫ e−a2λ2t |
|
∫ϕ(ξ)cos(λξ −λx)dξ dλ = |
||||
2π |
|
|||||||
|
|
−∞ |
|
|
|
|
||
|
|
|
|
−∞ |
|
|||
|
1 |
|
+∞ |
|
|
+∞ |
|
|
= |
|
∫ e−a2λ2t |
∫ϕ(ξ)cosλ(ξ − x)dξ dλ = |
|||||
|
2π |
|
||||||
|
|
|
−∞ |
|
−∞ |
|
||
|
|
|
|
|
|
|||
►Поменяем порядок интегрирования.◄
= 21π +∫∞ϕ(ξ) +∫∞ e−a2λ2t cosλ(ξ − x)dλ dξ =
−∞ −∞
►Используем, что
79
