Summary
In this paper the problem of smoothing a given data set by cubicC 2-splines is discussed. The spline may required to be convex in some parts of the domain and concave in other parts. Application of splines has the advantage that the smoothing problem is easily discretized. Moreover, the special structure of the arising finite dimensional convex program allows a dualization such that the resulting concave dual program is unconstrained. Therefore the latter program is treated numerically much more easier than the original program. Further, the validity of a return-formula is of importance by which a minimizer of the orginal program is obtained from a maximizer of the dual program.
The theoretical background of this general approach is discussed and, above all, the details for applying the strategy to the present smoothing problem are elaborated. Also some numerical tests are presented.
Similar content being viewed by others
References
Burmeister, W., Heß, W., Schmidt, J.W.: Convex spline interpolants with minimal curvature. Computing35, 219–229 (1985)
de Boor, C.: A Practical Guide to Splines. Berlin Heidelberg New York: Springer 1978
Dietze, S., Schmidt, J.W.: Determination of shape preserving spline interpolants with minimal curvature via dual programs. J. Approximation Theory52, 43–57 (1988)
Dietze, S., Schmidt, J.W.: Unconstrained dual programs to partially separable constrained programming problems. TU Dresden Preprint 07-10-89 (1989)
Dott, S.L., McAllister, D.F.: Algorithms for computing shpae preserving approximation to data. Numer. Math.46, 159–174 (1985)
Elfving, T., Andersson, L.-E.: An algorithm for computing constrained smoothing spline functions. Numer. Math.52, 583–595 (1988)
Elster, K.-H., Reinhardt, R., Schäuble, M., Donath, G.: Einführung in die nichtlineare Optimierung. Leipzig: Teubner 1977
Irvine, L.D., Martin, S.P., Smith, P.W.: Constrained interpolation and smoothing. Constr. Approximation2, 129–151 (1986)
Reinsch, C.: Smoothing by spline functions. Numer. Math.10, 177–183 (1967)
Rockafellar, R.T.: Convex Analysis. Princeton, N.J.: Princeton Univ. Press 1970
Schmidt, J.W.: On tridiagonal linear complementarity problems. Numer. Math.51, 11–21 (1987)
Schmidt, J.W.: An unconstrained dual program for computing convexC 1-spline approximants. Computing39, 133–140 (1987)
Schmidt, J.W., Heß, W.: Positivity of cubic polynomials on intervals and positive spline interpolation. BIT28, 340–352 (1988)
Schoenberg, I.J.: Spline functions and the problem of graduation. Proc. Nat. Acad. Sci. USA52, 947–950 (1964)
Schwetlick, H., Kunert, V.: Constrained least square approximation byB-splines of arbitrary order. (in preparation)
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Schmidt, J.W., Scholz, I. A dual algorithm for convex-concave data smoothing by cubicC 2-splines. Numer. Math. 57, 333–350 (1990). https://doi.org/10.1007/BF01386414
Received:
Issue Date:
DOI: https://doi.org/10.1007/BF01386414