Презентации и конспект лекции Синицын / Презентации_ВЭиМОвТС / ВычЭксперим / Лекция3_ДифУравнения
.pdfПример точного решения ДУ (продолжение3)
u |
|
1; |
u |
|
x 1 0; |
|
g ( x ) 1; |
f ( x ) 1 . |
|
0 |
0; q |
0 |
1; |
1 |
1; q |
1 |
0 |
|||||
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|||||||||||||||
x |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
x 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
u |
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
• |
Имеем сразу |
c1 |
g ( 0 ) |
|
|
1; |
p ( x ) |
|
f ( x ) d x x |
|
|
|
|
|||||||||
x |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
p (1) |
|
1 |
|
|
|
1 |
|
• Из |
q |
1 |
( |
c1 |
) |
1 |
( c2 |
|
|||
|
|
||||||||||
|
|
|
|||||||||
|
|
|
|
|
|||||||
|
|
g (1) |
|
|
|||||||
|
|
|
|
|
g (1) |
|
|
0 |
|||
|
|
|
|
|
|
|
|
|
|
1
p ( x ) d x c1 g ( x )
0
1 d x ) 1 ; g ( x )
• |
Получаем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
0 ( |
1 |
1 |
1 |
) 1( c2 |
1 |
x |
d x 1 |
1 |
|
1 |
d x ) 0; c2 1 .5 |
|
|
|
|||||||
|
|
|
|
1 |
|
1 |
|
|
|
|||||||||||||
|
1 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
0 |
|
|
|
0 |
|
|
|
|
|
|
|
|
|
• |
Окончательно |
|
x |
p ( x ) |
|
|
|
|
|
x |
1 |
|
x |
2 |
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
u ( x ) c2 |
|
|
|
d x c1 |
|
d x 1 .5 |
|
x |
||||||||||||
|
|
|
g ( x ) |
|
g ( x ) |
2 |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
05.01.2011 |
11 |
Постановка задач для ДУ в частных производных (ДУЧП)
•Большинство окружающих нас физических полей описывается ДУЧП второго порядка, или системами таких ДУЧП.
•Если процессы не слишком интенсивны, то можно ограничиться линейными дифференциальными уравнениями второго порядка
aij |
2 u |
bi |
u |
cu f |
x1 x , |
x 2 y , |
x3 z , |
x 4 t . |
xi x j |
|
|||||||
ij |
i |
xi |
|
|
|
|
•Свойства решений существенно зависят от коэффициентов, стоящих при старших производных. С помощью соответствующей замены независимых переменных уравнение может быть приведено к одному из трех типов
Параболические
Гиперболические
Эллиптические
05.01.2011 |
12 |
Параболические
• Канонический вид |
|
u |
|
2 |
u |
|
2 |
u |
|
2 |
u |
f |
|
|
|
x 2 |
y 2 |
z 2 |
|||||||
|
t |
|
|
|
• |
Уравнение |
|
|
|
|
|
|
|
|
|
• |
теплопроводности |
|
u |
|
|
|
u |
|
|
|
a |
|
|
g |
|
|
|
||||
|
|
|
|
|
|
|||||
|
|
|
t |
|
x |
x |
|
y |
u
g y
|
|
u |
f |
|
|
g |
|
|
|
|
z |
|||
z |
|
|
•распространение возмущения
|
|
|
u |
|
u |
|
2 u |
|
t=0 |
|
|
|
||
t |
x 2 |
|
|
|
|
|
t=t1>0 |
||
|
|
|
|
|
|
|
|
|
t=t2>t1 |
|
|
|
0 |
|
|
|
|
0 |
z |
|
|
|
|
05.01.2011 |
13 |
|
|
|
|
|
Гиперболические |
|
|
|
|
|
|||||||||||||
• |
Канонический вид |
|
2 u |
|
|
2 u |
|
2 u |
|
2 u |
f |
|
|
|
|
|
|||||||
|
|
|
|
|
|
t 2 |
x 2 |
y 2 |
z 2 |
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
• |
Волновое |
|
2 u |
|
|
|
|
u |
|
|
|
|
|
u |
|
|
u |
|
|
||||
• |
уравнение |
|
a |
2 |
|
g |
|
|
y |
g |
|
g |
|
f |
|||||||||
• |
Д’Аламбера |
t |
|
|
x |
|
x |
|
|
|
|
y |
|
z |
z |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
• |
Распространение одномерного возмущения |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
u |
|
|
|
|
|
|
|
|
|
t=0 |
|
|
|
|
|
|
2 u |
|
2 u |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
t |
2 |
z |
2 |
|
|
|
|
|
t=t2 |
|
t=t1 |
|
|
|
|
t=t1 |
t=t2 |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
|
|
|
Точное решение |
|
u |
|
f ( t |
|
|
x ) |
|
|
|
|
|
|||||||
|
05.01.2011 |
|
|
|
|
|
|
|
|
|
|
|
|
|
14 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Эллиптические (стационарные процессы)
• Канонический вид |
2 u |
|
2 u |
|
2 u |
f |
|
|
|
|
|||
|
x 2 |
y 2 |
z 2 |
|||
|
|
|
|
•Уравнение
•Пуассона
•(Лапласса)
|
|
u |
|
|
u |
|
|
|
u |
f |
||||
|
g |
|
|
|
g |
|
|
|
g |
|
|
|||
|
|
|
|
|
|
|||||||||
x |
x |
|
y |
y |
|
|
z |
z |
|
•При постановке задач для ДУ в частных производных обычно
требуется найти распределение , удовлетворяющее ДУ в некоторой области пространства с границей Г.
y
Г
x
05.01.2011 |
15 |
Граничные условия
•Общее решение ДУЧП содержит произвольные дифференцируемые функции, например
• |
Решением |
|
2 |
u |
|
2 |
u |
является |
u f (t z ) |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
t |
2 |
z |
2 |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
• |
Начальные условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
u |
|
t 0 U |
0 |
( x , y , z ) |
|
u |
|
|
U 0 ( x , y , z ) |
|
||||||||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
t |
|
t 0 |
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
• |
Граничные условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
первого рода |
|
|
|
второго рода |
|
|
|
третьего рода |
|
|||||||||||||||||
|
Дирихле |
|
|
|
|
|
Неймана |
|
|
|
|
Ньютона |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
u |
|
Г ( Г ) |
|
|
|
|
|
u |
|
|
( Г ) |
|
|
u |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u |
|
|
( Г ) |
|||||||||
|
|
|
|
|
|
|
n |
|
|
|
n |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
Г |
|
|
|
|
|
|
|
|
Г |
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
05.01.2011 |
16 |
Пример решения |
|
|
|
u |
|
2 u |
1 |
|||||
|
|
|
|
|
||||||||
u |
|
|
|
|
|
|
|
|
t |
z 2 |
||
|
|
|
|
|
|
|
|
|
|
|||
|
1; |
u |
|
x 1 0; |
u |
|
t 0 0 |
|
|
|
||
|
|
|
|
|
|
|
||||||
x |
|
|
|
|||||||||
|
x 0 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
05.01.2011 |
17 |
Как получают дифференциальные уравнения
•ДУ являются следствием фундаментальных законов природы и применения интегральных теорем
•Например, Французский математик и физик Жан Батист Жозеф Фурье в 1822 г. установил закон теплопроводности (закон Фурье):
количество тепла, проходящее через единицу площади в единицу времени, пропорционально проекции градиента температуры на нормаль к поверхности. Математически этот закон выражается следующим простым соотношением
|
|
|
|
|
|
|
|
q T n |
T |
||
|
|
|
|
|
|
|
|
|
|||
q |
|
|
|
|
|
|
|
|
|
T гр а д |
dS |
В т |
|
2 |
|
В т |
|
n |
|||||
м |
|
|
|
|
|||||||
|
|
|
|
|
|
м г р а д |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
05.01.2011 |
18 |
Пример получения стационарного двумерного уравнения теплопроводности
используя закон Фурье
Выберем элементарный |
|
y/2 |
|
объем V= x y |
|
|
|
в центре которого находится |
|
|
|
точка x,y . |
- x/2 |
|
x/2 |
|
|
||
|
|
|
|
|
|
- y/2 |
Допустим, что в таком кубике в |
|
|
единицу времени выделяется |
Q V |
qV ( x , y ) x y |
количество теплоты, равное |
05.01.2011 |
19 |
Составим уравнение баланса для элементарного объема
|
|
|
|
T |
|
qV |
x y |
|
|||
|
|
||||
x |
|||||
|
|
|
|
||
|
|
|
|
|
x x
2
T
x
x x
2
|
|
|
|
|
y |
|
|
||
|
||||
|
|
|
||
|
|
|||
|
|
|
|
T |
|
||
|
|
|
|
y |
y |
||
y |
|||
|
|
||
|
|
2 |
|
|
T |
|
||
|
|
|
|
|
|
y |
y |
||||
|
|
y |
|||
|
|
|
|
||
|
|
|
|
2 |
x
Разделим на x y:
|
|
|
|
T |
|
|
|
T |
|
|
|
|
|
|
T |
|
|
|
|
T |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
y |
|
y |
y |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
x x |
x |
|
|
x x |
x |
|
|
|
y |
y |
|
|
|
y |
||||||
qV |
|
|
|
|
|
2 |
|
|
|
|
2 |
|
|
|
|
|
|
2 |
|
|
|
|
|
2 |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Перейдем к пределу 0
T
xx
|
|
|
|
|
T |
qV ( x , y ) |
|
|
|
|
|
|
|
|
|||
|
y |
|
||||||
|
|
y |
|
|
|
|
05.01.2011 |
20 |