Table of Content

Open Access iconOpen Access

ARTICLE

A Highly Accurate Technique for Interpolations Using Very High-Order Polynomials, and Its Applications to Some Ill-Posed Linear Problems

by Chein-Shan Liu1, Satya N. Atluri2

Department of Civil Engineering, National Taiwan University, Taipei, Taiwan. E-mail: liuchei.nshan@msa.hinet.net
Center for Aerospace Research & Education, University of California, Irvine

Computer Modeling in Engineering & Sciences 2009, 43(3), 253-276. https://doi.org/10.3970/cmes.2009.043.253

Abstract

Since the works of Newton and Lagrange, interpolation had been a mature technique in the numerical mathematics. Among the many interpolation methods, global or piecewise, the polynomial interpolation p(x) = a0 + a1x + ... + anxn expanded by the monomials is the simplest one, which is easy to handle mathematically. For higher accuracy, one always attempts to use a higher-order polynomial as an interpolant. But, Runge gave a counterexample, demonstrating that the polynomial interpolation problem may be ill-posed. Very high-order polynomial interpolation is very hard to realize by numerical computations. In this paper we propose a new polynomial interpolation by p(x) = ā0 + ā1x/ R0 + ... + ānxn/ R0n, where R0 is a characteristic length used as a parameter, and chosen by the user. The resulting linear equations system to solve the coefficients āα is well-conditioned, if a suitable R0 is chosen. We define a non-dimensional parameter, R0* = R0/ (b-a) [where a and b are the end-points of the interval for x]. The range of values for R0* for numerical stability is identified, and one can overcome the difficulty due to Runge, as well as increase the accuracy and stability in the numerical interpolation by very-high-order polynomials, for these values of R0*. Numerical results which validate the current theory are presented for (i) the first and higher-order derivatives of noisy numerical data [such as determining da/dN in fatigue mechanics], (ii) the solution of the Abel integral equation under noisy data, and (iii) the numerical determination of an inverse Laplace transform under noisy data. These results confirm the validity of the present approach for very high-order polynomial interpolation.

Keywords


Cite This Article

APA Style
Liu, C., Atluri, S.N. (2009). A highly accurate technique for interpolations using very high-order polynomials, and its applications to some ill-posed linear problems. Computer Modeling in Engineering & Sciences, 43(3), 253-276. https://doi.org/10.3970/cmes.2009.043.253
Vancouver Style
Liu C, Atluri SN. A highly accurate technique for interpolations using very high-order polynomials, and its applications to some ill-posed linear problems. Comput Model Eng Sci. 2009;43(3):253-276 https://doi.org/10.3970/cmes.2009.043.253
IEEE Style
C. Liu and S. N. Atluri, “A Highly Accurate Technique for Interpolations Using Very High-Order Polynomials, and Its Applications to Some Ill-Posed Linear Problems,” Comput. Model. Eng. Sci., vol. 43, no. 3, pp. 253-276, 2009. https://doi.org/10.3970/cmes.2009.043.253



cc Copyright © 2009 The Author(s). Published by Tech Science Press.
This work is licensed under a Creative Commons Attribution 4.0 International License , which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • 1441

    View

  • 1100

    Download

  • 0

    Like

Share Link