Research Article  Open Access
Guanbin Gao, Huaishan Zhang, Xing Wu, Yu Guo, "Structural Parameter Identification of Articulated Arm Coordinate Measuring Machines", Mathematical Problems in Engineering, vol. 2016, Article ID 4063046, 10 pages, 2016. https://doi.org/10.1155/2016/4063046
Structural Parameter Identification of Articulated Arm Coordinate Measuring Machines
Abstract
Precise structural parameter identification of a robotic articulated arm coordinate measuring machine (AACMM) is essential for improving its measuring accuracy, particularly in robotic applications. This paper presents a constructive parameter identification approach for robotic AACMMs. We first develop a mathematical kinematic model of the AACMM based on the DenavitHartenberg (DH) approach established for robotic systems. This model is further calibrated and verified via the practical test data. Based on the difference between the calculated coordinates of the AACMM probe via the kinematic model and the given reference coordinates, a parameter identification approach is proposed to estimate the structural parameters in terms of the test data set. The Jacobian matrix is further analyzed to determine the solvability of the identification model. It shows that there are two coupling parameters, which can be removed in the regressor. Finally, a parameter identification algorithm taking the leastsquare solution of the identification model as the structural parameters by using the obtained poses data is suggested. Practical experiments based on a robotic AACMM test rig are carried out, and the results reveal the effectiveness and robustness of the proposed identification approach.
1. Introduction
The coordinate measuring machine (CMM) is a universal measuring instrument which can transform various geometric measurements into coordinate measurements [1]. This instrument has been widely used in the calibration and modeling of robotic systems. In particular, the articulated arm CMM (AACMM) is a new type of robotlike CMM with multiple degrees of freedom (DOF), which generally consists of a series of linkages connected by joints in series [2, 3]. The AACMM obtains the angles of joints by means of the angle encoder installed on the rotary joints. And the angles can be transformed into the three coordinates through the kinematic model. The AACMM possesses some specific and essential characteristics and advantages, for example, simple mechanical structure, small size, light weight, large measurement range, and flexible measurement in field [4].
However, the measuring accuracy of the AACMM is much lower than that of the orthogonal CMM [5, 6], which may greatly limit its applications. A potential strategy to improve the measuring accuracy of the AACMM is to choose high precision hardware components and to apply high requirements for manufacturing and assembly [3, 7]. However, the cost of the measuring machine increases dramatically by using more precise components. Moreover, in some specific applications, increasing the precision of individual components may not be able to ensure overall increased measuring precision and accuracy. In particular, small errors in some kinematic parameters may accumulate and thus influence the measuring accuracy of the AACMM greatly.
Another essential and economic way to eliminate errors of the structural parameters and to improve the measuring accuracy is to identify the robot’s structural parameters [8] by using appropriate parameter identification approaches. Generally, the structural parameter identification includes four steps [9]: (1) modeling: to establish a mathematical model describing the geometrical characteristics and kinematics of the robot; (2) measurement: to measure the coordinates of the end effector in the real world coordinate system; (3) identification: to identify the structural parameters of the obtained model by means of mathematical calculation of the data set; (4) compensation: to modify the parameters in the control system according to the identification results. These four steps are also applicable for the AACMM, and an improved identification method will be studied in this paper.
For the parameter identification of AACMM, Kovač and Frank [10] developed a new high precision device for the AACMM testing and calibration with the laser interferometer measurements along a line gauge beam. Santolaria et al. [11, 12] reported a method to calibrate an AACMM based on the DenavitHartenberg (DH) kinematic model parameters. These parameters are optimized by measuring a calibrated ball bar gauge located at different orientations and positions in the AACMM working space. Hamana et al. [13] presented a method, where the kinematic parameters of AACMM were calibrated using spherical center coordinates as the artifact. However, only part of the measuring space can be calibrated with the above methods. Thus, the AACMM cannot be calibrated by directly using these available results. In particular, the coupling relationships between the structural parameters and their effect on the measuring uncertainty of the AACMM were not considered in the above researches. Therefore, the robustness and efficiency of the identifications were affected due to those invalid calculations [14] of the redundant couplings.
In this paper, we propose an improved modeling and parameter identification method for AACMM robotic system. First, the kinematic model and structural identification matrix were established based on the DH method, and the coupling relationship between the structural parameters was obtained through further analysis of the structural identification matrix. Then the identification model of the AACMM was constructed, and a parameter estimation approach developed based on the LS method is proposed to identify the structural parameters. Practically collected data of the joint angles and coordinates of the probe are used to validate the model and identification approach. The redundancy embedded in the parameter matrix is further analyzed and eliminated to address the coupling effects and identifiability. Finally, practical experiments are conducted to verify the efficiency of the proposed identification method.
The advantages and the distinctive features of this proposed identification method in comparison to some other identification methods for AACMM (e.g., [12, 15, 16]) are as follows:(1)We do not need precise initial parameters, and even we do not need initial parameters (we can assign the initial parameters arbitrarily as long as they are not too exaggerated) in the identification. The identified values of the structural parameters can be solved through (18). However, in some available results, for example, [12, 15, 17], the initial identification parameters should be appropriately selected to achieve good identification results because an iteration calculation method is adopted.(2)In this paper, we conducted the coupling analysis such that those linearly dependant parameters are detected and removed from the parameters to be identified. Consequently, the calculation costs and the identification efficacy can be significantly improved. In fact, in our case study, only one time iteration calculation can provide fairly good results.(3)The time consumed by the proposed identification calculation is relatively shorter than the widely used iteration identification methods such as PSO, GA, and LS; that is, we can get the results just after one time iteration calculation.
The paper is organized as follows. Section 2 presents the kinematic modeling and validation; the parameter identification and the analysis are introduced in Section 3; and experimental results are given in Section 4. Section 5 provides conclusions.
2. Kinematic Modeling and Verification
2.1. Kinematic Modeling
As shown in Figure 1, the structure of the AACMM is similar to an articulated robot. Therefore, the AACMM model can be established by using available modeling methods developed for robotics. For the modeling of robotic kinematics, the most influential method is the DenavitHartenberg model (DH model) which has been widely used due to its clear physical meaning [18]. A homogeneous transform matrix is used to represent the spatial relations of adjacent joints coordinate systems [19]. Because all the adjacent joints of the AACMM are perpendicular [20], there is no nominally parallel problem, and we can use the DH method to establish the kinematic model.
In this paper, a 6DOF AACMM is studied. The schematic structure of this system is shown in Figure 1. For this 6DOF AACMM, the coordinate values of the probe in the reference coordinate system can be derived through 7 successive transformations.
According to the DH method, there are four groups of structural parameters in the AACMM, for example, linkage length , joint length , torsion angle , and joint angle .
The coordinates of the AACMM were established followed by the DH method, as shown in Figure 2. Some of the four groups of structural parameters were achieved via direct measurement and some of them can only be determined by using estimation method because they cannot be measured directly, for example, , as shown in Table 1. The measurement datum of the structural parameters is all axes of the joints and linkages, which are virtual lines only, and the positions of them are not accurate. Therefore, the structural parameters are not accurate as they will result in significant movement or measurement uncertainty of the AACMM.
 
