Презентации и конспект лекции Синицын / Презентации_ВЭиМОвТС / ВычЭксперим / Установочная лекция
.pdf
Граничные условия
•Общее решение ДУЧП содержит произвольные дифференцируемые функции, например
• |
Решением |
|
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 |
11 |
Суть метода сеток
•Суть метода сеток в том, что решение ДУ получают в виде достаточно подробной таблицы значений искомого решения в узлах сетки, покрывающей область определения решения.
•Получаемая таблица должна обладать свойством аппроксимации, т.е. возможностью восстановления всех значений искомого точного решения с заданной погрешностью.
•Будем иллюстрировать реализацию метода сеток на решении простейшей одномерной краевой задачи Дирихле
|
|
|
u |
|
f x , u ; |
|
|
u (0 ) ; |
u (b ) . |
|
g ( x , u ) |
|
|
|
|||||
|
|
|
|||||||
|
x |
x |
|
|
|
|
|
||
• В общем случае |
L u f ; |
u |
|
Г ( Г ), |
|
||||
|
|
||||||||
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
||
•Г - граница многомерной области , внутри которой необходимо получить
решение. В рассматриваемом частном случае представляет собой отрезок
[ 0 , b ]
05.01.2011 |
12 |
Результат решения по методу сеток
u |
Искомое |
|
решение |
|
|
|
|
u |
Решение в виде |
|
таблицы |

ui

