Numerical integration by the matrix method and evaluation of the approximation order of difference boundary value problems for non-homogeneous linear ordinary differential equations of the fourth order with variable coefficients

Abstract


The use of the second degree Taylor polynomial in approximation of derivatives by finite difference method leads to the second order approximation of the traditional grid method for numerical integration of boundary value problems for non-homogeneous linear ordinary differential equations of the second order with variable coefficients.The study considers a previously proposed method of numerical integration using matrix calculus which didn’t include the approximation of derivatives by finite difference method for boundary value problems of non-homogeneous fourth-order linear ordinary differential equations with variable coefficients. According to this method, when creating a system of difference equations, an arbitrary degree of the Taylor polynomial can be chosen in the expansion of the sought-for solution of the problem into a Taylor series.In this paper, the possible boundary conditions of a differential boundary value problem are written both in the form of derived degrees from zero to three, and in the form of linear combinations of these degrees. The boundary problem is called symmetric if the numbers of the boundary conditions in the left and right boundaries coincide and are equal to two, otherwise it is asymmetric.For a differential boundary value problem, an approximate difference boundary value problem in the form of two subsystems has been built. The first subsystem includes equations for which the boundary conditions of the boundary value problem were not used; the second one includes four equations in the construction of which the boundary conditions of the problem were used.Theoretically, the patterns between the order of approximation and the degree of the Taylor polynomial were identified.The results are as follows:a) the approximation order of the first and second subsystems is proportional to the degree of the Taylor polynomial used;b) the approximation order of the first subsystem is two units less than the degree Taylor polynomial with its even value and three units less with its odd value;c) the approximation order of the second subsystem is three units less than the degree Taylor polynomial regardless of both even-parity or odd-parity of this degree, and the degree of the highest derivative in the boundary conditions of the boundary value problem.The approximation order of the difference boundary value problem with all possible combinations of boundary conditions is calculated.The theoretical conclusions are confirmed by numerical experiments.

Full Text

