A new approach to the formation of systems of linear algebraic equations for solving ordinary differential equations by the collocation method; [Новый подход к формированию систем линейных алгебраических уравнений для решения обыкновенных дифференциальных уравнений методом коллокаций]

A new algorithm for the numerical solution of one-dimensional Cauchy problems and Poisson equations is implemented. The algorithm is based on the collocation method and representation of the solution as an expansion in Chebyshev polynomials. It is proposed instead of the usual approach, which consists in combining all known conditions — differential (the equation itself) and initial / boundary — into one system of approximate linear algebraic equations, to go to the method of solving the problem in several separate stages. First, spectral coefficients are identified that determine the “general” solution of the original problem. The collocation method determines the interpolation coefficients of the derivative of the solution, and thus the expansion coefficients of the solution itself (except for the initial ones). At this stage, the choice of a good basis with discrete orthogonality makes it possible to use very efficient algorithms for finding the desired coefficients. The complexity of reducing the matrix of a system of linear algebraic equations to a diagonal form becomes equivalent to the complexity of multiplying the Chebyshev matrix of coefficients by the vector of the right side of the system. Then the expansion coefficients of the solution itself (except for the first one or two) are obtained by multiplying the known tridiagonal integration matrix (inverse to the Chebyshev differentiation matrix) by the vector of interpolation coefficients of the derivative. At the last stage, considering the initial/boundary conditions select a “particular” desired solution, unambiguously redefining the missing coefficients of the desired expansion. © 2023 Vestnik Novosibirskogo Gosudarstvennogo Universiteta, Seriya: Istoriya, Filologiya. All rights reserved.

Авторы
Sevastianov L.A. , Lovetskiy K.P. , Kulyabov D.S.
Издательство
Saratov National Research State University
Номер выпуска
1
Язык
Русский
Страницы
36-47
Статус
Опубликовано
Том
23
Год
2023
Организации
  • 1 Peoples’ Friendship University of Russia (RUDN University), 6 Miklukho-Maklaya St., Moscow, 117198, Russian Federation
  • 2 Joint Institute for Nuclear Research, 6 Joliot-Curie St., Moscow Region, Dubna, 141980, Russian Federation
Ключевые слова
Chebyshev polynomials; collocation method; discrete orthogonality; Gauss–Lobatto sets; initial boundary value problems; numerical stability
Цитировать
Поделиться

Другие записи