Advances in Numerical Analysis
Volume 2016 (2016), Article ID 1397849, 11 pages
http://dx.doi.org/10.1155/2016/1397849

Research Article

Selection of an Interval for Variable Shape Parameter in Approximation by Radial Basis Functions

Jafar Biazar and Mohammad Hosami

Department of Applied Mathematics, Faculty of Mathematical Sciences, University of Guilan, P.O. Box 41335-1914, Rasht 41938 22697, Iran

Received 29 April 2016; Revised 12 July 2016; Accepted 20 July 2016

Academic Editor: Zhangxin Chen

Copyright © 2016 Jafar Biazar and Mohammad Hosami. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In radial basis function approximation, the shape parameter can be variable. The values of the variable shape parameter strategies are selected from an interval which is usually determined by trial and error. As yet there is not any algorithm for determining an appropriate interval, although there are some recipes for optimal values. In this paper, a novel algorithm for determining an interval is proposed. Different variable shape parameter strategies are examined. The results show that the determined interval significantly improved the accuracy and is suitable enough to count on in variable shape parameter strategies.

1. Introduction

Radial basis function methods are the tools for interpolating a multivariate data set, approximating a function, and solving partial differential equations [1–5]. A radial basis function, say , is a continuous univariate function that has been realized by composition with the Euclidean norm in . Radial basis functions (RBFs) may contain a parameter called the shape parameter. Most well-known RBFs have a shape parameter, such as Multiquadric (MQ), Inverse Multiquadric (IMQ), and Gaussian parameters. A well-known RBF without a shape parameter is thin-plate spline (TPS). This class of RBFs approximates with polynomial convergence rate, but three mentioned RBFs approximate with spectral convergence rate. This point leads to applying RBFs with a shape parameter, even though we know that the value of the shape parameter is an issue. In these radial basis functions, the value of the shape parameter plays an important role for the accuracy of the procedure. An important problem is determining the optimal value for the shape parameter. Several articles have been dedicated to introducing different algorithms to compute a constant value, as an optimal or a good value, for the shape parameter [6–15]. Kansa and Carlson showed that instead of a fixed shape parameter variable shape parameters are also useful [16]. They distribute some values in an interval to use them in variable shape parameter. Moreover, there are some articles applying the variable shape parameter in various problems [17–20]. But in many of them there is not any criterion to select the interval of the shape parameters. In fact the selected interval to distribute the shape parameters is user-defined as the selection of the value in constant shape parameter.

In this paper, we focus on the selection of a valid interval in variable shape parameter. An algorithm is suggested to obtain this interval. To show the efficiency of the proposed interval, it is investigated in different variable shape parameter strategies.

2. Shape Parameter in MQ-RBF Approximation

In the RBF methods, a function is approximated by as in the following:

By interpolation conditions

a system of linear equations is obtained as the matrix form , where and If the Multiquadric (MQ) radial basis function is used, the method is called MQ-RBF approximation. In this paper, the MQ is applied throughout the methods and examples.

In MQ-RBF approximation, is presented as the following form:

where and is the shape parameter. Sometimes, to show the dependence of in the accuracy of , the approximate function is written as

Now consider a boundary value problem defined by

where is a linear differential operator and is a linear boundary operator. is the domain of the problem and is the boundary of the domain. Suppose that are a set of distinct centre nodes in the domain such that are interior centre nodes and are boundary centres.

By substituting the approximation of given by RBF approximation (1), problem (4) can be written as

By substituting the centres and in (5) and (6), respectively, the system of equations is obtained, where

And the entries are defined as follows:

By solving the system of equations , the unknown variables are determined.

2.1. Variable Shape Parameter Strategies

The shape parameter in radial basis function approximation is usually a predefined constant parameter. It is acknowledged that in many cases using variable shape parameter strategies results in more accurate approximation [16, 17, 19]. A variable shape parameter strategy uses different values of the shape parameter at different centre points:

where is the shape parameter corresponding to the data point . Applying interpolation conditions (2) to results in the following system:

where in , the matrix of interpolation, each column has the same shape parameter, corresponding to each centre. In MQ-RBF, the interpolation matrix is as follows:

Kansa and Hon [21] showed that the coefficient matrices in global radial basis function approximation tend to become progressively more ill-conditioned, but the variable shape parameter can improve the conditioning of the coefficient matrix. In fact, using variable shape parameter leads to more distinct entries in the RBF matrices which in turn lead to lower condition numbers.

Although variable shape parameter strategy is very accurate, the coefficient matrix is no longer symmetric. This is the cost that should be paid for improving the accuracy‎.

In variable shape parameter, should be varied in a predetermined interval . There are different approaches to distribute the values of the variable shape parameter over an interval. Kansa and Carlson presented some variable shape parameter strategies with different shape parameters at different centres [16]. Their distributions of parameters are in the forms of increasing linearly, decreasing linearly, and exponentially varying shape parameter, respectively, as follows:

Their empirical results demonstrated that in comparison with the constant shape parameter the variable shape parameter can significantly reduce the RMS errors of MQ interpolations. They had shown that variable form of the MQ improves the accuracy of the interpolant in the interior regions where the fitted function varies relatively rapidly [16].

Sarra and Sturgill introduced a random variable shape parameter strategy in interpolations and two-dimensional linear elliptic boundary value problems as follows [19]:

where is a Matlab function that returns matrix of pseudorandom values induced from a uniform distribution on the unit interval. They showed that the random variable shape parameter produces the most accurate results, if the centres are uniformly spaced.

Recently, some other strategies have been applied to generalized Multiquadric RBF, . In [20] trigonometric variable shape parameter and exponent are presented as follows:

Another strategy for variable shape parameter has been introduced by Sanyasiraju and Satyanarayana [18]. They introduced a meshless scheme with a variable shape parameter for steady convection-diffusion equation. They minimized an error function to generate a variable shape parameter for each centre node.

3. Selection of an Interval in Variable Shape Parameter

In addition to selecting a suitable strategy, there is another important problem in variable shape parameter approach: distribution of the values in an appropriate interval. In [20], by using the classical variable MQ as , an interval such as is introduced, where and (the shape parameter of classical form of the MQ is shown by ). is the total number of nodes in one-dimensional approximations and is the node number of a row or a column in two-dimensional approximations. Sarra and Sturgill emphasize the importance of the length of the interval, , they take , and they study the errors by trial and errors, while and vary in a meaningful range [19]. Kansa and Carlson [16] consider a set of shape parameters which minimizes an error function over ‎some evaluation points. To find these values, the following objective function is minimized:‎

where

and are evaluation points to compute the RMS error. They used two approaches to find the optimal set for the shape parameters. In the first approach, a nonlinear system of normal equations is obtained by least square approach:

In the second one, the sum of squares over the evaluation points is minimized. For more details refer to [16].

In [19, 20], it is seen that the proposed values for and are empirical. In Kansa’s approach a set of optimal shape parameters has been obtained, but the shortage of this approach is that the exact function is needed. Kansa stated that the purpose of his report was some empirical observations to developments in the theory of the variable shape parameters. Therefore his approach is not applicable in practice. In this paper, we apply a new algorithm to find a reliable interval for the shape parameters.

3.1. Proposed Algorithm

In this section the algorithm is introduced.

Suppose that approximates the function , over a domain . In Section 2, it was noted that the accuracy of depends on . For all valid values of , keeping fixed, the corresponding values of should be convergent to . Therefore

If the curve is the plot of as a function of the shape parameter versus the shape parameter, then will be horizontal over an interval, say . This interval is a suitable range for the valid values of the shape parameter. According to this idea, the proposed algorithm to find an appropriate interval for the shape parameter can be summarized as follows.

Algorithm

(1)Let .(2)For to .(3)Plot the curve of versus . Do evaluate the value of .(4)Select the subinterval , over which the curve is almost parallel to axis.

