Convergence of the matrix method of numerical integration of the boundary value problems for linear nonhomogeneous ordinary differential second order equations with variable coefficients

Abstract


The problems of stability and convergence of previously proposed matrix method of numerical integration of boundary value problems with boundary conditions of the first, second and third kinds of nonhomogeneous linear ordinary differential second order equations with variable coeffcients are considered. Using of the Taylor polynomials of arbitrary degrees allowed to increase the approximation order of the method to an arbitrary natural number and to refuse from the approximation of derivatives by finite differences. When choosing the second degree Taylor polynomials the equation of the method coincided with the known equations of the traditional method of numerical integration of the boundary value problems where the derivatives are approximated by finite differences. It was shown that a suffcient criterion of stability when used in the method of Taylor polynomials of the third degree and more coincides with the suffcient criterion of stability of the traditional grid method for the numerical integration of boundary value problems with boundary conditions of the first, second and third kind. Theoretically, it is established that the degree of convergence of the matrix method for integration of boundary value problems with boundary conditions of the first kind is proportional to the degree of the used Taylor polynomials in the case, when the degree is even, and is proportional to the number that is one less than the degree if it is odd; when integrating the boundary value problems with boundary conditions of the second and third kind the degree of convergence of the method is proportional to the degree of the used Taylor polynomials regardless of its parity and one less than it. The obtained theoretical results are confirmed by numerical experiments.

Full Text

