 IV.

ЧИСЛЕННЫЕ МЕТОДЫ РЕШЕНИЯ КРАЕВЫХ ЗАДАЧ ДЛЯ

ОБЫКНОВЕННЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
ВТОРОГО ПОРЯДКА

Численные методы решения краевых задач для обыкновенных

дифференциальных уравнений рассматриваются на примере следующей
задачи для линейного уравнения с малым параметром при старшей
производной:
  u ( x )  a ( x )  u ( x )  b( x )  u ( x )  f ( x ), x  (0,1);

 0  u (0)  0    u (0)   0 , (IV.1)
  u (1)      u (1)   .
 1 1 1

В (IV.1):   0 - малый параметр; a(x), b(x), f(x) - непрерывные функции

на отрезке [0,1], причем b(x)  0 для x  [0,1] ; для параметров, входящих в
краевые условия предположим выполненными неравенства:

 0  0, 1  0, 0  0, 1  0;  0  0  0, 1  1  0.

Если к выше перечисленным условиям добавить еще одно:

1

 b(x )dx   0  1  0,
0

то можно гарантировать существование и единственность классического

решения задачи (IV.1) (см. [14]).
Задачи вида (IV.1) являются простейшими математическими
моделями диффузионно-конвективных процессов и родственных
физических явлений [15]. Наличие относительно малых подобластей с
большими градиентами решения, которые возникают при малых
значениях параметра  , делают такие задачи сложными для численной
реализации и требует использования разностных схем, учитывающих их
специфику [16].
Л А Б О Р А Т О Р Н А Я Р А Б О Т А IV.1

МОНОТОННЫЕ РАЗНОСТНЫЕ СХЕМЫ ДЛЯ ЛИНЕЙНЫХ

УРАВНЕНИЙ

Основная цель. Познакомиться с монотонными разностными

схемами для решения простейших краевых задач с краевыми условиями
общего вида. В том числе, рассмотреть задачи, решения которых
формируют пограничные слои.

Теория и основные формулы. На отрезке [0,1] рассмотрим

сетку
0  x1  x 2  ...  x n  1
с равноотстоящими друг от друга узлами:
i 1
xi  , i = 1,2,...,n;
n 1
1
и шагом h  . Примем также следующие обозначения:
n 1
ai  h
a i  a (x i ) , b i  b(x i ) , f i  f ( x i ) , R i  , (i=1,2,…,n);
2

u ih1  u ih u ih  u ih1
D u  h
i , i  1,2,...,n  1; D u i 
h
, i  2,3,...,n .
h h

Для аппроксимации задачи (IV.1) рассматривается следующее семейство

разностных схем, полученное с использованием простейшего метода
дискретизации – метода конечных разностей:

D  u ih  D  u ih D u h  D  u ih
  i   ai   i  b i  u ih  f i , i  2,3,...n  1; (IV.1.1)
h 2

 0 u1h  0  D  u1h  0 ,
(IV.1.2)
1u hn  1  D  u hn  1 .
В (IV.1.1)  i (i  1,2,...,n) - параметры аппроксимационной вязкости, их
выбор и определяет разностную схему в семействе (IV.1.1); неизвестными
в (IV.1.1), (IV.1.2) являются величины u ih ,i  1,2,...,n, аппроксимирующие
значения u(xi), (i=1,2,…,n) точного решения задачи (IV.1). Разностные
краевые условия (IV.1.2) могут быть заменены на более точные:
1  1 1  1
(  0  b1  h  )  u 1h  0   ( 1  R 1 )  D  u1h   0  0  f1  h  ,
2 2 (IV.1.3)
1  n 1  n
( 1  b n  h  )  u n  1  (  n  R n )  D  u n  1  1  f n  h 
h h
.
2 2
Здесь:
i 1
i  , i = 1,2,...,n.
Ri

Для выбора параметров i (i  1,2,...,n) можно предложить следующие

известные из литературы варианты :
1)  i  1 , (i = 1,2,...,n) - схема с центральной разностью [15];

2)  i  1 R i , (i = 1,2,...,n) - схема с направленной разностью [15];

R i2
3)  i  1  , (i = 1,2,...,n) - схема Самарского А.А. [17];
1 Ri

R 3i
4)  i  1  , (i = 1,2,...,n) - схема Булеева Н.И. и
1  R i  R i2
ТимухинаТ.И. [18];
1 2 Ri
5)  i  1  R i2  , (i = 1,2,...,n)-схема Булеева Н.И. [19];
3  3  R i  2  R i2

6)  i  R i  cthRi , (i = 1,2,...,n) - схема Ильина А.М. [16,20].

Следующая схема (El-Mistikawy & Werle [21]) не содержится в семействе

(IV.1.1), для ее записи примем обозначения:

a ( x i )  a ( x i1 ) b( x i )  b( x i1 ) f ( x i )  f ( x i1 )

a i1 2  , b i1 2  , f i1 2  ;
2 2 2
 a i1 2  h 1