mm. 
According to the principle of homogeneous transformation, the transform from the coordinate system to is equivalent to a process in which conducts rotation and translation and then completely coincides with . The whole transform procedure can be expressed mathematically as
The probe coordinates in the reference coordinate system can be expressed as follows:
There is one group of parameters which are the variables in (2). It is the joint angle , the value of which can be acquired by the joint angular sensor.
2.2. Verification of the Kinematic Model
To validate the proposed kinematic model (2), we collect realistic operation data via sensors, for example, joint angle of the AACMM, which are then used to calculate the coordinate system. For this purpose, a software according to the kinematic model (2) was developed to acquire the joint angles, which can also be used to verify the kinematic model of the AACMM. The interface of the built software with C++ is shown in Figure 3. As we find all the structural parameters can be imported into the software through the interface. Then with the joint angles acquired from the AACMM and the structural parameters, the coordinates of the AACMM probe can be calculated by the software.
By comparing the coordinates of the probe calculated by the software and the reference coordinates (which can be treated as the true values of the coordinates of the probe), we found that the calculated coordinates of the probe are very close to the reference coordinates. The differences between the two group coordinates are always less than 1.5 mm, which means that the movement uncertainty of the AACMM is about 1.5 mm. In this sense, the kinematic model (2) of the AACMM is correct because a wrong kinematic model will lead to uncontrollable movement uncertainty. However, it should be noted that although such movement uncertainty is acceptable for industrial robots, it is not applicable for AACMM as this is a much more precision instrument. Thus further analysis will be conducted in the following Section 3 to identify the exact structural parameters by introducing a discretetime parameter identification method.
3. Identification for AACMM Uncertain Parameters
3.1. Modeling Error Analysis
It is known that the accuracy of the AACMM is affected greatly by the manufacture, assembly, and component selection [21, 22], which can all result in inconsistence between the structural parameters in the kinematic model and the actual ones. However, the proposed identification method in this paper can achieve high accuracy for AACMM modeling without increasing the cost to construct the parts and buy extra sensors. This section will further address the improvement of accuracy by analyzing the error dynamics and then introducing an associated identification scheme.
From (2), we know that the uncertainty of the AACMM is affected by the error of structural parameters including the linkage length error , the torsion angle error , the joint length error , the probe length error , and the joint angle deviation at the zero position [23]. Thus the effects of these uncertainties should be identified and then compensated. For this purpose, (2) can be rewritten as a function group of the structural parameters, which is shown in
To show the effect of uncertainties, we calculate the differential of (3) as
Thus (4) can be rewritten in a compact form as given in
It is shown that (5) is in a linear form with perturbing modeling uncertainties. Since we consider sufficiently small sampling interval in the discretetime implementation, the structural parameters of the AACMM are relatively small; that is, and . Then (5) can be rewritten as
We denote as the Jacobian matrix, as the difference between the coordinates of the probe calculated by the kinematic model and the measured reference coordinates, and as the uncertainties perturbing the structure parameters. Then (6) can also be simplified aswhere denotes the difference of the regressor within each iteration step .
Equation (10) is in a strictly linearly parameterized form within each iteration step, and thus it can be used for the purpose of structure parameter identification to estimate .
3.2. Identification of Structural Parameters
From (6), we know that (10) can be used to solve the error of structural parameters . It should be noted that there are 25 unknown parameters in (10), for example, , , , , and . However, there are only 3 measurements (e.g., ). Therefore, we need to collect at least 9 groups of reference coordinates and poses (joint angles) in various operation regimes to creates equations, which are then used to identify these 25 unknown parameters. To get more robust results, we need even more than 9 groups of reference coordinates and joint angles. Thus, for , we can write augmented (10) as
It is clear that the uncertain structural parameters can be calculated based on (11) when the regressor or matrix is nonsingular. This is the wellknown persistent excitation condition and can be fulfilled in our case study by operating the robotic in extensive scenarios. To solve in (11) conveniently, (11) can be further rewritten in a more compact formwhere and denote the augmented regressor and measurement matrices.
3.3. Analysis of Couplings of the Structural Parameters
In (12), if the Jacobian matrix is not a column full rank matrix, that is, there are linear correlation structural parameters, we cannot calculate its solution directly. Therefore, we should check if there are linearly correlated structural parameters in . Based on the matrix theory, the linearly related row can be found by applying singular value decomposition elementary row transformation of via the orthogonal matrix decomposition [24, 25]. Premultiplying on both sides of (12), we can obtain
We denote ; then by using the matrix singular value decomposition, we havewhere and are orthogonal matrices of , and r is the rank of the matrix and the Jacobian matrix . Therefore, the number of linearly related parameters is . From (13) and (14), we can obtain
Since is a symmetric matrix, , and are rotation matrices, then is equivalent to . In this case, the linearly related structural parameters in can be found by the elementary row transformation of the last 4 lines, as shown in where is a matrix of , is a vector of , is the matrix of , and is the rank of . The calculation result shows that the rank of is 23; that is, ; thus there are two linearly related structural parameters. Specifically, we can obtain linearly related parameters by carrying out the singular value decomposition and the elementary row transform. As shown in (17), we find that is linearly related with and is linearly related with . Thus in the identification calculation, and are treated as the redundant parameters and they do not need to be identified. Also, and can be determined as the redundant parameters and are not identified in the identification calculation, which can also reduce the dimension of the regressor matrices and thus the required calculation costs.
3.4. Implementation of Identification Calculation
According to the analysis presented in last section, we know that there are two structural parameters which are related to others and do not need to be identified in terms of (12). Instead, they can be calculated based on (17) provided other 23 parameters are identified. Therefore, the two corresponding columns of the Jacobian matrix can be removed and the new identification equations can be obtained as shown in where denotes the change between each iteration interval , which is the vector of the structural parameters of the AACMM after the removal of and , and is the Jacobian matrix after removal of the two columns of and .
The overall identification algorithm has been carried out to solve (18) in Matlab, and the flowchart shown in Figure 4 provides the practical implementation of the suggested algorithm.
In the literature, there are many different identification methods, for example, PSO and GA. However, these algorithms need much more computational costs, which may limit their applicability for AACMM. The time consumed by the proposed identification calculation is relatively shorter than the widely used iteration identification methods such as PSO and GA; that is, we can get the results just after one time iteration calculation.
The stepbystep implementation procedure can be given as follows:(1)initialize the structural parameters of with the data in Table 1 and set iteration variable as ;(2)calculate and with (2), (5), , and the poses acquired from the AACMM;(3)calculate with the reference coordinates;(4)calculate the least squares solution of ;(5)set and correct the parameters as by using the composition .
4. Experiments
To show the effectiveness of the suggested modeling and identification algorithms, an experimental study was carried out based on an AACMM test rig, which is shown in Figure 5. The nominal point repeatability and volumetric accuracy of the AACMM are 0.05 mm and 0.06 mm, respectively. However, the structural parameters and working principle of the AACMM are not clear for the users. This is the motivation for our current work, which is dedicated to estimate the parameters to achieve an identification model of the AACMM. We conducted our experiments on the normal room temperature (in the range of 17°C to 23°C) by a same operator. On the other hand, it should be noted that the main aim of this paper is to provide a parameter identification method for the AACMM and other similar link mechanisms. Thus, this paper does not consider the temperature compensation and calibration test according to ASME B.89.4.22, VDI/VDE 2617 Part 9 or ISO 10360 Part 12.
In the data acquisition procedure of the AACMM, many poses are used to obtain enough information about the AACMM. The information is very important for the following identification. And acquiring poses is commonly used in the calibration of AACMM and robots [15–17]. In this paper, we not only obtain many poses but also get the information of the structural parameters except measuring them directly. We operate the AACMM robot in extensive scenarios, and the joint angles and reference coordinates of 200 poses were acquired to test the uncertainty of the AACMM. The differences between the coordinates calculated joint angles based on the kinematic model and the reference coordinates are the movement uncertainty of the AACMM, which is shown in Figure 6. The maximum movement uncertainty is 1.485 mm, 1.366 mm, and −0.978 mm in the directions of , , and , respectively, and the average of the absolute movement uncertainties are 0.524 mm, 0.545 mm, and 0.421 mm for , , and , respectively, as shown in Table 2. Although the values of the structural parameters were estimated after repeated tests and measurements, the uncertainties of the AACMM were still relatively large. According to the analysis in Section 3, it is known that such movement uncertainty may be caused by the uncertainties of the structural parameters.

