Research Article | Open Access

Volume 2010 |Article ID 541809 | https://doi.org/10.1155/2010/541809

Nenggen Ding, Saied Taheri, "Application of Recursive Least Square Algorithm on Estimation of Vehicle Sideslip Angle and Road Friction", Mathematical Problems in Engineering, vol. 2010, Article ID 541809, 18 pages, 2010. https://doi.org/10.1155/2010/541809

# Application of Recursive Least Square Algorithm on Estimation of Vehicle Sideslip Angle and Road Friction

Revised05 Dec 2009
Accepted10 Feb 2010
Published22 Apr 2010

#### Abstract

A recursive least square (RLS) algorithm for estimation of vehicle sideslip angle and road friction coefficient is proposed. The algorithm uses the information from sensors onboard vehicle and control inputs from the control logic and is intended to provide the essential information for active safety systems such as active steering, direct yaw moment control, or their combination. Based on a simple two-degree-of-freedom (DOF) vehicle model, the algorithm minimizes the squared errors between estimated lateral acceleration and yaw acceleration of the vehicle and their measured values. The algorithm also utilizes available control inputs such as active steering angle and wheel brake torques. The proposed algorithm is evaluated using an 8-DOF full vehicle simulation model including all essential nonlinearities and an integrated active front steering and direct yaw moment control on dry and slippery roads.

#### 1. Introduction

The performance of a vehicle active safety system depends on not only the control algorithm, but also on the estimation of some key states if they can not directly be measured. Among these states to be estimated online, vehicle sideslip angle and tire-road friction coefficient have been extensively studied in the literature. It is noted that road friction is also used for determination of target response to driver’s steering inputs in the entire range of operation.

There are several strategies  for estimation of sideslip angle and road friction, such as Kalman filter (KF), RLS algorithms, closed-loop feedback observers  and sliding mode observers . Gustafsson and many others  designed observers based on tire-road force models and single-track vehicle models. In , Yi et al. designed an observer of tire-road friction coefficient using RLS methods for vehicle collision warning/avoidance system. Wenzel et al. , and Baffet et al. , independently proposed a dual extended Kalman filter (DEKF) for estimation of vehicle states and parameters which is intended for various active chassis control systems. From the viewpoint of online estimation, the DEKF is too complex to be used with the current systems, due to the limited computational authority of microprocessors and availability of only a few sensors.

A common feature of most state observer and KF/RLS based algorithms for estimation of sideslip angle is that they rely heavily on an accurate tire model, which may vary during vehicle operation. To overcome the limitation, Hac and Simpson , combined the state estimation method based on vehicle dynamic model with a closed-loop nonlinear observer to estimate yaw rate and sideslip angle of the vehicle. This resulted in good estimates for maneuvers on high-friction roads using “pseudomeasurement” of yaw rate as preliminary estimate to supply additional feedback to the observer, which gets rid of the sensor for direct measurement of yaw rate. However, sideslip angle cannot be estimated with enough accuracy on very low-friction roads. Cheli et al. , estimated sideslip angle as a weighted mean of the results provided by a kinematic formulation and those obtained through a state observer based on vehicle single-track model. The basic idea is to make use of the information provided by the kinematic formulation during a transient maneuver to update the single-track model parameters (tire cornering stiffness).

In , tire-road friction estimation (TRFE) methodologies are classified as four types. The first approach is slip-based and the other three use dedicated optical, acoustic [15, 16], or strain gage type sensors [7, 1518]. Slip-based approaches use wheel slip calculated based on the difference between the wheel velocities of driven and nondriven wheels at normal driving conditions [3, 14, 19]. The slip-based approach can be further extended to model-based friction estimation [1, 8, 12, 14], in which wheel dynamics and/or brake pressure model, and vehicle longitudinal and/or lateral dynamics, are employed.

The main idea of most slip-based friction estimation approaches is to predict the maximum friction based on the collected low-slip and low-friction data at normal driving, where normally acceleration/deceleration is less than 0.2 g  and slip ratios are rarely greater than 5%. In these cases, maximum tire-road friction is estimated according to the slope at the low values of slip/friction curve of the driven wheels, which is mainly determined by the tire carcass stiffness rather than the road condition and thus quite sensitive to tire type, inflation pressure, tire wear, and possibly vehicle configuration .