Введение. Предложенный в работе [1] метод, использующий средства матричного исчисления и численного интегрирования краевых задач для неоднородных линейных обыкновенных дифференциальных уравнений второго порядка (ОДУ2) с переменными коэффициентами позволяет удерживать произвольное число членов разложения в ряд Тейлора искомого решения задачи, отказавшись при этом от аппроксимации производных конечными разностями. Известно, что использование конечных разностей приводит ко второму порядку аппроксимации (ПА) при численном интегрировании как краевых задач для ОДУ2 [2–7], так и ряда краевых задач для дифференциальных уравнений в частных производных [5–11]. Последнее обусловлено тем, что при аппроксимации производных удерживается всего три члена в разложении в ряд Тейлора искомого решения задачи или, что то же самое, был использован многочлен Тейлора второй степени. 1. Обозначения, терминология, некоторые предварительные результаты и постановка задачи. Далее будем придерживаться принятых в [5] обозначений: а) — область интегрирования, ограниченная отрезком [ , ], — узлы сетки, определяемые значениями = 0 + , = 1, 2, . . . , , 0 = , ( ) = , = / , + 1 — число узлов сетки; б) ( ) — непрерывная функция, являющаяся точным решением краевой задачи для ОДУ4; в) [ ] — сеточная функция, совпадающая с точным решением в узлах сетки ; г) ( ) — искомая сеточная функция. Для краткости примем для любой функции обозначение ( ) = , где — узел сетки . В дальнейшем опустим индекс в наименованиях сеточных функций [ ] , ( ) и будем оговаривать особо случаи, в которых будет использоваться непре138 рывная функция ( ), являющаяся точным решением, при сохранении обозначений ( ) = для нее в узлах сетки. Примем следующую терминологию. 1. Перечислим виды граничных условий (ГУ) в левой границе 0 = сетки , которые будут использоваться ниже в краевых задачах для ОДУ4 (9): ( ) = 0 , (1) ( ) = 0 , (2) ( ) = (3) ( ) = 0 , 0, (4) 0 ( ) + 0 ( ) + 0 ( ) + 0 ( ) = 0 , (5) 0 , 0 , 0 , 0 , 0 , 0 — заданные числа; числа 0 , 0 , 0 , 0 где 0 , 0 , 0, не обращаются в нуль одновременно. Условия (1)–(4) назовем граничными условиями той или иной степени в соответствии с используемой степенью производной в них; условие (5) назовем смешанным ГУ. Аналогично могут быть перечислены ГУ в правой границе = сетки . 2. Дифференциальную краевую задачу (ДКЗ) для ОДУ4 назовем симметричной, если количество ГУ (1)–(5) в левой и правой границах области интегрирования совпадает и равно двум; в противном случае задачу назовем несимметричной. 3. Дифференциальную краевую задачу для ОДУ4 назовем смешанной, если она содержит одно или несколько ГУ (5). Смешанная ДКЗ может оказаться как симметричной, так и несимметричной. Например, задача ( ) (4) ( ) + ( ) ( ) + ( ) ( )+ + ( ) ( ) + ( ) ( ) = ( ), [ , ], (6) 0 ( ) + 0 ( ) + 0 ( ) + 0 ( ) = 0 , ( ) = , ( ) = , ( ) = является смешанной несимметричной ДКЗ. Согласно [5], разностная краевая задача (РКЗ), аппроксимирующая дифференциальную задачу, может быть записана в компактной символической форме как = , (7) где есть степень используемого многочлена Тейлора в разложении искомого решения задачи в ряд Тейлора. Подробный вид равенства (7) будет дан ниже. Сеточная функция , = 0, 1, . . . , , являющаяся решением некоторой РКЗ, записанной в форме (7), при подстановке в уравнения этой РКЗ обратит их в верные равенства. В [5] показано, что подстановка в уравнения задачи (7) значений сеточной функции [ ], отличающихся от , приведет к некоторому отличию от верных равенств. Эти отличия и характеризует невязка . Иными словами, подстановка [ ] в задачу (7) приведет к равенству [ ] = + . (8) Согласно [5, 7], РКЗ (8) аппроксимирует краевую за дифференциальную дачу на точном решении ( ), если норма > 0 при > 0. Если при этом 139 имеет место неравенство 6 , > 0, > 0 — некоторые постоянные, не зависящие от , то говорят, что имеет место аппроксимация порядка относительно величины . В соответствии с [5] норму невязки задачи (8) можно записать как ( ) = max | 0 |, | 1 |, | 2 |, . . . , | | , где второй нижний индекс означает принадлежность компоненты к номеру уравнения задачи (8) при начале нумерации с нуля. Поставим целью исследование возможности численного интегрирования матричным методом краевых задач c различными ГУ для неоднородных линейных ОДУ4 с переменными коэффициентами ( ) (4) ( ) + ( ) ( ) + ( ) ( ) + ( ) ( ) + ( ) ( ) = ( ), (9) где ( ), ( ), ( ), ( ), ( ), ( ) — заданные функции, дифференцируемые нужное число раз, ( ) — искомая функция, являющаяся точным решением задачи; с последующим вычислением ПА в зависимости от степени используемого многочлена Тейлора. 2. Численное интегрирование краевых задач для ОДУ4. Для ДКЗ (6) составим ее аппроксимирующую РКЗ. В соответствии с матричным методом численного интегрирования составим систему линейных алгебраических уравнений (СЛАУ) для пятиточечного шаблона 2 , 1 , , +1 , +2 , = 2, 3, . . . , 2, в которую при фиксированном > 4 внесем: а) четыре многочлена Тейлора степени , полученных из четырех разложений в ряд Тейлора искомого точного решения ( ) в окрестностях слева и справа от некоторого внутреннего узла (центрального узла шаблона) , = 2, 3, . . . , 2, сетки ; б) уравнения ( (4) ) ( ) ( ) + + + = , = 0, 1, . . . , 4, + полученные дифференцированием раз обеих частей ОДУ4 (9) и записанные в узле . В итоге получим замкнутую СЛАУ 3 4 2 + 4 8 + 16 (4) + · · · + ( 2) ( ) = 2 2 , 2! 3! 4! ! 2 3 4 + + (4) + · · · + ( 1) ( ) = 1 , 2! 3! 4! ! 2 3 4 (4) ( ) + + + + · · · + = +1 , + 2! 3! 4! ! 2 3 4 (4) (10) +4 + 8 + 16 ( ) = + 2 + · · · + 2 +2 , 2! 3! 4! ! (4) + + + + = , ( ) ( ) ( ) ( ) (4) (5) + + + + + + + + + = , ... ( 4) + · · · + ( ) = ( 4) . 140 В матричной форме система уравнений (10) имеет вид = в обозначениях 3 4 2 1 2 4 2! 8 3! 16 4! . . . ( 2) ! 3 2 4 1 3! . . . ( 1) ! 2! 4! 2 3 4 1 ... 2! 3! 4! ! = 1 , (11) 2 3 4 2 4 2! 8 3! 16 4! . . . 2 ! ... 0 ... ... ... ... ... ... 0 ( 4) ... ... ... ... ... [ ] (4) ( ) = , . . . (12) [ = 2 1 +1 +2 . . . ( 4) ] , (13) где — символ транспонирования. Здесь и ниже первый верхний индекс означает степень используемого многочлена Тейлора, если речь не идет о показателях алгебраических степеней, степенях производных, символов обратных матриц и транспонирования; второй из пары верхних индексов в наименованиях матриц и их элементов, если таковой присутствует, означает номер центрального узла пятиточечного шаблона, в котором записана матрица. Матрицы , как и ранее в [12–14], будем называть локальными матрицами. Предполагая существование обратной матрицы = ( ) 1 от локальной матрицы , найдем = (14) или в координатной форме = 11 2 + 12 1 + 13 +1 + 14 +2 + 15 + +1 ( 5) , (15) ( 5) , (16) 1 =6 = 21 2 + 22 1 + 23 +1 + 24 +2 + 25 + +1 2 =6 = 31 2 + 32 1 + 33 +1 + 34 +2 + 35 + +1 ( 5) 3 , (17) =6 ... ( ) = ( +1)1 2 + ( +1)2 1 + ( +1)3 +1 + + ( +1)4 +2 + ( +1)5 + +1 ( 5) ( +1) , (18) =6 где , = 1, 2, . . . , +1, = 1, 2, . . . , +1 — элементы матрицы в узле . При = 4 последние суммы в соотношениях (15)–(18) отсутствуют. 141 Из равенств (15), являющихся разностными уравнениями четвертого порядка [5] для пятиточечного шаблона 2 , 1 , , +1 , +2 , = 2, 3, . . . , 2, составим СЛАУ 11 12 13 14 + = 2 1 +1 +2 15 15 15 15 15 (19) = +1 1 ( 5) = + , = 2, 3, . . . , 2. =6 15 Обратную матрицу от локальной матрицы и ее элементы в случае, когда при построении разностного уравнения не было использовано никакого ГУ, будем ниже обозначать как и соответственно. Система (19) содержит 4 уравнения с неизвестными , = 0, 1, . . . , , т.е. необходимо составить дополнительно еще четыре уравнения. Однако заметим, что еще не были использованы четыре ГУ задачи (6). Исследуем смешанное ГУ (5). Имеем следующие многочлены при фиксированном > 4: 2 3 ( ) 2 8 + · · · + ( 2) = 0 , 2! 3! 2 ! 2 1 3 (4) 2 ( ) 1 8 + · · · + ( 2) = 0 , 2 2 2 + 4 2! 2 3! 2 ( 1)! 2 2 (4) 3 (5) 2 ( ) 2 2 2 8 2 + · · · + ( 2) 2 = 0 , 2 +4 2! 3! ( 2)! 2 2 (5) 3 (6) 3 ( ) (4) 2 8 2 + · · · + ( 2) 3 = 2 2 2 + 4 0. 2! 3! ( 3)! 2 2 2 2 + 4 (20) (21) (22) (23) Обе части равенства (20) умножим на 0 , равенства (21) — на 0 , равенства (22) — на 0 , равенства (23) — на 0 и сложим, в итоге получим ) ( 2 0 2 0 + 0 2 + 0 2 + ( 2 0 + 0 ) 2 + 4 2! ( 3 ) 0 2 0 + 8 +4 2 0 + 0 2 + ···+ 3! 2! ( 0 1 0 2 0 3 ) ( ) 0 + ( 2) 3 8 +4 2 + = 0 ! ( 1)! ( 2)! ( 3)! 2 или ( ) 0 [ 2 ] + 1 [ 2 ] + 2 [ 2 ] + 3 [ 0 , 2 ] + · · · + [ 2 ] = (24) где 142 0 2 0 = 0 , 1 = 2 0 + 0 , 2 = 4 2 0 + 0 , 2! ( 0 1 0 2 0 3 ) 0 = ( 2) 3 8 +4 2 + , = 3, 4, . . . , . ! ( 1)! ( 2)! ( 3)! Сохраним приведенные выше выкладки при составлении СЛАУ (19), а для учета смешанного ГУ (5) при построении первого дополнительного разностного уравнения в СЛАУ (10) при = 2 вместо первого приближенного равенства используем многочлен (24); в этом случае вместо (11)–(13) получим 2 = 0 1 1 1 2 ... ( 4) 2 1 2 2 2! 2 2! 2 2 4 2! 2 2 ... ... ... ... 3 3 3! 3 3! 3 8 3! 2 ... ... 4 4 4! 4 4! 4 16 4! 2 ... ... ... . . . ( 1) ! ... ! ... 2 ! , ... 0 ... 0 ... 2 (25) [ ] (4) ( ) 2 = 2 2 2 , . . . 2 2 2 ] [ 2 = 0 1 3 4 2 2 2 . . . 2( 4) . Предполагая существование обратной матрицы 2 = ( 2 ) 1 , вместо (14) найдем 2 = 2 2 . (26) Первая строка матричного равенства (26) после преобразований примет вид соотношения +1 2 2 2 2 2 2 13 14 1 ( 5) 11 12 + , 2 1 + 2 2 3 2 4 = 2 + 2 0 15 15 15 15 2 2 15 =6 15 которое примем в качестве первого дополнительного разностного уравнения, при построении которого использовано смешанное ГУ (5). Напомним, что при = 4 сумма в правой части последнего соотношения отсутствует и при этом задача фактически сводится к классическому методу сеток. При построении трех оставшихся дополнительных уравнений, учитывающих три оставшихся ГУ задачи (6), воспользуемся описанной выше процедурой, в которой в СЛАУ (10) при = 2 вместо четвертого приближенного равенства используем последовательно один из трех следующих многочленов: 2 3 (4) 1 ( ) 2 + 8 2 + · · · + 2 1 = , (27) 2! 3! ( 1)! 2 2 (4) 3 (5) 2 ( ) 2 + 2 2 + 8 2 + · · · + 2 2 = , 2 + 4 2! 3! ( 2)! 2 2 (5) 3 (6) 3 ( ) (4) 2 + 8 2 + · · · + 2 3 = 2 + 2 2 + 4 . 2! 3! ( 3)! 2 2 + 2 2 + 4 143 Объединим четыре дополнительных уравнения в систему +1 2 2 2 2 2 2 13 14 1 ( 5) 11 12 + = + + , 1 3 4 2 2 2 2 2 2 2 2 0 15 15 15 15 15 15 =6 ( 2) ( 2) ( 2) 12 2 13 11 + = 4 ( 2) 3 ( 2) ( 2) 1 ( 2) 15 15 15 15 +1 ( 2) ( 2) 14 ( 5) 1 = + + , 2 ( 2) 2 ( 2) =6 15 15 ( 2) ( 2) ( 2) 2 12 11 4 ( 2) 3 + ( 2) 13 = ( 2) ( 2) 1 15 15 15 15 +1 ( 2) ( 2) 14 ( 5) 1 = 2 + + ( 2) , ( 2) 2 15 =6 15 ( 2) ( 2) ( 2) 2 11 4 12 3 + ( 2) 13 = ( 2) ( 2) ( 2) 1 15 15 15 15 +1 ( 2) ( 2) 1 14 ( 5) = + + . 2 ( 2) 2 ( 2) 15 =6 15 (28) Здесь и ниже матрицы и коэффициенты разностного уравнения, при построении которого не были использованы ГУ, что уже было выше отмечено, или было использовано ГУ в форме значений искомой непрерывной функции ( ) в границе сетки , обозначены как и соответственно; было использовано ГУ в форме значений производной первой степени — как и соответственно; было использовано ГУ в форме значений производной второй степени — как и соответственно; было использовано ГУ в форме значений производной третьей степени — как и соответственно; бы ло использовано смешанное ГУ — как и соответственно. Указанные обозначения введены для внесения ясности в силу того, что во всех уравнениях в СЛАУ (28) центральным узлом шаблона является либо узел 2 , либо узел 2 . Последнее приводит к совпадению номеров пар индексов в наименованиях ряда коэффициентов в уравнениях системы (28). Вследствие этого зафиксировать отличия одного уравнения от другого оказалось возможным лишь различными наименованиями входящих в них коэффициентов. Совокупность СЛАУ (19), (28) есть РКЗ, аппроксимирующая ДКЗ (6), решение которой позволит найти искомые значения сеточной функции , = 0, 1, . . . , . Как только числа , = 0, 1, . . . , , станут известными, в центральных узлах , = 2, 3, . . . , 2, пятиточечных шаблонов по формулам (16)–(18) вычислим производные искомой функции ( ) вплоть до степени . Недостающие значения первой производной 0 , 1 , 1 , , если они отсутствуют в условии ДКЗ, найдем по формулам (21), (27) и с использованием разложений 2 3 (4) 1 ( ) 2 2 + · · · + ( 1) 1 , 2! 3! ( 1)! 2 1 ( ) 2 3 (4) = 2 + 2 + 2 + · · · + . 2 + 2! 3! ( 1)! 2 1 = 2 2 + 1 144 (29) (30) По формулам, аналогичным (21), (27), (29), (30), можно найти недостаю( ) ( ) ( ) ( ) щие значения производных 0 , 1 , 1 , , = 2, 3, . . . , . 3. Вычисление порядка аппроксимации разностных краевых задач для ОДУ4. При вычислении ПА, с которым разностная краевая задача аппроксимирует дифференциальную, в [5] обоснована целесообразность разбиения РКЗ на подзадачи (подсистемы); поэтому в отдельные подзадачи будем выделять СЛАУ, в уравнения которой не входят постоянные значения из ГУ, и СЛАУ, в уравнения которой входят все постоянные значения из ГУ, а затем отдельно будем вычислять ПА этих двух подзадач. Отметим, что выше разбиение на эти две подзадачи уже осуществлено в форме СЛАУ (19) и (28). Первую подзадачу (19) запишем в компактной символической форме в обозначениях [5]: , (31) , = , вторую подзадачу (28) — в виде , = , , (32) где второй нижний индекс означает, что рассматривается смешанная ДКЗ. или = будем понимать Ниже под обозначениями , = , первую подзадачу. эту же подзадачу для краткости Наряду с обозначением , = , будем ниже обозначать и как , , возможно, без второго нижнего индекса в случае, когда будет рассматриваться несмешанная ДКЗ. Рассмотрим первую и вторую подзадачи (31) и (32). Подстановка [ ] в подзадачи (31) и (32) приведет к и + , , [ ] = , (33) , [ ] = , + , . (34) Вычислим ПА разностных подзадач (33) и (34). Для подзадачи (33) в соответствии с [5] в качестве оценки величины невязки примем норму ( ) , = max | 2 |, | 3 |, . . . , | ( 2) | , (35) где компонента, второй нижний индекс которой совпадает с номером центрального узла шаблона, характеризует меру отличий, появление которой обусловлено уравнением подзадачи (33) c номером на единицу меньше этого второго нижнего индекса, и норму для подзадачи (34) ( ) , = max | 1 |, | 2 |, | 3 |, | 4 | , (36) где второй нижний индекс означает принадлежность компоненты к номеру уравнения подзадачи (34) при начале нумерации с единицы. Норму всей РКЗ (33), (34) в соответствии с [5] запишем как ( ) = max , , , . (37) 145 Исследуем подзадачу (33). При фиксированном > 4 при построении всех разностных уравнений сохраним все приведенные выше выкладки для задачи (19) с тем лишь отличием, что в СЛАУ (10) вместо четырех первых приближенных равенств используем следующие точные равенства: 3 2 ( ) [ ] 8 [ [ ] = [ 2 ] 2 , ] + · · · + ( 2) 2! 3! ! 3 2 ( ) [ ] [ ] + [ ] [ [ ] = [ 1 ] 1 , ] + · · · + ( 1) 2! 3! ! ( ) 3 2 [ ] = [ +1 ] +1 , [ ] + [ ] + [ ] + [ ] + ··· + 2! 3! ! 2 3 ( ) [ ] + 2 [ ] + 4 [ ] + 8 [ [ ] = [ +2 ] +2 , ] + ··· + 2 2! 3! ! [ ] 2 [ ] + 4 +1 (2 ) , , , ( +1) ( ) = ( +1 ), ( , где 2 +2 ) — допол 1 +1 +2 = ( +1)! нительные члены разложений в ряд Тейлора в форме Лагранжа [15]. Локальная матрица (11) сохранит свой вид; вместо (12)–(14) получим [ ] (4) ( ) [ ] = [ ] [ ] [ ] [ , ] [ ] . . . [ ] [ 2 ] 2 [ 1 ] 1 [ +1 ] +1 [ +2 ] +2 [ ] = , ... ( 4) [ ] = [ ]. (38) Выпишем первую строку матричного равенства (38), являющуюся разностным уравнением: ( ) ( ) ( ) 11 [ 2 ] 2 + 12 [ 1 ] 1 + 13 [ +1 ] +1 + + 14 ( [ +2 ] +2 ) + 15 2 + +1 ( 5) 1 2 = [ 2 ] =6 или [ ] 12 13 14 11 [ ] [ ] + [ ] [ +2 ] = 2 1 +1 15 15 15 15 15 = 2 + 146 +1 1 ( 5) 11 2 + 12 1 + 13 +1 + 14 +2 . (39) 2 15 15 =6 Сравнивая уравнения (19) и (39), заключаем, что последняя дробь в (39) . Это позволяет записать есть, в соответствии с (33), величина невязки систему невязок задачи так: { + + + { , 11 2 12 1 13 +1 14 +2 , (40) , = = 15 = 2, 3, . . . , 2 = 2, 3, . . . , 2. есть величина невязки, появляющаяся в уравнении с номером 1 Здесь в подзадаче (33). подзадачи (33). При вычислении алгебВычислим норму невязки , раических дополнений транспонированных локальных матриц ( ) в нижние индексы их обозначений будем вносить для ясности наименования используемых обратных матриц от локальных матриц . элемента матрицы ( ) ; Рассмотрим алгебраическое дополнение 11, 11 здесь определяется соотношением (11). Используя известные свойства определителя [16], имеем 2 . . . 2 2 2 4 2! . . . 2! 2! 3 3 3 8 3! . . . 3! 3!4 4 4 11, = 16 4! . . . = 4! 4! ... ... ... . . . . . . . . . . . . 1 1 1 ( 1) 1 1 0 . . . ( 1)! ( 1)! 2 ( 1)! ( 1) ! 2 ( )! 0 . . . 0 ! 2 . . . 2 2 4 2! . . . 2! 3 3 8 3! . . . = 3!4 + 4 ! 4! 16 . . . 4! ... ... . . . . . . . . . . . . 1 1 1 0 . . . ( 1)! 2 ( 1)! 2 . . . 2 2 4 2! . . . 2! 3 3 3! 8 3! . . . + ( 1) + 4 4 ! 16 . . . 4! 4! ... ... . . . . . . . . . . . . 1 1 1 ( 1) 0 . . . ( 1)! 2 ( 1)! . . . 2 2 . . . 2! 3 2!3 3! . . . 3!4 + ( 2) +1 + 4 ! . . . 4! 4! ... ... . . . . . . . . . . . . 1 ( 1) 1 0 . . . ( 1)! ( 1)! 147 2 2! 3 3! + 4 4! ... 1 ( 1) 1 ( 1)! 2 2! 3 3! 4 4! ... 2 2 4 2! 3 8 3! 4 16 4! ... 1 ( 1)! 2 1 ( 1)! ... 0 1 . . . . . . . . . ... . . . . . . . . . 3 3 { } ( 1) ( 1) 3 + · · · + 1 2 + 11 + 11 ; = +3 здесь и ниже в фигурных скобках приведены лишь суммы степеней, а сомножители, не зависящие от , опущены; — коэффициенты, не зависящие от ; , , , , , , , , — некоторые функции от , , , , и их производных. Другими словами, 11, 3 3 ( 1) + 11, . (41) = +3 Имеем ( 1) 11, 2 2! 3 3! = 4 4! ... 1 ( 1) 1 ( 1)! 2 2! 3 3! 4 4! ... 2 2 4 2! 3 8 3! 4 16 4! ... 1 ( 1)! 2 1 ( 1)! 1 ... 0 . . . . . . . . . ... . . . . . . . . . { } 3 6 6 + · · · + 1 2 3 , (42) =6 где — коэффициенты, не зависящие от . Сравнивая (41) и (42), заключаем, что при > 5 (при = 4 первая полученная ниже формула теряет смысл) в соотношении (41) младшая степень первого слагаемого в первой сумме превосходит младшую степень второго слагаемого как минимум на две единицы, следовательно, можно положить ( 1) 11, 11, . (43) Несколько повторных использований последней формулы приводит к 4 11, 4 11, . (44) Аналогичным образом могут быть получены приближенные равенства 4 1 , 4 1 , , = 2, 3, 4, 5. (45) Отсутствие в настоящей работе полного вывода формул (45) обусловлено лишь громоздкостью выкладок. При исследовании краевых задач для ОДУ2 вывод формул, аналогичных (44), (45), дан в [12]. 148 На основании оценок (44), (45) и очевидных равенств 1 15 имеем 1 15 = 4 , Вычислим оценки 1 , имеем 2 2 2 4 2 4 2!3 11, = 2! 3 2!3 8 43! 3!4 3! 16 4 4! 4 12, 4! 4! 1 , 15, 4 1 , 4 15, , = 1, 2, 3, 4, , = 1, 2, 3, 4. (46) = 1, 2, . . . , 5. Пренебрегая старшими степенями, ( 2 2 2 2 3 ) 6 = 6 + + 3! 4! 3!4! 3!4! ( 6 6 ( 2 ) ) = 6 , 3! 4 2 ( ) 6 6 ( 2 ) 4 14, + = 6 + , 3! 4 2 ( ( ) 8 2 4 3 ) = 24 6 + + 6 4 , 3! 4! 3!4! 3!4! ( ) ( ) 4 13, 24 6 = 6 4 + , 3! 4! 4 15, = 10 . (47) (48) (49) (50) (51) Отметим закономерность в парах формул (47), (48) и (49), (50): знаки первых слагаемых совпадают, вторых — противоположны. Для рядов Тейлора 2 3 + · · · = 2! 3! ( ) ( ) = ( 1) + ( 1) = 1 + 1 , ! ! 1 = + =0 +1 = + + = +1 2 3 + + · · · = 2! 3! ( ) ( ) = + = +1 + +1 , ! ! =0 = +1 , — многочлены Тейлора степени , , — дополнительные где 1 +1 1 +1 члены, в [12] показана справедливость следующих формул: +1 1 = ( +2 ), +1 + 1 = ( +1 ) (52) 149 для нечетного и +1 1 = ( +1 ), +1 + 1 = ( +2 ) (53) для четного . Аналогично тому, как было сделано в [12], можно показать справедливость следующих формул: +2 2 = ( +2 ), +2 + 2 = ( +1 ) (54) +2 + 2 = ( +2 ) (55) для нечетного и +2 2 = ( +1 ), для четного . С учетом соотношений (46)–(51), из (40) при = 2, 3, . . . , 2 имеем 11 2 + 12 1 + 13 +1 + 14 +2 15 [ ( ) ( ) ] 6 2 2 + + 2 +2 (4 ) 1 (4 + ) +1 = 10 ) ( 2 ( 2 + +2 2 +2 ) = 2 4 ) 2 ( 8 ( 1 + +1 1 +1 ) . (56) 2 4 = Подстановка оценок (52), (54) в (56) дает 2 ( +1 ) ( +2 ) 8 ( +1 ) + 2 ( +2 ) 2 4 ( 3 ) + ( 1 ) + ( 3 ) + ( 1 ) = ( 3 ) или ( ) 3 , (57) , = ( 3 ) (58) откуда в соответствии с (35) для нечетного . Подстановка оценок (53), (55) в (56) дает или 150 2 ( +2 ) ( +1 ) 8 ( +2 ) + 2 ( +1 ) 2 4 ( 2 ) + ( 2 ) + ( 2 ) + ( 2 ) = ( 2 ) ( 2 ), (59) откуда , = ( 2 ) (60) для четного . Из оценок (58), (60) следует, что при > 2, — натуральное число, зада2 +1 чи 2 , и , имеют одинаковый ПА. Аналогичная ситуация имела место и при исследовании краевых задач для ОДУ2 [12] и для систем ОДУ2 [13] при > 1 и не имела места для краевых задач для ОДУ3 [14] при использовании пятиточечного шаблона, где ситуация оказалась противоположной: одинаковый ПА имели задачи 2 1 и 2 при > 3; при использовании четырехточечного шаблона в задаче для ОДУ3 закономерности, связанные с зависимостью ПА от четности или нечетности , отсутствовали. Исследуем подзадачу (34). Используем описанную выше процедуру, выполненную при построении первого уравнения системы (28) с заменой приближенных равенств (20)–(23) на следующие точные: 3 2 ( ) [ 2 ] 8 [ ] + · · · + ( 2) [ ] = [ 0 ] 0 , 2! 3! 2 ! 2 1 2 ( ) 1 [ 2 ] 2 [ 2 ] + 4 [ [ ] = [ 0 ] 0 1 , 2 ] + · · · + ( 2) 2! ( 1)! 2 2 2 (4) ( ) 2 [ 2 ] 2 [ ] + 4 [ ] + · · · + ( 2) [ ] = [ 0 ] 0 2 , 2 2! 2 ( 2)! 2 3 2 (5) (4) ( ) 3 3 [ ] 2 [ ] + 4 [ ] + · · · + ( 2) [ 2 ] = [ 2 0 ] 0 . 2 2 2! ( 3)! [ 2 ] 2 [ 2 ] + 4 В частности, вместо (24) получим ( ) 0 [ 2 ] + 1 [ 2 ] + 2 [ 2 ] + 3 [ 2 ] + · · · + [ 2 ] = = 0 0 0 0 0 1 0 0 2 0 0 3 . (61) В итоге локальная матрица 2 сохранит свой вид (25), что приводит к матричному равенству [ 2 ] = 2 [ 2 ] (62) в обозначениях [ ] (4) ( ) [ 2 ] = [ 2 ] [ 2 ] [ 2 ] [ , ] [ ] . . . [ ] 2 2 2 ( ) 1 2 = 2 , 0 0 0 0 0 1 0 0 2 0 0 3 [ 1 ] 1 [ 3 ] 3 [ 4 ] 4 [ 2 ] = . 2 2 2 ... ( 4) 2 151 Выпишем первую строку матричного равенства (62), являющуюся разностным уравнением, при построении которого было использовано смешанное ГУ (5): ( ) ( ) 2 2 11 0 0 0 0 0 1 0 0 2 0 0 3 + 12 [ 1 ] 1 + +1 ( ) ( ) 2 2 2 2 ( 5) 2 = [ 2 ] + 13 [ 3 ] 3 + 14 [ 4 ] 4 + 15 2 + 1 =6 или +1 2 2 2 2 2 [ 2 ] 13 14 1 ( 5) 11 12 [ ] + [ ] [ ] = + + 1 3 4 2 2 2 2 2 2 0 15 15 15 15 2 2 15 =6 15 2 ( + 1 + 2 + 3 ) + 2 + 2 + 2 11 0 0 0 0 0 0 0 0 12 1 13 3 14 4 , (63) 2 15 , в наименовании которой где последняя дробь есть величина невязки второй нижний индекс указывает на использование элементов матрицы 42 при построении разностного уравнения (63). подзадачи (34). Для матрицы 2 , как Вычислим норму невязки , это было сделано выше для матрицы , можно показать справедливость оценок 42 2 1 , 1 42 , = 1, 2, 3, 4. 2 15, 15 42 , = 1, 2, . . . , 5. Отметим, что Вычислим алгебраические дополнения 1 , 42 не является вычисление точных значений алгебраических дополнений 1 , особо трудоемкой процедурой; тем не менее нет строгой необходимости в нахождении точных значений в силу того, что для вычисления ПА (здесь речь не идет о выполнении численных экспериментов) необходимы лишь главные части этих алгебраических дополнений в их разложениях по степеням ; поэтому лишь для сокращения объема выкладок допустим пренебрежение старшими степенями, прежде всего в разложениях по степеням коэффициентов 1 , 2 , 3 , 4 , при нахождении алгебраических дополнений элементов первого столбца локальной матрицы 42 (см. (24), (25)). 42 Анализ матриц (11), (25) указывает на независимость значения 11, от чисел 0 , 0 , 0 , 0 , тогда как оставшиеся оценки алгебраических дополнений зависят от значений этих чисел [16]; действительно, пренебрегая старшими степенями, имеем 42 11, = 6 6 ( 3! ( 2 2 2 2 3 ) 2 ) + + 6 2 2 6 , (64) 4! 3!4! 3!4! 2 ( 5 6 7 ) ( 7 2 5 2 7 ) 2 2 2 + + 2 2 4 + 3 4 18 12 18 ( ( 7 2 5 2 6 ) 2 5 ) 3 2 3 + + 4 2 3 2 4 + 12 4 3 42 12, = 1 152 ( 2 1 2 ) 3 2 4 2 3 + 2 2 , (65) 3 ( ) 42 13, 3 3 2 4 3 2 3 + 2 2 + 2 1 2 , ( 2 2 2 2 1 2 ) 42 14, 3 2 4 + 2 3 + , 12 6 ( 3 2 2 1 3 ) 42 + . 15, = 6 4 2 12 12 (66) (67) (68) разобьем на две составляющие: Невязку = 1 2 ( + 1 + 2 + 3 ) 11 0 0 0 0 0 0 0 0 , 2 15 (69) 2 + 2 + 2 12 1 13 3 14 4 2 15 (70) = 2 и исследуем их отдельно. Число комбинаций одновременно ненулевых коэффициентов 0 , 0 , 0 , 0 смешанного ГУ определяется, очевидно, как 41 + 42 + 43 + 44 = 15, где — число сочетаний из по . Главные части 1 , 2 , 3 , 4 коэффициентов левой части равенства (61) в зависимости от значений 0 , 0 , 0 , 0 приведены в табл. 1. Здесь единица означает присутствие соответствующего коэффициента в смешанном ГУ, нуль — его отсутствие. Отметим, что четыре последние строки в табл. 1 соответствуют ГУ той или иной степени в форме одного слагаемого (1)–(4). В главных частях постоянные сомножители, не зависящие от , опущены. Таблица 1 Главные части коэффициентов левой части равенства (61) [The main parts of coefficients of the left side equality (61)] № 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 0 1 1 1 0 0 0 1 1 1 0 1 1 0 0 0 0 1 0 0 1 1 0 1 1 0 1 1 0 1 0 0 0 0 1 0 1 0 1 1 0 1 1 1 0 0 1 0 0 0 0 1 0 1 1 0 1 1 1 1 0 0 0 1 1 0 0 0 0 0 0 2 0 0 0 0 0 0 0 0 3 0 2 2 0 0 0 2 0 0 0 0 0 0 0 0 0 2 0 0 0 2 3 4 2 0 0 0 0 0 0 0 0 0 0 0 0 3 0 2 0 0 0 3 0 2 0 0 2 0 0 0 2 0 0 0 0 0 4 0 3 0 2 0 Если опустить не зависящие от постоянные сомножители и пренебречь старшими степенями, то соотношения (64)–(68) и данные табл. 1 позволяют 153 записать следующие предварительные оценки: 42 11, 6 , 42 12, + 42 13, 42 15, 42 + 14, 4 6 , (71) 3 3 , (72) (73) 3 = 1 , 4 на основании которых и очевидных равенств 42 1 42 15 = 42 1 , 42 15, , (74) = 1, 2, 3, 4, из (69) получим = 1 2 ( + 1 + 2 + 3 ) 11 0 0 0 0 0 0 0 0 2 15 42 ( + 1 + 2 + 3 ) 11, 0 0 0 0 0 0 0 0 42 15, 6 ( 0 3 + 0 2 + 0 + 0 ) ( 2 ) = 4 6 ( 0 3 + 0 2 + 0 + 0 ) ( 2 ) = , (75) 4 а из (70) — = 2 2 + 2 + 2 12 1 13 3 14 4 2 15 42 + 42 + 42 ) ( +1 ) ( 12, 13, 14, 42 15, 3 3 ( +1 ) 1 ( +1 ) = = ( 3 ). (76) 6 4 3 Оценка (75) для каждого набора 0 , 0 , 0 , 0 требует дальнейшего уточнения результата с использованием данных табл. 1, тогда как оценка (76) дает окончательный результат. Пусть 0 = 0. В этом случае, что следует из табл. 1, 4 = 0 . Тогда, пренебрегая старшими степенями, из (75) получим оценку ( 0 3 + 0 2 + 0 + 0 ) ( 2 ) 0 ( 2 ) = ( 3 ), 4 0 справедливую для всех строк табл. 1 с номерами 3, 5, 6, 8–11, 15. Пусть 0 =0, 0 = 0. В этом случае, что следует из табл. 1, 4 = 0 2 . Тогда, пренебрегая старшими степенями, из (75) получим оценку 1 1 154 ( 0 3 + 0 2 + 0 + 0 ) ( 2 ) 0 ( 2 ) = ( 3 ), 4 0 2 справедливую для всех строк табл. 1 с номерами 2, 4, 7, 14. Пусть 0 =0, 0 = 0, 0 = 0. В этом случае, что следует из таблицы 1, 4 = 0 3 . Тогда, пренебрегая старшими степенями, из (75) получим оценку ( 0 3 + 0 2 + 0 + 0 ) ( 2 ) 0 2 ( 2 ) = ( 3 ), 4 0 3 справедливую для всех строк табл. 1 с номерами 1, 13. для строк 13, 14 и 15 табл. 1 указыПолученные выше оценки невязок 1 вают на независимость ПА от степени старшей производной в ГУ, записанных в форме одного слагаемого. Особое внимание уделим оставшейся строке 12 табл. 1 ( 0 = 1, 0 = 0, 0 = 0, 0 = 0), соответствующей ГУ ( ) = 0 . В этом случае 1 42 42 , (77) а невязка примет вид = 1 + 2 = 2 2 2 2 11 0 + 12 1 + 13 3 + 14 4 . 2 15 (78) 42 , = 1, 2, 3, 4, что позволяет 42 = 1 , Тождество (77) приводит к 1 , воспользоваться оценками (47)–(50), наличие закономерности в которых в парах формул (47), (48) и (49), (50) (а именно: знаки первых слагаемых совпадают, вторых — противоположны) и использование соотношений (52)–(55) приводит к зависимости значения невязки (78) от четности или нечетности , т.е. здесь правомерны оценки вида (57) для нечетного и (59) для четного . Для любой другой строки табл. 1 не следует ожидать повторения аналогич будет содержать как минимум один ной ситуации в силу того, что невязка дополнительный член степени меньше , поэтому надежда на использование соотношений (52)–(55) отсутствует. Следовательно, для смешанного ГУ, в котором 0 = 1, 0 = 0, 0 = 0, 0 = 0, имеем следующие оценки: ( 2 ) (79) для четного ; ( 3 ) для нечетного и во всех остальных случаях при любых сочетаниях 0 , 0 , 0 , 0 независимо от четности . Заметим, что вывод формул для вычисления алгебраических дополнений 42 , = 1, 2, . . . , 5, в узле 1 , 2 , что соответствует правой границе сетки , приведет к формулам (64)–(68), в которых будут заменены на 2 индексы в функциях , , , , значения которых, как показано выше, отсутствуют в оценках (71)–(74) и поэтому не использованы в (75), (76) при вычислении ПА; последнее означает, что на значение невязки не оказывает влияние граница области интегрирования, в которой ГУ записано. Несмотря на то, что при четном имеет место оценка (79), норма невязки второй подзадачи в соответствии с (36) независимо от четности или нечетности окажется всегда равной , = ( 3 ) (80) 155 за счет оставшихся ГУ задачи (6) или, что то же самое, за счет второго, третьего и четвертого уравнений второй подзадачи (28), записанных, с учетом замечания, в узле 2 . Норму всей РКЗ (33), (34), используя (58), (60), (80), в соответствии с (37) запишем как ( ) = max , , , 6 3 , (81) где — не зависящий от коэффициент, откуда следует ПА всей РКЗ (33), (34), равный 3. Отметим, что оценка (81) справедлива для любой краевой задачи, как симметричной, так и несимметричной, для ОДУ4 при использовании матричного метода в силу того, что среди ГУ задачи обязательно найдутся как минимум два ГУ в форме производной той или иной степени больше нуля или в форме смешанного ГУ, в котором хотя бы два коэффициента из чисел 0 , 0 , 0 , 0 отличны от нуля. При вычислении ПА были использованы главные части величин , = 1, 2, 3, 4, тогда как при выполнении далее численных экспериментов будут использоваться точные значения , = 0, 1, . . . , . 5. Оценка погрешностей. При выполнении численных экспериментов использованы следующие нормы: – в качестве суммарной оценки относительной погрешности v ( ) 2 =0 [ ] · 100%, (82) = [ ] =0 которую можно трактовать как некий аналог коэффициента вариации в статистике, характеризующего меру разброса в процентах [17]; – в качестве оценки абсолютной погрешности [5, 6] (83) = max [ ] , = 0, 1, . . . , . При апробации численного метода и оценке его погрешности использовано ОДУ4 ( ) ( ) + 2 (4) + + 2 + 6 + 4 + = 5 , (84) которое вместе со своим общим решением ( ) 1 + 2 = + 1 + 2 + 3 2 + 4 3 24 взято из [18]. При выполнении расчетов использованы следующие наборы ГУ, [2, 6]: { (2) = 2.285, (2) = 0.480, (85) (6) = 0.746, (6) = 0.159; { (2) = 2.285, (2) = 0.480, (86) (6) = 0.403, (6) = 0.159; { (2) = 0.135, (2) = 0.480, (87) (6) = 0.403, (6) = 2.18 · 10 4 ; 156 { { (2) = 2.285, (2) = 0.135, (2) = 0.480, (6) = 2.18 · 10 4 ; (2) = 2.285, (2) = 0.135, (2) = 0.480, (6) + 2 (6) + 3 (6) + 4 (6) = 0.415. (88) (89) Было принято = 20, = 0.20. Результаты численных экспериментов для решений ( ) и первых производных ( ) приведены в табл. 2, в которой нормы , , характеризующие оценки относительных и абсолютных погрешностей, вычислены по формулам (82), (83) для производных ( ). В рассмотренных задачах по данным анализа табл. 2 с увеличением степени используемого многочлена Тейлора относительная и абсолютная погрешности: а) уменьшаются, как это имело место при исследовании краевых задач для ОДУ2 [1, 12], систем ОДУ2 [13] и ОДУ3 [14]; б) уменьшаются довольно «плавно», что свидетельствует о независимости ПА от четности или нечетности . Отметим, что при = 4 имеем классический метод сеток. Выводы. Сформулируем основные выводы по работе. 1. Для дифференциальной краевой задачи составлена аппроксимирующая ее разностная краевая задача в виде двух подсистем: в первую подсистему вошли уравнения, при построении которых не были использованы граничные условия (ГУ) краевой задачи; во вторую подсистему вошли четыре уравнения, при построении которых были использованы ГУ задачи. 2. Теоретически выявлены закономерности между порядком аппроксимации (ПА) и степенью используемого многочлена Тейлора. Установлено следующее: а) ПА первой и второй подсистем пропорционален степени используемого многочлена Тейлора; б) ПА первой подсистемы меньше степени многочлена Тейлора на две единицы при ее четном значении и меньше на три единицы при ее нечетном значении; в) ПА второй подсистемы меньше степени многочлена Тейлора на три единицы независимо как от четности или нечетности этой степени, так и от степени старшей производной в ГУ краевой задачи. 3. Вычислен ПА разностной краевой задачи со всеми возможными комбинациями ГУ, включая смешанные граничные условия.

About the authors

Vladimir Nikolaevich Maklakov

Samara State Technical University


Candidate of physico-mathematical sciences, Associate professor

Mariya Aleksandrovna Ilicheva

Samara State Technical University


References

  1. Радченко В. П., Усов А. А., "Модификация сеточных методов решения линейных дифференциальных уравнений с переменными коэффициентами на основе тейлоровских разложений", Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2008, № 2(17), 60-65
  2. Keller H. B., "Accurate difference methods for nonlinear two-point boundary value problems", SIAM J. Numer. Anal., 11:2 (1974), 305-320
  3. Lentini M., Pereyra V., "A variable order finite difference method for nonlinear multipoint boundary value problems", Math. Comp., 28:128 (1974), 981-1003
  4. 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, Part I: Survey Lectures, ed. A. K. Aziz, Academic Press, New York, 1975, 27-88
  5. Годунов С. К., Рябенький В. С., Разностные схемы. Введение в теорию, Наука, М., 1977, 439 с.
  6. Формалеев В. Ф., Ревизников Д. Л., Численные методы, Физматлит, М., 2004, 400 с.
  7. Самарский А. А., Теория разностных схем, Наука, М., 1977, 656 с.
  8. Самарский А. А., Гулин А. В., Численные методы, Наука, М., 1973, 432 с.
  9. Самарский А. А., Гулин А. В., Устойчивость разностных схем, Наука, М., 1973, 416 с.
  10. Boutayeb A., Chetouani A., "Global Extrapolations Of Numerical Methods For Solving A Parabolic Problem With Non Local Boundary Conditions", Intern. J. Comp. Math., 80:6 (2003), 789-797
  11. Boutayeb A., Chetouani A., "A numerical comparison of different methods applied to the solution of problems with non local boundary conditions", Appl. Math. Sci., 1:44 (2007), 2173-2185
  12. Маклаков В. Н., "Оценка порядка аппроксимации матричного метода численного интегрирования краевых задач для линейных неоднородных обыкновенных дифференциальных уравнений второго порядка", Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2014, № 3(36), 143-160
  13. Маклаков В. Н., "Оценка порядка аппроксимации матричного метода численного интегрирования краевых задач для систем линейных неоднородных обыкновенных дифференциальных уравнений второго порядка с переменными коэффициентами. Сообщение 2. Краевые задачи с граничными условиями второго и третьего рода", Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 21:1 (2017), 55-79
  14. Маклаков В. Н., Стельмах Я. Г., "Численное интегрирование матричным методом краевых задач для линейных неоднородных обыкновенных дифференциальных уравнений третьего порядка с переменными коэффициентами", Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 22:1 (2018), 153-183
  15. Фихтенгольц Г. М., Курс дифференциального и интегрального исчисления, т. 1, Наука, М., 1970, 608 с.
  16. Курош А. Г., Курс высшей алгебры, Наука, М., 1971, 431 с.
  17. Закс Л., Статистическое оценивание, Статистика, М., 1976, 598 с.
  18. Камке Э., Справочник по обыкновенным дифференциальным уравнениям, Наука, М., 1976, 576 с.

Statistics

Views

Abstract - 18

PDF (Russian) - 4

Cited-By


PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2020 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