The significant uncertainties shown in Table 2 are unacceptable in practical application as this may severely deteriorate the overall measurement accuracy. Thus, the identification of uncertain structural parameters should be further conducted, which are then used to compensate the modeling errors in the DH model and then to improve the performance.
From the experimental studies, it is found that the movement uncertainty becomes relatively stable when the number of identification poses is greater than 30, and it is almost invariant when the number of identification poses is larger than 50. Therefore, the joint angles and reference coordinates of 50 poses were acquired to identify the structural parameters of the AACMM, and the identification algorithm presented in Section 3.3 is used. The identification results are shown in Table 3. And the results can be reached just after one time iteration. Further experiments show that there is no further benefit to increase iteration times.
 
mm. 
Compared to Table 1, one may find that the identified parameters in Table 3 are very close to the nominal values. This result indirectly validates the efficacy of the suggested identification algorithm. There are also some structural parameters (e.g., , , , and ) which are not changed for the reason that they are redundant and do not need to be identified. To further verify the identification results, another 200 groups of data were acquired to test the movement uncertainty of the AACMM. Figure 7 shows that the maximum movement uncertainties are −0.112 mm, 0.116 mm, and 0.117 mm in the directions of , , and , respectively, and the averages of the absolute movement uncertainties are 0.036 mm, 0.038 mm, and 0.044 mm, which are summarized in Table 4. It is found from Tables 2 and 4 that the overall accuracy can be significantly improved by using the proposed identification and compensation method. The comparisons between the movement uncertainties after identification and the ones before identification are shown in Figure 8. As it is shown, all the movement uncertainty indexes are greatly improved after identification.