The interval or a subset of this interval can be applied in the variable shape parameter strategies as . The approximate function is recomputed by using the proposed interval.

4. Numerical Examples

In this section to demonstrate the reliability and the efficiency of the algorithm in various types of variable shape parameter strategies, three examples, for interpolation in one- and two-dimensional spaces, are presented. Also two boundary value problems in two-dimensional spaces are investigated. Through all the examples we use the uniform centre nodes. To validate the numerical accuracy, two error functions such as Max error and RMS error are calculated:

where are some evaluation points.

In random variable shape parameter strategy, due to the nature of random selection, the results may be different. In all of the numerical results, to reduce this impact, the average error of ten time successive runs is provided, when the random variable shape parameter strategy is used.

4.1. One- and Two-Dimensional Interpolation

The test functions in one-dimensional interpolation are

First, the algorithm is applied to find an interval for each test function and then the approximate function is recomputed by using the variable shape parameter strategy in the proposed interval. To compare the results, the obtained accuracy is compared with the accuracy of approximate functions computed by using variable shape parameter in different intervals.

Figures 1(a), 1(b), and 1(c) show the curve of versus shape parameter in the interval for three different , using 60 centre nodes. The curve is almost parallel to axis, over the interval . The interval can be selected as , in variable shape parameter. By rerunning the “plot command” on subinterval , a more accurate interval is obtained by Figures 1(d), 1(e), and 1(f). By considering and , the error functions in interpolation of by constant, exponential, linear, and random shape parameter strategies are plotted in Figure 3(a). 40 evaluation points and 60 centre points are applied. In this figure, the impact of different distribution strategies is shown. To illustrate the influence of the interval , which is the aim of this study, the results of applying the proposed interval for two strategies (random and exponential) are shown in Figure 4.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 1: The plots of for , , and (a, b, and c) versus the shape parameters and (d, e, and f).

The plots in Figure 4 show that the proposed interval obtained by the algorithm is reliable and accurate for variable shape parameter strategies. By selecting the shape parameters in the proposed range the accuracy significantly improved.

By visual inspection, in Figures 2(a), 2(b), and 2(c), as a function of the shape parameter has the same value in the interval . By rerunning the plots in , Figures 2(d), 2(e), and 2(f) demonstrate that the suitable interval can be selected as . Using the interval in variable shape parameter, the error functions obtained by constant, exponential, linear, and random shape parameter strategies are shown in Figure 3(b). To demonstrate the reliability of this interval, the error functions obtained by exponential and random variable shape parameter in the intervals , , , and are plotted in Figure 4. The errors demonstrate that using the proposed interval , obtained by the proposed algorithm, results in better accuracy.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: The plots of for , , and (a, b, and c) versus the shape parameters and (d, e, and f).

(a)

(b)

Figure 3: The error functions for approximation of (a) and (b) by constant, exponential, linear, and random shape parameters in the proposed intervals.

(a)

(b)

(c)

(d)

Figure 4: The error functions by random (a, c) and exponential (b, d) variable shape parameters in different intervals for approximating (a, b) and (c, d).

It is noted that, for various test points , the optimum range may be a little different. The experimental results show that the proposed interval may not be the optimal one to approximate in the whole domain, but it is a valid interval. Therefore, in next examples, the results of only one test point are plotted.

In two-dimensional interpolation consider the test function

in the rectangle . Using 144 centre points, the algorithm is applied for an arbitrary point . Figure 5(a) shows the curve obtained in step of the algorithm. The plot demonstrates that the appropriate interval can be chosen as . The error function of approximate function by using linear variable shape parameter in this interval is shown in Figure 6. Also some numerical results are summarized in Table 1.

Table 1: The comparison of the RMS error and Max error for different intervals.

⁢
Exponential
Random

RMS error
Max error
RMS error
Max error

Proposed interval
0.5, 1.5
4.79e − 6
2.32e − 5
7.15e − 7
1.72e − 6

