- •В.П. Первадчук, Е.М. Кадырова, В.Ю. Соколов
- •УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ: МЕТОДЫ РЕШЕНИЯ ЗАДАЧ
- •Глава 1. КЛАССИФИКАЦИЯ УРАВНЕНИЙ ВТОРОГО ПОРЯДКА. ОБЩЕЕ РЕШЕНИЕ УРАВНЕНИЯ
- •1.2. Классификация квазилинейных уравнений второго порядка с двумя независимыми переменными
- •1.3. Приведение уравнений к каноническому виду
- •1.4. Упрощение уравнения в каноническом виде
- •2.1. Нахождение общего решения
- •{мдydy = jo dy,
- •3.1. Распространение волн в бесконечной струне. Задача Коши
- •Рекомендуем решить:
- •3.2. Полуограниченная прямая. Метод ограничений
- •4.1. Принцип суперпозиции
- •7.2. Метод Фурье для уравнений гиперболического типа. Неоднородные задачи
- •7.3. Метод Фурье для уравнений параболического типа
- •7.5. Понятие функций Бесселя
- •(xVXf
- •7.6. Метод Фурье для многомерных задач
- •УРАВНЕНИЯ МАТЕМАТИЧЕСКОЙ ФИЗИКИ: МЕТОДЫ РЕШЕНИЯ ЗАДАЧ
u(x,t) = |
* h M x* - M |
x + ' ( t m ± m ± M |
£ x , (dT. |
|||
|
2IT |
T |
JJ |
Tp |
|
^ |
|
+ a 0 +/3Qe+ jgo{rXt-T)dr |
|
|
|||
|
|
|
о |
|
|
|
+ |
cos— |
+ p n sin— |
+ hn{t) |
^ |
Tthx |
|
cos---- |
l
Здесь a„,fin,hn определяются по формулам (7.2.4)
7.3. Метод Фурье для уравнений параболического типа
Задача 7.9. Тонкий однородный стержень, на боковой поверхности кото рого происходит теплообмен с окружающей средой, а конец х = / - теп лоизолирован. Найти распределение температур, если начальная темпера тура V - const.
ul =a2uxx-flu, 0 <х<1, t> О, ux(0,t)-hu(Q,t) = 0,
ы(х,0) =V, h> 0.
Решение ищем в виде u(x,t) = X (x)r(t). Для нахождения Х(х) поставим задачу Штурма - Лиувилля.
ГХ" + ЛХ = 0,
\ аг'(о) - аяг(о) = А"(/)=0.
Собственные значения Лп ={р„)2, где р п - положительные корни урав нения Яtg Я/ = h , а собственные функции Х п(х) = р пcos /лпх +hsin р пх ,
00
далее подставим u(x,t)= ^ T n(t)Xп(х) в исходное уравнение и граничное
Л=1
условие. Y jTnX n - ° 2Y^TnX n - P 'E jn X n ~ Таким образом, получаем задачу для нахождения Тп (t).
|
т:+а2м Х - /Я „ = 0 , |
|
|
Т„ = |
|
1 ^ |
cosf l n x + Asm fi„ x )d x = - r - |
2V |
т„(о)= A J = -— -Г |
-----rr— . |
|
Wt . |
l(h2 + K ) + h |
|
Итак, u(x,t) = 2FAУ — г t |
1—ГЛ— .е-("2Л!-/<>(и cos finx + Asm ft x). |
Z l M „W ! + ^ J + AJ
Задача 7.10. Найти распределение температур в стержне, концы которого поддерживаются при нулевой температуре, с боковой поверхности про исходит лучеиспускание в окружающую среду, в стержне имеется тепло-
7DC |
|
|
|
|
|
|
вой источник sin — . Начальная температура равна нулю. |
||||||
и( = а |
2 |
|
|
п |
■ Я* |
|
|
ихх - ри +sm— , |
|||||
и\ п =и\ |
, =0, |
«I |
= 0, |
|||
!дг=0 |
|
I |
х=1 |
’ |
1/=0 |
’ |
х я + л х = о,
х(о)= *(/)= о,
|
я"=(т ) |
^ W =sin^ f ’ п=1’2»---» |
|
Е З Д , - 5 > 2М + Е / В Д . = sinf • |
|
Таким образом, |
+ (а2Лп - /з)гп = | ^ ” * |
|
|
2 |
2 |
Ti +\ ^ - ~ P I7! =1> Г,(0)=0.
л„г
ап
Отсюда |
Ш = А е |
~Р |
|
|
J + |
„ 2 —2 |
|
||
|
|
|
а к |
- Р |
|
|
|
|
|
Из начального условия А{ = ------- - |
- Д 2 |
|
|
|
|
a V |
|
|
|
|
|
|
„2„2 |
|
|
|
|
а к |
|
|
/ 2 |
1 -е |
-Р |
. лх |
|
а2я 2 -/Я 2 |
|
sin— . |
|
|
|
|
/ |
7.4. Метод Фурье для уравнений параболического типа (неоднородные задачи)
Задача 7.11.
Uf^a^Uxx, 0 <х<1, О О,
M* L ~ At>
их\х=1 = Т '
“1,-0 = Т
Решение исходной задачи ищем в виде u(x,t) = v(x, t)+ w(x, t) , где w(x,t) удовлетворяет неоднородным граничным условиям
w (x,t)= y\Х2 + у 2х + у 3^ + {з^х2 +S2x + S3)r,
тогда
w, 1 ^ = (2У\х + у2 )Ж +(2^iх + 82)г|^ = y2At + S2T = A t,
wx\x=i =(2r\x + r2)At + (2$\x + 52)r \x=i = (2y{l + y2)At + (28J + S2)r = T
Наиболее удобное решение с нашей точки зрения
w(x,0 = - x 2(T-A t)+ Axt
Получаем первую задачу для определения v(x,t) с однородными гранич ными условиями.
2 |
( |
|
|
2 |
\ |
2 |
|
|
Г |
|
1 2 |
л |
а |
|
2 ^ ~ А?''i |
|
vt = a |
|
|
|
wxx)=a |
|
Vxx+ \^Y 1X A + A x~ a |
— j ~ у |
|||||||||
|
v ,= a |
2 |
|
A |
2 |
|
J |
|
|
2 T ~ At |
|
|
|
|
||
|
|
|
v „ + — x |
|
-A x + a |
— -— , |
|
|
|
|||||||
|
* |
|
|
” |
21 |
|
|
|
|
|
|
l |
|
|
|
|
|
vJ |
|
= 0, |
vJ |
|
. = 0, |
|
|
|
|
|
|
|
|||
|
*1*=0 |
|
|
5 |
x \x=l |
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
f |
2 \ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1- |
— ' |
|
|
|
Составим задачу Штурма - Лиувилля: |
|
|
|
4 |
2 lJ |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
fX" + AX = 0, |
|
|
|
|
|
|
|
|||||
|
|
|
|
tr (o )=*'(/)=о. |
|
|
|
|
|
|||||||
|
|
|
|
|
rmS1 |
П = 0,1,2,.. |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
1 , |
ПФ О, |
|
|
|
|
|
X n = c o s ? f, |
|
\\ХпГ=< I |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
I, |
п = 0. |
|
|
|
|
|
|
|
|
v (x ,» )= 2 r A |
|
' |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
п=1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Г |
|
„2 'N |
а2 {Г -A t) |
|
||||
|
|
|
|
|
|
|
|
Х-- |
2/ |
|
|
|
|
|
||
при и>1 |
|
|
|
~ / Г |
|
|
|
2 ^ |
2 |
(T - A tj |
|
7ШХ |
, |
|||
7’; + а 2Я„7’„ = у | |
|
|
|
21 |
|
a |
|
|
cos---- |
ах, |
||||||
|
|
|
|
* |
о . |
|
|
|
|
|
|
|
|
I |
|
|
|
I f |
|
|
2 \ |
|
т х |
, |
|
l ^ |
|
|
|
|
|||
|
|
|
|
|
cos---- |
ах = |
|
|
|
|
|
|
||||
|
i l |
|
21 j |
|
|
l |
|
|
\ m ) |
|
|
|
|
a2(T -A t) l( 7 nx j |
|
— ------ |
Jcos-y- dx = 0, |
при |
п = О |
|
|
i ' j |
|
|
, a2(T -A t)] dx = |
||||
|
|
|
|
2 0J l |
2'J |
|
l |
|
|
||
|
|
|
|
|
|
|
|
||||
|
1 Ах2 |
|
Ах3 |
а2 (Г - A t) |
|
^ Л / |
Л/2 |
| а2(Г -A t) |
|||
|
----------------------------+ |
— |
-------------- |
|
|
|
|
|
|||
|
/ 22 |
6 |
6 |
|
/ |
|
|
|
|
|
|
Из начальных условий |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
mix |
, |
/ ,\n+i 2/ |
||
|
|
|
|
|
-----х cos— |
dx = ( - l f |
|||||
|
» - |
|
т ] |
Н |
2/ |
|
|
1 |
|
|
[m) |
при |
n = 0 |
|
|
|
|
X3^ |
|
= 1 - 1 . |
|
|
|
|
r0(oM X ------ |
|
|
|
|||||||
|
|
|
|
|
|
6/ ' |
|
|
6 |
|
|
Таким образом, получаем две задачи: |
|
|
|
|
|
||||||
|
|
|
п, |
|
I/ та |
\ 2 |
|
|
|
п>1, |
|
|
|
|
Г» Ч / |
Г. = - ^ г , |
|||||||
|
|
|
" “ W 2 |
|
|
||||||
|
|
|
|
|
|
i«+i |
2/ |
|
|
|
|
|
|
|
r„ (o )-(-i) |
|
;2 * |
|
|
||||
|
|
|
|
|
|
м |
|
|
|
|
|
|
|
|
|
|
-4/ |
Л/2 |
. а2{Г-A t) |
|
|||
|
|
|
г» = т |
' “ б” + |
|
; |
|
■ |
|||
|
|
|
г0(о ) - ‘ - т - |
|
|
|
|
|
|||
Далее достаточно решить две задачи Коши |
|
|
|||||||||
|
|
|
|
|
|
(яигЛ2 |
|
|
|
|
|
|
Г .Й - |
^ |
|
i- Д —J* + 2(-l)”+1 |
Tl |
||||||
|
|
M V |
|
|
|
|
|
|
W * |
||
|
|
|
|
|
Al |
a2T , |
a h 2 ,, |
l |
|||
|
|
|
r , « - T '+ — « - — • |
|
|||||||
|
|
w(x, t) = w(x, t) + TQ{ t ) x 0(x) + ^ |
T„ |
( 4 |
n=1