(a) In the direction of
(b) In the direction of
(c) In the direction of
5. Conclusions
A constructive parameter identification approach for articulated arm coordinate measuring machines has been presented in this paper. A structural kinematic model is established based on the DH method and verified through experiments. Based on the difference between the coordinates of the probe calculated by the kinematic model and the reference coordinates, a mathematical parameter identification model is further developed to decrease the uncertainties in the DH model. The analysis of the Jacobian matrix in the identification model shows that there are two structural parameters which are related to others in our case studies. Therefore, these structural parameters are removed from the identification model. Then the structural parameter identification with the aim to get the leastsquare solution of the identification model can be carried out by using the obtained poses and the reference coordinates of the AACMM. To facilitate practical implementations, experimental studies have been conducted. These experimental results have revealed the effectiveness of the proposed structural parameter identification for AACMM.
The proposed modeling and identification method can be extended to the calibration of serial robots, where the identification method is required due to their dynamic operation environments. This will be further studied in our future work.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant no. 51465027).
References
 Y. Li and P. R. Nomula, “Surfaceopening feature measurement using coordinatemeasuring machines,” International Journal of Advanced Manufacturing Technology, vol. 79, no. 9–12, pp. 1915–1929, 2015. View at: Publisher Site  Google Scholar
 W. Wang, G. Gao, Y. Wu, and Z. Chen, “Analysis and compensation of installation errors for circular grating angle sensors,” Advanced Science Letters, vol. 4, no. 67, pp. 2446–2451, 2011. View at: Publisher Site  Google Scholar
 X. H. Li, B. Chen, and Z. R. Qiu, “The calibration and error compensation techniques for an Articulated Arm CMM with two parallel rotational axes,” Measurement: Journal of the International Measurement Confederation, vol. 46, no. 1, pp. 603–609, 2013. View at: Publisher Site  Google Scholar
 D. Zheng, Z. Xiao, and X. Xia, “Multiple measurement models of articulated arm coordinate measuring machines,” Chinese Journal of Mechanical Engineering (English Edition), vol. 28, no. 5, pp. 994–998, 2015. View at: Publisher Site  Google Scholar
 F. Romdhani, F. Hennebelle, M. Ge, P. Juillion, R. Coquet, and J. F. Fontaine, “Methodology for the assessment of measuring uncertainties of articulated arm coordinate measuring machines,” Measurement Science and Technology, vol. 25, no. 12, Article ID 125008, 2014. View at: Publisher Site  Google Scholar
 A. PiratelliFilho, F. H. T. Fernandes, and R. V. Arencibia, “Application of virtual spheres plate for AACMMs evaluation,” Precision Engineering, vol. 36, no. 2, pp. 349–355, 2012. View at: Publisher Site  Google Scholar
 J. Sładek, K. Ostrowska, and A. Gaęska, “Modeling and identification of errors of coordinate measuring arms with the use of a metrological model,” Measurement, vol. 46, no. 1, pp. 667–679, 2013. View at: Publisher Site  Google Scholar
 Y. Meng and H. Zhuang, “Autonomous robot calibration using vision technology,” Robotics and ComputerIntegrated Manufacturing, vol. 23, no. 4, pp. 436–446, 2007. View at: Publisher Site  Google Scholar
 J. U. Dolinsky, I. D. Jenkinson, and G. J. Colquhoun, “Application of genetic programming to the calibration of industrial robots,” Computers in Industry, vol. 58, no. 3, pp. 255–264, 2007. View at: Publisher Site  Google Scholar
 I. Kovač and A. Frank, “Testing and calibration of coordinate measuring arms,” Precision Engineering, vol. 25, no. 2, pp. 90–99, 2001. View at: Publisher Site  Google Scholar
 J. Santolaria, J.J. Aguilar, D. Guillomía, and C. Cajal, “A crenellatedtargetbased calibration method for laser triangulation sensors integration in articulated measurement arms,” Robotics and ComputerIntegrated Manufacturing, vol. 27, no. 2, pp. 282–291, 2011. View at: Publisher Site  Google Scholar
 J. Santolaria, A. C. Majarena, D. Samper, A. Brau, and J. Velázquez, “Articulated arm coordinate measuring machine calibration by laser tracker multilateration,” The Scientific World Journal, vol. 2014, Article ID 681853, 11 pages, 2014. View at: Publisher Site  Google Scholar
 H. Hamana, M. Tominaga, M. Ozaki, and R. Furutani, “Calibration of articulated arm coordinate measuring machine considering measuring posture,” International Journal of Automation Technology, vol. 5, no. 2, pp. 109–114, 2011. View at: Publisher Site  Google Scholar
 G. Gao, J. Lu, and H. Yang, “Study on the structrual parameter coupling of articulated arm coordinate measuring machines,” Telkomnika, vol. 11, no. 5, pp. 2454–2459, 2013. View at: Publisher Site  Google Scholar
 L. Zhu, W. Li, Z. Pan, Y. Guo, and Q. Chen, “Research on parameter selfcalibration method for partlybonded articulated arm coordinate measuring machine,” Chinese Journal of Scientific Instrument, vol. 35, no. 3, pp. 572–579, 2014. View at: Google Scholar
 D. Zheng, Z. Xiao, and Y. Zhou, “Key technologies to improve articulated arm coordinate measuring machines' measurement accuracy,” Journal of Hebei University of Science and Technology, vol. 35, no. 1, pp. 20–23, 2014. View at: Google Scholar
 W. Cheng, Y. Fei, L. Yu, and R. Yang, “Probe parameters calibration for articulated arm coordinate measuring machine,” in Proceedings of the 6th International Symposium on Precision Engineering Measurements and Instrumentation, vol. 7544 of Proceedings of SPIE, Hangzhou, China, August 2010. View at: Publisher Site  Google Scholar
 J. Denavit and R. S. Hartenberg, “A kinematic notation for lowerpair mechanisms based on matrices,” ASME Journal of Applied Mechnics, vol. 22, no. 6, pp. 215–221, 1955. View at: Google Scholar  MathSciNet
 L. Q. Zhu, W. X. Li, Z. K. Pan, Y. K. Guo, and Q. S. Chen, “Coordinate measuring method with two operation modes based on the adjustable articulated arms,” Optical Engineering, vol. 53, no. 12, Article ID 122407, 2014. View at: Publisher Site  Google Scholar
 A. Brau, M. Valenzuela, J. Santolaria, and J. J. Aguilar, “Evaluation of different probing systems used in articulated arm coordinate measuring machines,” Metrology and Measurement Systems, vol. 21, no. 2, pp. 233–246, 2014. View at: Publisher Site  Google Scholar
 R. J. Lapeer, S. J. Jeffrey, J. T. Dao et al., “Using a passive coordinate measurement arm for motion tracking of a rigid endoscope for augmentedreality imageguided surgery,” International Journal of Medical Robotics and Computer Assisted Surgery, vol. 10, no. 1, pp. 65–77, 2014. View at: Publisher Site  Google Scholar
 G.M. Daniel, B. Joaquín, C. Eduardo, and M.P. Susana, “Influence of human factor in the AACMM performance: a new evaluation methodology,” International Journal of Precision Engineering and Manufacturing, vol. 15, no. 7, pp. 1283–1291, 2014. View at: Publisher Site  Google Scholar
 G. Gao, W. Wang, and J. Zhou, “Study on the error transfer of articulated arm coordinate measuring machines,” TELKOMNIKA Indonesian Journal of Electrical Engineering, vol. 11, no. 2, pp. 637–641, 2013. View at: Publisher Site  Google Scholar
 S. Sagara and Y. Taira, “Cooperative manipulation of a floating object by some space robots: application of a tracking control method using the transpose of the generalized Jacobian matrix,” Artificial Life and Robotics, vol. 12, no. 12, pp. 138–141, 2008. View at: Publisher Site  Google Scholar
 G. Pond and J. A. Carretero, “Formulating Jacobian matrices for the dexterity analysis of parallel manipulators,” Mechanism and Machine Theory, vol. 41, no. 12, pp. 1505–1519, 2006. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2016 Guanbin Gao et al. 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.