Estimation of vehicle sideslip angle relies heavily on an accurate tire model, but the computing power required in such detailed models easily exceeds the control cycle time. For example, the famous Magic Formula tire model can very accurately represent the force and moment properties. However, due to trigonometric and exponential functions associated with such formulation with several associated coefficients, the time required for online calculation, far exceeds the control cycle time. On the other hand, if tire and vehicle dynamic models do not include the required details, estimation accuracy will be lost. Consequently, several efficient and accurate tire models with simpler expressions of tire forces have been investigated in the literature. In other words, a compromise between accuracy and complexity of tire model is required for online implementation.

In this paper, an RLS algorithm is proposed to estimate sideslip angle and road friction for online application during activation of active front steering and direct yaw moment control. Two main means are adopted to reduce the computational time of the algorithm. The first one is use of a modified Dugoff tire model, for which simple expressions of tire forces are used and parametric differences with respect to tire normal forces can be easily functionalized using polynomials. The second one is online linearization of the model and iterative computation for the proposed algorithm which are distributed to different control cycles without sacrificing estimation accuracy. A significant merit of the algorithm lies in the fact that it can provide estimates with reasonable accuracy without additional sensors. It makes adequate use of the data available through the control logic for correction steering angle and wheel brake torques. Comparison between estimated results and simulation data using Matlab/Simulink and an 8-DOF full vehicle model shows that the proposed algorithm is promising for practical use in active safety systems.

#### 2. Algorithm Description

The algorithm considered in this paper is intended for real-time online estimation of sideslip angle and road friction for the vehicle stability control systems using active front steering and direct yaw moment control. Due to the limited computational authority of microprocessors used in vehicle control systems, the algorithm should not be too complex. A simple two-degrees-of-freedom (DOF) vehicle model is used to develop the estimation algorithm. Shown in Figure 1 is the scheme of the proposed algorithm for sideslip angle and road friction estimation.

The estimator in Figure 1 consists of the vehicle model and the RLS algorithm. All the inputs of the estimator are from the controller (i.e., electric control unit, or ECU), either directly measured by sensors or estimated by the control algorithm. Here, it is assumed that such typical sensors as those for lateral acceleration and yaw rate of the vehicle, steering wheel angle, and wheel speeds are used. Vehicle speed is estimated using the wheel speeds.

In the vehicle model, tire forces are computed according to the estimated vehicle states and , estimated road friction and , and control inputs of steering angle and brake torque “vector” . The lateral and yaw accelerations are then estimated and used in the RLS algorithm. Errors of lateral acceleration and yaw acceleration of the vehicle are computed according to inputs from the vehicle model and the controller, and the RLS algorithm uses these errors to recursively estimate the sideslip angle and road friction.

##### 2.1. A 2-DOF Vehicle Model

Shown in Figure 2 is the 2-DOF vehicle model with 4 wheels. The model is intended for online computation of lateral and yaw accelerations based on vehicle states and road friction conditions. Equations governing the lateral and yaw motion of the vehicle are as follows

Note that the vehicle velocity at the center of gravity, , is the resultant vector of longitudinal speed and lateral speed . Tire slip angles can be determined according to kinematic relationships shown in Figure 2 as follows:

Wheel load transfer is included in calculation of tire normal force as follows: where and need to be properly selected for each specific tire, and is the estimated longitudinal acceleration according to the estimated longitudinal speed. Although lateral load transfer should be distributed between the two axles according to suspension roll stiffness and some other factors, here only the average lateral inertia force is included and distributed.

During a typical intervention of AFS/DYC, the tires often operate at or near the friction limit and combined-slip conditions may arise. Therefore, a nonlinear tire model capable of simulating the friction ellipse phenomena is required. A modified Dugoff tire model is used here for the estimator. First, lateral tire forces at pure-slip conditions are calculated using the modified Dugoff model and longitudinal forces are determined from the brake torques. Then, the lateral forces are further amended according to the magnitude of the longitudinal force. Variation of tire-road friction with respect to slip is included in the calculation of the lateral forces.

The pure-slip lateral force is first calculated for dry asphalt road with a nominal tire-road friction coefficient = 1.0 and then is adjusted for the estimated friction coefficient ( or ). Calculation of can be summarized as follows: where the coefficients ’s () in (2.4) can be determined according to tire test data or drawn from other tire models which have a high accuracy but are not suitable for online computation, and can be determined from (2.2) as follows: where has been set equal to due to the fact that, during AFS/DYC intervention, the total steering angle at front wheels is not likely to exceed 20° (the relative error at 20° is only approximately 4.1%).

