MATH 590: Meshfree Methods
|
|
- Clare McDaniel
- 5 years ago
- Views:
Transcription
1 MATH 590: Meshfree Methods Chapter 2: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 MATH 590 Chapter 2 1
2 Outline 1 Radial (Basis) Functions 2 Radial Basis Function Interpolation fasshauer@iit.edu MATH 590 Chapter 2 2
3 Outline Radial (Basis) Functions 1 Radial (Basis) Functions 2 Radial Basis Function Interpolation fasshauer@iit.edu MATH 590 Chapter 2 3
4 Next goals Radial (Basis) Functions Get a feel for RBFs before we study them in detail Want to overcome limitations of distance matrix interpolation, but keep overall structure: combine distance matrix interpolants with good basic functions fasshauer@iit.edu MATH 590 Chapter 2 4
5 Radial (Basis) Functions Example The Gaussian ϕ(r) = e (εr)2, r R, has a shape parameter ε related to the variance σ 2 of the normal distribution: ε 2 = 1/(2σ 2 ). fasshauer@iit.edu MATH 590 Chapter 2 5
6 Radial (Basis) Functions Example The Gaussian ϕ(r) = e (εr)2, r R, has a shape parameter ε related to the variance σ 2 of the normal distribution: ε 2 = 1/(2σ 2 ). Compose the Gaussian with Euclidean distance function 2 and get Φ k (x) = e ε2 x x k 2 2, x R s, where x k R s is some fixed center or knot. fasshauer@iit.edu MATH 590 Chapter 2 5
7 Radial (Basis) Functions Example The Gaussian ϕ(r) = e (εr)2, r R, has a shape parameter ε related to the variance σ 2 of the normal distribution: ε 2 = 1/(2σ 2 ). Compose the Gaussian with Euclidean distance function 2 and get Φ k (x) = e ε2 x x k 2 2, x R s, where x k R s is some fixed center or knot. Remark Note that ϕ is univariate, but Φ k is a multivariate function such that Φ k (x) = ϕ( x x k 2 ) radial basis function (RBF). fasshauer@iit.edu MATH 590 Chapter 2 5
8 Radial (Basis) Functions Definition A function Φ : R s R is called radial provided there exists a univariate function ϕ : [0, ) R such that Φ(x) = ϕ(r), where r = x, and is some norm on R s usually the Euclidean norm. fasshauer@iit.edu MATH 590 Chapter 2 6
9 Radial (Basis) Functions Definition A function Φ : R s R is called radial provided there exists a univariate function ϕ : [0, ) R such that Φ(x) = ϕ(r), where r = x, and is some norm on R s usually the Euclidean norm. Thus, for a radial function x 1 = x 2 = Φ(x 1 ) = Φ(x 2 ), x 1, x 2 R s. fasshauer@iit.edu MATH 590 Chapter 2 6
10 Radial (Basis) Functions Definition A function Φ : R s R is called radial provided there exists a univariate function ϕ : [0, ) R such that Φ(x) = ϕ(r), where r = x, and is some norm on R s usually the Euclidean norm. Thus, for a radial function Example x 1 = x 2 = Φ(x 1 ) = Φ(x 2 ), x 1, x 2 R s. The Euclidean distance function Φ(x) = x 2 (or ϕ(r) = r) is a special case of a radial (basis) function. Terminology: ϕ: basic function, Φ k ( 2 ) (centered at x k ): radial basis function fasshauer@iit.edu MATH 590 Chapter 2 6
11 Radial (Basis) Functions Figure: Gaussian with ε = 1 (left) and ε = 3 (right) centered at the origin. A smaller value of ε (i.e., larger variance) causes the function to become flatter, so it s like an inverse length scale. Increasing ε leads to a more peaked RBF. fasshauer@iit.edu MATH 590 Chapter 2 7
12 Radial (Basis) Functions Figure: Gaussian with ε = 1 (left) and ε = 3 (right) centered at the origin. A smaller value of ε (i.e., larger variance) causes the function to become flatter, so it s like an inverse length scale. Increasing ε leads to a more peaked RBF. The choice of ε influences both accuracy and numerical stability of the solution to our interpolation problem (details later). fasshauer@iit.edu MATH 590 Chapter 2 7
13 Radial (Basis) Functions Some properties of radial function interpolants Invariance under all Euclidean transformations (i.e., translations, rotations, and reflections). It doesn t matter whether we first compute the RBF interpolant and then apply a Euclidean transformation, or if we first transform the data and then compute the interpolant. True since Euclidean transformations are characterized by orthogonal transformation matrices and are therefore 2-norm-invariant. Invariance under translation, rotation and reflection is often desirable in applications. fasshauer@iit.edu MATH 590 Chapter 2 8
14 Radial (Basis) Functions Some properties of radial function interpolants Invariance under all Euclidean transformations (i.e., translations, rotations, and reflections). It doesn t matter whether we first compute the RBF interpolant and then apply a Euclidean transformation, or if we first transform the data and then compute the interpolant. True since Euclidean transformations are characterized by orthogonal transformation matrices and are therefore 2-norm-invariant. Invariance under translation, rotation and reflection is often desirable in applications. Insensitive to the dimension s of the space in which the data sites lie. Instead of having to deal with a multivariate function Φ we can work with the same univariate function ϕ for all choices of s. fasshauer@iit.edu MATH 590 Chapter 2 8
15 The Kernel Notation Radial (Basis) Functions Sometimes we use more general (reproducing) kernels K. A kernel is a function of two variables K : R s R s R, K : (x, y) K (x, y) We will discuss the reproducing concept later. fasshauer@iit.edu MATH 590 Chapter 2 9
16 Radial (Basis) Functions Kernels can be radial: (also called isotropic) K (x, y) = κ( x y ), x, y R s, κ : R + 0 R fasshauer@iit.edu MATH 590 Chapter 2 10
17 Radial (Basis) Functions Kernels can be radial: (also called isotropic) K (x, y) = κ( x y ), x, y R s, κ : R + 0 R (only) translation invariant: (also called stationary) K (x, y) = K (x y), x, y R s, K : R s R fasshauer@iit.edu MATH 590 Chapter 2 10
18 Radial (Basis) Functions Kernels can be radial: (also called isotropic) K (x, y) = κ( x y ), x, y R s, κ : R + 0 R (only) translation invariant: (also called stationary) K (x, y) = K (x y), x, y R s, K : R s R some other special form, e.g., zonal: K (x, y) = κ(x y), x, y S 2, κ : [ 1, 1] R anisotropic (but stationary): K (x, y) = s K d (x d y d ), x d, y d R, K d : R R d=1 fasshauer@iit.edu MATH 590 Chapter 2 10
19 Radial (Basis) Functions Kernels can be radial: (also called isotropic) K (x, y) = κ( x y ), x, y R s, κ : R + 0 R (only) translation invariant: (also called stationary) K (x, y) = K (x y), x, y R s, K : R s R some other special form, e.g., zonal: K (x, y) = κ(x y), x, y S 2, κ : [ 1, 1] R anisotropic (but stationary): K (x, y) = s K d (x d y d ), x d, y d R, K d : R R d=1 very general: K (x, y), x, y R s, K : R s R s R A nice paper that discusses different classes of kernels from a statistical perspective is [Genton (2001)]. fasshauer@iit.edu MATH 590 Chapter 2 10
20 Examples Radial (Basis) Functions Radial kernel: K (x, y) = κ( x y ) = e ε2 x y 2, κ(r) = e ε2 r 2 fasshauer@iit.edu MATH 590 Chapter 2 11
21 Examples Radial (Basis) Functions Radial kernel: K (x, y) = κ( x y ) = e ε2 x y 2, κ(r) = e ε2 r 2 Anisotropic translation invariant kernel: s K (x, y) = e ε2 d (x d y d ) 2 = e s d=1 ε2 d (x d y d ) 2 d=1 becomes radial if ε d = ε for all d = 1,..., s fasshauer@iit.edu MATH 590 Chapter 2 11
22 Examples Radial (Basis) Functions Radial kernel: K (x, y) = κ( x y ) = e ε2 x y 2, κ(r) = e ε2 r 2 Anisotropic translation invariant kernel: s K (x, y) = e ε2 d (x d y d ) 2 = e s d=1 ε2 d (x d y d ) 2 d=1 becomes radial if ε d = ε for all d = 1,..., s Zonal kernel: 1 K (x, y) = κ(x y) = 1 + ε 2 2εx y, κ(t) = ε 2 2εt fasshauer@iit.edu MATH 590 Chapter 2 11
23 Examples Radial (Basis) Functions Radial kernel: K (x, y) = κ( x y ) = e ε2 x y 2, κ(r) = e ε2 r 2 Anisotropic translation invariant kernel: s K (x, y) = e ε2 d (x d y d ) 2 = e s d=1 ε2 d (x d y d ) 2 d=1 becomes radial if ε d = ε for all d = 1,..., s Zonal kernel: 1 K (x, y) = κ(x y) = 1 + ε 2 2εx y, κ(t) = ε 2 2εt General kernel: K (x, y) = min{x, y} xy, s = 1 (Brownian bridge) fasshauer@iit.edu MATH 590 Chapter 2 11
24 Outline Radial Basis Function Interpolation 1 Radial (Basis) Functions 2 Radial Basis Function Interpolation fasshauer@iit.edu MATH 590 Chapter 2 12
25 Radial Basis Function Interpolation Instead of distance matrices we now use a radial basis function expansion to solve the scattered data interpolation problem by assuming N P f (x) = c k ϕ ( x x k 2 ), x R s. (1) k=1 fasshauer@iit.edu MATH 590 Chapter 2 13
26 Radial Basis Function Interpolation Instead of distance matrices we now use a radial basis function expansion to solve the scattered data interpolation problem by assuming N P f (x) = c k ϕ ( x x k 2 ), x R s. (1) k=1 Using the interpolation conditions P f (x i ) = f (x i ), i = 1,..., N, we get c k from ϕ ( x 1 x 1 2 ) ϕ ( x 1 x 2 2 )... ϕ ( x 1 x N 2 ) c 1 ϕ ( x 2 x 1 2 ) ϕ ( x 2 x 2 2 )... ϕ ( x 2 x N 2 ) c = ϕ ( x N x 1 2 ) ϕ ( x N x 2 2 )... ϕ ( x N x N 2 ) c N f (x 1 ) f (x 2 ). f (x N ) fasshauer@iit.edu MATH 590 Chapter 2 13
27 Radial Basis Function Interpolation Instead of distance matrices we now use a radial basis function expansion to solve the scattered data interpolation problem by assuming N P f (x) = c k ϕ ( x x k 2 ), x R s. (1) k=1 Using the interpolation conditions P f (x i ) = f (x i ), i = 1,..., N, we get c k from ϕ ( x 1 x 1 2 ) ϕ ( x 1 x 2 2 )... ϕ ( x 1 x N 2 ) c 1 ϕ ( x 2 x 1 2 ) ϕ ( x 2 x 2 2 )... ϕ ( x 2 x N 2 ) c = ϕ ( x N x 1 2 ) ϕ ( x N x 2 2 )... ϕ ( x N x N 2 ) c N f (x 1 ) f (x 2 ). f (x N ) Question: For what type of basic functions ϕ is the system matrix non-singular? fasshauer@iit.edu MATH 590 Chapter 2 13
28 Radial Basis Function Interpolation Setup for Example Use s = 2 Basic function: Gaussians or linear function ϕ(r) = r Code in RBFInterpolation2D.m written for 2D (can easily be generalized to sd), and uses DistanceMatrix.m. Use Franke s function f (x, y) = 3 4 e 1/4((9x 2)2 +(9y 2) 2) e (1/49)(9x+1)2 (1/10)(9y+1) e 1/4((9x 7)2 +(9y 3) 2) 1 5 e (9x 4)2 (9y 7) 2 Use Halton data sites (get others easily from CreatePoints) Compute errors on uniform grid fasshauer@iit.edu MATH 590 Chapter 2 14
29 Radial Basis Function Interpolation Figure: Franke s test function. fasshauer@iit.edu MATH 590 Chapter 2 15
30 Radial Basis Function Interpolation Program (RBFInterpolation2D.m) 1 rbf exp(-(e*r).^2); ep = 21.1; 2 f1 0.75*exp(-((9*x-2).^2+(9*y-2).^2)/4); 3 f2 0.75*exp(-((9*x+1).^2/49+(9*y+1).^2/10)); 4 f3 0.5*exp(-((9*x-7).^2+(9*y-3).^2)/4); 5 f4 0.2*exp(-((9*x-4).^2+(9*y-7).^2)); 6 testfunction f1(x,y)+f2(x,y)+f3(x,y)-f4(x,y); 7 N = 1089; gridtype = h ; 8 dsites = CreatePoints(N,2,gridtype); 9 ctrs = dsites; 10 neval = 40; M = neval^2; 11 epoints = CreatePoints(M,2, u ); 12 rhs = testfunction(dsites(:,1),dsites(:,2)); 13 DM_data = DistanceMatrix(dsites,ctrs); 14 IM = rbf(ep,dm_data); 15 DM_eval = DistanceMatrix(epoints,ctrs); 16 EM = rbf(ep,dm_eval); 17 Pf = EM * (IM\rhs); 18 exact = testfunction(epoints(:,1),epoints(:,2)); 19 maxerr = norm(pf-exact,inf) 20 rms_err = norm(pf-exact)/neval fasshauer@iit.edu MATH 590 Chapter 2 16
31 Radial Basis Function Interpolation Stationary and Non-Stationary Approximation Remark Note the conflicting use of the term stationary: here vs. statistics = translation-invariant fasshauer@iit.edu MATH 590 Chapter 2 17
32 Radial Basis Function Interpolation Stationary and Non-Stationary Approximation Remark Note the conflicting use of the term stationary: here vs. statistics = translation-invariant We usually verify convergence of a method numerically by looking at a sequence of experiments with increasingly larger sets of data. fasshauer@iit.edu MATH 590 Chapter 2 17
33 Radial Basis Function Interpolation Stationary and Non-Stationary Approximation Remark Note the conflicting use of the term stationary: here vs. statistics = translation-invariant We usually verify convergence of a method numerically by looking at a sequence of experiments with increasingly larger sets of data. Non-stationary approximation: Use one fixed value of ε for all of the experiments. fasshauer@iit.edu MATH 590 Chapter 2 17
34 Radial Basis Function Interpolation Stationary and Non-Stationary Approximation Remark Note the conflicting use of the term stationary: here vs. statistics = translation-invariant We usually verify convergence of a method numerically by looking at a sequence of experiments with increasingly larger sets of data. Non-stationary approximation: Use one fixed value of ε for all of the experiments. Stationary approximation: Scale the shape parameter ε according to the fill distance (or meshsize) h so that we end up using peaked basis functions for densely spaced data and flat basis functions for coarsely spaced data. More details later fasshauer@iit.edu MATH 590 Chapter 2 17
35 Fill distance Radial Basis Function Interpolation The fill distance (or covering radius) is usually defined as h = h X,Ω = sup min x x j 2. (2) x j X x Ω fasshauer@iit.edu MATH 590 Chapter 2 18
36 Fill distance Radial Basis Function Interpolation The fill distance (or covering radius) is usually defined as h = h X,Ω = sup min x x j 2. (2) x j X x Ω It indicates how well the data in the set X fill out the domain Ω. fasshauer@iit.edu MATH 590 Chapter 2 18
37 Fill distance Radial Basis Function Interpolation The fill distance (or covering radius) is usually defined as h = h X,Ω = sup min x x j 2. (2) x j X x Ω It indicates how well the data in the set X fill out the domain Ω. In MATLAB: hx = max(min(dm_eval )) (3) fasshauer@iit.edu MATH 590 Chapter 2 18
38 Fill distance Radial Basis Function Interpolation The fill distance (or covering radius) is usually defined as h = h X,Ω = sup min x x j 2. (2) x j X x Ω It indicates how well the data in the set X fill out the domain Ω. In MATLAB: hx = max(min(dm_eval )) (3) Remark Since min works along columns of a matrix, transposition of the non-symmetric evaluation matrix corresponds to finding for each evaluation point the distance to the corresponding closest data site, and then setting h X,Ω as the worst of those distances. fasshauer@iit.edu MATH 590 Chapter 2 18
39 Radial Basis Function Interpolation y h X x Figure: The fill distance for N = 25 Halton points (h X,Ω ). It s the radius of the largest possible empty ball that can be placed among the data locations inside Ω. fasshauer@iit.edu MATH 590 Chapter 2 19
40 Radial Basis Function Interpolation Test of non-stationary interpolation Fix large ε to prevent severe ill-conditioning with large N (if ε = 1, then N = 25 already very ill-conditioned) Consequence: very localized basis functions that don t capture enough information on small point sets (see left plot) Figure: Gaussian RBF interpolant with ε = 21.1 at N = 289 (left) and at N = 1089 Halton points (right). fasshauer@iit.edu MATH 590 Chapter 2 20
41 Radial Basis Function Interpolation Gaussian Distance matrix k N RMS-error max-error RMS-error max-error e e e e e e e e e e e e e e e e e e e e e e e e-003 Table: Non-stationary RBF interpolation to Franke s function using Gaussians (ε = 21.1) and Euclidean distance matrices. fasshauer@iit.edu MATH 590 Chapter 2 21
42 Radial Basis Function Interpolation Gaussian Distance matrix k N RMS-error max-error RMS-error max-error e e e e e e e e e e e e e e e e e e e e e e e e-003 Table: Non-stationary RBF interpolation to Franke s function using Gaussians (ε = 21.1) and Euclidean distance matrices. Remark Surprising observation in above example: distance matrix fit more accurate than Gaussian. fasshauer@iit.edu MATH 590 Chapter 2 21
43 Radial Basis Function Interpolation Gaussian Distance matrix k N RMS-error max-error RMS-error max-error e e e e e e e e e e e e e e e e e e e e e e e e-003 Table: Non-stationary RBF interpolation to Franke s function using Gaussians (ε = 21.1) and Euclidean distance matrices. Remark Surprising observation in above example: distance matrix fit more accurate than Gaussian. Would expect Gaussian to be more accurate find a better ε fasshauer@iit.edu MATH 590 Chapter 2 21
44 Radial Basis Function Interpolation Same test function and Gaussians as before. Figure: Maximum (blue) and RMS (red) errors vs. ε for 81 (top left), 289 (top right), 1089 (bottom left), and 4225 Halton points (bottom right). MATH 590 Chapter 2 22
45 Radial Basis Function Interpolation What do we learn from the ε error curves? Errors decrease with decreasing ε (not the same as convergence for h 0). fasshauer@iit.edu MATH 590 Chapter 2 23
46 Radial Basis Function Interpolation What do we learn from the ε error curves? Errors decrease with decreasing ε (not the same as convergence for h 0). Error curves are not monotonic. fasshauer@iit.edu MATH 590 Chapter 2 23
47 Radial Basis Function Interpolation What do we learn from the ε error curves? Errors decrease with decreasing ε (not the same as convergence for h 0). Error curves are not monotonic. Optimal value of ε seems to exist which minimizes errors. fasshauer@iit.edu MATH 590 Chapter 2 23
48 Radial Basis Function Interpolation What do we learn from the ε error curves? Errors decrease with decreasing ε (not the same as convergence for h 0). Error curves are not monotonic. Optimal value of ε seems to exist which minimizes errors. For small enough values of ε the computational results become unpredictable, and the error curves become erratic (ill-conditioning) fasshauer@iit.edu MATH 590 Chapter 2 23
49 Radial Basis Function Interpolation What do we learn from the ε error curves? Errors decrease with decreasing ε (not the same as convergence for h 0). Error curves are not monotonic. Optimal value of ε seems to exist which minimizes errors. For small enough values of ε the computational results become unpredictable, and the error curves become erratic (ill-conditioning) We ll consider ε to be safe if we don t get a MATLAB warning about ill-conditioning. Note that for small N optimal ε is safe, but for larger N not. We obtain highly accurate solutions from severely ill-conditioned linear systems! This is known as the uncertainty or trade-off principle. fasshauer@iit.edu MATH 590 Chapter 2 23
50 Radial Basis Function Interpolation What do we learn from the ε error curves? Errors decrease with decreasing ε (not the same as convergence for h 0). Error curves are not monotonic. Optimal value of ε seems to exist which minimizes errors. For small enough values of ε the computational results become unpredictable, and the error curves become erratic (ill-conditioning) We ll consider ε to be safe if we don t get a MATLAB warning about ill-conditioning. Note that for small N optimal ε is safe, but for larger N not. We obtain highly accurate solutions from severely ill-conditioned linear systems! This is known as the uncertainty or trade-off principle. Interesting problem: how to compute optimal solution in a stable way. One possible approach is [Fornberg & Wright (2004)] with details later, also ongoing work with Mike McCourt (see next slide). fasshauer@iit.edu MATH 590 Chapter 2 23
51 Radial Basis Function Interpolation Figure: Preliminary results for stable Gaussian evaluation. QR corresponds to new method, Direct to standard method. MATH 590 Chapter 2 24
52 Radial Basis Function Interpolation Best possible errors for Gaussian interpolation smallest safe ε smallest RMS-error k N ε RMS-error max-error ε RMS-error max-error e e e e e e e e e e e e e e e e e e e e e e e e-007 Table: Optimal RBF interpolation to Franke s function using Gaussians. fasshauer@iit.edu MATH 590 Chapter 2 25
53 Radial Basis Function Interpolation Best possible errors for Gaussian interpolation smallest safe ε smallest RMS-error k N ε RMS-error max-error ε RMS-error max-error e e e e e e e e e e e e e e e e e e e e e e e e-007 Table: Optimal RBF interpolation to Franke s function using Gaussians. Remark Errors for safe ε now comparable (or smaller) than those for distance matrix interpolation. fasshauer@iit.edu MATH 590 Chapter 2 25
54 Radial Basis Function Interpolation Best possible errors for Gaussian interpolation smallest safe ε smallest RMS-error k N ε RMS-error max-error ε RMS-error max-error e e e e e e e e e e e e e e e e e e e e e e e e-007 Table: Optimal RBF interpolation to Franke s function using Gaussians. Remark Errors for safe ε now comparable (or smaller) than those for distance matrix interpolation. Much better for optimal ε. However, ε = 6.2 with N = 1089 Halton points yields RCOND = e-020. fasshauer@iit.edu MATH 590 Chapter 2 25
55 Radial Basis Function Interpolation Summary Remark If the data are not sampled from a known test function, then we will not be able to choose an optimal shape parameter by monitoring the RMS error. fasshauer@iit.edu MATH 590 Chapter 2 26
56 Radial Basis Function Interpolation Summary Remark If the data are not sampled from a known test function, then we will not be able to choose an optimal shape parameter by monitoring the RMS error. New challenge: how to find optimal ε based only on the data? fasshauer@iit.edu MATH 590 Chapter 2 26
57 Radial Basis Function Interpolation Summary Remark If the data are not sampled from a known test function, then we will not be able to choose an optimal shape parameter by monitoring the RMS error. New challenge: how to find optimal ε based only on the data? We will study later. ill-conditioning and preconditioning, optimal shape parameter selection, and alternate stable evaluation methods via a Contour-Padé algorithm or eigenfunction expansions fasshauer@iit.edu MATH 590 Chapter 2 26
58 References I Appendix References Buhmann, M. D. (2003). Radial Basis Functions: Theory and Implementations. Cambridge University Press. Fasshauer, G. E. (2007). Meshfree Approximation Methods with MATLAB. World Scientific Publishers. Higham, D. J. and Higham, N. J. (2005). MATLAB Guide. SIAM (2nd ed.), Philadelphia. Iske, A. (2004). Multiresolution Methods in Scattered Data Modelling. Lecture Notes in Computational Science and Engineering 37, Springer Verlag (Berlin). Wendland, H. (2005a). Scattered Data Approximation. Cambridge University Press (Cambridge). fasshauer@iit.edu MATH 590 Chapter 2 27
59 References II Appendix References Fornberg, B. and Wright, G. (2004). Stable computation of multiquadric interpolants for all values of the shape parameter. Comput. Math. Appl. 47, pp Genton, M. C. (2001). Classes of kernels for machine learning: a statistics perspective. J. Machine Learning Research 2, pp MATLAB Central File Exchange. available online at fasshauer@iit.edu MATH 590 Chapter 2 28
MATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 1 Part 3: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 6: Scattered Data Interpolation with Polynomial Precision Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 37: RBF Hermite Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 33: Adaptive Iteration Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 33 1 Outline 1 A
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 39: Non-Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 14: The Power Function and Native Space Error Estimates Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 40: Symmetric RBF Collocation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 5: Completely Monotone and Multiply Monotone Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 9: Conditionally Positive Definite Radial Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 1 Part 2: Scattered Data Interpolation in R d Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 7: Conditionally Positive Definite Functions Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2 Part 3: Native Space for Positive Definite Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH
More informationKernel-based Approximation. Methods using MATLAB. Gregory Fasshauer. Interdisciplinary Mathematical Sciences. Michael McCourt.
SINGAPORE SHANGHAI Vol TAIPEI - Interdisciplinary Mathematical Sciences 19 Kernel-based Approximation Methods using MATLAB Gregory Fasshauer Illinois Institute of Technology, USA Michael McCourt University
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 43: RBF-PS Methods in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 43 1 Outline
More informationMultivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness
Multivariate Interpolation with Increasingly Flat Radial Basis Functions of Finite Smoothness Guohui Song John Riddle Gregory E. Fasshauer Fred J. Hickernell Abstract In this paper, we consider multivariate
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 34: Improving the Condition Number of the Interpolation Matrix Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu
More informationRadial Basis Functions I
Radial Basis Functions I Tom Lyche Centre of Mathematics for Applications, Department of Informatics, University of Oslo November 14, 2008 Today Reformulation of natural cubic spline interpolation Scattered
More informationA new stable basis for radial basis function interpolation
A new stable basis for radial basis function interpolation Stefano De Marchi and Gabriele Santin Department of Mathematics University of Padua (Italy) Abstract It is well-known that radial basis function
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 43: RBF-PS Methods in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590 Chapter 43 1 Outline
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More informationMeshfree Approximation Methods with MATLAB
Interdisciplinary Mathematical Sc Meshfree Approximation Methods with MATLAB Gregory E. Fasshauer Illinois Institute of Technology, USA Y f? World Scientific NEW JERSEY LONDON SINGAPORE BEIJING SHANGHAI
More informationStable Parameterization Schemes for Gaussians
Stable Parameterization Schemes for Gaussians Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado Denver ICOSAHOM 2014 Salt Lake City June 24, 2014 michael.mccourt@ucdenver.edu
More informationStability of Kernel Based Interpolation
Stability of Kernel Based Interpolation Stefano De Marchi Department of Computer Science, University of Verona (Italy) Robert Schaback Institut für Numerische und Angewandte Mathematik, University of Göttingen
More informationRadial basis functions topics in slides
Radial basis functions topics in 40 +1 slides Stefano De Marchi Department of Mathematics Tullio Levi-Civita University of Padova Napoli, 22nd March 2018 Outline 1 From splines to RBF 2 Error estimates,
More informationGreen s Functions: Taking Another Look at Kernel Approximation, Radial Basis Functions and Splines
Green s Functions: Taking Another Look at Kernel Approximation, Radial Basis Functions and Splines Gregory E. Fasshauer Abstract The theories for radial basis functions (RBFs) as well as piecewise polynomial
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods The Connection to Green s Kernels Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction
More informationRBF Collocation Methods and Pseudospectral Methods
RBF Collocation Methods and Pseudospectral Methods G. E. Fasshauer Draft: November 19, 24 Abstract We show how the collocation framework that is prevalent in the radial basis function literature can be
More informationPositive Definite Kernels: Opportunities and Challenges
Positive Definite Kernels: Opportunities and Challenges Michael McCourt Department of Mathematical and Statistical Sciences University of Colorado, Denver CUNY Mathematics Seminar CUNY Graduate College
More informationSolving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels
Solving the 3D Laplace Equation by Meshless Collocation via Harmonic Kernels Y.C. Hon and R. Schaback April 9, Abstract This paper solves the Laplace equation u = on domains Ω R 3 by meshless collocation
More informationLeast Squares Approximation
Chapter 6 Least Squares Approximation As we saw in Chapter 5 we can interpret radial basis function interpolation as a constrained optimization problem. We now take this point of view again, but start
More informationInterpolation by Basis Functions of Different Scales and Shapes
Interpolation by Basis Functions of Different Scales and Shapes M. Bozzini, L. Lenarduzzi, M. Rossini and R. Schaback Abstract Under very mild additional assumptions, translates of conditionally positive
More informationDIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION
Meshless Methods in Science and Engineering - An International Conference Porto, 22 DIRECT ERROR BOUNDS FOR SYMMETRIC RBF COLLOCATION Robert Schaback Institut für Numerische und Angewandte Mathematik (NAM)
More informationToward Approximate Moving Least Squares Approximation with Irregularly Spaced Centers
Toward Approximate Moving Least Squares Approximation with Irregularly Spaced Centers Gregory E. Fasshauer Department of Applied Mathematics, Illinois Institute of Technology, Chicago, IL 6066, U.S.A.
More informationScattered Data Approximation of Noisy Data via Iterated Moving Least Squares
Scattered Data Approximation o Noisy Data via Iterated Moving Least Squares Gregory E. Fasshauer and Jack G. Zhang Abstract. In this paper we ocus on two methods or multivariate approximation problems
More informationStability constants for kernel-based interpolation processes
Dipartimento di Informatica Università degli Studi di Verona Rapporto di ricerca Research report 59 Stability constants for kernel-based interpolation processes Stefano De Marchi Robert Schaback Dipartimento
More informationRecent Results for Moving Least Squares Approximation
Recent Results for Moving Least Squares Approximation Gregory E. Fasshauer and Jack G. Zhang Abstract. We describe two experiments recently conducted with the approximate moving least squares (MLS) approximation
More informationTheoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions
Theoretical and computational aspects of multivariate interpolation with increasingly flat radial basis functions Elisabeth Larsson Bengt Fornberg June 0, 003 Abstract Multivariate interpolation of smooth
More informationA Posteriori Error Bounds for Meshless Methods
A Posteriori Error Bounds for Meshless Methods Abstract R. Schaback, Göttingen 1 We show how to provide safe a posteriori error bounds for numerical solutions of well-posed operator equations using kernel
More informationNumerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical analysis of heat conduction problems on 3D general-shaped domains by means of a RBF Collocation Meshless Method To cite this article: R
More informationApproximation by Conditionally Positive Definite Functions with Finitely Many Centers
Approximation by Conditionally Positive Definite Functions with Finitely Many Centers Jungho Yoon Abstract. The theory of interpolation by using conditionally positive definite function provides optimal
More informationSTABLE CALCULATION OF GAUSSIAN-BASED RBF-FD STENCILS
STABLE CALCULATION OF GAUSSIAN-BASED RBF-FD STENCILS BENGT FORNBERG, ERIK LEHTO, AND COLLIN POWELL Abstract. Traditional finite difference (FD) methods are designed to be exact for low degree polynomials.
More informationRBF-FD Approximation to Solve Poisson Equation in 3D
RBF-FD Approximation to Solve Poisson Equation in 3D Jagadeeswaran.R March 14, 2014 1 / 28 Overview Problem Setup Generalized finite difference method. Uses numerical differentiations generated by Gaussian
More informationKernel-Based Contrast Functions for Sufficient Dimension Reduction
Kernel-Based Contrast Functions for Sufficient Dimension Reduction Michael I. Jordan Departments of Statistics and EECS University of California, Berkeley Joint work with Kenji Fukumizu and Francis Bach
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter V: Least Squares Problems Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2011 fasshauer@iit.edu MATH
More informationNumerical analysis of heat conduction problems on irregular domains by means of a collocation meshless method
Journal of Physics: Conference Series PAPER OPEN ACCESS Numerical analysis of heat conduction problems on irregular domains by means of a collocation meshless method To cite this article: R Zamolo and
More informationSlide05 Haykin Chapter 5: Radial-Basis Function Networks
Slide5 Haykin Chapter 5: Radial-Basis Function Networks CPSC 636-6 Instructor: Yoonsuck Choe Spring Learning in MLP Supervised learning in multilayer perceptrons: Recursive technique of stochastic approximation,
More informationRKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets Class 22, 2004 Tomaso Poggio and Sayan Mukherjee
RKHS, Mercer s theorem, Unbounded domains, Frames and Wavelets 9.520 Class 22, 2004 Tomaso Poggio and Sayan Mukherjee About this class Goal To introduce an alternate perspective of RKHS via integral operators
More information1. Introduction. A radial basis function (RBF) interpolant of multivariate data (x k, y k ), k = 1, 2,..., n takes the form
A NEW CLASS OF OSCILLATORY RADIAL BASIS FUNCTIONS BENGT FORNBERG, ELISABETH LARSSON, AND GRADY WRIGHT Abstract Radial basis functions RBFs form a primary tool for multivariate interpolation, and they are
More informationCIS 520: Machine Learning Oct 09, Kernel Methods
CIS 520: Machine Learning Oct 09, 207 Kernel Methods Lecturer: Shivani Agarwal Disclaimer: These notes are designed to be a supplement to the lecture They may or may not cover all the material discussed
More informationReproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Differential Operators
Reproducing Kernels of Generalized Sobolev Spaces via a Green Function Approach with Differential Operators Qi Ye Abstract In this paper we introduce a generalization of the classical L 2 ( )-based Sobolev
More informationON DIMENSION-INDEPENDENT RATES OF CONVERGENCE FOR FUNCTION APPROXIMATION WITH GAUSSIAN KERNELS
ON DIMENSION-INDEPENDENT RATES OF CONVERGENCE FOR FUNCTION APPROXIMATION WITH GAUSSIAN KERNELS GREGORY E. FASSHAUER, FRED J. HICKERNELL, AND HENRYK WOŹNIAKOWSKI Abstract. This article studies the problem
More informationA orthonormal basis for Radial Basis Function approximation
A orthonormal basis for Radial Basis Function approximation 9th ISAAC Congress Krakow, August 5-9, 2013 Gabriele Santin, joint work with S. De Marchi Department of Mathematics. Doctoral School in Mathematical
More informationScattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions
Chapter 3 Scattered Data Interpolation with Polynomial Precision and Conditionally Positive Definite Functions 3.1 Scattered Data Interpolation with Polynomial Precision Sometimes the assumption on the
More informationNumerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method
Numerical solution of nonlinear sine-gordon equation with local RBF-based finite difference collocation method Y. Azari Keywords: Local RBF-based finite difference (LRBF-FD), Global RBF collocation, sine-gordon
More informationIn the Name of God. Lectures 15&16: Radial Basis Function Networks
1 In the Name of God Lectures 15&16: Radial Basis Function Networks Some Historical Notes Learning is equivalent to finding a surface in a multidimensional space that provides a best fit to the training
More informationApplied Linear Algebra in Geoscience Using MATLAB
Applied Linear Algebra in Geoscience Using MATLAB Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional Plots Programming in
More informationA scaling perspective on accuracy and convergence in RBF approximations
A scaling perspective on accuracy and convergence in RBF approximations Elisabeth Larsson With: Bengt Fornberg, Erik Lehto, and Natasha Flyer Division of Scientific Computing Department of Information
More informationD. Shepard, Shepard functions, late 1960s (application, surface modelling)
Chapter 1 Introduction 1.1 History and Outline Originally, the motivation for the basic meshfree approximation methods (radial basis functions and moving least squares methods) came from applications in
More informationAN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION
J. KSIAM Vol.19, No.4, 409 416, 2015 http://dx.doi.org/10.12941/jksiam.2015.19.409 AN ELEMENTARY PROOF OF THE OPTIMAL RECOVERY OF THE THIN PLATE SPLINE RADIAL BASIS FUNCTION MORAN KIM 1 AND CHOHONG MIN
More informationMultiscale RBF collocation for solving PDEs on spheres
Multiscale RBF collocation for solving PDEs on spheres Q. T. Le Gia I. H. Sloan H. Wendland Received: date / Accepted: date Abstract In this paper, we discuss multiscale radial basis function collocation
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Accuracy and Optimality of RKHS Methods Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2014 fasshauer@iit.edu MATH 590 1 Outline 1 Introduction
More informationRichard F. Bass Krzysztof Burdzy University of Washington
ON DOMAIN MONOTONICITY OF THE NEUMANN HEAT KERNEL Richard F. Bass Krzysztof Burdzy University of Washington Abstract. Some examples are given of convex domains for which domain monotonicity of the Neumann
More informationRBF-Partition of Unity Method: an overview of recent results
RBF-Partition of Unity Method: an overview of recent results Alessandra De Rossi a In collaboration with Roberto Cavoretto a, Emma Perracchione b a Department of Mathematics G. Peano University of Torino
More informationBreaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points
Breaking Computational Barriers: Multi-GPU High-Order RBF Kernel Problems with Millions of Points Michael Griebel Christian Rieger Peter Zaspel Institute for Numerical Simulation Rheinische Friedrich-Wilhelms-Universität
More informationMATH 350: Introduction to Computational Mathematics
MATH 350: Introduction to Computational Mathematics Chapter IV: Locating Roots of Equations Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Spring 2011 fasshauer@iit.edu
More informationA numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices
A numerical study of a technique for shifting eigenvalues of radial basis function differentiation matrices Scott A. Sarra Marshall University and Alfa R.H. Heryudono University of Massachusetts Dartmouth
More informationIntroduction to Gaussian Processes
Introduction to Gaussian Processes Neil D. Lawrence GPSS 10th June 2013 Book Rasmussen and Williams (2006) Outline The Gaussian Density Covariance from Basis Functions Basis Function Representations Constructing
More informationAMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning
AMS526: Numerical Analysis I (Numerical Linear Algebra) Lecture 23: GMRES and Other Krylov Subspace Methods; Preconditioning Xiangmin Jiao SUNY Stony Brook Xiangmin Jiao Numerical Analysis I 1 / 18 Outline
More informationLinear Models for Regression. Sargur Srihari
Linear Models for Regression Sargur srihari@cedar.buffalo.edu 1 Topics in Linear Regression What is regression? Polynomial Curve Fitting with Scalar input Linear Basis Function Models Maximum Likelihood
More informationStrictly Positive Definite Functions on a Real Inner Product Space
Strictly Positive Definite Functions on a Real Inner Product Space Allan Pinkus Abstract. If ft) = a kt k converges for all t IR with all coefficients a k 0, then the function f< x, y >) is positive definite
More informationSTABLE EVALUATION OF GAUSSIAN RBF INTERPOLANTS
STABLE EVALUATION OF GAUSSIAN RBF INTERPOLANTS GREGORY E. FASSHAUER AND MICHAEL J. MCCOURT Abstract. We provide a new way to compute and evaluate Gaussian radial basis function interpolants in a stable
More informationKernel Method: Data Analysis with Positive Definite Kernels
Kernel Method: Data Analysis with Positive Definite Kernels 2. Positive Definite Kernel and Reproducing Kernel Hilbert Space Kenji Fukumizu The Institute of Statistical Mathematics. Graduate University
More informationCOMS 4721: Machine Learning for Data Science Lecture 19, 4/6/2017
COMS 4721: Machine Learning for Data Science Lecture 19, 4/6/2017 Prof. John Paisley Department of Electrical Engineering & Data Science Institute Columbia University PRINCIPAL COMPONENT ANALYSIS DIMENSIONALITY
More informationA STABLE ALGORITHM FOR DIVERGENCE-FREE AND CURL-FREE RADIAL BASIS FUNCTIONS IN THE FLAT LIMIT
A STABLE ALGORITHM FOR DIVERGENCE-FREE AND CURL-FREE RADIAL BASIS FUNCTIONS IN THE FLAT LIMIT by Kathryn Primrose Drake A thesis submitted in partial fulfillment of the requirements for the degree of Master
More informationKey words. Radial basis function, scattered data interpolation, hierarchical matrices, datasparse approximation, adaptive cross approximation
HIERARCHICAL MATRIX APPROXIMATION FOR KERNEL-BASED SCATTERED DATA INTERPOLATION ARMIN ISKE, SABINE LE BORNE, AND MICHAEL WENDE Abstract. Scattered data interpolation by radial kernel functions leads to
More informationNumerical Methods I Non-Square and Sparse Linear Systems
Numerical Methods I Non-Square and Sparse Linear Systems Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 September 25th, 2014 A. Donev (Courant
More informationSolving Boundary Value Problems (with Gaussians)
What is a boundary value problem? Solving Boundary Value Problems (with Gaussians) Definition A differential equation with constraints on the boundary Michael McCourt Division Argonne National Laboratory
More informationFitting Linear Statistical Models to Data by Least Squares I: Introduction
Fitting Linear Statistical Models to Data by Least Squares I: Introduction Brian R. Hunt and C. David Levermore University of Maryland, College Park Math 420: Mathematical Modeling February 5, 2014 version
More informationNumerical solution of surface PDEs with Radial Basis Functions
Numerical solution of surface PDEs with Radial Basis Functions Andriy Sokolov Institut für Angewandte Mathematik (LS3) TU Dortmund andriy.sokolov@math.tu-dortmund.de TU Dortmund June 1, 2017 Outline 1
More informationWendland Functions A C++ code to compute them
Wendland Functions A C++ code to compute them Carlos Argáez 1, Sigurdur Hafstein 1 and Peter Giesl 2 1 Faculty of Physical Sciences, University of Iceland, 107 Reykjavík, Iceland 2 Department of Mathematics,
More information5.6 Nonparametric Logistic Regression
5.6 onparametric Logistic Regression Dmitri Dranishnikov University of Florida Statistical Learning onparametric Logistic Regression onparametric? Doesnt mean that there are no parameters. Just means that
More informationy(x) = x w + ε(x), (1)
Linear regression We are ready to consider our first machine-learning problem: linear regression. Suppose that e are interested in the values of a function y(x): R d R, here x is a d-dimensional vector-valued
More informationKernel Methods. Foundations of Data Analysis. Torsten Möller. Möller/Mori 1
Kernel Methods Foundations of Data Analysis Torsten Möller Möller/Mori 1 Reading Chapter 6 of Pattern Recognition and Machine Learning by Bishop Chapter 12 of The Elements of Statistical Learning by Hastie,
More information10-701/ Recitation : Kernels
10-701/15-781 Recitation : Kernels Manojit Nandi February 27, 2014 Outline Mathematical Theory Banach Space and Hilbert Spaces Kernels Commonly Used Kernels Kernel Theory One Weird Kernel Trick Representer
More informationFast and stable rational RBF-based Partition of Unity interpolation
Fast and stable rational RBF-based Partition of Unity interpolation S. De Marchi, A. Martínez, E. Perracchione Dipartimento di Matematica Tullio Levi-Civita, Università di Padova, Italia Abstract We perform
More informationJoint distribution optimal transportation for domain adaptation
Joint distribution optimal transportation for domain adaptation Changhuang Wan Mechanical and Aerospace Engineering Department The Ohio State University March 8 th, 2018 Joint distribution optimal transportation
More informationStability and Lebesgue constants in RBF interpolation
constants in RBF Stefano 1 1 Dept. of Computer Science, University of Verona http://www.sci.univr.it/~demarchi Göttingen, 20 September 2008 Good Outline Good Good Stability is very important in numerical
More informationNumerical cubature on scattered data by radial basis functions
Numerical cubature on scattered data by radial basis functions A. Sommariva, Sydney, and M. Vianello, Padua September 5, 25 Abstract We study cubature formulas on relatively small scattered samples in
More informationCompact Matérn Kernels and Piecewise Polynomial Splines Viewed From a Hilbert-Schmidt Perspective
Noname manuscript No (will be inserted by the editor) Compact Matérn Kernels and Piecewise Polynomial Splines Viewed From a Hilbert-Schmidt Perspective Roberto Cavoretto Gregory E Fasshauer Michael McCourt
More informationComplexity and regularization issues in kernel-based learning
Complexity and regularization issues in kernel-based learning Marcello Sanguineti Department of Communications, Computer, and System Sciences (DIST) University of Genoa - Via Opera Pia 13, 16145 Genova,
More informationLinear Regression (9/11/13)
STA561: Probabilistic machine learning Linear Regression (9/11/13) Lecturer: Barbara Engelhardt Scribes: Zachary Abzug, Mike Gloudemans, Zhuosheng Gu, Zhao Song 1 Why use linear regression? Figure 1: Scatter
More informationStable computations with flat radial basis functions using vector-valued rational approximations
Stable computations with flat radial basis functions using vector-valued rational approximations Grady B. Wright a,, Bengt Fornberg b a Department of Mathematics, Boise State University, Boise, ID 83725-1555,
More informationSupport Vector Machines.
Support Vector Machines www.cs.wisc.edu/~dpage 1 Goals for the lecture you should understand the following concepts the margin slack variables the linear support vector machine nonlinear SVMs the kernel
More informationThis is a postprint version of the following published document:
This is a postprint version of the following published document: Gonzalez-Rodriguez, P., Moscoso, M., & Kindelan, M. (05). Laurent expansion of the inverse of perturbed, singular matrices. Journal of Computational
More informationLinear Models for Regression
Linear Models for Regression Machine Learning Torsten Möller Möller/Mori 1 Reading Chapter 3 of Pattern Recognition and Machine Learning by Bishop Chapter 3+5+6+7 of The Elements of Statistical Learning
More informationNumerical Minimization of Potential Energies on Specific Manifolds
Numerical Minimization of Potential Energies on Specific Manifolds SIAM Conference on Applied Linear Algebra 2012 22 June 2012, Valencia Manuel Gra f 1 1 Chemnitz University of Technology, Germany, supported
More informationarxiv: v1 [math.na] 17 Oct 2016
Stable computations with flat radial basis functions using vector-valued rational approximations Grady B. Wright a,, Bengt Fornberg b a Department of Mathematics, Boise State University, Boise, ID 83725-1555,
More informationA MODIFIED TSVD METHOD FOR DISCRETE ILL-POSED PROBLEMS
A MODIFIED TSVD METHOD FOR DISCRETE ILL-POSED PROBLEMS SILVIA NOSCHESE AND LOTHAR REICHEL Abstract. Truncated singular value decomposition (TSVD) is a popular method for solving linear discrete ill-posed
More information