0 |
b |
x |
0 |
xi |
b |
x |
|
|
|
|
05.01.2011 |
13 |
Получение конечноразностной схемы
• Решение u(x) ищется в виде таблицы значений в узлах
выбранной сетки
u h u i u ( xi )
• дифференциальное уравнение L u f , заменяется системой алгебраических уравнений, связывающих между собой значения искомой функции в соседних узлах.
•Такая система алгебраических уравнений называется конечноразностной схемой
•Обозначим
, u h u1 , u 2 , ..., u n 1f hLh u h
• Имеется много способов получения конечно-разностной схемы
05.01.2011 |
14 |
Конечно-разностная схема
|
|
|
|
|
|
|
|
u |
|
|
f x , u ; |
|
u (0 ) ; |
|
u (b ) . |
|||||||||||||
|
g ( x , u ) |
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
x |
|
|
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
g |
i 1 / 2 |
|
|
|
|
|
g |
|
1 / 2 |
g |
i 1 / 2 |
|
|
|
|
g |
i 1 / 2 |
|
|
|
|
|
|
||||
|
|
|
|
|
i |
|
|
|
|
|
f i , i 2 ... n. |
|||||||||||||||||
|
|
u i 1 |
|
|
|
u i |
|
u i 1 |
||||||||||||||||||||
|
|
h 2 |
|
|
|
|
|
h 2 |
|
|
h 2 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
ai |
|
|
|
|
|
|
|
|
|
|
|
i 2 ... n. |
|
ui 1 bi ui ci ui 1 f i , |
|||||||||||||
|
|
|
|
|
|
||||||||
u1 , u n 1 , |
|
||||||||||||
05.01.2011 |
15 |
система конечно-разностных уравнений
• Стандартная система с трехдиагональной матрицей:
b1 |
c1 |
0 |
0 |
... |
0 |
0 |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
a 2 |
b2 |
c 2 |
0 |
... |
0 |
0 |
|||||
|
0 |
a 3 |
b3 |
c3 |
... |
0 |
0 |
||||
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||
... ... ... ... ... ... |
... |
||||||||||
0 |
0 |
0 |
0 |
... |
a |
n |
b |
n |
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
0 |
0 |
0 |
0 |
... |
0 |
a |
n 1 |
||||
|
|
|
|
|
|
|
|
|
|||
0
0
0
...
c n bn 1
|
|
|
|
|
||||||
u1 |
||||||||||
|
|
|
|
|
|
|||||
u |
2 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u |
3 |
|||||||||
|
|
|||||||||
|
|
|
|
|
|
|
|
|||
|
... |
|||||||||
|
|
|
|
|
|
|
||||
u |
|
n |
||||||||
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|||
|
|
u |
n 1 |
|||||||
|
|
|
|
|
|
|||||
|
d 1 |
|
d 2
|
|
|
d 3 |
|
|
|
|
|
|||
|
|
... |
|
||
|
|
d |
n |
|
|
|
|
|
|
|
|
|
|
d |
n 1 |
|
|
|
|
|
|
|
|
n1 n 1; b1 1; c1 0; d 1 ; a n 1 0; bn 1 1; d n 1 ;
a i g i 1 / 2 |
; |
ci |
g i 1 / 2 ; |
|
|
a |
|
c |
|
|
|
h |
2 |
|
; i 2 ... n. |
b |
|
|
; |
d |
|
f |
|||||||||
|
|
i |
i |
i |
|
||||||||||
05.01.2011 |
|
|
|
|
|
i |
|
|
|
|
i |
16 |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Погрешность аппроксимации
•При замене дифференциального уравнения системой алгебраических уравнений вносится так называемая погрешность аппроксимации
конечно-разностной схемой дифференциального уравнения
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
Lh |
|
|
|
|
вместо значения u h |
|||
• |
подставим в конечно-разностную схему |
|
u h f h , |
||||||||||||||||
|
точного решения u h . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
• |
Ввиду того, что |
|
|
|
, после такой подстановки получается невязка |
||||||||||||||
u h u h |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
• |
Погрешность аппроксимации |
|
h |
Lh u h |
f h 0 |
||||||||||||||
• |
Основное требование: |
при |
h 0 |
|
h 0 |
|
|
|
|||||||||||
• |
Ассимптотическая оценка |
|
|
|
C h p |
|
|
|
|||||||||||
|
h |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
•Порядок погрешности аппроксимации = p
05.01.2011 |
17 |
Оценка погрешности решения Понятие устойчивости
• |
Погрешность решения : |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
h |
u h |
u h |
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
• |
Для сходимости к точному решению |
|
h 0 п р и |
h 0 |
|
|
|||||||||||||||||||||
• |
кроме |
h 0 |
|
необходима |
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
• устойчивость к ошибкам округления |
|
|
|
|
|
|
|
|
|
C 0 |
|
f h f h |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u h u h |
|
|
|
|
|||||||||
• |
|
|
|
|
|
Основная теорема |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
L h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u h |
f h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
L h u h f h h f h |
|
|
|
h |
|
C 0 |
|
h |
|
C 0 C h |
p |
|||||||||||||||
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
h |
|
|
|
|
|
|
|
|
||||||||||||
|
L h ( |
u h u h ) L h h |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
05.01.2011 |
18 |
Проекционные методы решения краевых задач
L u f ; |
u |
( Г ) |
- Краевая задача |
|
|
Г |
|
|
|
|
|
1 ( x ), 2 ( x ), ..., N ( x ) |
- базис из функций |
||
N |
|
|
u N ( x ) a k k ( x ). |
- Представление решения |
|
k 1 |
|
|
N |
|
|
a k L k , i f , i ; |
i 1 ..N |
- Проекционное уравнение |
k 1
05.01.2011 |
19 |
Решение одномерной краевой задачи методом Галеркина
• |
Найти решение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
u |
|
f ; |
|
|
|
|
q 0 |
u |
|
0 u |
|
|
|
|
0 ; |
|
q 1 |
u |
1u |
|
1 . |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
g |
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
x |
x |
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
x 0 |
|
|
|
|
x |
|
|
x 1 |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
• Ищем решение в виде |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
u N ( x ) 0 ( x ) a k k ( x ) |
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k 1 |
|
|
|
|
|
|
|
|
|||
• |
Проекционное уравнение |
1 |
|
|
|
u |
N |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
g |
|
|
|
i d x |
|
f i d x ; |
i 1 ... N . |
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
0 |
x |
|
x |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
u N |
|
|
|
1 |
|
1 |
|
|
u N |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
• |
преобразуем |
|
g |
i |
|
|
g |
|
|
|
|
i |
|
d x |
|
|
f i d x |
|
|
|
|||||||||||||||||
|
x |
|
|
|
|
x |
|
|
|
x |
|
|
|
||||||||||||||||||||||||
05.01.2011 |
|
|
|
|
|
|
|
|
0 |
|
0 |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
20 |
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