The coefficient is used to compensate for overlimitation of tire force value at large slip rates in the original Dugoff tire model. Shown in Figure 3 are some examples of lateral forces using the modified Dugoff model and a Magic Formula model, where the coefficients in the former are drawn from the latter. When slip angles are small or very large, the two models are close. For a real vehicle, the slip angle is normally less than 10°, and thus the modified model is suitable for practical use. However, the computational power required for the modified Dugoff model using (2.4) and (2.5) is much less than that of Pacejka Magic Formula model.

When DYC is activated, brake torque is applied on some of the wheels. For simplification, tire longitudinal forces are calculated using the following equation where is the brake torque commanded by the controller. Since a wheel slip controller is usually incorporated in the AFS/DYC control system, the above simplification is made based on the following assumptions:()The target brake torque can always be realized without delay;()The longitudinal slip rate is constant during activation of DYC.

Further, is assumed to be 0.2 for determination of the final lateral forces at combined-slip conditions as follows In order to reduce the computational requirements, (2.8) was fitted to (2.7) where takes a fixed value of 0.2:

When slip angle is less than 10° (this is almost always true as mentioned above), or equivalently , the results using (2.7) and (2.8) are very close as shown in Figure 4.

##### 2.2. RLS Algorithm

The recursive least square algorithm introduced here was developed based on a simple two-degrees-of-freedom vehicle model. Due to the nonlinearities involved in the equations, online linearization becomes a dynamic part of the algorithm.

As shown in Figure 1, tire-road nominal friction coefficients for the left and right sides and the sideslip angle at the vehicle’s center of gravity need to be estimated. Since both vehicle lateral acceleration and yaw acceleration are nonlinear functions of the estimated parameters, linearization of these functions is required when employing the RLS algorithm. Let and define the functions at a given point as follows:

At this point, functions and can be linearized as follows where

Equation (2.11) can be rewritten as follows: where

Now the unknown parameter vector and the state vector can be defined using the following equations:

The state vector error can be expressed as Where

Define an index with a forgetting factor as follows:

By adopting the formulations given above and using the procedure in  for minimizing the index given by (2.18), the unknown parameters can be estimated as follows where is always set to zero (this is only specific for the problem stated here), is a large number, and is a unit matrix with the proper dimensions.

On certain situations such that when estimated road friction coefficients or sideslip angle are far from current operating point, a new operating point is needed and linearization of the functions and need to be renewed since (2.11) or (2.13) hold only when is small. Therefore, should be renewed using according to the criteria defined in the next section. Matrix must be recalculated when operating point is renewed.

#### 3. Algorithm Implementation

##### 3.1. Numerical Implementation of Partial Derivatives

The following equations are used to compute the partial derivatives at :

Computation of the partial derivatives using (3.1) may be time-consuming. Therefore, the alternative solution adopted in this paper is distributing the computation among different control cycles to help with the limited computational authority of the microprocessor.

##### 3.2. Criteria of Reparameterization

When any of the following conditions holds, substitute with and a new round of recursive process will be initiated: where is a preset small positive number, is the number of iterative computation using (2.19) since last re-parameterization, and is a known threshold.

If the value of estimated parameters changes too quickly, restrictions for their increment are applied. In this paper, maximum increment for , and are 0.1, 0.1, and 0.005 rad, respectively. For example, if and , the value of will be restricted to 0.46; if rad and rad, the value of will be restricted to −0.0027 rad. Whenever maximum increment is violated, the value for matrix is reset as .

##### 3.3. Summary of Procedure for the Estimator

To facilitate understanding of parameter estimation using the proposed RLS algorithm, the main steps of the procedure are outlined as follows.

Step 1 (determinaion of and at ). and are estimated using the 2-DOF vehicle model according to (3.3), which is rearranged from (2.1) with consideration of and using approximation for and under small angle assumption. and are used as the first group of parameters for a new round as in (2.19). In the following steps, (3.3) is always used if the 2-DOF vehicle model is involved.
Go to Step 2.

