Abstract
We have analyzed two algorithms, close in composition, for constructing a smoothing spline, which imply a change only in the last link of the spline when new experimental data arrive. The main feature of the N. D. Dicoussar algorithm is the form of a polynomial representation in order to describe a link of the spline. It is shown that a given polynomial is one of the hierarchical form of the Hermitian polynomial. We have proposed a modification to the D. A. Silaev algorithm for constructing a smoothing spline with different orders of smoothness: from zero to the second, aimed at enhancing the stability of this algorithm. To this end, we substantiated recommendations related to the form of polynomials representation, which describe the links of splines of the specified form. For this purpose, we estimated conditionality of matrices used in the algorithm. For the spline of zero-order smoothness, the most advisable is to apply a polynomial in the N. D. Dicoussar form, and for splines with higher orders of smoothness of joining the links, it is appropriate to use different forms of the Hermitian polynomials. Based on computational examples, a possibility was demonstrated to generalize the D. A. Silaev algorithm to construct a spline with links of various lengths, which is determined by the rate of change in the examined parameter. That makes it possible to reduce the volume of information that contains a description of the spline itself, and to prevent such a widespread shortcoming of approximation when using polynomials as parasitic oscillations. It was shown as well that in the presence of significant measurement errors in experimental data there may occur a need to decrease the length of the spline's link (compared to that derived by the D. A. Silaev rule) in order to provide the spline with a property of robustness
Highlights
The task on finding an effective analytical description of a time series is commonly in many application areas
In cases when the study relates to one link of the spline in order to reduce the number of indexes we shall equate the region for determining this link to section [xa; xb ]
Dicoussar polynomial coefficients (2) takes the form: Author of paper [9] obtianed an estimate of the magnitude of parameter q through the value of function f(x) and its derivatives at the ends of segment [xa; xb]. We show that it can be achieved in another way
Summary
The task on finding an effective analytical description of a time series is commonly in many application areas. The task on constructing an optimal smoothing spline implies sorting out all the possible cases of distribution of grid knots among the links of the spline. It is a relevant task to build algorithms for constructing the smoothing splines, close to optimal, provided there is the possibility to modify a small number of spline links in proportion to the arrival of new data. Such algorithms would make it possible to achieve the approximation of an experimental dependence with acceptable quality for much less time than when solving a given problem in classical statement
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have