R i1 2  ,  i1 2  cthR i1 2  .
2 R i1 2

В этих обозначениях разностные уравнения имеют вид:

 
 (1  R    R)i  D  u ih   (1  R    R)i D  u ih  1 1
h h 2 2

1  i  1  i
(bi   bi   )  u ih 
1 1

1
2
1
2
(IV.1.4)
2
2 2 2

1  i  1  i 
 fi   fi  , i  2,3,...,n  1;
1 1
2 2
1 1
2
2 2 2

1 
( 0  b  h  )  u1h  0    (1  R    R)  D  u1h 
1
2
1 1
2
2 2

1 
 0  0  f  h 
1

1
2
;
2
2
(IV.1.5)
1  n
(1  b n   h  )  u nh  1    (1  R    R) n   D  u nh 
1
2
1 1
2
2 2

1  n
 1  1  f n  h
1

1 . 2

2
2

Разностные задачи (IV.1.1), (IV.1.2); (IV.1.1), (IV.1.3) и (IV.1.4), (IV.1.5)

могут быть записаны в виде:

B1  u1h  C1  u 2h  F1 ,
Ai  u ih1  Bi  u ih  Ci  u ih1  Fi , i  2,3,...,n  1; (IV.1.6)
A n  u nh 1  Bn  u nh  Fn .

и решены при помощи метода прогонки [17] . «Внутренние» свойства

дискретной задачи (IV.1.6) определяются понятием монотонная
разностная схема.
О п р е д е л е н и е IV.1.1 (см. [17]). Дискретная задача (IV.1.6)
называется монотонной, если для нее выполнены условия:
1) B1  0 , C1  0 , Bn  0 , A n  0 ;
2) A i  0 , C i  0 , i =2,3,...,n-1;
 3) B1  C1 , Bi  A i  C i , (i = 2,3,...,n-1), Bn  A n .
и хотя бы одно из этих неравенств - строгое.
Разностные схемы (IV.1.1), (IV.1.2); (IV.1.1), (IV.1.3) и (IV.1.4),
(IV.1.5) называются монотонными, если соответствующие им разностные
задачи (IV.1.6) монотонны для любых положительных значений
параметров h и  
Известно, что монотонность разностной схемы гарантирует
однозначную разрешимость соответствующей задачи (IV.1.6) (это
доказано, например, в [17]), кроме того, для такой разностной схемы
обеспечено выполнение достаточных условий корректности и
устойчивости монотонной прогонки.
Предположим, что u  ( x ) - решение дифференциальной задачи
(IV.1), в обозначения включена зависимость от параметра   (0,  0 ] , пусть

u h  {u h,i }in1 - сеточная функция являющаяся решением некоторой

разностной задачи, аппроксимирующей задачу (IV.1), например задачи
(IV.1.1), (IV.1.2); (IV.1.1), (IV.1.3) или (IV.1.4), (IV.1.5). Для произвольной
сеточной функции v h  {v ih }in1 введем «сильную» сеточную норму по

формуле:
vh  max v ih ,
 1i n

напомним также, что (u ) h - проекция непрерывной функции u(x) на сетку.

О п р е д е л е н и е IV.1.2 (см. [16]). Решение u h разностной

(дискретной) задачи сходится к решению u  ( x ) дифференциальной задачи
с порядком p > 0, если найдется такая константа C > 0, возможно
зависящая от  , но не зависящая от h («классическая» сходимость), что:
(u  ) h  u h  С  h p .

Если константа «C» в последней оценке не зависит от h и от  , то говорят

о "равномерной по  " сходимости с порядком «p» 
 Из литературы известно (например, [16]), что разностные схемы 1) -
4) не обладают равномерной по  сходимостью. Тем не менее, порядки
классической сходимости для этих схем следующие: для схемы 2) –
первый, для схем 1), 3), 4) - второй.
Схемы 5) (Ильина А.М.) и 6) (El_Mistikawy & Werle) обе обладают
вторым порядком классической сходимости, кроме того, первая из этих
схем гарантирует первый порядок равномерной по  сходимости, а вторая
схема - второй порядок равномерной по  сходимости.
Все схемы, за исключением схемы с центральной разностью 1),
являются монотонными в смысле Определения IV.1.1, схема 1)
гарантирует монотонность только в случае выполнения ограничения на
выбор параметров h и  :
h  (a ) h  2 ,


поэтому эту схему обычно классифицируют как условно-монотонную.

Тестовые задачи.
Задача 1. «Уравнение с постоянными коэффициентами», считаем,
что коэффициенты уравнения в (IV.1) постоянные:
a (x)  a, b(x)  b  0, f (x)  f ,
тогда решение задачи (IV.1) можно выписать в явном виде, сделаем это,
вводя обозначение для следующей функции:
e-x , x  50,
E(x)=  (IV.1.7)
0, x  50.
и рассматривая три различных случая:
a) b  0;

a  a 2  4b a  a 2  4b
s1  , s2  ,
2 2
t 01  0  0s1, t 02   0  0s2 ,
t11  1  1s1 , t12  1  1s 2 ,
 d  t 01  t12  E(s 2 )  t 02  t11  E(s1 ),