Step 2 (calculation of partial derivatives with respect to ). and () are computed using the 2-DOF vehicle model. Then and are determined according to the first formula in (3.1).

Go to Step 3.

Step 3 (calculation of partial derivatives with respect to ). Similarly, and () are computed using the 2-DOF vehicle model. Then and are determined according to the second formulae in (3.1).

Go to Step 4.

Step 4 (calculation of function values of ). and are computed using the 2-DOF vehicle model.

Go to Step 5.

Step 5 (calculation of function values of ). and are computed using the 2-DOF vehicle model. Then and are determined according to the third formulae in (3.1).

Go to Step 6.

Step 6 (ierative calculation using the RLS algorithm). Matrix and vector are calculated using (2.19). Estimations of the unknown parameters are available in this step.

If any of the conditions in (3.2) holds, substitute with and then go to Step 1; otherwise, continue with Step 6.

From the viewpoint of online application, each of the steps is intended to be executed within one control cycle.

#### 4. Evaluation of the RLS Algorithm by Simulation

The algorithm is evaluated using the data from simulation of an AFS/DYC-based integrated control system. Simulation of a double lane change maneuver is conducted using Matlab/Simulink. A nonlinear 8-DOF vehicle model along with a combined-slip tire model and a single-point preview driver model is used. Control commands are executed through correction steering angle on front wheels and brake torque applied on one of the four wheels.

The data for the steering angles at front wheels, brake torques on the four wheels, yaw rates, lateral acceleration, and vehicle speeds are used as inputs to the RLS based estimator. Estimated results of vehicle sideslip angles and road friction coefficients are compared with those from the simulation of double lane change maneuver using Matlab/Simulink. This enables the reader to evaluate whether the results are sufficiently precise to be used in control.

Two scenarios of double lane change maneuvers are involved: one is on high friction road surface and the other is on low friction road surface, and the target vehicle speeds for the two scenarios are 110 km/h and 40 km/h, respectively.

The initial sideslip angle and nominal tire-road friction coefficients on both sides are assumed to be 0, 0.8, and 0.8, respectively. The forgetting factor taking a value of 0.975 and in (2.19) is set to 1. The results are shown in Figures 5 and 6.

In each figure, the first one or two diagrams illustrate the inputs to the estimator. Estimated results are plotted in the second diagram, together with the actual data for comparison. When integrated control quits from intervention, the estimated sideslip angle and nominal tire-road friction coefficients are reset to their initial values. This is due to the fact that, with the current sensors onboard vehicles equipped with active safety systems, it is not possible to determine the surface coefficient of adhesion as long as vehicle remains within the linear range of operation .

For the double lane change maneuver performed on dry road with and at 110 km/h, comparison of estimated and actual data in Figure 5 shows that the estimates of the yaw rate track the actual values with reasonable accuracy, and that the estimates of road friction coefficients are, on average, less than the actual values. However, the estimates of road friction coefficients can still provide useful information of road adhesion for the control algorithm and are acceptable in the sense of road conditions in terms of slipperiness: normal (), slippery (), and very slippery () .

Figure 6 illustrates a double lane change maneuver performed on slippery road with and at 40 km/h. This is a difficult maneuver from the estimation viewpoint, because of extremely slippery surface and low speed. Nevertheless, the estimates of the yaw rate and road friction coefficients track the actual values quite well during the maneuver. It is shown that the estimates of road friction after about 11s, when the double lane change maneuver has been completed, are not quite accurate. However, this inaccuracy has no adverse effect on control because information about the road friction within the linear range of operation is not required.

For evaluating the accuracy of the above estimated sideslip angle using the RLS algorithm, some results cited from  for a vehicle performing the same maneuver but using a methodology that combines a kinematic formulation and a state observer based on a single track vehicle model are shown in Figure 7 for comparison with the results shown in Figures 5(b) and 6(b). The results of Figure 7 show that the methodology used in  yields high accuracy of estimation. It is found from Figure 5(b) that when sideslip angle changes abruptly such as those from 4.2 s to 5.8 s, the estimate can not catch up with its actual value fast enough and thus a relatively large error arises. However, this error can be corrected by combining the proposed RLS algorithm with kinematic formulation, just like the methodology used in .