Some other intervals
1, 2
1.00e − 4
4.61e − 4
2.89e − 5
6.84e − 5

3, 4
4.10e − 4
2.29e − 3
1.56e − 3
5.44e − 3

4, 5
6.81e − 4
3.74e − 3
6.35e − 4
3.32e − 3

Problem (12)
Proposed interval
0.4, 1
2.69e − 5
7.82e − 5
5.13e − 5
9.22e − 5

Some other intervals
1, 2
1.50e − 3
2.84e − 3
9.28e − 4
2.02e − 3

2, 3
2.01e − 3
3.57e − 3
4.69e − 3
9.07e − 3

4, 5
4.93e − 3
7.58e − 3
4.73e − 3
7.02e − 3

Problem (13)
Proposed interval
0.2, 1
2.40e − 9
6.86e − 9
2.25e − 12
5.09e − 12

Some other intervals
1, 2
1.17e − 4
2.04e − 4
1.04e − 6
1.76e − 6

2, 3
1.67e − 4
3.03e − 4
7.23e − 4
2.56e − 3

3, 4
4.18e − 4
6.70e − 4
1.17e − 3
2.25e − 3

(a)

(b)

(c)

Figure 5: Plot of (a) and the solution of problems (12) and (13), (resp., (b) and (c)), versus the shape parameter .

(a)

(b)

Figure 6: The error function of the approximation of by linear variable shape parameter in the intervals (a) and (b).

4.2. Boundary Value Problems

To illustrate the accuracy and the efficiency of the algorithm, two linear boundary value problems are considered. First consider the following BVP:

where and is the boundary of the domain. The exact solution is . Using 100 centre points, the algorithm is applied to approximate the solution of problem (24). The plot of , for a predetermined , versus shape parameter is shown in Figure 5(b). The curve determines the appropriate interval . The comparison of the error functions by applying this interval to various shape parameter strategies with those obtained by some other intervals is shown in Figure 7 and Table 1. Numerical results demonstrated that one can achieve good numerical results by choosing the interval determined by the algorithm.

(a)

(b)

Figure 7: The error function of the approximation of problem (24) by random variable shape parameters in intervals (a) and (b).

The algorithm is also performed for the following problem:

where and are specified such that the exact solution is

Using 100 centre points and an arbitrary point , Figure 5(c) illustrates that the interval can be selected for applying in the variable shape parameter strategies. The absolute error of the approximate solution by applying this interval, which has been shown in Figure 8, indicates that the proposed interval significantly increases the accuracy of the results. The results of Table 1 illustrate the impact of applying a suitable interval in different variable shape parameter strategies.

(a)

(b)

Figure 8: The error function of the approximation of problem (25) by random variable shape parameter in intervals (a) and (b).

5. Conclusion

In this paper, a new algorithm is suggested to determine an interval for variable shape parameter strategies. By distributing the values of the shape parameters in the proposed interval, it was shown that the new numerical results will be more accurate. To show the efficiency of the proposed interval for different strategies, some famous variable shape parameter strategies were investigated. The results demonstrate that in one- and two-dimensional interpolation problems and in two-dimensional boundary value problems the algorithm works well.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of the paper.

References

G. E. Fasshauer, Meshfree Approximation Methods with MATLAB, vol. 6, World Scientific, River Edge, NJ, USA, 2007.

A. J. M. Ferreira, L. M. S. Castro, and S. Bertoluzza, “A high order collocation method for the static and vibration analysis of composite plates using a first-order theory,” Composite Structures, vol. 89, no. 3, pp. 424–432, 2009.

Q. Shen, “A meshless method of lines for the numerical solution of KdV equation using radial basis functions,” Engineering Analysis with Boundary Elements, vol. 33, no. 10, pp. 1171–1180, 2009.

J. Lin, W. Chen, and K. Y. Sze, “A new radial basis function for Helmholtz problems,” Engineering Analysis with Boundary Elements, vol. 36, no. 12, pp. 1923–1930, 2012.

