کد مقاله | کد نشریه | سال انتشار | مقاله انگلیسی | نسخه تمام متن |
---|---|---|---|---|
4640562 | 1341279 | 2010 | 7 صفحه PDF | دانلود رایگان |
![عکس صفحه اول مقاله: Error saturation in Gaussian radial basis functions on a finite interval Error saturation in Gaussian radial basis functions on a finite interval](/preview/png/4640562.png)
Radial basis function (RBF) interpolation is a “meshless” strategy with great promise for adaptive approximation. One restriction is “error saturation” which occurs for many types of RBFs including Gaussian RBFs of the form ϕ(x;α,h)=exp(−α2(x/h)2)ϕ(x;α,h)=exp(−α2(x/h)2): in the limit h→0h→0 for fixed αα, the error does not converge to zero, but rather to ES(α)ES(α). Previous studies have theoretically determined the saturation error for Gaussian RBF on an infinite, uniform interval and for the same with a single point omitted. (The gap enormously increases ES(α)ES(α).) We show experimentally that the saturation error on the unit interval, x∈[−1,1]x∈[−1,1], is about 0.06exp(−0.47/α2)‖f‖∞0.06exp(−0.47/α2)‖f‖∞ — huge compared to the O(2π/α2)exp(−π2/[4α2])O(2π/α2)exp(−π2/[4α2]) saturation error for a grid with one point omitted. We show that the reason the saturation is so large on a finite interval is that it is equivalent to an infinite grid which is uniform except for a gap of many points. The saturation error can be avoided by choosing α≪1α≪1, the “flat limit”, but the condition number of the interpolation matrix explodes as O(exp(π2/[4α2]))O(exp(π2/[4α2])). The best strategy is to choose the largest αα which yields an acceptably small saturation error: If the user chooses an error tolerance δδ, then αoptimum(δ)=1/−2log(δ/0.06).
Journal: Journal of Computational and Applied Mathematics - Volume 234, Issue 5, 1 July 2010, Pages 1435–1441