To evaluate robustness of the proposed RLS algorithm with repect to certain variantons that may occur during vehicle operation (mass, moment of inertia, tire cornering stiffness etc.), more simulation was performed. As an example, Figure 8 illustrates the results obtained in a double lane change maneuver performed with different vehicle inertia properties. The vehicle inertia parameters are designated as = 1535 kg, = 1318 kg, = 445 kgm2, = 2513 kgm2 for the estimator, while those for the double lane change maneuver are = 1997 kg, = 1780 kg, = 601 kgm2, = 3269 kgm2. Though partly deviated from the actual states, the estimates are generally acceptable. In Figure 8a, the average estimated road friction coefficients deviate more from their actual values than those in Figure 5(b). For both cases shown in Figure 8, there are certain periods of time when the estimates of sideslip angle have a large error and lag, which indicates that the parameters for the algorithm should be further tuned to improve its robustness. Again, these errors appear during abrupt change of sideslip angle and can be reduced by combining the proposed RLS algorithm with kinematic formulation.

#### 5. Conclusion

A model-based recursive least square algorithm for estimation of sideslip angle and road friction using data from the active front steering and dynamic yaw control logic is proposed. The estimates are evaluated through simulation of double lane change maneuvers using Matlab/Simulink. The results indicate that the strategy of estimation is valid and successful without using additional sensors, on both high and low friction road surfaces. Robustness of the algorithm is evaluated through more simulation with variation of vehicle inertia properties, and results show that the estimates are generally acceptable but the parameters for the algorithm need to be further tuned.

Though not yet included in our investigation, we propose that the RLS algorithm developed in this research be combined with kinematic formulation to enhance estimation accuracy during abrupt change of sideslip angle.

Future work of the research may include evaluation of the methodology through hardware-in-the-loop and road tests and implementation of the estimation algorithm on a vehicle stability enhancement system for online applications.

#### Nomenclature

Subscript, Abbreviation, and Symbol
 fl: Front left fr: Front right rl: Rear left rr: Rear right f: Front r: Rear COG: Center of gravity : Indicator for estimated value : Indicator for error between measured and estimated values.
Parameters and Variables
 : Horizontal distance between vehicle COG and front axle : Longitudinal acceleration of vehicle : Lateral acceleration of vehicle B: Matrix for RLS algorithm : Horizontal distance between vehicle COG and rear axle : Tire longitudinal slip stiffness : Tire cornering stiffness : Peak value of lateral force of tire : Longitudinal tire force in tire -axis (of wheel plane) : Lateral tire force in tire -axis (of wheel plane) at pure-slip condition : Lateral tire force in tire -axis (of wheel plane) at combined-slip condition : Peak value of lateral tire force in tire -axis : Lateral tire force at pure lateral sliding in tire -axis : Vertical force on tire : Vertical static force on tire : COG height of total vehicle mass with respect to ground : Roll moment of inertia (about vehicle -axis) Yaw moment of inertia (about vehicle -axis) : Wheel base : Total vehicle mass : Sprung mass of vehicle : Yaw rate : Tire static loaded radius : Tire lateral slip rate : Brake torque vector for all the four wheels, defined as : Brake torque on a single wheel : Time : Wheel track : Longitudinal velocity : Velocity vector at vehicle COG : Lateral velocity x: State vector : Tire sideslip angle : Vehicle sideslip angle at COG : Total applied steer angle at wheels : Applied steer angle at wheels, result of driver’s input : Correction steer angle at wheels supplied by AFS : Longitudinal slip rate : Angle between velocity vector and vehicle -axis : Forgetting factor : Parameter vector to be estimated : Index : Tire-road nominal friction coefficient.

#### Acknowledgement

This work is supported by National Science Fund of China, with an approval number of 50475003.