В классическом методе численного интегрирования краевых задач для неоднородных обыкновенных дифференциальных уравнений второго порядка (ОДУ2) с переменными коэффициентами u + p(x)u + q(x)u = f (x), u0 = u0 , un = un , (1) где p(x), q(x), f (x) - заданные функции, дифференцируемые нужное число раз, [a, b] - область интегрирования, u0 , un - заданные числа, конечные разности при аппроксимации производных дают второй порядок точности [1-6]. Такой же порядок аппроксимации имеют методы численного интегрирования ряда краевых задач для уравнений в частных производных [4-10]. Второй порядок обусловлен тем, что при аппроксимации производных конечными разностями удерживается всего три члена разложения в ряд Тейлора искомого решения задачи. Метод, использующий средства матричного исчисления и численного интегрирования разностных краевых задач для неоднородных ОДУ2, позволяющий увеличить количество удержанных членов до произвольного натурального числа в разложении искомой функции в ряд Тейлора, предложен в [11], при этом аппроксимация производных конечными разностями не использовалась. Оценка порядка аппроксимации предложенного в [11] матричного метода интегрирования для различных значений удержанного числа членов в рядах Тейлора дана в [12] для краевых задач с граничными условиями первого, второго и третьего рода. Поставим целью исследовать вопросы устойчивости и, как следствие, сходимости матричного метода численного интегрирования разностных краевых задач для ОДУ2 с различными граничными условиями. 1. Обозначения и некоторые предшествующие результаты. Остановимся кратко на полученных в [11] результатах, придерживаясь следующих обозначений, принятых в [4]: 1) D - область интегрирования, ограниченная отрезком [a, b], Dh - узлы сетки, определяемые значениями xi = x0 + ih, i = 1, 2, . . . , n, x0 = a, xn = b, h = (b - a)/n, n + 1 - число узлов сетки; 2) u(x) - непрерывная функция, являющаяся точным решением краевой задачи (1); 3) [u]h - сеточная функция, совпадающая с точным решением задачи (1) в узлах сетки Dh ; 4) u(h) - искомая сеточная функция. Для краткости примем для любой функции обозначение ϕ(xi ) = ϕi , где xi - узел сетки Dh . В дальнейшем опустим индекс h в наименованиях сеточных функций [u]h , (h) и будем оговаривать особо случаи, в которых будет использоваться непреu рывная функция u(x), являющаяся точным решением задачи (1). Для некоторого фиксированного k 2, k - степень используемого многочлена Тейлора, в узлах сетки с номерами i = 1, 2, . . . , n-1 составлена система уравнений, состоящая из (k + 1)-го уравнения [11]. В систему внесены первые k +1 члена разложения в ряд Тейлора точной функции u(x) в узлах ui-1 , ui+1 соответственно и k - 1 уравнение, которые получены дифференцированием 560 Сходимость матричного метода численного интегрирования краевых задач . . . по x обеих частей уравнения (1), так что qi ui + pi ui + ui (r) (r) = fi , где r = 0, 1, . . . , k - 2. Получено  3 k 2  u - hu + h u - h u + · · · + (-1)k h u(k) = u ,  i  i-1 i  2! i 3! i k! i    h2 h3 hk (k)  ui + hu + u + u + · · · + u = ui+1 ,   i 2! i 3! i k! i qi ui + pi ui + ui = fi ,  q ui + (qi + p ) u + pi u + u = f ,  i i i i i i   . . .     (k-2) (k) (k-2) q ui + · · · + ui = fi . i (2) В матричной форме система уравнений (2) имеет вид Ak Vik = Gk в обоi i значениях  3 k h2 - h . . . (-1)k h 1 -h 2! 3! k!  1  h2 h3 hk h ...   2! 3! k!  q  pi 1 0 ... 0 k  i , Ai =  (3)  qi qi + p i p i 1 ... 0    ... ... ... ... ... ...  (k-2) qi ... ... ... ... 1     ui ui-1  ui   ui+1  u   f   i   i  k k V i =  u  , Gi =  f  .  i   i   ...   ...  (k) (k-2) ui fi В дальнейшем верхний индекс будет означает степень многочлена Тейлора k, если речь не идёт о показателях алгебраических степеней, степенях производных, символов обратных матриц и транспонирования. При фиксированном номере i узла сетки Dh все элементы матрицы Ak могут быть вычисi лены, следовательно, (2) является системой линейных алгебраических уравнений (СЛАУ). В предположении существования обратной матрицы (Ak )-1 i от матрицы Ak найдено i Vik = (Ak )-1 Gk , i i или в координатной форме (k-2) , (4) (k-2) , (5) (k-2) , (6) ui = aki ui-1 + aki ui+1 + aki fi + aki fi + · · · + aki fi 11 12 13 14 1,k+1 ui = aki ui-1 + aki ui+1 + aki fi + aki fi + · · · + aki fi 21 22 23 24 2,k+1 ui = aki ui-1 + aki ui+1 + aki fi + aki fi + · · · + aki fi 31 32 33 34 3,k+1 561 М а к л а к о в В. Н. ... (k) ui = aki ui-1 + aki ui+1 + aki fi + k+1,1 k+1,2 k+1,3 (k-2) +aki fi + · · · + aki k+1,4 k+1,k+1 fi , (7) где aki - соответствующие элементы обратной матрицы (Ak )-1 в узле с ноi lm мером i. Из уравнений (4), каждое из которых при i = 1, 2, . . . , n-1, есть разностное уравнение второго порядка, связывающее значения искомой сеточной функции ui в узлах xi-1 , xi , xi+1 , при построении которого не использовались конечные разности для аппроксимации производных, составлена СЛАУ, решение которой для граничных условий u0 = u0 , un = un даёт искомое приближённое решение ui , i = 1, 2, . . . , n - 1, задачи (1). В [11, 12] показано, что увеличение числа удержанных членов в разложении в ряд Тейлора функции u(x), являющейся точным решением задачи (1), приводит к уменьшению погрешности между точным и найденным приближённым решениями задачи в узлах сетки Dh . После нахождения решений ui , i = 1, 2, . . . , n - 1, производные искомой сеточной функции вплоть до k-того порядка включительно могут быть найдены по формулам (5)-(7) во внутренних узлах сетки Dh . Значения производных в граничных узлах могут быть вычислены, например, в x0 , как hk-1 (k) h2 u . u0 = u1 - hu1 + u1 + · · · + (-1)k-1 2! (k - 1)! 1 2. Некоторые предварительные оценки. В дальнейшем квадратные матрицы Ak будем называть локальными матрицами, размерность которых равна i (k + 1). Уравнение (4) преобразуем во всех внутренних узлах сетки Dh как - k+1 ki a1j (j-3) aki ui aki 11 ui-1 + ki - 12 ui+1 = fi + f , ki ki a13 a13 a13 aki i 13 j=4 (8) и, введя соответствующие обозначения в узле xi :  ki Mk  k  a = - a11 = - (11),i ,  i  k  aki M(13),i  13   k   k b = 1 = ∆i ,  i  k  aki M(13),i 13 k ki  k c = - a12 = - M(12),i ,  i  k   aki M(13),i  13    Mk  aki  k d = 1j = (1j),i , j = 4, 5, . . . , k + 1,  ij  aki Mk 13 (9) (13),i k где ∆k - определитель матрицы Ak ; M(1j),i - алгебраическое дополнение элеi i мента, расположенного на пересечении первой строки со столбцом с номером 562 Сходимость матричного метода численного интегрирования краевых задач . . . j, j = 1, 2, . . . , k + 1, матрицы (Ak ) , - символ транспонирования, разностi ную краевую задачу, аппроксимирующую первую дифференциальную краевую задачу (1), запишем как  ak ui-1 + bk ui + ck ui+1 = fi +  i i i u0 = u0 ,  u =u , n n k+1 k (j-3) , j=4 dij fi i = 1, 2, . . . , n - 1, (10) или, по аналогии с [4], в компактной символьной форме k Lk u = fh , h причём наряду с введённым обозначением эту же задачу для краткости будем обозначать и как Lk . Дифференциальную краевую задачу, по аналогии с [4], h будем обозначить как Lu = f . При k = 2 из (8), (9) следует разностное уравнение 2 - h2 qi 2 + hpi 2 - hpi ui-1 - ui + ui+1 = fi , 2 2 2h h 2h2 (11) которое совпадает с уравнением традиционного метода численного интегрирования [4], использующего конечные разности при аппроксимации производных задачи (1). Исходя их определения матрицы Ak в форме (3) для любого натурального i числа k 3 имеем 1 ∆k = i -h 1 qi ... h pi ... ri ti (k-3) qi (k-2) qi h2 2! h2 2! ... 1 ... ... ... si . . . vi . . . = + ... hk k! (-1)k+1 hk k! (-1)k-1 hk-1 (k-1)! hk-1 (k-1)! (-1)k hk k! hk k! 0 ... 1 wi 0 ... 0 1 1 qi ... (k-3) qi (k-2) qi 1 qi ... (k-3) qi (k-2) qi 2 h h 2! pi 1 ... ... ri si t i vi ... ... ... ... ... 2 -h h 2! pi 1 ... ... ri si ti vi ... ... ... ... ... = hk-1 (k-1)! 0 ... 1 wi + (-1)k-1 hk-1 (k-1)! 0 ... 1 wi + 563 М а к л а к о в В. Н. 1 + -h 1 qi ... h pi ... ri (k-3) qi = hk k! k-1 am hm + m=0 h2 2! h2 2! ... ... 1 ... ... ... si . . . (-1)k+1 hk k! (-1)k-1 hk-1 (k-1)! hk-1 (k-1)! 0 ... 1 = k-1 bm hm + ∆k-1 = O(hk ) + ∆k-1 , i i m=0 где ri , si , ti , vi , wi - некоторые записанные в узле xi функции от p(x), q(x) и их производных; am , bm - независящие от h величины; или k-1 ∆k = ∆i + O(hk ). i (12) Уменьшая степень используемого многочлена Тейлора в равенстве (12) на единицу, получим k-2 ∆k-1 = ∆i + O(hk-1 ). (13) i Подставляя в (12) равенство (13), получим, пренебрегая старшими степенями, ∆k = ∆k-2 + O(hk-1 ). Выполнение аналогичной процедуры ещё k - 4 i i раза приводит к выражению ∆k = ∆2 + O(h3 ) i i (14) в узле xi , i = 1, 2, . . . , n - 1, сетки Dh . Для первой краевой задачи (1) в [12] получены следующие оценки в узле с номером i: k-1 k M(1j),i = M(1j),i + O(hk ), j = 1, 2. (15) Выполняя описанную выше процедуру для равенства (15), найдём k 2 M(1j),i = M(1j),i + O(h3 ), j = 1, 2. (16) В узлах сетки Dh с номерами i = 2, 3, . . . , n - 2 второй и третьей краевых задач, граничные условия которых имеют, соответственно, вид u0 = u0 , un = un и α0 u0 + β0 u0 = z0 , α1 un + β1 un = zn , где u0 , un , α0 , β0 , z0 , αn , βn , zn - заданные числа, ситуация для оценок (14), (16) полностью совпадает с изложенной выше для первой краевой задачи. k k k k Найдём недостающие оценки ∆k , M(11),1 , M(12),1 и ∆k , M(11),n-1 , M(12),n-1 1 n-1 в узлах x1 и xn-1 для второй и третьей краевых задач Lk соответственно. h Для третьей краевой задачи в узле x1 при фиксированном k запишем разложения u1 - hu1 + 564 h3 hk (k) h2 u1 - u1 + · · · + (-1)k u1 = u0 , 2! 3! k! (17) Сходимость матричного метода численного интегрирования краевых задач . . . u1 - hu1 + h2 hk-1 (k) u1 + · · · + (-1)k-1 u = u0 . 2! (k - 1)! 1 (18) Умножив обе части равенства (17) на α0 , а (18) - на β0 и сложив, вместо первого уравнения системы (2), записанного при i = 1, получим α0 u1 + (-α0 h + β0 ) u1 + · · · + (-1)k-1 -α0 hk-1 hk (k) + β0 u = k! (k - 1)! 1 = α0 u0 + β0 u0 = z0 . Тогда локальная матрица Ak третьей краевой задачи примет вид 1  hk hk-1 2 α0 -α0 h + β0 α0 h - β0 h . . . (-1)k-1 -α0 + β0  2! k! (k - 1)!   1 h2 hk h ...  2! k! Ak =  q1 p1 1 ... 0 1   q q1 + p 1 p1 ... 0  1  ... ... ... ... ... (k-2) q1 ... ... ... 1      .     Имеем ∆k = (-1)2k+1 -α0 1 hk-1 hk + β0 k! (k - 1)! + (-1)k+1 k-1 am hm + m=0 k-1 hk k! bm hm + ∆k-1 = O(hk-1 ) + ∆k-1 . (19) 1 1 m=0 Очевидно, полученное равенство (19) не нарушится при α0 = 0, что соответствует второй краевой задаче. В работе [12] в узле x1 показана справедливость оценки (16) при j = 1 для второй и третьей краевых задач и справедливость оценки k-1 k M(12),1 = M(12),1 + ak hk-1 + bk hk , где ak , bk - не зависящие от h величины для третьей краевой задачи, откуда непосредственно следует k-1 k M(12),1 = M(12),1 + O(hk-1 ) (20) как для третьей, так и для второй краевой задачи. Из соотношений (19), (16) при j = 1 и (20) имеем следующие оценки для второй и третьей краевых задач в узле x1 : ∆k = ∆2 + O(h2 ), 1 1 k M(11),1 k M(12),1 = = 2 M(11),1 2 M(12),1 (21) 3 + O(h ), (22) + O(h2 ). (23) 565 М а к л а к о в В. Н. k k Оценки ∆k , M(11),n-1 , M(12),n-1 вычисляются аналогично изложенным n-1 выше способом. 3. Устойчивость матричного метода. Разностную краевую задачу ai ui-1 + bi ui + ci ui+1 = fi , u0 = u0 , un = un , i = 1, 2, . . . , n - 1, (24) называют устойчивой [4], если она имеет единственное решение при произвольных u0 , un , f1 , f2 , . . . , fn-1 и имеет место оценка u C f , (25) где C - не зависящее от h число, w - некоторая норма вектора w = [w0 w1 . . . wn ] , включая и выбранную в [4] норму w = max (|w0 |, |w1 |, . . . , |wn |) . (26) Правомерность выбора нормы в форме (26) для линейной разностной краевой задачи L2 обоснована в [4]. h При выборе нормы в форме (26) в неравенстве (25) следует положить u = [u0 u1 . . . un ] , f = [u0 un f1 f2 . . . fn-1 ] . Разностную краевую задачу (24) называют хорошо обусловленной [4], если она при всех достаточно больших n имеет единственное решение при произвольных u0 , un , f1 , f2 , . . . , fn-1 и если значения u0 , u1 , . . . , un , образующие решение, удовлетворяют неравенству |ui | P max (|u0 |, |un |, |f1 |, . . . , |fn-1 |) , i = 0, 1, . . . , n, (27) где P - не зависящее от h число. Выполнение неравенств (27) означает, что погрешности, допущенные при задании u0 , un , f1 , f2 , . . . , fn-1 , не накапливаются и не приводят к возрастающим ошибкам с увеличением n [4]. Неравенство (25) и совокупность неравенств (27) эквивалентны в смысле совпадения оценок ввиду выбора нормы в форме (26) и в силу совпадения их правых частей, следовательно, хорошая обусловленность разностной краевой задачи влечёт за собой её устойчивость. Достаточным условием (критерием) хорошей обусловленности разностной краевой задачи (24) является выполнение неравенств |bi | - |ai + ci | |bi | + |ai | + |ci | Θ > 0, i = 1, 2, . . . , n - 1, (28) при условии гладкости коэффициентов |ak - am | |ck - cm | D D k-m r , n k-m r , n |bk - bm | D D > 0, r > 0 k-m r n , (29) в случае их ограниченности |ai | 566 Q, |bi | Q, |ci | Q, i = 1, 2, . . . , n - 1 (30) Сходимость матричного метода численного интегрирования краевых задач . . . и в предположении немалости одновременно gi = max (|ai |, |bi |, |ci |) B > 0, i = 1, 2, . . . , n - 1, (31) где числа Θ, Q, B не зависят от номера узла i и значения n [4]. Для первой разностной краевой задачи Lk , k 3, с учётом (9), (14), (16) h имеем bk - ak + ck i i i = bk + ak + ck i i i = ∆k i k M(13),i k M(11),i k M(13),i = ≈ Mk - - M(11),i - k (13),i + ∆k i k M(13),i + k M(12),i k M(13),i k M(12),i = k M(13),i k k ∆k - M(11),i + M(12),i i k k M(11),i + ∆k + M(12),i i 2 2 ∆2 + O(h3 ) - M(11),i + M(12),i + O(h3 ) i 2 2 M(11),i + O(h3 ) + ∆2 + O(h3 ) + M(12),i + O(h3 ) i 2 2 ∆2 - M(11),i + M(12),i i 2 2 M(11),i + ∆2 + M(12),i i ∆2 i 2 M(13),i = 2 M(11),i 2 M(13),i M2 - - M(11),i - k (13),i + ∆2 i 2 M(13),i = + = ≈ 2 M(12),i 2 M(13),i 2 M(12),i = k M(13),i b2 - a2 + c2 i i i , (32) 2 + a2 + c2 bi i i откуда следует, что если устойчивой является краевая задача L2 , т.е. наряду h с выполнением (29)-(31) имеет место неравенство |b2 | - |a2 + c2 | i i i |b2 | + |a2 | + |c2 | i i i Θ > 0, i = 1, 2, . . . , n - 1, (33) то устойчивой является и задача Lk при условии справедливости оценок (29)- h (31) для неё. Выявим достаточное условие устойчивости для задачи L2 . В силу того, h что, как видно из (11), a2 + c2 = i i имеем |b2 | = i 1 pi 1 pi 2 - + + = 2 > 0, h2 2h h2 2h h ∆2 2 - h2 qi 2 i = = 2 - qi = a2 + c2 - qi . i i 2 2 h h M(13),i (34) Тогда, по крайней мере при q(x) < 0, x ∈ [a, b], (35) 567 М а к л а к о в В. Н. из равенства (34) следует |b2 | = |a2 + c2 | + |qi | или i i i |b2 | - |a2 + c2 | = |qi |, i i i i = 1, 2, . . . , n - 1. (36) При условии справедливости оценок (29)-(31) для первых разностных краевых задач Lk , k 2, из (33) с учётом (36) имеем оценку h |qi | |b2 | - |a2 + c2 | i i i 2 | + |a2 | + |c2 | = |b2 | + |a2 | + |c2 | |bi i i i i i Θ > 0, i = 1, 2, . . . , n - 1, где Θ= 1 min|q(x)|, 3B [a, b] что обеспечивает устойчивость задачи L2 и, как следствие, задачи Lk при h h выполнении условия (35). 2, для узлов Для второй и третьей разностных краевых задач Lk , k h сетки xi , i = 2, 3, . . . , n - 2, ситуация полностью совпадает с изложенной выше. Для узлов x1 , xn подстановка в левую часть (32) оценок (21)-(23) вместо (14), (16) приводит к аналогичному (32) соотношению. Следовательно, методы оценки устойчивости, разработанные в классическом методе численного интегрирования разностной краевой задачи L2 [4, 6, 8], могут быть испольh зованы при исследовании матричным методом задачи Lk , k 3. h В [4] доказана эквивалентность достаточного критерия устойчивости (28)- (31) задачи (24) следующему условию: |bi | |ai | + |ci | + δ, i = 1, 2, . . . , n - 1, где δ > 0 не зависит от номера узла i и числа n [4]. При i = 1, i = n - 1 разностное уравнение (24) содержит только две неизвестные: u1 , u2 при i = 1 и un-2 , un-1 при i = n - 1; вследствие чего критерий устойчивости в [5] представлен как |bi | |b1 | |c1 |, |ai | + |ci |, i = 2, 3, . . . , n - 2, |bn-1 | |an-1 |, (37) (38) (39) причём строгое неравенство должно иметь место хотя бы в одном из выражений (37)-(39). 4. Сходимость матричного метода. Решение u0 , u1 , . . . , un разностной краk евой задачи Lk u = fh при измельчении сетки Dh сходится к решению дифh ференциальной краевой задачи Lu = f , если [u] - u → 0 при h → 0. Если при этом имеет место неравенство [u] - u 568 P hk , Сходимость матричного метода численного интегрирования краевых задач . . . где P , k - не зависящие от h числа, то в этом случае говорят, что разностная краевая задача имеет k-тый порядок точности или что имеет место сходимость порядка hk (задача сходится со скоростью O(hk )) [4, 6]. k Разностная краевая задача Lk u = fh аппроксимирует дифференциальную h краевую задачу Lu = f , если k δfh → 0 при h → 0, k где δfih - невязка в узле с номером i сетки Dh , характеризующая отличия, возникающие при подстановке сеточной функции [ui ] в разностную краевую k задачу; иными словами, подстановка [u] в Lk u = fh приводит к выражению h k k Lk [u] = fh + δfh . h Если при этом выполнено неравенство k δfh Chk , где C, k - не зависящие от h числа, то говорят, что имеет место аппроксимация порядка hk , или порядка k относительно величины h [4, 6]. Согласно доказанной в [4] теореме, если разностная краевая задача Lk u = h k = fh аппроксимирует дифференциальную краевую задачу Lu = f c порядком hk и является устойчивой, то решение разностной краевой задачи ui сходится к [ui ], i = 1, 2, . . . , n - 1, причём порядок точности совпадает с порядком аппроксимации k относительно величины h. Аналогичный вывод сделан и в работе [6]. Для первой разностной краевой задачи (10) в работе [12] получены следующие оценки: k Chk для чётного k; δfh k δfh Chk-1 для нечётного k. Откуда имеем, что задачи L2m и L2m+1 для любого натурального числа m 1 h h имеют одинаковый порядок аппроксимации. Следовательно, при условии выполнения критерия устойчивости получаем, что первая разностная краевая задача сходится со скоростью O(hk ) для чётного k и со скоростью O(hk-1 ) для нечётного k. В работе [12] показано, что порядок аппроксимации второй и третьей краk евых задач Lk u = fh на единицу меньше коэффициента k независимо от h его чётности, следовательно, в этом случае задача Lk сходится со скоростью h O(hk-1 ). Использование предложенного в [12] метода повышения порядка аппроксимации на единицу приведёт к скорости сходимости O(hk ) второй и третьей краевых задач Lk . h 5. Численный эксперимент. При выполнении численного эксперимента для разностной задачи Lk в качестве меры суммарной оценки относительh ной погрешности в узлах сетки Dh принята величина k E = n-1 i=1 ([ui ] - n-1 i=1 [ui ] ui )2 · 100%, 569 М а к л а к о в В. Н. которую можно трактовать как некий аналог коэффициента вариации в статистике, являющегося безразмерной мерой рассеяния (разброса) и характеризующего эту меру в процентах [13]. При выполнении численного эксперимента полагалось, что n = 20. Рассмотрим ОДУ2 2 4 sin x 4 cos x u - 2u = + , x x x x2 имеющее аналитическое решение u - u = C1 x 4 + (40) C2 - sin x , x для которого при любых граничных условиях выполнено обеспечивающее устойчивость условие (35). В табл. 1-3 приведены результаты исследования краевых задач для ОДУ2 (40) с граничными условиями первого рода u(1) = -1.827, u(5) = 8.742, (41) с граничными условиями третьего рода 2u(1) - u (1) = -5.012, u(5) + 2u (5) = -22.631 (42) и с граничными условиями смешанного типа 2u(1) - u (1) = -5.012, u (5) = 6.945. (43) В третьей строке табл. 1 (и ниже в табл. 4) указано число узлов n сетки Dh , при котором была достигнута погрешность E 9 , соответствующая задаче L9 , а в четвёртой строке для числа узлов n сетки Dh указано суммарное число h арифметических операций Q при вычислении обратных матриц от локальных матриц Ak , i = 1, 2, . . . , n - 1, когда каждое из них вычислялось по формуле i [14]: 8 1 1 Q = (k + 1)3 - (k + 1)2 - (k + 1) + 1. 3 2 6 Отметим, что в задачах для уравнения (40) при всех k = 2, 3, . . . , 9 в узлах xi , i = 2, 3, . . . , n - 2, сетки Dh с граничными условиями первого рода (41), условиями третьего рода (42) и условиями смешанного типа (43) устойчивость, а вместе с ней и сходимость, обеспечивались за счёт выполнения неравенства (35). В узле x1 устойчивость обеспечена тем, что |c1 | |b1 | 6.089 · 10-1 для всех упомянутых выше граничных условий, что повлекло за собой выполнение (37). В узле xn-1 выполнение условия |an-1 | |bn-1 | 570 9.176 · 10-1 Сходимость матричного метода численного интегрирования краевых задач . . . повлекло за собой выполнение обеспечивающего устойчивость неравенства (39). Аналогичная ситуация наблюдалась при исследовании краевых задач для ряда других ОДУ2, в частности для u + 1 1 u - 2 u = 1, x x для которого условия устойчивости (35), (37), (39) оказались не нарушенными. Рассмотрим ОДУ2 u - 2 2 u + 2 u = x cos x, x x (44) имеющее аналитическое решение u = C1 x + C1 x2 - x cos x, для которого при любых граничных условиях нарушено обеспечивающее устойчивость условие (35). В табл. 4-6 приведены результаты исследования краевых задач для ОДУ2 (44) с граничными условиями первого рода u(5) = 22.331, u(13) = 23.953 (45) и с граничными условиями смешанного типа: u(5) + 3u (5) = 17.597, u(13) = 23.953, (46) u(5) + 3u (5) = 17.597, u (13) = 4.054. (47) Приведённые в табл. 4-6 результаты свидетельствуют о наличии сходимости краевых задач для уравнения (44) с граничными условиями (45)-(47) для всех k = 2, 3, . . . , 9, однако неравенство (35), которое обеспечивает выполнение достаточного условия устойчивости в узлах сетки xi , i = 2, 3, . . . , n - 2, было нарушено для всех типов используемых граничных условий. Неравенства (37), (39) были нарушены для задач с граничными условиями смешанного типа, в частности |c1 | 1.278 |b1 | для условий (46), (47) и |an-1 | |bn-1 | 1.003 для граничного условия (47). Полученные результаты объясняются лишь достаточностью условий устойчивости (35), (37), (39). Аналогичная ситуация наблюдалась при исследовании краевых задач для ряда других ОДУ2, в частности для u + 2 u + u = 1. x Очевидно, для ОДУ2 u - 2 4 u - 2 u = x ln x x x (48) 571 572 k E ,% k k E ,% k k E ,% n Q k -2 -1 1.050 · 10 2 4.503 · 10 91 29029 4 -4 2.306 · 10 79 43924 5 -4 -6 9.877 · 10 32 28416 6 5.148 · 10 29 38599 7 -6 7.678 · 10 29 38020 8 2.691 · 10 3 -2 4.042 · 10 4 -3 1.069 · 10 5 -3 -4 1.504 · 10 6 3.481 · 10 7 -5 6.211 · 10 8 2.862 · 10 3 -2 4.183 · 10 4 -3 1.107 · 10 5 -3 -4 1.551 · 10 6 3.567 · 10 7 -5 5.979 · 10 8 Значения погрешностей для решения граничной задачи (40), (43) [Estimated values of the errors for the solution of the boundary problem (40), (43)] -1 1.055 · 10 2 -3 7.975 · 10 865 139265 3 Значения погрешностей для решения граничной задачи (40), (42) [Estimated values of the errors for the solution of the boundary problem (40), (42)] 3.572 · 10 1580 104280 2 Значения погрешностей для решения граничной задачи (40), (41) [Estimated values of the errors for the solution of the boundary problem (40), (41)] -6 -6 -7 5.297 · 10-6 9 Таблица 3 5.149 · 10-6 9 Таблица 2 7.648 · 10-7 20 52280 9 Таблица 1 М а к л а к о в В. Н. k E ,% k k E ,% k k E ,% n Q k -1 8.370 · 10 2 -2 4.842 · 10 94 29986 4 -4 4.176 · 10 85 47260 5 -4 -6 1.957 · 10 52 46176 6 1.506 · 10 48 63888 7 -6 5.884 · 10 29 38020 8 1.187 · 10 3 -1 9.661 · 10 4 -3 5.250 · 10 5 -4 -5 6.136 · 10 6 1.457 · 10 7 -6 2.599 · 10 8 9.320 · 10 3 -2 1.641 · 10 4 -2 6.713 · 10 5 -4 -4 1.153 · 10 6 3.545 · 10 7 -6 4.686 · 10 8 Значения погрешностей для решения граничной задачи (44), (47) [Estimated values of the errors for the solution of the boundary problem (44), (47)] -1 6.076 · 10 2 7.293 · 10 2300 37000 3 Значения погрешностей для решения граничной задачи (44), (46) [Estimated values of the errors for the solution of the boundary problem (44), (46)] -2 8.106 · 10 2400 15800 2 Значения погрешностей для решения граничной задачи (44), (45) [Estimated values of the errors for the solution of the boundary problem (44), (45)] -7 -7 -7 1.552 · 10-7 9 Таблица 6 1.672 · 10-7 9 Таблица 5 5.896 · 10-7 20 52280 9 Таблица 4 Сходимость матричного метода численного интегрирования краевых задач . . . 573 М а к л а к о в В. Н. условие (35) выполнено для всех k = 2, 3, . . . , 9 в узлах xi , i = 2, 3, . . . , n - 2 любого отрезка интегрирования [a, b]. Численный эксперимент для ОДУ2 (48) со смешанными граничными условиями u(-5) + u (-5) = -33.617, u(-1) = -1.050 показал отсутствие сходимости для всех k = 2, 3, . . . , 9. Оказалось, что |c1 | |b1 | 1.237, причём с ростом k погрешность увеличивается от E 2 = 1215.2 до E 9 = = 70796.6. Для того же ОДУ2 со смешанными граничными условиями u(1) + u (1) = -2.594, u(5) = 151.979 оказалось, что |c1 | |b1 | 0.802, а сама разностная краевая задача оказалась сходящейся, чего и следовало ожидать. Результаты численного эксперимента свидетельствуют, что выполнение достаточных условий устойчивости (35), (37), (39) обеспечивает устойчивость и, как следствие, сходимость разностных краевых задач для ОДУ2. При невыполнении хотя бы одного из условий (35), (37), (39) вопрос о сходимости остаётся открытым. Анализ четвёртых строк в табл. 1, 4 ещё раз подтверждает сделанный в [12] вывод о целесообразности использования на практике задач Lk с чётным h номером k. Выводы. В настоящей работе теоретически установлены следующие факты. 1. Достаточный критерий устойчивости матричного метода при использовании многочленов Тейлора степени три и выше эквивалентен известному достаточному критерию устойчивости традиционного метода сеток для численного интегрирования разностных краевых задач для ОДУ2 с граничными условиями первого, второго и третьего рода. 2. Скорость сходимости матричного метода интегрирования разностных краевых задач с граничными условиями первого рода пропорциональна степени используемого многочлена Тейлора в случае, когда эта степень является чётной, и пропорциональна числу, меньшему степени на единицу, когда эта степень является нечётной. 3. При интегрировании разностных краевых задач с граничными условиями второго и третьего рода скорость сходимости метода пропорциональна степени используемого многочлена Тейлора независимо от его чётности и меньше него на единицу.