t 12  E(s1  x  s 2 )  t 11  E(s 2  x  s1 )
z 0 (x)  ,
d
t  E(s 2  x )  t 02  E(s1  x )
z1 ( x )  01 ,
d
Точное решение:
 f  0   f 1  f
u ( x)  0    z0 ( x)  1    z1 ( x)  .
 b   b  b

b) b  0, a  0;
Точное решение:
   f  0   f  (1  1 )  f x
u (x )   0    z 0 ( x )   1    z1 ( x )  .
 a   a  a
c) b  0, a  0;
d   0  1    ( 0  1  1  0 ),
1  (1  x )    1   x    0
z 0 (x)  , z1 ( x )  0 ,
d d
Точное решение:
 f  (1  21 )  f  x2
u ( x )  0  z 0 ( x )   1  
 1 z ( x )  .
 2   2 

Задача 2.

Положим в (IV.1):

2 2 1 2  x
a ( x)   , b( x)  0, f ( x)   [  (1  x)  1]   ,
1  x (1  x) 2
(1  x) 2

1

 0   1  1, 0  1  0, 0  0.5  e  , 1  0.

Точное решение:
1 2
 
 (1 x )
x  (1  x) e 
e
u ( x)   1
.
2 
2(1  e ) 
 Задача 3.

Положим в (IV.1):

2 3  3  (1  x )
a ( x )  3  (1  x ) 2 , b( x )  0, f ( x )   ,
1 x 2  (1  x ) 2
2
1  1 ln 2
 0  1  1, 0  , 1  0, 0   7
, 1  1  .
3 6 1 e  2

Точное решение :
1(1 x )3

1 e 
ln(1  x )
u(x)  7
 .
1 e  2

Задача 4.

Положим в (IV.1):

2 4 1  2x e  
1
2
a(x)  , b( x )  , f (x)    1
 ,
(1  x ) 2
(1  x ) 3
 1  x 1 e    (1  x )
3

 1
 1  e 
 0  2, 1  2, 0  1, 1  4, 0  1    1
, 1  1   
 1
.
 1 e   1 e 

Точное решение:
1 2x
 
(1 x )
x e 
e
u(x)   .
1 x 
1

2(1  e ) 

Задача 5.

Положим в (IV.1):

2 2 4
a (x)   , b( x )  0, f ( x )  ,
1  x (1  x ) 2
(1  x ) 4
 0  1  1, 0  1  0, 0  1, 1  1.

Точное решение:
 1 2x
 
(1 x )
2x e 
e
u(x)   .
1 x 
1

(1  e ) 

Задача 6.

Положим в (IV.1):

a (x)  0, b(x)  (1  x) 2 ,
 0  1  1, 0  1  0, 0  1  0.

Точное решение:

u(x) 

Ex  
  E (1  x )    [cos(x)] ,
 
2

1 E 1 

где функция E(x) определена формулой (IV.1.7), а правая часть f(x) в

уравнении (IV.1) определяется после подстановки решения u(x) в левую
часть этого уравнения.

Требования к программе. Программа должна включать:

1) Разностную схему (IV.1.1), (IV.1.2), или (IV.1.1), (IV.1.3), с
возможностью выбора параметров i (i  2,3...,n-1) , или при их
фиксированных значениях, или схему (IV.1.4), (IV.1.5) (по
усмотрению преподавателя).
2) Тестовые задачи (по выбору преподавателя).
3) Вывод относительных погрешностей по формуле:
max u ( x i )  u ih
Err  1i  n
 100(%) .
max u ( x i )
1i  n

4) Возможность выбора числа узлов сетки «n», а также параметра «  »

(если он есть по условию тестовой задачи).
 5) Графика: одновременная отрисовка точного и приближенных
решений (проинтерполированных кусочно-линейным сплайном).

Задание для работы с программой. В процессе работы с программой

и при написании Отчета необходимо ответить на следующие вопросы:
1. Как ведут себя решения дифференциальных задач при 
стремящемся к нулю?
2. Как влияет увеличение параметра n (при фиксированном  ) на
точность разностных схем? Сопоставить численные результаты с
теоретическими оценками сходимости.
3. Влияет ли уменьшение параметра  (при фиксированном n) на
точность разностных схем (по каждой схеме отдельно)?
4. Сравнить точность разностных схем между собой и сопоставить
экспериментально полученные результаты с теорией.
6. Исследовать монотонность используемых разностных схем.
Проявляется ли наличие или отсутствие свойства монотонности в работе
разностной схемы?
7. Задача 1 предполагает возможность выбора постоянных
коэффициентов в уравнении и граничных условиях. При работе с этой
задачей необходимо рассмотреть различные варианты пограничных слоев:
а) пограничный слой у левой границы (a > 0, b  0),
б) пограничный слой у правой границы (a < 0, b  0),
в) пограничные слои у обеих границ (a = 0, b > 0).