Все выпуски

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

Получено альтернативное представление для невязки, представляющее собой сумму двух положительно определенных квадратичных форм от новых конечных разностей, заданных на границе. Матрица первой формы имеет второй порядок и очевидный спектр. Элементы второй матрицы порядка N + 1 выражаются через многочлены Чебышева, матрица обратима и такова, что обратная матрица имеет трехдиагональный вид. Эта особенность позволяет получить для спектра матрицы верхние и нижние оценки, не зависящие от размерности N. Данное обстоятельство позволяет провести исследование на качество аппроксимации для разных размерностей N и весовых коэффициентов ω ∈ [−1, 1]. Показано, что наилучшее приближение дает параметр ω = 0, а невязка стремится к нулю с ростом N.

In the previous paper of the authors the parameter family of finite-dimensional spaces of special quadratic splines of Lagrange’s type has been defined. In each space, as a solution to the initial-boundary problem for the simplest heat conduction equation, we have proposed the optimal spline, which gives the smallest residual. We have obtained exact formulas for coefficients of this spline and its residual. The formula for coefficients of this spline is a linear form of initial finite differences. The formula for the residual is a positive definite quadratic form of these quantities, but because of its bulkiness it is ill-suited for analyzing of the approximation quality of the input problem at the variation with the parameters.

For the purposes of the present paper, we have obtained an alternative representation for the residual, which is the sum of two positive definite quadratic forms of the new finite differences defined on the boundary. The matrix of the first form has second order and the apparent spectrum. The elements of the second matrix of order N + 1 are expressed in terms of Chebyshev’s polynomials, the matrix is invertible and the inverse matrix has a tridiagonal form. This feature allows us to obtain, for the spectrum of the matrix, upper and lower bounds that are independent of the dimension N. Said fact allows us to make a study of the quality of approximation for different dimensions N and weights ω ∈ [−1, 1]. It is shown that the parameter ω = 0 gives the best approximation and the residual tends to zero as N increasing.