Z.-J. Fu, W. Chen, and L. Ling, “Method of approximate particular solutions for constant- and variable-order fractional diffusion models,” Engineering Analysis with Boundary Elements, vol. 57, pp. 37–46, 2015.

R. E. Carlson and T. A. Foley, “The parameter R2 in multiquadric interpolation,” Computers and Mathematics with Applications, vol. 21, no. 9, pp. 29–42, 1991.

S. Rippa, “An algorithm for selecting a good value for the parameter c in radial basis function interpolation,” Advances in Computational Mathematics, vol. 11, no. 2-3, pp. 193–210, 1999.

V. Bayona, M. Moscoso, and M. Kindelan, “Optimal constant shape parameter for multiquadric based RBF-FD method,” Journal of Computational Physics, vol. 230, no. 19, pp. 7384–7399, 2011.

A. H. D. Cheng, “Multiquadric and its shape parameter—a numerical investigation of error estimate, condition number, and round-off error by arbitrary precision computation,” Engineering Analysis with Boundary Elements, vol. 36, no. 2, pp. 220–239, 2012.

G. E. Fasshauer and J. G. Zhang, “On choosing ‘optimal’ shape parameters for RBF approximation,” Numerical Algorithms, vol. 45, no. 1–4, pp. 345–368, 2007.

B. Fornberg and G. Wright, “Stable computation of multiquadric interpolants for all values of the shape parameter,” Computers & Mathematics with Applications, vol. 48, no. 5-6, pp. 853–867, 2004.

C.-S. Huang, C.-F. Lee, and A. H.-D. Cheng, “Error estimate, optimal shape factor, and high precision computation of multiquadric collocation method,” Engineering Analysis with Boundary Elements, vol. 31, no. 7, pp. 614–623, 2007.

C. M. Roque and A. J. Ferreira, “Numerical experiments on optimal shape parameters for radial basis functions,” Numerical Methods for Partial Differential Equations, vol. 26, no. 3, pp. 675–689, 2010.

C. H. Tsai, J. Kolibal, and M. Li, “The golden section search algorithm for finding a good shape parameter for meshless collocation methods,” Engineering Analysis with Boundary Elements, vol. 34, no. 8, pp. 738–746, 2010.

M. Uddin, “On the selection of a good value of shape parameter in solving time-dependent partial differential equations using RBF approximation method,” Applied Mathematical Modelling, vol. 38, no. 1, pp. 135–144, 2014.

E. J. Kansa and R. E. Carlson, “Improved accuracy of multiquadric interpolation using variable shape parameters,” Computers & Mathematics with Applications, vol. 24, no. 12, pp. 99–120, 1992.

B. Fornberg and J. Zuev, “The Runge phenomenon and spatially variable shape parameters in RBF interpolation,” Computers & Mathematics with Applications, vol. 54, no. 3, pp. 379–398, 2007.

Y. Sanyasiraju and C. Satyanarayana, “Local Hermite-RBF based grid-free scheme with a variable (optimal) shape parameter for steady convection-diffusion equations,” International Journal of Numerical Analysis and Modeling, Series B, vol. 4, pp. 382–393, 2013.

S. A. Sarra and D. Sturgill, “A random variable shape parameter strategy for radial basis function approximation methods,” Engineering Analysis with Boundary Elements, vol. 33, no. 11, pp. 1239–1245, 2009.

S. Xiang, K.-M. Wang, Y.-T. Ai, Y.-D. Sha, and H. Shi, “Trigonometric variable shape parameter and exponent strategy for generalized multiquadric radial basis function approximation,” Applied Mathematical Modelling, vol. 36, no. 5, pp. 1931–1938, 2012.

E. J. Kansa and Y. C. Hon, “Circumventing the Ill-conditioning problem with multiquadric radial basis functions: applications to elliptic partial differential equations,” Computers & Mathematics with Applications, vol. 39, no. 7-8, pp. 123–137, 2000.