Mathematical Construction of Interpolation and Extrapolation Function by Taylor Polynomials Arxiv:2002.11438V1 [Math.NA] 26 Fe
Total Page:16
File Type:pdf, Size:1020Kb
Load more
Recommended publications
-
Stable Extrapolation of Analytic Functions
STABLE EXTRAPOLATION OF ANALYTIC FUNCTIONS LAURENT DEMANET AND ALEX TOWNSEND∗ Abstract. This paper examines the problem of extrapolation of an analytic function for x > 1 given perturbed samples from an equally spaced grid on [−1; 1]. Mathematical folklore states that extrapolation is in general hopelessly ill-conditioned, but we show that a more precise statement carries an interesting nuance. For a function f on [−1; 1] that is analytic in a Bernstein ellipse with parameter ρ > 1, and for a uniform perturbation level " on the function samples, we construct an asymptotically best extrapolant e(x) as a least squares polynomial approximant of degree M ∗ given explicitly. We show that the extrapolant e(x) converges to f(x) pointwise in the interval −1 Iρ 2 [1; (ρ+ρ )=2) as " ! 0, at a rate given by a x-dependent fractional power of ". More precisely, for each x 2 Iρ we have p x + x2 − 1 jf(x) − e(x)j = O "− log r(x)= log ρ ; r(x) = ; ρ up to log factors, provided that the oversampling conditioning is satisfied. That is, 1 p M ∗ ≤ N; 2 which is known to be needed from approximation theory. In short, extrapolation enjoys a weak form of stability, up to a fraction of the characteristic smoothness length. The number of function samples, N + 1, does not bear on the size of the extrapolation error provided that it obeys the oversampling condition. We also show that one cannot construct an asymptotically more accurate extrapolant from N + 1 equally spaced samples than e(x), using any other linear or nonlinear procedure. -
On Multivariate Interpolation
On Multivariate Interpolation Peter J. Olver† School of Mathematics University of Minnesota Minneapolis, MN 55455 U.S.A. [email protected] http://www.math.umn.edu/∼olver Abstract. A new approach to interpolation theory for functions of several variables is proposed. We develop a multivariate divided difference calculus based on the theory of non-commutative quasi-determinants. In addition, intriguing explicit formulae that connect the classical finite difference interpolation coefficients for univariate curves with multivariate interpolation coefficients for higher dimensional submanifolds are established. † Supported in part by NSF Grant DMS 11–08894. April 6, 2016 1 1. Introduction. Interpolation theory for functions of a single variable has a long and distinguished his- tory, dating back to Newton’s fundamental interpolation formula and the classical calculus of finite differences, [7, 47, 58, 64]. Standard numerical approximations to derivatives and many numerical integration methods for differential equations are based on the finite dif- ference calculus. However, historically, no comparable calculus was developed for functions of more than one variable. If one looks up multivariate interpolation in the classical books, one is essentially restricted to rectangular, or, slightly more generally, separable grids, over which the formulae are a simple adaptation of the univariate divided difference calculus. See [19] for historical details. Starting with G. Birkhoff, [2] (who was, coincidentally, my thesis advisor), recent years have seen a renewed level of interest in multivariate interpolation among both pure and applied researchers; see [18] for a fairly recent survey containing an extensive bibli- ography. De Boor and Ron, [8, 12, 13], and Sauer and Xu, [61, 10, 65], have systemati- cally studied the polynomial case. -
Easy Way to Find Multivariate Interpolation
IJETST- Vol.||04||Issue||05||Pages 5189-5193||May||ISSN 2348-9480 2017 International Journal of Emerging Trends in Science and Technology IC Value: 76.89 (Index Copernicus) Impact Factor: 4.219 DOI: https://dx.doi.org/10.18535/ijetst/v4i5.11 Easy way to Find Multivariate Interpolation Author Yimesgen Mehari Faculty of Natural and Computational Science, Department of Mathematics Adigrat University, Ethiopia Email: [email protected] Abstract We derive explicit interpolation formula using non-singular vandermonde matrix for constructing multi dimensional function which interpolates at a set of distinct abscissas. We also provide examples to show how the formula is used in practices. Introduction Engineers and scientists commonly assume that past and currently. But there is no theoretical relationships between variables in physical difficulty in setting up a frame work for problem can be approximately reproduced from discussing interpolation of multivariate function f data given by the problem. The ultimate goal whose values are known.[1] might be to determine the values at intermediate Interpolation function of more than one variable points, to approximate the integral or to simply has become increasingly important in the past few give a smooth or continuous representation of the years. These days, application ranges over many variables in the problem different field of pure and applied mathematics. Interpolation is the method of estimating unknown For example interpolation finds applications in the values with the help of given set of observations. numerical integrations of differential equations, According to Theile Interpolation is, “The art of topography, and the computer aided geometric reading between the lines of the table” and design of cars, ships, airplanes.[1] According to W.M. -
Spatial Interpolation Methods
Page | 0 of 0 SPATIAL INTERPOLATION METHODS 2018 Page | 1 of 1 1. Introduction Spatial interpolation is the procedure to predict the value of attributes at unobserved points within a study region using existing observations (Waters, 1989). Lam (1983) definition of spatial interpolation is “given a set of spatial data either in the form of discrete points or for subareas, find the function that will best represent the whole surface and that will predict values at points or for other subareas”. Predicting the values of a variable at points outside the region covered by existing observations is called extrapolation (Burrough and McDonnell, 1998). All spatial interpolation methods can be used to generate an extrapolation (Li and Heap 2008). Spatial Interpolation is the process of using points with known values to estimate values at other points. Through Spatial Interpolation, We can estimate the precipitation value at a location with no recorded data by using known precipitation readings at nearby weather stations. Rationale behind spatial interpolation is the observation that points close together in space are more likely to have similar values than points far apart (Tobler’s Law of Geography). Spatial Interpolation covers a variety of method including trend surface models, thiessen polygons, kernel density estimation, inverse distance weighted, splines, and Kriging. Spatial Interpolation requires two basic inputs: · Sample Points · Spatial Interpolation Method Sample Points Sample Points are points with known values. Sample points provide the data necessary for the development of interpolator for spatial interpolation. The number and distribution of sample points can greatly influence the accuracy of spatial interpolation. -
On Sharp Extrapolation Theorems
ON SHARP EXTRAPOLATION THEOREMS by Dariusz Panek M.Sc., Mathematics, Jagiellonian University in Krak¶ow,Poland, 1995 M.Sc., Applied Mathematics, University of New Mexico, USA, 2004 DISSERTATION Submitted in Partial Ful¯llment of the Requirements for the Degree of Doctor of Philosophy Mathematics The University of New Mexico Albuquerque, New Mexico December, 2008 °c 2008, Dariusz Panek iii Acknowledgments I would like ¯rst to express my gratitude for M. Cristina Pereyra, my advisor, for her unconditional and compact support; her unbounded patience and constant inspi- ration. Also, I would like to thank my committee members Dr. Pedro Embid, Dr Dimiter Vassilev, Dr Jens Lorens, and Dr Wilfredo Urbina for their time and positive feedback. iv ON SHARP EXTRAPOLATION THEOREMS by Dariusz Panek ABSTRACT OF DISSERTATION Submitted in Partial Ful¯llment of the Requirements for the Degree of Doctor of Philosophy Mathematics The University of New Mexico Albuquerque, New Mexico December, 2008 ON SHARP EXTRAPOLATION THEOREMS by Dariusz Panek M.Sc., Mathematics, Jagiellonian University in Krak¶ow,Poland, 1995 M.Sc., Applied Mathematics, University of New Mexico, USA, 2004 Ph.D., Mathematics, University of New Mexico, 2008 Abstract Extrapolation is one of the most signi¯cant and powerful properties of the weighted theory. It basically states that an estimate on a weighted Lpo space for a single expo- p nent po ¸ 1 and all weights in the Muckenhoupt class Apo implies a corresponding L estimate for all p; 1 < p < 1; and all weights in Ap. Sharp Extrapolation Theorems track down the dependence on the Ap characteristic of the weight. -
3.2 Rational Function Interpolation and Extrapolation
104 Chapter 3. Interpolation and Extrapolation do 11 i=1,n Here we find the index ns of the closest table entry, dift=abs(x-xa(i)) if (dift.lt.dif) then ns=i dif=dift endif c(i)=ya(i) and initialize the tableau of c’s and d’s. d(i)=ya(i) http://www.nr.com or call 1-800-872-7423 (North America only), or send email to [email protected] (outside North Amer readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, v Permission is granted for internet users to make one paper copy their own personal use. Further reproduction, or any copyin Copyright (C) 1986-1992 by Cambridge University Press. Programs Copyright (C) 1986-1992 by Numerical Recipes Software. Sample page from NUMERICAL RECIPES IN FORTRAN 77: THE ART OF SCIENTIFIC COMPUTING (ISBN 0-521-43064-X) enddo 11 y=ya(ns) This is the initial approximation to y. ns=ns-1 do 13 m=1,n-1 For each column of the tableau, do 12 i=1,n-m we loop over the current c’s and d’s and update them. ho=xa(i)-x hp=xa(i+m)-x w=c(i+1)-d(i) den=ho-hp if(den.eq.0.)pause ’failure in polint’ This error can occur only if two input xa’s are (to within roundoff)identical. den=w/den d(i)=hp*den Here the c’s and d’s are updated. c(i)=ho*den enddo 12 if (2*ns.lt.n-m)then After each column in the tableau is completed, we decide dy=c(ns+1) which correction, c or d, we want to add to our accu- else mulating value of y, i.e., which path to take through dy=d(ns) the tableau—forking up or down. -
Comparison of Gap Interpolation Methodologies for Water Level Time Series Using Perl/Pdl
Revista de Matematica:´ Teor´ıa y Aplicaciones 2005 12(1 & 2) : 157–164 cimpa – ucr – ccss issn: 1409-2433 comparison of gap interpolation methodologies for water level time series using perl/pdl Aimee Mostella∗ Alexey Sadovksi† Scott Duff‡ Patrick Michaud§ Philippe Tissot¶ Carl Steidleyk Received/Recibido: 13 Feb 2004 Abstract Extensive time series of measurements are often essential to evaluate long term changes and averages such as tidal datums and sea level rises. As such, gaps in time series data restrict the type and extent of modeling and research which may be accomplished. The Texas A&M University Corpus Christi Division of Nearshore Research (TAMUCC-DNR) has developed and compared various methods based on forward and backward linear regression to interpolate gaps in time series of water level data. We have developed a software system that retrieves actual and harmonic water level data based upon user provided parameters. The actual water level data is searched for missing data points and the location of these gaps are recorded. Forward and backward linear regression are applied in relation to the location of missing data or gaps in the remaining data. After this process is complete, one of three combinations of the forward and backward regression is used to fit the results. Finally, the harmonic component is added back into the newly supplemented time series and the results are graphed. The software created to implement this process of linear regression is written in Perl along with a Perl module called PDL (Perl Data Language). Generally, this process has demonstrated excellent results in filling gaps in our water level time series. -
Bayesian Interpolation of Unequally Spaced Time Series
Bayesian interpolation of unequally spaced time series Luis E. Nieto-Barajas and Tapen Sinha Department of Statistics and Department of Actuarial Sciences, ITAM, Mexico [email protected] and [email protected] Abstract A comparative analysis of time series is not feasible if the observation times are different. Not even a simple dispersion diagram is possible. In this article we propose a Gaussian process model to interpolate an unequally spaced time series and produce predictions for equally spaced observation times. The dependence between two obser- vations is assumed a function of the time differences. The novelty of the proposal relies on parametrizing the correlation function in terms of Weibull and Log-logistic survival functions. We further allow the correlation to be positive or negative. Inference on the model is made under a Bayesian approach and interpolation is done via the posterior predictive conditional distributions given the closest m observed times. Performance of the model is illustrated via a simulation study as well as with real data sets of temperature and CO2 observed over 800,000 years before the present. Key words: Bayesian inference, EPICA, Gaussian process, kriging, survival functions. 1 Introduction The majority of the literature on time series analysis assumes that the variable of interest, say Xt, is observed on a regular or equally spaced sequence of times, say t 2 f1; 2;:::g (e.g. Box et al. 2004; Chatfield 1989). Often, time series are observed at uneven times. Unequally spaced (also called unevenly or irregularly spaced) time series data occur naturally in many fields. For instance, natural disasters, such as earthquakes, floods and volcanic eruptions, occur at uneven intervals. -
Multivariate Lagrange Interpolation 1 Introduction 2 Polynomial
Alex Lewis Final Project Multivariate Lagrange Interpolation Abstract. Explain how the standard linear Lagrange interpolation can be generalized to construct a formula that interpolates a set of points in . We will also provide examples to show how the formula is used in practice. 1 Introduction Interpolation is a fundamental topic in Numerical Analysis. It is the practice of creating a function that fits a finite sampled data set. These sampled values can be used to construct an interpolant, which must fit the interpolated function at each date point. 2 Polynomial Interpolation “The Problem of determining a polynomial of degree one that passes through the distinct points and is the same as approximating a function for which and by means of a first degree polynomial interpolating, or agreeing with, the values of f at the given points. We define the functions and , And then define We can generalize this concept of linear interpolation and consider a polynomial of degree n that passes through n+1 points. In this case we can construct, for each , a function with the property that when and . To satisfy for each requires that the numerator of contains the term To satisfy , the denominator of must be equal to this term evaluated at . Thus, This is called the Th Lagrange Interpolating Polynomial The polynomial is given once again as, Example 1 If and , then is the polynomial that agrees with at , and ; that is, 3 Multivariate Interpolation First we define a function be an -variable multinomial function of degree . In other words, the , could, in our case be ; in this case it would be a 3 variable function of some degree . -
Interpolation, Extrapolation, Etc
Interpolation, extrapolation, etc. Prof. Ramin Zabih http://cs100r.cs.cornell.edu Administrivia Assignment 5 is out, due Wednesday A6 is likely to involve dogs again Monday sections are now optional – Will be extended office hours/tutorial Prelim 2 is Thursday – Coverage through today – Gurmeet will do a review tomorrow 2 Hillclimbing is usually slow m Lines with Lines with errors < 30 errors < 40 b 3 Extrapolation Suppose you only know the values of the distance function f(1), f(2), f(3), f(4) – What is f(5)? This problem is called extrapolation – Figuring out a value outside the range where you have data – The process is quite prone to error – For the particular case of temporal data, extrapolation is called prediction – If you have a good model, this can work well 4 Interpolation Suppose you only know the values of the distance function f(1), f(2), f(3), f(4) – What is f(2.5)? This problem is called interpolation – Figuring out a value inside the range where you have data • But at a point where you don’t have data – Less error-prone than extrapolation Still requires some kind of model – Otherwise, crazy things could happen between the data points that you know 5 Analyzing noisy temporal data Suppose that your data is temporal – Radar tracking of airplanes, or DJIA There are 3 things you might wish to do – Prediction: given data through time T, figure out what will happen at time T+1 – Estimation: given data through time T, figure out what actually happened at time T • Sometimes called filtering – Smoothing: given data through -
Comparison of Spatial Interpolation Methods for the Estimation of Air Quality Data
Journal of Exposure Analysis and Environmental Epidemiology (2004) 14, 404–415 r 2004 Nature Publishing Group All rights reserved 1053-4245/04/$30.00 www.nature.com/jea Comparison of spatial interpolation methods for the estimation of air quality data DAVID W. WONG,a LESTER YUANb AND SUSAN A. PERLINb aSchool of Computational Sciences, George Mason University, Fairfax, VA, USA bNational Center for Environmental Assessment, Washington, US Environmental Protection Agency, USA We recognized that many health outcomes are associated with air pollution, but in this project launched by the US EPA, the intent was to assess the role of exposure to ambient air pollutants as risk factors only for respiratory effects in children. The NHANES-III database is a valuable resource for assessing children’s respiratory health and certain risk factors, but lacks monitoring data to estimate subjects’ exposures to ambient air pollutants. Since the 1970s, EPA has regularly monitored levels of several ambient air pollutants across the country and these data may be used to estimate NHANES subject’s exposure to ambient air pollutants. The first stage of the project eventually evolved into assessing different estimation methods before adopting the estimates to evaluate respiratory health. Specifically, this paper describes an effort using EPA’s AIRS monitoring data to estimate ozone and PM10 levels at census block groups. We limited those block groups to counties visited by NHANES-III to make the project more manageable and apply four different interpolation methods to the monitoring data to derive air concentration levels. Then we examine method-specific differences in concentration levels and determine conditions under which different methods produce significantly different concentration values. -
3.1 Polynomial Interpolation and Extrapolation
108 Chapter 3. Interpolation and Extrapolation f(x, y, z). Multidimensional interpolation is often accomplished by a sequence of one-dimensional interpolations. We discuss this in 3.6. § CITED REFERENCES AND FURTHER READING: Abramowitz, M., and Stegun, I.A. 1964, Handbook of Mathematical Functions, Applied Mathe- visit website http://www.nr.com or call 1-800-872-7423 (North America only),or send email to [email protected] (outside North America). readable files (including this one) to any servercomputer, is strictly prohibited. To order Numerical Recipes books,diskettes, or CDROMs Permission is granted for internet users to make one paper copy their own personal use. Further reproduction, or any copying of machine- Copyright (C) 1988-1992 by Cambridge University Press.Programs Numerical Recipes Software. Sample page from NUMERICAL RECIPES IN C: THE ART OF SCIENTIFIC COMPUTING (ISBN 0-521-43108-5) matics Series, Volume 55 (Washington: National Bureau of Standards; reprinted 1968 by Dover Publications, New York), 25.2. § Stoer, J., and Bulirsch, R. 1980, Introduction to Numerical Analysis (New York: Springer-Verlag), Chapter 2. Acton, F.S. 1970, Numerical Methods That Work; 1990, corrected edition (Washington: Mathe- matical Association of America), Chapter 3. Kahaner, D., Moler, C., and Nash, S. 1989, Numerical Methods and Software (Englewood Cliffs, NJ: Prentice Hall), Chapter 4. Johnson, L.W., and Riess, R.D. 1982, Numerical Analysis, 2nd ed. (Reading, MA: Addison- Wesley), Chapter 5. Ralston, A., and Rabinowitz, P. 1978, A First Course in Numerical Analysis, 2nd ed. (New York: McGraw-Hill), Chapter 3. Isaacson, E., and Keller, H.B. 1966, Analysis of Numerical Methods (New York: Wiley), Chapter 6.