2024
ВЕСТНИК ТОМСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА
Математика и механика Tomsk State University Journal of Mathematics and Mechanics
№ 88
Научная статья
УДК 517.927.4 MSC: 65Q10, 76M20, 41A50
doi: 10.17223/19988621/88/2
Об интегральном подходе при использовании метода коллокации
Оксана Владимировна Гермидер1, Василий Николаевич Попов2
12 Северный (Арктический) федеральный университет им. М.В. Ломоносова,
Архангельск, Россия 1 [email protected] 2 v.popov@narfu. ru
Аннотация. Предложен матричный вариант реализации метода коллокации с использованием интегрального подхода для построения решения неоднородного дифференциального уравнения четвертого порядка. Метод основан на полиномиальной аппроксимации Чебы-шева производной четвертого порядка искомой функции. В качестве базисных функций использованы многочлены Чебышева первого рода. Проведен анализ сходимости решений, полученных методом коллокации с применением интегрального подхода с использованием в качестве узлов точек экстремумов и нулей многочленов Чебышева первого рода. Представлены результаты сравнения с аналогичными результатами, представленными в открытой печати.
Ключевые слова: метод коллокации, многочлены Чебышева первого рода, неоднородные дифференциальные уравнения
Благодарности: Исследование поддержано грантом РНФ № 24-21-00381 «Развитие методов полиномиальной аппроксимации Чебышева для решения нелинейных задач математической физики».
Для цитирования: Гермидер О.В., Попов В.Н. Об интегральном подходе при использовании метода коллокации // Вестник Томского государственного университета. Математика и механика. 2024. № 88. С. 14-25. doi: 10.17223/19988621/88/2
Original article
About the integral approach using the collocation method
Oksana V. Germider1, Vasilii N. Popov2
12 Northern (Arctic) Federal University named after M. V. Lomonosov, Arkhangelsk, Russian Federation 1 [email protected] 2 v.popov@narfu. ru
Abstract. The article describes a matrix method of polynomial Chebyshev approximation using an integral approach to construct a solution to a nonhomogeneous fourth-order differential equation
© О.В. Гермидер, В.Н. Попов, 2024
with mixed boundary conditions of the first kind. The proposed method is based on the expansion of the fourth-order derivative of the desired function into a series in terms of Chebyshev polynomials of the first kind and the representation of the partial sum of this series as a product of matrices whose elements are, respectively, the Chebyshev polynomials and the coefficients in this expansion. In this paper, using analytical formulas for calculating integrals of Chebyshev polynomials, we obtain a representation of the desired function in terms of the product of the matrices defined above. The use of points of extrema and zeros of Chebyshev polynomials of the first kind as nodes, as well as the properties of the sums of products of Chebyshev polynomials at these points, made it possible to reduce the boundary value problem by the collocation method to a system of inhomogeneous linear algebraic equations with a sparse matrix of this system. It is shown that the solution constructed in this way satisfies the differential equation at all nodes, including the boundary ones, in contrast to the approximate solution obtained by approximating the exact solution in the form of a finite sum of the Chebyshev series. The effectiveness of the proposed method is demonstrated by considering a boundary value problem with a known analytical solution. The convergence analysis of the constructed solution is carried out. Keywords: collocation method, Chebyshev polynomials of the first kind, inhomogeneous differential equations
Acknowledgments: The research was supported by the Russian Academy of Sciences grant No. 24-21-00381 "Development of Chebyshev polynomial approximation methods for solving nonlinear problems of mathematical physics".
For citation: Germider, O.V., Popov, V.N. (2024) About the integral approach using the collocation method. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 88. pp. 14-25. doi: 10.17223/19988621/88/2
Введение
Хорошо известно, что численное решение обыкновенных дифференциальных уравнений четвертого и более высокого порядка, встречающихся в различных разделах гидродинамики, в области биотехнологий, геологии, конструирования авиационно-космической техники и т.д. [1-5], вызывает ряд трудностей. Это связано с наличием в уравнениях производных высокого порядка, оказывающих существенное влияние на обусловленность исходной краевой задачи. В связи с этим значительное внимание уделяется разработке высокоточных численных методов.
В представленной работе предлагается реализация метода полиномиальной аппроксимации Чебышева в матричной форме с применением интегрального подхода для решения неоднородного дифференциального уравнения четвертого порядка со смешанными граничными условиями. Производная четвертого порядка искомой функции записывается в виде усеченного ряда по многочленам Чебышева первого рода в матричной нотации. Для восстановления самой функции используются точные формулы для вычисления интегралов от многочленов Чебышева. Далее для нахождения коэффициентов полученного разложения производной поставленная задача сводится к решению линейной системы алгебраических уравнений. В данной работе происходит развитие спектрального метода колло-кации, описанного в [6, 7]. В отличие от этих публикаций производная третьего и более низкого порядка искомой функции, как и сама функция, записывается с использованием произведения матриц, введенного для представления производной четвертого порядка в виде усеченного ряда по многочленам Чебышева
первого рода; проводится анализ сходимости полученных решении при использовании как точек экстремумов, так и нулей многочленов Чебышева при выборе их в качестве узлов интерполирования. Здесь необходимо заметить, что скорость сходимости спектральных методов ограничена только регулярностью интерполируемой функции в отличие от методов конечных разностей и конечных элементов [8], а в случае использования полиномов Чебышева в качестве базисных функций и выбора точек коллокации в нулях или точках экстремумов этих полиномов наблюдается устойчивость к ошибкам округления [9].
Рассмотрим обыкновенное дифференциальное уравнение четвертого порядка и(4) + £ (х,и,и(1) ,и(2) ,и(3)) = /2(х) , -1 < х < 1, (1)
с граничными условиями
«(("!)') = 8п, «(1)((-1)"') = Я,-2, " = 1,2 . (2)
Предполагаем, что функция £ имеет в области определения непрерывные ограниченные частные производные по переменным х, и, и(1), и(2) и и(3), функция /2 тождественно не равна нулю и имеет непрерывную ограниченную производную по х, а сама краевая задача (1), (2) имеет на отрезке [-1; 1] единственное решение.
Метод решения краевой задачи
Представляем функцию и 4) в виде частичной суммы ряда полиномов Чебы-
шева первого рода
|Tj(x) = cos(jarccosx), (j = 0,n)j [10]:
uf (x) = ^ ajTj (x) = T (x) о ISA, (3)
j=0
где T(x) - матрица-строка размером 1 x«' («'=« + 5):
T(x) = (T0 (x), 7j (x),... Тп+Ъ (x), Tn+4 (x)) , знаком о обозначено произведение Адамара двух матриц [11], Is - матрица-строка размером 1х П с ненулевыми элементами /j0i = 1 (i = 0, n ). Здесь и ниже нумерация строк и столбцов начинается с нуля. Матрица-столбец A имеет размер п'х 1, ее элементами являются коэффициенты a (j = 0, n) в разложении (3)
и со > С1 > С2 > сз:
A = (a0,al,...,an_l,an,c0,cl,c2,c3)T . Выберем в качестве узлов интерполяции в (1) точки экстремума многочлена Tn (x) [10]:
x* = cos ———, k = 0, n. (4)
n
При таком выборе узлов интерполяции происходит медленный рост константы Лебега Ли на отрезке [-1;1], и при фиксированном числе узлов она стремится к своему минимальному значению, уменьшая погрешность алгебраического ин-
n
терполирования и обеспечивая меньшую чувствительность решения по отношению к ошибкам округления [12]. В случае, когда n - нечетное число, точное значение Лп можно определить по формуле [13]
n
(2к - 1)л
Лn = - У ctg( ) . (5)
n ¿-^ 4n
к=1
Вблизи границ области x0 = -1 и хп = 1 происходит сгущение узловых точек (4), расположение которых симметрично относительно нуля. Узлы x0 и xn принимают значения -1 и 1 соответственно, что хорошо согласуется с постановкой краевой задачи (1) и (2). И последнее, при увеличении числа интерполяционных узлов (4) в два раза полученные значения функций в точках экстремума до удвоения сохраняются.
Из равенства — (x) = cos( j arccos x) при подстановке в него узлов (4) получаем
T (x) = cos Uj(n - к) , j, к = М . (6)
n
Интегрируя (3) по переменной x, получаем
n
Mf(X) = У aj j Tj (x)dx + c0. (7)
j=0
Для нахождения интегралов от многочленов Чебышева первого рода в (7) учитывая, что — (x) = 1 и — (x) = x [10], имеем с точностью до константы
г г x2
j — (x) dx = x , j — (x) dx = — , (8)
для четных j > 2 согласно [10] получаем
с T,, (x) T , (x)
2j Tj (x) dx = - ^, (9)
j j J +1 J -1
для нечетных j > 3
, (x) T , (x) 2 j(-1)(j+1)/2
2 j -. (x) dx = tJ+1(x) - tJ-1(x) - -. (10)
j j j +1 j -1 j2 -1
Постоянная в (10) получена с использованием представления [12]:
[J/2]
УУУ
- < x)=yv. *=• <»>
(12)
где [' / 2] - целая часть числа 1 /2 .
В случае нечетных ' > 3 из (11) находим эту постоянную:
^(1+1)/2 ^_ 1)( 1 +1)/2 1 + 1 1 - 1 / -1 . Далее, пользуясь (8)—(10), запишем (7) в виде:
и^ (х) = (Т(х)(С °1) + Р1 (х))А , (13)
где С - квадратная матрица размером п 'х п', в которой последний столбец нулевой, ненулевые элементы первой строки
G0,1 = 1 , G0,2 j+1 = (-1)
2 j +1
(2j +1)2 -1
j = 1, [(n'-1) / 2] -1,
(14)
ненулевые элементы второй строки С10 = 1, С12 = -1/2, одиночные ненулевые элементы предпоследней и последней строк Си,_2 п,_3 = 1 / (2п'- 4), Оп,_1п,_2 = = 1 / (2п'- 2), парные ненулевые элементы остальных строк
(-1)''+1
G.
j = 2,n'-3 , i = 1,2 ;
(15)
1,1+(-1)' 21
где I - квадратная матрица размером п'хп' со сроками 18, (х) - матрица-строка размером 1 х п', в которой один ненулевой элемент Р10п+1 = 1. Например, при п = 5 квадратная матрица С имеет вид:
G =
_5 24
— 0--0 0
7 48
Г 1 3
0 1 0 -3 0 48
10 -1 0 0 0 0 0 00 2
Л
0 -4
0 0
0 — 0 0 0 4
0 0 0
0 0 0 0 0 1
000-0 — 0 0 00 8
0 0 0 0
1 10
— 0--
1 10
0 0 0 0 0 12
0 0 0 0 0 0
— 0--
0 0 0 0 0
1 12
J_ 14
1
0--0
14
000000 0 — 00
16
000000 0 0 —0
18 /
Последовательно интегрируя (13) по переменной х, получаем и(2), и(1) и ^
(х) = ^ cij || Tj (х) d2x + c0x + cl =(t(x)(G2°i) + P2(x))a, (16)
j=0
и 2
и®(x) = ^a.|||7;.(x)J3x + c0— + с,х + с2 =(t(x)(G3 °l) + P3(x)U, (17)
j=о 2
u(x) = у a ЦЦ—(x) d4x+c0 +c ^+c2x+С = (T(x)(G4ol) + P4(x))A,
j =0
р (х) - матрица-строка размером 1 х п', в которой ненулевые элементы
Рдл+,- = 1, =^т—^, У = 11 -1, 1 = 2,4 .
Подставляя (3), (4), (6), (13), (16)-(18) в (1) и используя (2), приходим к системе линейных п' -уравнений в матричной форме:
3
ВА = Р , В = £ В , (19)
1=1
где В (I = 1 , 3) - квадратные матрицы размером п 'х п', первые две из них получены из уравнения (1) в узлах (4), последние четыре строки В и В2 нулевые, ^-я ненулевая строка матрицы В определяется с использованием (3): и(4) (х^) = (В \ А, т.е. (В1 )к = Т(х/;) о 15 (к = 0,п), строка В2 с таким же индексом восстанавливается на основе (16)-(18): (х,и(хк),и(1)(х),и(2)(х),и(3)(х)) = (В)* А ; матрица В строится из граничных условий (4), ее ненулевые строки - последние четыре: (Вз1=Т((-1),)(С4о1) + РДН)')' (Вз)и+2+; =т((-1/)(С3°1) + Р3((-1/), / = 1,2.
= (Л (хо X /2 <Л),..., /2 (хи ), Ян, §21, §12,, § 22, ) . Для приведения матрицы В в (19) к разреженной матрице можно воспользоваться свойствами конечных сумм многочленов Чебышева в выбранных узлах интерполяции (4). Левые и правые части уравнения (19) умножаем на матрицу 8
размером п 'х п', в которой ^ = 2п- (Б\) (к = 0, п), J = 1 (У = 1,4). Здесь
верхний индекс Т у В обозначает операцию транспонирования, а штрихом обозначено деление строк и столбцов с индексами 0 и п на 2. В частности, в результате умножения имеем V = 8В , V = SJ, где J и V (' = 1,2) - квадратные
матрицы размером п'хп', для J ^й ненулевой строкой является Т( хк) (к = 0, п). У матрицы V на главной диагонали п +1 элемент равен 1: V** = 1 (к = 0, п), а остальные равны нулю, что следует из свойств конечных сумм полиномов Чебышева, вычисленных в точках экстремума Тп (х) [10, 14]. У матрицы V2 те же элементы равны 1 и дополнительно к ним элементы К2я 1я+( = 1 (/ = 1,4 ). Обозначая через О матрицу размером п'хп', в которой первые к строк равны соответственно Р4 (хк) (к = 0, п), а остальные четыре строки нулевые, в результате умножения 8О имеем матрицу V , в которой отличны от нуля шесть элементов: ^з,2,и+1 = 1/24, У,Хп+1 = 1/8, ¥з 0 п+2= УЗАп+2 = 1/4, ¥зА„+ 4 = 3 = 1. Матрица B, построенная из граничных условий (2), при умножении на 8 слева сохраняет свой вид: 8В3 = В. Таким образом, уравнение (19) приводим к виду:
УА = 8Ж , (20)
V = V + ОВ + В.
Решение уравнения (20) находим Z^-методом в системе компьютерной алгебры Maple [15] с числом значащих цифр 22. Число обусловленности матрицы v системы (20) определяем согласно [16, 17]:
v(V) = ||v ||-||v-1||,
n+4
где норму матрицы V находим как V = max У | V I [16].
II II 0<i<и+4 1
J=0
Функцию u(x) получаем, используя (18) и найденные значения элементов матрицы A.
Аналогично коэффициенты в (18) могут быть найдены, если в качестве узлов в уравнении (1) для переменной х выбрать нули многочлена —и+1 (x) [10]:
. n(2n - 2к +1) —
x* = cos —--1, к = 0, n . (21)
2(n +1)
Результаты вычислений и их анализ
Рассмотрим уравнение [7]
2 2
u(4) + U(2) = л4 sin(ny) - л2 sin(ny), - 2 < y < 2 , (22)
с граничными условиями
u[(-1)i 2] = (-1)^ , uУ((-1)i 2) = -f , i = 1,2 . (23)
Аналитическое решение задачи (22) и (23) имеет вид: ue (y) = sin(ny) [7].
Произведем замену переменных в уравнении (22) и граничных условиях (23):
2 x
y = у, x е[-1;1]. (24)
В результате приходим к уравнению
|]4u^ +(3 Juf =n4 sin(^ j-n2sin(^ j , -1 < x < 1, (25) с граничными условиями
и((-1)')= С-!)'^3 , и?((-1)')=-|, ' = 1,2 . (26)
Явный вид приближенного аналитического решения, полученного методом коллокации с использованием интегрального подхода, для задачи (22), (23) при п = 5 в случае узловых точек (4) записывается следующим образом:
и(у) = 0.06826у9 - 0.5892у7 + 2.548у5 -5.168у3 + 3.142у , а в случае узловых точек (21):
и(у) = 0.07013у9 - 0.5922у7 + 2.549у5 -5.168у3 + 3.142у .
В таблице представлены результаты вычислений с использованием интегрального подхода (1А) на основе (3), (7)-(18) для случая узловых точек, которые являются точками экстремума многочлена Чебышева степени п (4) и нулями (х) (21), где для расчета среднеквадратичного отклонения полученных значе-
ний в контрольных точках построенного решения от соответствующих значений точного ие, вычисленных в этих же точках, как и в [18, 19], применено 100 равно-
мерно распределенных контрольных точек на (-2 / 3; 2 / 3] (, = 1; 100). В этом случае во внутренних точках отрезка у е [-2 / 3; 2 / 3] отклонение определяем согласно [6]:
£ (и. (У,) - и(у1 ))2
е
г=1
Nе,п =
" £и.2(у,)
1 ,=1
где пе - число внутренних точек отрезка [-2 / 3; 2 / 3]. Число обусловленности матрицы системы (20) при п = 5 составляет 46.0 для узловых точек (4) и (21), при п = 11 равно 50.8, при п = 15 равно 52.9.
Порядок сходимости находим, как и работах [18-20]:
N....
Ne.2n+1
Значения отклонения N „ и порядка сходимости гв 1
г =
е.п
п N.,п
Л, (4) ChP, (4) 1Л, (21) ар, (21) [7]
4 2.0 -10-3 7.3 -10-2 7.3 -10-4 7.3 -10-2 1.0 -10-1
5 3.4 -10-5 1.4 -10-2 1.6 -10-5 1.6 -10-2 -
6 2.3 -10-5 3.9-10-3 4.6 -10-6 7.6 -10-3 2.3 -10-3
7 3.0 -10-7 2.7 -10-4 6.110-8 5.7 -10-4 -
8 8.4 -10-8 1.0 -10-4 1.9 -10-8 2.7 -10-4 1.2 -10-5
9 6.7 -10-10 3.7 -10-6 2.0 -10-10 1.0 -10-5 -
10 2.6 -10-10 1.4 -10-6 7.8 -10-11 4.8 -10-6 2.4 -10-7
11 1.3 -10-12 3.2 -10-8 3.9 -10-13 1.2 -10-7 -
12 7.8 -10-13 1.3 -10-8 2.7 -10-13 5.3 -10-8 4.9 -10-10
13 2.7 -10-15 2.0 -10-10 9.3 -10-16 8.8-10-10 -
п г е,п
4 3.0-105 2.0 -104 3.7-106 7.3-103 -
5 2.6-107 4.4 -105 4.1-107 1.3-105 -
6 8.5-109 2.0 -107 4.9-109 8.6-106 -
В таблице приведены значения N „ для численного решения краевой задачи
(25) и (26), представленного в виде усеченного ряда по многочленам Чебышева первого рода. Результаты, полученные в этом случае без использования интегрального подхода, в таблице имеют аббревиатуру ChP. Здесь приходим к системе п +1 линейных уравнений, полученных при использовании точек экстремума многочлена Чебышева степени п (4) и нулей Ти+1 (х) (21) в методе коллокации для (25), и осуществляем замену уравнений согласно граничным условиям (26) в точках, для которых х = х0 и х = хп, соответственно на уравнения
<((-1У)= (-1)'^ , (27)
а в точках x = x и x = xBl - на уравнения
u?((-1У )=-П , i = 1,2 . (28)
Из таблицы видно, что высокая точность полученного решения на основе интегрального подхода (IA) достигается при сравнительно малых значениях n как при использовании точек экстремумов многочленов Чебышева в качестве узлов интерполирования, так и нулей этих полиномов; наблюдается быстрая сходимость. Используя (5), получаем, что для всех значений n, указанных в таблице, отношение Ли / Л2и+1 не превышает 2 для узлов (4). В частности, при n = 5 значение константы Лебега составляет Л5 = 1.98854, а при n = 11 оно равно Лп = 2.900825. Учитывая, что для нечетного n имеет место равенство Ли = Ли_1 [12], где ЛиЧ -значение константы Лебега для узлов (21), получаем, что Л4 = 1.988854 и Л10 = 2.900825 . Таким образом, медленный рост константы Лебега уменьшает погрешность алгебраического интерполирования.
Заключение
В работе предложен матричный метод полиномиальной аппроксимации Чебы-шева с использованием интегрального подхода построения решения неоднородного дифференциального уравнения четвертого порядка со смешанными граничными условиями. Построена матрица, позволяющая получить значения интегралов от полиномов Чебышева при умножении на нее слева матрицы, составленной из этих полиномов. Краевая задача при использовании точек экстремумов и нулей многочленов Чебышева сведена к системе неоднородных линейных уравнений с разреженной матрицей этой системы. Показано, что полученное таким образом решение при интегральном подходе сходится с повышенным порядком при использовании как точек экстремумов, так и нулей многочленов Чебышева в качестве точек коллокации.
Список источников
1. Tac V., Rausch M.K., Costabal F.S., Tepole A.B. Data-driven anisotropic finite viscoelasticity
using neural ordinary differential equations // Computer Methods in Applied Mechanics and Engineering. 2023. V. 411. Art. 116046. doi: 10.1016/j.cma.2023.116046
2. Nanshan M., Zhang N., Xun X., Cao J. Dynamical modeling for non-Gaussian data with high-
dimensional sparse ordinary differential equations // Computational Statistics & Data Analysis. 2022. V. 173. Art. 107483. doi: 10.1016/j.csda.2022.107483
3. Khader M.M., Mahdy A.M.S., Shehata M.M. An Integral Collocation Approach Based on Le-
gendre Polynomials for Solving Riccati, Logistic and Delay Differential Equations // Applied Mathematics. 2014. V. 5. P. 2360-2369. doi: 10.4236/am.2014.515228
4. Gimeno J., Jorba A., Jorba-Cusco M., Miguel N., Zou M. Numerical integration of high-order
variational equations of ODEs // Applied Mathematics and Computation. 2023. V. 442. Art. 127743. doi: j.amc.2022.127743
5. Лун-Фу А.В., Бубенчиков М.А., Жамбаа С., Цыдыпов С.Г. Определение частот поперечных
колебаний переходников и тупиковых ответвлений газопроводов // Вестник Томского
государственного университета. Математика и механика. 2020. № 68. С. 95-105. doi: 10.17223/19988621/68/9
6. Mai-Duy N., Tanner R.I. A spectral collocation method based on integrated Chebyshev poly-
nomials for two-dimensional biharmonic boundary-value problems // Journal of Computational and Applied Mathematics. 2007. V. 201. P. 30-47. doi: 10.1016/j.cam.2006.01.030
7. Mai-Duy N., See H., Tran-Cong T. A spectral collocation technique based on integrated Che-
byshev polynomials for biharmonic problems in irregular domains // Appl. Math. Model. 2009. V. 33 (1). P. 284-299. doi: 10.1016/j.apm.2007.11.002
8. Shao W., Wu X. An effective Chebyshev tau meshless domain decomposition method based
on the integration-differentiation for solving fourth order equations // Appl. Math. Model. 2015. V. 39 (9). P. 2554-2569. doi: 10.1016/j.apm.2014.10.048
9. Baseri A., Abbasbandy S., Babolian E. A collocation method for fractional diffusion equation
in a long time with Chebyshev functions // Applied Mathematics and Computation. 2018. V. 322. P. 55-65. doi: 10.1016/j.amc.2017.11.048
10. Mason J., Handscomb D. Chebyshev polynomials. Florida : CRC Press, 2003.
11. Liu S., Trenkler G. Hadamard, Khatri-Rao, Kronecker and other matrix products // International Journal of Information and Systems Sciences. 2008. V. 4 (1). P. 160-177. Corpus ID: 36767622
12. Ibrahimoglu B.A. Lebesgue functions and Lebesgue constants in polynomial interpolation // Journal of Inequalities and Applications. 2016. V. 93. P. 1-15. doi: 10.1186/s13660-016-1030-3
13. McCabe J.H., Phillips G.M. On a certain class of Lebesgue constants // BIT. 1973. Vol. 13. P. 434-442. doi: 10.1007/BF01933407
14. Гермидер О.В., Попов В.Н. О решении модельного кинетического уравнения ES // Чебышевский сборник. 2022. Т. 23, № 3. С. 37-49. doi: 10.22405/2226-8383-2022-23-337-49
15. Corless R.M., Jeffrey D.J. The Turing factorization of a rectangular matrix // ACM SIGSAM Bulletin. 1997. V. 31 (3). P. 20-30. doi: 10.1145/271130.271135
16. Бахвалов Н.С., Жидков Н.П., Кобельков Г.М. Численные методы. М. : Наука, 1987. 598 с.
17. Антюфеев В.С. Вероятностная оценка числа обусловленности матрицы // Сибирский журнал чистой и прикладной математики. 2018. Т. 18, № 1. С. 28-34. doi: 10.17377/ PAM.2018.18.3
18. Голушко С.К., Идимешев С.В., Шапеев В.П. Метод коллокаций и наименьших невязок в приложении к задачам механики изотропных пластин // Вычислительные технологии. 2013. Т. 18, № 6. С. 31-43.
19. Шапеев В.П., Брындин Л.С., Беляев В.А. hp-Вариант метода коллокации и наименьших квадратов с интегральными коллокациями решения бигармонического уравнения // Вестник Самарского государственного технического университета. Сер. Физико-математические науки. 2022. Т. 26, № 3. С. 556-572. doi: 10.14498/vsgtu1936
20. Chen G., Li Zh., Lin P. A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow // Adv. Comput. Math. 2008. V. 29. P. 113-133. doi: 10.1007/s10444-007-9043-6
References
1. Tac V., Rausch M.K., Costabal F.S., Tepole A.B. (2023) Data-driven anisotropic finite viscoelas-
ticity using neural ordinary differential equations. Computer Methods in Applied Mechanics and Engineering. 411(116046). DOI: 10.1016/j.cma.2023.116046.
2. Nanshan M., Zhang N., Xun X., Cao J. (2022) Dynamical modeling for non-Gaussian data
with high-dimensional sparse ordinary differential equations. Computational Statistics & Data Analysis. 173(107483). DOI: 10.1016/j.csda.2022.107483.
3. Khader M.M., Mahdy A.M.S., Shehata M.M. (2014) An integral collocation approach based
on Legendre polynomials for solving Riccati, logistic and delay differential equations. Applied Mathematics. 5. pp. 2360-2369. DOI: 10.4236/am.2014.515228.
4. Gimeno J., Jorba A., Jorba-Cusco M., Miguel N., Zou M. (2023) Numerical integration
of high-order variational equations of ODEs. Applied Mathematics and Computation. 442 (127743). DOI: j.amc.2022.127743.
5. Lun-Fu A.V., Bubenchikov M.A., Zhambaa S., Tsydypov S.G. (2020) Opredeleniye chastot
poperechnykh kolebaniy perekhodnikov i tupikovykh otvetvleniy gazoprovodov [Determining of transverse vibration frequencies for crossovers and dead ends of gas pipelines]. Vestnik Tomskogo gosudarstvennogo universiteta. Matematika i mekhanika - Tomsk State University Journal of Mathematics and Mechanics. 68. pp. 95-105. DOI: 10.17223/19988621/68/9.
6. Mai-Duy N., Tanner R.I. (2007) A spectral collocation method based on integrated Chebyshev
polynomials for two-dimensional biharmonic boundary-value problems. Journal of Computational and Applied Mathematics. 201. pp. 30-47. DOI: 10.1016/j.cam.2006.01.030.
7. Mai-Duy N., See H., Tran-Cong T. (2009) A spectral collocation technique based on integrated
Chebyshev polynomials for biharmonic problems in irregular domains. Applied Mathematical Modelling. 33(1). pp. 284-299. DOI: 10.1016/j.apm.2007.11.002.
8. Shao W., Wu X. (2015) An effective Chebyshev tau meshless domain decomposition method
based on the integration-differentiation for solving fourth order equations. Applied Mathe-maticalModelling. 39(9). pp. 2554-2569. DOI: 10.1016/j.apm.2014.10.048.
9. Baseri A., Abbasbandy S., Babolian E. (2018) A collocation method for fractional diffusion
equation in a long time with Chebyshev functions. Applied Mathematics and Computation. 322. pp. 55-65. DOI: 10.1016/j.amc.2017.11.048.
10. Mason J., Handscomb D. (2003) Chebyshev Polynomials. Boca Raton: CRC Press. 2003.
11. Liu S., Trenkler G. (2008) Hadamard, Khatri-Rao, Kronecker and other matrix products. International Journal of Information and Systems Sciences. 4(1). pp. 160-177.
12. Ibrahimoglu B.A. (2016) Lebesgue functions and Lebesgue constants in polynomial interpolation. Journal of Inequalities and Applications. 93. pp. 1-15. DOI: 10.1186/s13660-016-1030-3.
13. McCabe J.H., Phillips G.M. (1973) On a certain class of Lebesgue constants BIT Numerical Mathematics. 13. pp. 434-442. DOI: 10.1007/BF01933407.
14. Germider O.V., Popov V.N. (2022) O reshenii model'nogo kineticheskogo uravneniya ES [On the solution of the model kinetic equation ES]. ChebyshevskiiSbornik. 23(3). pp. 37-49. DOI: 10.22405/2226-8383-2022-23-3-37-49.
15. Corless R.M., Jeffrey D.J. (1997) The Turing factorization of a rectangular matrix. ACM SIGSAMBulletin. 31(3). pp. 20-30. DOI:10.1145/271130.271135.
16. Bakhvalov N.S., Zhidkov N.P., Kobelkov G.M. (1987) Chislennyye metody [Numerical methods]. Moscow: Nauka.
17. Antyufeev V. S. (2018) Probabilistic estimation of matrix condition number. Siberian Journal of Pure and Applied Mathematics. 18(1). pp. 28-34. DOI:10.17377/PAM.2018.18.3.
18. Golushko S.K., Idimeshev S.V., Shapeev V.P. (2013) Metod kollokatsiy i naimen'shikh nevyazok v prilozhenii k zadacham mekhaniki izotropnykh plastin [Application of collocations and least residuals method to problems of the isotropic plates theory] Vychislitel'nyye. Tekhnologii. - Computational Technologies. 18(6). pp. 31-43.
19. Shapeyev V.P., Bryndin L.S., Belyayev V.A. (2022) hp-variant metoda kollokatsii i naimen'shikh kvadratov s integral'nymi kollokatsiyami resheniya bigarmonicheskogo uravneniya [The hp-version of the least-squares collocation method with integral collocation for solving a biharmonic equation]. Vestnik Samarskogo gosudarstvennogo tekhnicheskogo universiteta. Seriya Fiziko-matematicheskiye nauki — Journal of Samara State Technical University, Seriya Physical and Mathematical Sciences. 26(3). pp. 556-572. DOI: 10.14498/vsgtu1936.
20. Chen G., Li Zh., Lin P. (2008) A fast finite difference method for biharmonic equations on irregular domains and its application to an incompressible Stokes flow. Advances in Computational Mathematics. 29. pp. 113-133. DOI: 10.1007/s10444-007-9043-6.
Сведения об авторах:
Гермидер Оксана Владимировна - кандидат физико-математических наук, доцент кафедры инженерных конструкций, архитектуры и графики Северного (Арктического) федерального университета им. М.В. Ломоносова, Архангельск, Россия. E-mail: [email protected] Попов Василий Николаевич - доктор физико-математических наук, профессор, профессор кафедры высшей и прикладной математики Северного (Арктического) федерального университета им. М.В. Ломоносова, Архангельск, Россия. E-mail: [email protected]
Information about the authors:
Germider Oksana V. (Candidate of Physics and Mathematics, Northern (Arctic) Federal University named after M.V. Lomonosov, Arkhangelsk, Russian Federation). E-mail: [email protected] Popov Vasilii N. (Doctor of Physics and Mathematics, Northern (Arctic) Federal University named after M.V. Lomonosov, Arkhangelsk, Russian Federation). E-mail: [email protected]
Статья поступила в редакцию 26.04.2023; принята к публикации 10.04.2024
The article was submitted 26.04.2023; accepted for publication 10.04.2024