About the authors

Vladimir N Maklakov

Samara State Technical University

Email: makvo63@yandex.ru
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
(Cand. Phys. & Math. Sci.), Associate Professor, Dept. of Higher Mathematics and Applied Informatics

References

  1. Keller H. B. Accurate difference methods for nonlinear two-point boundary value problems // SIAM J. Numer. Anal., 1974. vol. 11, no. 2. pp. 305-320. doi: 10.1137/0711028.
  2. Lentini M., Pereyra V. A variable order finite difference method for nonlinear multipoint boundary value problems // Math. Comp., 1974. vol. 28, no. 128. pp. 981-1003. doi: 10.1090/s0025-5718-1974-0386281-4.
  3. Keller H. B. Numerical solution of boundary value problems for ordinary differential equations: Survey and some resent results on difference methods / Numerical solutions of boundary value problems for ordinary differential equations; ed. A. K. Aziz. New York: Academic Press, 1975. pp. 27-88. doi: 10.1016/b978-0-12-068660-5.50007-7.
  4. Годунов С. К., Рябенький В. С. Разностные схемы. М.: Наука, 1977. 439 с.
  5. Формалеев В. Ф., Ревизников Д. Л. Численные методы. М.: Физматлит, 2004. 400 с.
  6. Самарский А. А. Теория разностных схем. М.: Наука, 1977. 656 с.
  7. Самарский А. А., Гулин А. В. Численные методы. М.: Наука, 1989. 432 с.
  8. Самарский А. А., Гулин А. В. Устойчивость разностных схем. М.: Наука, 1973. 416 с.
  9. Boutayeb A., Chetouani A. Global extrapolations of numerical methods for a parabolic problem with nonlocal boundary conditions // International Journal of Computer Mathematics, 2003. vol. 80, no. 6. pp. 789-797. doi: 10.1080/0020716021000039209.
  10. Boutayeb A., Chetouani A. A numerical comparison of different methods applied to the solution of problems with non local boundary conditions // Applied Mathematical Sciences, 2007. vol. 1, no. 44. pp. 2173-2185, http://www.m-hikari.com/ams/ams-password-2007/ams-password41-44-2007/boutayebAMS41-44-2007.pdf.
  11. Радченко В. П., Усов А. А. Модификация сеточных методов решения линейных дифференциальных уравнений с переменными коэффициентами на основе тейлоровских разложений // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2008. № 2(17). С. 60-65. doi: 10.14498/vsgtu646.
  12. Маклаков В. Н. Оценка порядка аппроксимации матричного метода численного интегрирования краевых задач для линейных неоднородных обыкновенных дифференциальных уравнений второго порядка // Вестн. Сам. гос. техн. ун-та. Сер. Физ.мат. науки, 2014. № 3(36). С. 143-160. doi: 10.14498/vsgtu1364.
  13. Закс Л. Статистическое оценивание. М.: Статистика, 1976. 598 с.
  14. Турчак Л. И. Основы численных методов. М.: Наука, 1987. 320 с.

Statistics

Views

Abstract - 42

PDF (Russian) - 6

Cited-By


Article Metrics

Metrics Loading ...

PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2015 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies