## 1. Introduction

So far today, numerous types of modern robotic platforms that perform complex tasks are composed of underactuated mechanisms that were deployed by ancient civilizations. Underactuated mechanisms prevail as the most efficient systems because they take more advantage applying the law of conservation of energy than redundantly kinematic systems. Unlike redundant systems, the underactuated systems pose a reduced number of actuators and less independent control variables and naturally take more advantage of the inertial and gravity forces. The redundant systems have a larger number of control variables than variables in the working space. As a consequence, each controlled actuator somehow has to counteract gravity and inertial forces to establish own kinematic behaviors. A major aspect in robotics engineering concerns physical modeling and control of mobility for a robot to provide autonomous navigation. Similar to any biological entity, the ability to purposely navigate is fundamental for an intelligent robot. Over 2000 years ago, the Chinese invented the south-seeking chariot (SSC), which was used to maintain an absolute orientation along very long trips of hundreds of kilometers. Thus, they basically created one of the first absolute direction compass devices that did not require any other element to function, but its inner mechanisms only. The invention’s compass was used to adjust toward a desired orientation at the beginning of a trip, and it was invented nearly 800 years in advance when the magnetic compass was invented. The SSC is basically a differential gearing system with a pointing-out silhouette above (**Figure 1**). The gearing system compensates the chariot turns by gear transmission relations, keeping the statue arm pointing out always to the same direction. **Figure 1a** shows the SSC design, and **Figure 1b** depicts the SSC prototype with its gearing system made of straight wooden gears. However, this system has largely been studied in modern times by numerous authors [1–7].

Any type of navigational and path-tracking task depends on steering systems. And, in order to infer navigation references, mobile robots use a diversity of exteroceptive sensing devices, such as ultrasound sonar, infrared range detectors, cameras, GPS, and so forth. Nevertheless, a common disadvantage using these types of devices is that they have to obtain external measurements with respect to (w.r.t.) the robot’s fixed Cartesian coordinate system and use them to estimate orientation through geometry models with cumulative errors.

Moreover, some types of proprioceptive sensing devices offer inner measurements that are relative to global inertial system (e.g., magnetic compass, GPS). Unlike local exteroceptive sensors, global proprioceptive sensors yield noncumulative overtime measurement errors. For instance, a magnetic compass implemented as a global orientation system measures angle w.r.t. the earth magnetic axis. And then, such measurement angle is arithmetically used to infer a global destination angle likely described in another external inertial system. However, magnetic compasses are sensitive and are affected by other nearby magnetic fields or tided to suffer errors from rotations that set the device perpendicular to the earth’s magnetic axis. Therefore, other types of sensing devices must be combined to recover from missed observations. Furthermore, the global positioning system (GPS) is another type of global orientation measurement device, which requires at least two successive position observations overtime to provide an instantaneous robot’s global angle w.r.t. the earth’s geographic north pole. As a major disadvantage, GPS cannot provide measurements nearly of inside buildings, forests, or in cloudy days because electromagnetic signals are blocked missing observations during arbitrary periods of time.

A diversity of works has reported navigation systems that combine numerous types of measuring orientation devices such as GPS and magnetic compasses, which provide high precision. This type of technology produces sensing measurements tided to established references (e.g., magnetic/north pole). Hence, useful global-specific references have to be inferred by different methods, which may imply in these calculations inverse/direct solutions of geometric triangulations or algebraic models [8–12].

Reliable local/global navigation for a wheeled mobile robot has fundamentals on controllability and maneuverability. And, both robot’s abilities, respectively, must depend on robust driving control models and the steering kinematic designs [13, 14]. These models are the bases for planning and motion control [15, 16] in navigation, and most of them relay on different numerical mathematical solutions in robotics [17]. There are some complex works on robot’s navigation seeking absolute orientation references [19, 20]. There are other navigation works with major emphasis on collision avoidance relaying on kinematic approaches [18, 21].

In this work, the mechanical design, the physical model, and a control system for a tricycle-wheeled robot with fundamentals on underactuated mechanical functions are proposed. The self-steer robot design proposed in this work has been inspired on the south-seeking chariot, in part, to take advantage of the underactuated mechanical compass with absolute direction to the robot’s goal. Therefore, a specific and complex gearing mechanism system was designed to self-steer commuting between an obstacle-avoiding compass and a leading-to-the-goal compass (**Figure 2**). The global goal-reference compass permanently maintains orientation information available in a direct manner for the robot without sensors, despite experimenting multiple evasive maneuvers.

A kinematic model for the local/global compass and commutation mechanisms are deduced. Moreover, a kinematic control law deduced for a three-wheeled structure with one drive-steer frontal wheel and a pair of lateral passive wheels is analyzed and disclosed. The proposed control law estimates the robot’s posture and combines the interactive switching between goal-leading and obstacle avoidance navigation control. One underactuated mechanism directly leads the robot to the goal. Another underactuated mechanism leads the robot toward free-collision routes.

This manuscript presents simulation results to validate this novel approach that combines ancient underactuated automaton type, with a modern-wheeled robot focus. So far today, the authors are not conscious about other similar approaches reported in the scientific literature. This work does not pretend to introduce a comparative analysis nor efficiency with other state-of-the-art robotic trends. The authors of this research believe that the preliminary results presented in this manuscript will evolve into an efficient technological approach in the near future. Its application will establish a novel approach because it allows to directly have absolute angle observation overtime, allowing the mechanisms to lead the robot with naturally global navigation, reducing computational efforts to other algorithmic tasks, and complementing other sensing devices to improve control and perception.

Sections of this chapter have been organized in the following manner. Chapter 2 discloses the kinematic models for the absolute/local mechanisms. In Chapter 3, the commutation system kinematic and time delays are discussed. In Chapter 4, the control model and simulation results are presented. Finally, Chapter 5 discusses some conclusive remarks about this work.

## 2. Absolute/local compass mechanism

### 2.1. Global underactuated compass

The proposed system sets the robot’s desired Cartesian goal as the absolute reference orientation. The absolute compass directly compensates its bearing by the differential angle provided by the two lateral wheels (**Figure 3a**). The mechanical compass directly self-steer the robot toward the global reference without any actuator (**Figure 3b**). When the robot is not leading along the global compass angle, the synchronization chain gradually reorientates the robot’s steering wheel until matching the goal’s absolute angle. **Figure 3c** illustrates the self-steering mechanism parts: (1) passive lateral wheels, (2) frontal driving wheel with passive steering, (3) wheel differential gearing mechanisms, (4) transmission gears coupling the absolute compass, (5) global goal-angle compass, (6) driving actuator/motor, (7) self-steer synchronization chain, and (8) lateral wheel shaft.

In **Figure 4**, a more detailed depiction of the differential mechanism is illustrated. The absolute compass (5) is composed of three differential systems (3) and (5), which transmit rotary motion from the lateral wheels (1) up to the steering wheel (2) (see **Figure 4b**). The compass (5) is compensated by the gearing relations (**Figure 4b**), and the differential rotary motion is transmitted to the gear (5), namely, *e _{abs}*. It is worth noting that the differential motion (3) of the lateral wheels’ shaft (8) is fixed to the main axis (4). The lateral wheels’ angular speeds, dubbed

*w*

_{1}and

*w*

_{2}, are asynchronous w.r.t. the main shaft (4), similarly differential gears

*e*

_{1}and

*e*

_{2}. The differential shafts of gears

*dif*

_{1}and

*dif*

_{2}(

**Figure 4b**) are perpendicularly joined to the main axis (

**Figure 4a**), and the rotary speed is basically the same for both differentials (4). Nevertheless, the angular velocity for gears (

*dif*1

_{1},

*dif*1

_{2},

*dif*2

_{1}, and

*dif*2

_{2}) might be different for each differential mechanism (3); it would depend on each wheel’s (1) instantaneous velocity.

Finally, the compass, differential mechanism (5) that is composed of the gears *e*_{1}, *e*_{2}, and *e _{abs}* is the mechanism providing the absolute orientation reference toward the robot’s goal.

Without loss of generality, it follows that the deduction of the differential angular velocity model transmitted between *e*_{1} and *e*_{2} is provided next. Thus, the angular velocity for *ω*_{e1} is an averaged differential value:

where both lateral angular speeds have clockwise (left-sided wheel) and counterclockwise (right-sided wheel) signs, respectively:

The differential angular speed is equivalent to the main shaft rotary speed, such that

and, in general, it is assumed that the gears’ (*dif*1 and *dif*2) angular velocities are averaged values:

as well as

In such a manner, the angular velocity transmitted to the compass *e*_*abs* poses the following relation:

In the proposed design, it is assumed that the gear *e*_*abs* doubles its radius w.r.t. gears *e*_{1} and *e*_{2}. Therefore, by rewriting Eq. (6)

and by substituting in the previous expression, the differential rotary speed model *ω*_{e1} since ultimately depends on velocities *ω*_{1} and *ω*_{2}:

Therefore, from the previous expression, let us assume that the robot moves along a straight trajectory line, then *ω*_{1} = *ω*_{2}, and no compass lateral motion is yielded; hence, *ω _{eabs}* = 0. Following

**Figure 2b**notation for the gears, the steering wheel with gear

*E*

_{7}is synchronized with

*ω*by the mechanism system illustrated in

_{eabs}**Figure 5a**.

A first gear connection *E*0 − *E*1 has the following speed relation:

Similarly, the gear connection *E*_{1} − *E*_{2} is algebraically simplified:

The parallel connection *E*_{5} with *E*_{2} has the same relation:

Finally, the angular velocity for *E*_{7} w.r.t. *E*_{0}

Thus, the relation (12) means that *E*_{7} poses opposite rotary movement w.r.t. the absolute compass sense of rotation and half of its angular speed as well.

### 2.2. Local evasive compass

Another compass mechanism with local reference frame steers the robot for obstacle avoidance leading the robot along a safe instantaneous angle (**Figure 5b**). When the robot detects near obstacles, the absolute compass is suspended, and a commutator device switches to the evasive local compass mechanism, activating the gear *E*_{8}. **Figure 5b** describes the gearing mechanism that steers the robot for evasive navigation.

In the kinematic model for evasive steering through *E*_{7}, w.r.t. *E*_{8} is deduced next. The relation between *E*_{8} and *E*_{6} is

and the angular velocity for *E*_{4} is

For the perpendicular connection with *E*_{3}, the angular speed model is

For *E*_{3} and *E*_{5} connected in parallel

Finally, the model for *E*_{7} is

Let us highlight that for the evasive local compass *E*_{8}, the escaping orientation is instantaneously set up by an actuator only when any nearby obstacle is detected. And, such local orientation transmits motion to the steering gear (*E*_{7}) through the mechanism of **Figure 5b**. Thus, the local/global compass mechanisms are physical controllers that substitute models and computer algorithms.

## 3. Absolute/local commuting mechanism

The commutation mechanism interactively couples and uncouples either the absolute compass or the local compass (**Figure 6a**). The commuting device switches into the local compass immediately where any nearby obstacle is detected. Alternatively, it switches to the global compass retaking orientation toward the goal as soon as obstacles are no longer detected. **Figure 7a** depicts the gearing transmission that commutates the different steering mechanisms, the local compass (the front), and the global compass (the back).

According to **Figure 6b**, the commutation mechanism yields linear motion (the commutator device of **Figure 6a**); it activates the local compass by rotating a servomotor and shrinking a sliding crank link *L*. Inversely, the absolute compass is activated by turning off the servomotor, and then a spring-mass-damper system stretches the sliding crank link *L*.

The motion transmission system (**Figure 6a**) is inspired by the model of a modern vehicle’s speed transmission box. For the present case, a common shaft works for two asynchronous speed gearing systems. In addition, **Figure 7a** and **7b** illustrates the motion transmission flow by the red color arrowed lines. When either the local compass mechanism or the global compass mechanism is commutated, the flow transmission motion only concerns with its own gearing mechanism featuring a physical continuous controller.

The angular speed model *ω _{l}* [

*s*

^{−1}] for the commutation rotary servomotor (

**Figure 6b**) to switch into the local compass was obtained through fitting empirical observations with the next theoretical model:

where *κ*_{1} [rad/s] is an amplitude factor and *κ*_{2} [1/s] is a slope rate factor. Subsequently, to estimate the commutation time *t _{l}* [s] required to mechanically coupling into the local compass mechanisms, the next model is deduced:

Moreover, when no obstacles are detected, the commuter device switches into its initial state by coupling the absolute compass. In this case, an underactuated system commutates the state by a spring-mass-damper mechanical system with critically dumped configuration, modeled by second-order linear differential equation. Where *m* [kg] is the mass pulling the spring, *c* [kg/s] is the damping effect constant, *k* [kg/s ^{2}] is the spring elasticity constant, and *x*,

Thus, it may be solved as a first-order linear equation by temporally omitting the second-order term

reorganizing and completing the integrals

Thus, by solving the improper integrals and multiplying both sides of the equation by the Euler number e,

and a solution is obtained; to simplify let us define *λ* ≔ − *k*/*c*,

Assuming the integration constant *C*_{1} = 0 and developing the higher order solutions

In addition, by substituting such functions in the next expression

In order to decrease the commutation time of the sliding linear mechanism, a critically damped system is assumed and modeled by (*c*^{2} − 4*mk*) = 0, with terms *λ*_{1} = *λ*_{2}:

It is desired to speed up as much as possible the commutation time. Then, the mechanical device linear displacement is modeled as a function of time by

And, solving for the commutation time for switching to the absolute compass *t _{a}*,

**Figure 8** depicts a sequence of transition times taken by the commuting mechanism. Times *t _{l}* are the times behavior of the slider crank servomotor when shrinking the link. Times

*t*are the times behavior of the spring-mass-damper device when stretching the link. Such time magnitudes are reached due to mechanical motions. However, such commutation times are physically fast enough, and a robot’s trajectory is not affected in the precise points when commutations occur.

_{a}## 4. Kinematic control law

In this section, a kinematic control law is deduced and analyzed. The proposed controller simultaneously controls driving and self-steering velocities and keeps track of the robot’s posture. The robotic platform is a tricycle-type kinematic structure with two lateral passive wheels at the back (**Figure 9a**) and a central active-drive and passive self-steering wheel at the front (**Figure 9b**). Each wheel kinematic is described by three parameters, *α*, *ℓ*, and *r*, where *α* [rad] is the angle of a wheel w.r.t. the *x*-axis of the robot’s fixed inertial system, *ℓ* [m] is the distance between the robot’s centroid and each wheel’s contact point with the ground, and *r* [m] is an ideal wheel’s radius. In addition, the robot’s kinematic structure considers two controlled kinematic variables *β _{t}* and

*ϕ*, where

_{t}*β*[rad] is the instantaneous steering angle of the front-sided wheel and

_{t}*ϕ*[rad] is the instantaneous driving wheel’s angle.

The kinematic parameters and variables describing the robot platform are summarized in **Table 1** according to each type of wheel (fixed passive, steerable, and drivable active).

The robot’s wheel kinematic parameter is modeled by the following constraint equation, which is stated from the wheel plane,

Likewise, the kinematic constraint is modeled in the wheel’s orthogonal plane by

where **J**_{1}, **C**_{1} ∈ ^{N × 3} basically describes the wheel kinematic constraints. Matrix **J**_{2} poses the ideal wheel’s radius, and **R**(*θ*) is the Euler rotation matrix. Vector *ξ** _{t}* describes the state space system, and its components are the robot’s posture

**= (**

*ξ**x*,

*y*,

**,**

*θ***)**

*β*^{Τ}, as well as the front-sided wheel driving speed

*ϕ*

*, and its product with*

_{t}**J**

_{2}establishes a diagonal matrix of the wheels’ tangential velocities.

Moreover, for the particular case of fixed wheels and central orientable wheels, the following kinematic constraints apply for individual wheels.

In the wheel’s plane

In the wheel’s orthogonal plane

where matrix **J**_{1} is the fixed passive back-sided wheels matrix of kinematic constraints, which is obtained by substituting the parameters of **Table 1** in expression (32). Such that, **J*** _{f}* for the fixed wheels is

and by simplifying

Taking into account the kinematic constraints that are similar to the central orientable wheel, the vector model **J*** _{o}* is obtained:

For the front-sided central orientable wheel, since the angle *β _{t}* is nonstationary, and in order to easy its algebraic solution the following trigonometric identities are substituted in (36) and obtain a simplified form for the vector

**J**

*. Thus,*

_{o}Likewise,

and the following row vector is produced:

Therefore, for a tricycle-type robotic structure as in the present context, and with parameters of **Table 1**, the matrix **J**_{1} is consequently defined as

where **J*** _{f}* represents the fixed wheel kinematic constraints and

**J**

*represents the central orientable wheel constraints. The constraints denoted in the general form (31) implicate that vector*

_{o}being **C*** _{f}* and

**C**

*as vectors of the fixed and the central orientable wheel constraint, respectively, with notation*

_{β}Hence, the robot’s posture first-order derivative

and the null space of any matrix **A**, written as a **A** · **x** = 0 and stated in set notation by

Thus, for the particular case of Eq. (41), the null space is

A manner to solve for the null space matrix condition is to reduce it to the echelon form Ψ:

Thus, a vector **x** that satisfies Eq. (43) must be defined then:

The previous expression (45) is a set of linear systems and is rewritten as an augmented matrix [**C**_{1}| **0**]:

The augmented matrix (46) is algebraically solved by reducing it to echelon form, starting by making zeros the first column, but one the first element:

Then, sum up Ψ_{1, i}Ψ_{3, 1} where Ψ_{3, 1} = sin(*β _{t}*):

Since the pivot element Ψ_{2, 2} = 0, rows 2 and 3 are exchanged:

Now, Ψ_{i, 2} = **0** except its pivot element Ψ_{2, 2}, thus let us divide Ψ_{2, i}/ cos(*β _{t}*) to set to 1 such pivoting element:

Thus, by having Ψ in the reduced form, now the system is solved as the system (45):

Hence, rewriting such a solution in the matrix form

Being the vector **x** of this null space of matrix **C**_{1}, such that

For our application purpose, *x*_{3} = cos *β*, then the null space solution

Therefore, Eq. (41) is rewritten because for each instant time *t*, a temporal variable *v _{t}* exist .

and

where *v _{t}* is the instantaneous robot’s absolute velocity and

*w*is instantaneous robot’s yaw rate. As a matter of fact

_{t}**Σ**(

*β*) is the null space of the kinematic constraint matrix in the wheel’s orthogonal plane (Eq. (54)), defined by

_{t}This state control law is expressed in a more compact form:

where

Thus, by substituting expression (57) in the model (58), the kinematic control law is rewritten:

The robot’s posture is essentially modeled by Eq. (60) and validates by the Cartesian trajectory depicted in **Figure 10a**. For this case, *v* was set at constant value overtime, and *w* gradually increased according to the range {0, …, *π*/4}. Then, the robot turns into its original angle and then continues straight.

**Figure 10b** and **10c** shows *β* and *θ* behaviors w.r.t. orientation error *θ _{E}*, when

*θ*= 0 the robot is bearing toward the desired goal. Actually,

_{E}*β*depicts how the steering wheel performs from its starting angle up to the final angle. Because of the absence of obstacles, the steering wheel’s angle experiences no commutation toward the local mechanism. However, the absolute compass controls the wheel from the beginning and gradually steers until orientation is aligned to the global destination. Its behavior along the navigation task is illustrated in plot 10b.

Moreover, in order to instantaneously estimate the deviation angle *θ _{E}*, let us define the global destination

*θ*[rad]. And, such desired orientation is perturbed by escaping angle

_{G}*θ*[rad] due to relative obstacle locations (

_{R}**Figure 11a**). Likewise,

*θ*is established every time the local compass mechanism is switched to. Obstacles are detected within a territorial radius

_{R}*d*, or ignored when the robot moves farther away from the scope

_{l}*d*(

_{l}**Figure 11b**).

For instance, sensor 1 detects a far obstacle, and its repulsive acceleration influence is low. However, sensor 6 yields a greater repulsive acceleration influence since it is nearer to the robot’s location.

Thus, let us postulate:

From the control law, the state vector **z** is now estimated by directly using the instantaneous error angle *θ _{E}*:

and by obtaining the inverse solution for **u**, now the input system involves error angle:

The robot’s controlled orientation is fed back by *θ _{E}* and calculated by numerical successive approximations. This controller is a recursive system that involves both, the direct and the inverse solutions, to reduce the numerical error

*θ*.

_{E}The instantaneous escaping angle *θ _{R}* is obtained within a local coordinate system by directional derivatives, which approach a cosine function, similarly to the work reported in Ref. [20]. A gradient operator of a cosine navigation function approximates a repulsive partial differential equation to evade any observable obstacle. The robot’s repulsive acceleration vector function

*δ*[m]:

_{α}where *k*_{2} [m] is the obstacle diametrical territory and *k _{α}* [

*ms*

^{−2}] is an adjustment constant factor of the physical acceleration amplitude. Thus, the Cartesian geometric model (not measurement)

*r*describes the distance between the robot

_{α}**x**

*= (*

_{r}*x*,

_{r}*y*)

_{r}^{┬}and any obstacle

**x**

*= (*

_{o}*x*,

_{o}*y*)

_{o}^{┬}, which is geometrically modeled by

The angle limits {0 ≤ *φ _{α}* ≤

*π*/2} is a transformation relationship of the obstacle range

*δ*, and the condition

_{α}*δ*<

_{α}*d*limits the repulsive acceleration effects:

_{l}Substituting the functional form of *δ _{α}* in (64) and temporally considering

*k*= 1 for analysis purpose

_{α}and algebraically expanding

Thus, applying the gradient operator ∇_{x, y} w.r.t. *x* and *y* components and algebraically simplifying

and

Therefore, the obstacle-repulsive directional vector expressed in terms of Cartesian components *XY* is

Thus, the model for multiple obstacles *α* produces a controlled escaping direction in robot’s local coordinate framework:

It follows that the evasive acceleration magnitude is defined by

and the instantaneous escaping angle *θ _{R}*

Finally, the proposed kinematic controller has the following scheme described in Algorithm 1.

**Algorithm 1.** Local/global robot’s underactuated system controller

Hereafter, **Figure 12** validates our kinematic model approach by depicting how the robot’s trajectory reaches the final Cartesian goal (−7, −8) in global coordinates. In addition, **Figure 12a** shows the Cartesian trajectory among multiple obstacles and the local/global commuting steering modalities. Likewise, as a manner to show validation, **Figure 12b** illustrates only the component *X* versus the robot’s angle during the same navigational task.

## 5. Conclusion

The interest of this chapter was to introduce the analysis of an alternative kinematic steering controller by using underactuated mechanical compasses and to demonstrate its feasibility, controllability approach, and natural efficiency. Despite the complexity in its implementation, a compass mechanism is proposed because it allows to directly set up the global goal as an absolute reference. Moreover, the compass mechanism itself directly steers the robot to the goal simultaneously avoiding obstacles. These functional features are important advantages for a robotic platform w.r.t. other traditional orientation systems, because the global orientation depends neither on feedback from sensor devices nor on computational complexity expenses to instantaneously estimate the global destination orientation. Unlike kinematic redundant structures, the proposed approach limits electric energy use for one driving actuator, not for steering actuators. Steering controllers no longer spend algorithmic computational resources, reserving such resources either for other robotic algorithmic tasks or for increasing additional hardware devices. Mechanical controllers are slow if compared with software algorithms, even though mechanisms are fast enough w.r.t. the available robot’s mobility speeds. In this regard, a critical issue is the commutation linear actuator, which is passive/active based on a slider crank combined with a spring-mass-damper system. The commutator device was fast enough, proving that sophisticated active linear actuators were not needed. It was found out that commutation times between local/global compasses did not negatively affect the system performance or the trajectory stability. Local/global compasses commuting activity warranted the robot to reach the global destination while avoiding collisions. The controlled trajectories yielded were concatenations of inter-switching segments with no discontinuities found. We concluded that from the implementation perspective, the proposed approach is neither necessarily better nor worst in effectiveness than a traditional redundantly actuated and multisensor approach. As a matter of fact, a traditional approach is easier to physically implement. Nevertheless, a redundant traditional approach is in disadvantage, if driving wheels grow in number, not to mention incrementing steering actuators. As a consequence, the more redundant is a system, the more discretized are the controlled motions and trajectories lose stability, which is an inherent difference from the naturally continuous motions produced by underactuated systems. A compact linear kinematic control law to switch inter-compass usage was deduced, with direct and inverse solutions. It elegantly combined the underactuated mechanism control with a very short algorithm to detect obstacles and to estimate the instantaneous escaping orientation. It was found that the robot’s trajectory continuity may be altered, if kinematic evasion parameters are readjusted consequently changing the mechanism commutation response. Through simulation results, the inter-mechanism interactions and functions were validated. For the specific case of the proposed kinematic structure, the local/global compass showed as much efficiency as any similar redundant system. The future work will focus on fault recovery from slips, sliding, or collision dynamics that get the compass orientation uncalibrated. In addition, not only to further orientation analysis but also positioning sensors to supplement and improve this underactuated approach.