Скачиваний:
43
Добавлен:
15.06.2014
Размер:
440.62 Кб
Скачать

Граничные условия

Общее решение ДУЧП содержит произвольные дифференцируемые функции, например

Решением

 

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