1. K. Yi, K. Hedrick, and S.-C. Lee, “Estimation of tire-road friction using observer based identifiers,” Vehicle System Dynamics, vol. 31, no. 4, pp. 233–261, 1999. View at: Publisher Site | Google Scholar
2. A. Hac and M. D. Simpson, “Estimation of vehicle sideslip angle and yaw rate,” SAE Technical Paper Series 2000-01-0696, 2000. View at: Google Scholar
3. K. Li, J. A. Misener, and K. Hedrick, “On-board road condition monitoring system using slip-based tyre-road friction estimation and wheel speed signal analysis,” Proceedings of the Institution of Mechanical Engineers, Part K, vol. 221, no. 1, pp. 129–146, 2007. View at: Publisher Site | Google Scholar
4. U. Kiencke and A. Daiss, “Observation of lateral vehicle dynamic,” Control Engineering Practice, vol. 5, no. 8, pp. 1045–1050, 1997. View at: Google Scholar
5. E. Esmailzadeh, A. Goodarzi, and G. R. Vossoughi, “Optimal yaw moment control law for improved vehicle handling,” Mechatronics, vol. 13, no. 7, pp. 659–675, 2003. View at: Publisher Site | Google Scholar
6. C. Canudas de Wit, H. Olsson, K. J. Åström, and P. Lischinsky, “A new model for control of systems with friction,” IEEE Transactions on Automatic Control, vol. 40, no. 3, pp. 419–425, 1995.
7. F. Gustafsson, “Slip-based tire-road friction estimation,” Automatica, vol. 33, no. 6, pp. 1087–1099, 1997. View at: Publisher Site | Google Scholar | MathSciNet
8. L. Alvarez, J. Yi, R. Horowitz, and L. Olmos, “Dynamic friction model-based tyre-road friction estimation and emergency braking control,” ASME Journal of Dynamic Systems, Measurement, and Control, vol. 127, no. 1, pp. 22–32, 2005. View at: Publisher Site | Google Scholar
9. L. R. Ray, “Nonlinear tire force estimation and road friction identification: simulation and experiments,” Automatica, vol. 33, no. 10, pp. 1819–1833, 1997.
10. J. Yi, L. Alvarez, X. Claeys, and R. Horowitz, “Emergency braking control with an observer-based dynamic tire/road friction model and wheel angular velocity measurement,” Vehicle System Dynamics, vol. 39, no. 2, pp. 81–97, 2003. View at: Publisher Site | Google Scholar
11. T. A. Wenzel, K. J. Burnham, M. V. Blundell, and R. A. Williams, “Dual extended Kalman filter for vehicle state and parameter estimation,” Vehicle System Dynamics, vol. 44, no. 2, pp. 153–171, 2006. View at: Publisher Site | Google Scholar
12. G. Baffet, A. Charara, D. Lechner, and D. Thomas, “Experimental evaluation of observers for tire-road forces, sideslip angle and wheel cornering stiffness,” Vehicle System Dynamics, vol. 46, no. 6, pp. 501–520, 2008. View at: Publisher Site | Google Scholar
13. F. Cheli, E. Sabbioni, M. Pesce, and S. Melzi, “A methodology for vehicle sideslip angle identification: comparison with experimental data,” Vehicle System Dynamics, vol. 45, no. 6, pp. 549–563, 2007. View at: Publisher Site | Google Scholar
14. T. Shim and D. Margolis, “Model-based road friction estimation,” Vehicle System Dynamics, vol. 41, no. 4, pp. 249–276, 2004. View at: Publisher Site | Google Scholar
15. U. Eichhorn and J. Roth, “Prediction and monitoring of tyre-road friction,” in Proceedings of the 24th Congress on Safety, the Vehicle, and the Road (FISTA '92), vol. 2, pp. 67–74, 1992. View at: Google Scholar
16. B. Breuler, U. Eichhorn, and J. Roth, “Measurement of tyre-road friction ahead of the car and inside the tyre,” in Proceedings of the International Symposium on Advanced Vehicle Control (AVEC '92), pp. 347–353, Yokohama, Japan, September 1992. View at: Google Scholar
17. T. Bachmann, “The importance of the integration of road, tyre, and vehicle technologies,” in Proceedings of 20th World Congress on Federation of Societies of Automobile Engineering (FISITA '95), Montreal, Canada, 1995. View at: Google Scholar
18. B. Breuer, M. Bartz, B. Karlheinz et al., “The mechatronic vehicle corner of Darmstadt University of technology-interaction and cooperation of a sensor tire, new low-energy disc brake and smart wheel suspension,” in Proceedings of the International Federation of Societies of Automobile Engineering (FISITA '00), Seoul, Korea, June 2000, Paper F2000G281. View at: Google Scholar
19. S. Müller, M. Uchanski, and K. Hedrick, “Estimation of the maximum tyre-road friction coefficient,” Journal of Dynamic Systems, Measurement, and Control, vol. 125, pp. 607–617, 2003. View at: Publisher Site | Google Scholar