## Optimization-based Control for Bearing-only Target Search with a Mobile Vehicle

Zhuo Li, Keyou You, Shiji Song, Anke Xue

Abstract -This work aims to design an optimization-based controller for a discrete-time Dubins vehicle to approach a target with unknown position as fast as possible by only using bearing measurements. To this end, we propose a bi-objective optimization problem, which jointly considers the performance of estimating the unknown target position and controlling the mobile vehicle to a known position, and then adopt a weighted sum method with normalization to solve it. The controller is given based on the solution of the optimization problem in ties with a least-square estimate of the target position. Moreover, the controller does not need the vehicle's global position information. Finally, simulation results are included to validate the effectiveness of the proposed controller.

Index Terms -target search, controller design, bi-objective optimization, Dubins vehicle, estimation, local.

task where sensor measurements are limited. Clearly, it is not necessary to control the vehicle to such an unreliable estimated position. In fact, the controller should depend on the quality of the target position estimate, which is again significantly affected by the controlled trajectory of the mobile vehicle. That is, the control and the estimation objectives should be jointly considered. To achieve it, we propose a bi-objective optimization problem subject to the dynamical constraint of a discrete-time Dubins vehicle, where the first objective is to control the target as fast as possible to the target (assuming it is known) and the second objective is to optimize the estimation performance of the unknown target position.

## I. INTRODUCTION

T ARGET search is a problem of steering single or multiple mobile vehicles to a target with unknown position. It has a broad range of applications, such as search and rescue in disaster response settings [1], planetary and undersea exploration [2], [3], search for lost targets [4], and environmental monitoring [5]. In these scenarios, it is significant and challenging for mobile vehicles to reach target positions as fast as possible, which is the motivation of this work.

If the target position is known, quite a few approaches have been proposed to minimize the searching time, see [6]-[9] and references therein. Clearly, these approaches cannot be directly applied here. To estimate unknown target positions, various sensor measurements have been utilized, such as bearing angles [10], [11], ranges [12], and time difference of arrivals [13], [14]. However, most of existing works only consider either controlling a mobile vehicle to a known target position as fast as possible [6]-[9] or maximizing the estimation performance of sensor measurements [10]-[17]. In this work, they are called the control objective and the estimation objective, respectively, and are both essential for the target search problem. Based on these observations, a naive idea is to separately study these two objectives [18]-[20], i.e., directly control the vehicle to the estimated target position.

Unfortunately, a reliable estimate cannot always be obtained, particularly in the initial stage of the target search

This work was supported by National Key Research and Development Program of China under Grant No.2016YFC0300801, and National Natural Science Foundation of China under Grant No.41576101.

Z. Li, K. You and S. Song are with Department of Automation and BNRist, Tsinghua University, Beijing 100084, China (youky@tsinghua.edu.cn). A. Xue is with the Key Laboratory for IOT and Information Fusion Technology of Zhejiang, Institute of Information and Control, Hangzhou Dianzi University, Hangzhou 310018, China (e-mail: akxue@hdu.edu.cn).

Usually, the vehicle cannot achieve the control and the estimation objectives simultaneously. Note that the trajectory for the control objective is a line segment from its initial position pointing at the true target, which potentially leads to a poor estimate of the target position [21]. However, the trajectory for the estimation objective significantly deviates from the line segment. Thus, it is essential to balance the two objectives. This work adopts a weighted sum method [22] to transform the bi-objective optimization problem into a composite optimization problem with a weighting parameter representing the relative importance between the two objectives. Although this idea has been widely adopted in [23]-[26], it is still difficult to select an appropriate weighting parameter due to different physical meanings of the objective functions. For example, the objective functions in [26] represent the angle difference and the distance difference, which are also of different magnitudes. To this end, we rewrite the composite objective problem by normalizing the two objective functions and introducing a new weighting parameter. It proves to be much easier to select an appropriate value of the new weighting parameter, since both objective functions are normalized to be dimensionless and between zero and one.

However, the normalization always results in a complicated form of the composite objective function. Furthermore, the discrete-time Dubins vehicle [27] adopted in this work brings a challenging constraint in the composite optimization problem, which is different from the single integrators in [13], [14], [18], [19]. In the literature, evolutionary algorithms are commonly-used, such as the genetic algorithm in [23], pigeon inspired optimization method in [24], and value iterationlike algorithm in [28]. Moreover, the graph-based searching method is thoroughly investigated in [25], [26]. Different from the aforementioned approaches, we attack the complicated optimization problem by decomposing the forward velocity of the vehicle and then transforming the optimization problem into a solvable form. Particularly, an explicit solution can

be directly obtained for some weighting parameters with this approach.

Then, the controller for target search can be given in terms of the solution to the optimization problem and the vehicle's global position information. Since the solution to the optimization problem depends on the unknown target position, the controller cannot be directly implemented. To solve this problem, we design a recursive nonlinear least-square estimate [29] of the target position for the optimization-based controller by utilizing the certainty-equivalence principle. This substantially distinguishes our controller from those in [18]-[20].

In addition, the vehicle's global position is not always available in GPS-denied environments, although it is used in [13]-[20]. In this work, we further propose a controller without the vehicle's global position. Observe that the vehicle only aims to reach the target position, for which the relative position information between the vehicle and the target is sufficient. Therefore, we utilize the bearing angles relative to the target in the vehicle's own coordinate, and derive an optimizationbased controller for the global position system (GPS) -denied target search problem together with the vehicle's orientation angles from a compass. Similarly, the GPS information is not required in [30], [31] to estimate the target position, whereas they employ a group of vehicles to measure the bearing angles relative to their neighbors and to communicate individual estimates with their neighbors. In particular, the communication graphs have to be connected.

The rest of this paper is organized as follows. In Section II, the target search problem is formulated as a bi-objective optimization problem for a discrete-time Dubins vehicle. In Section III, we adopt the weighted sum method to solve the biobjective optimization problem, which is rewritten by transformation and normalization. The optimization-based controller design is given in Section IV, with a recursive estimator and an optimization solver. Section V considers the target search problem in GPS-denied environments, where the controller is based on the bearing angles in the vehicle's coordinate. Simulations are conducted in Section VI to illustrate the effectiveness of the proposed controllers. Finally, some conclusion remarks are drawn in Section VII.

## II. PROBLEM FORMULATION

## A. The target search problem with a mobile vehicle

Consider the scenario in Fig. 1, where a target is located at an unknown position p T = [ x T , y T ] ′ ∈ R 2 and the vehicle's position is p ( k ) = [ x k , y ( ) ( k )] ′ ∈ R 2 at the k -th time step in the global coordinate Σ . We first assume that the GPS information p ( k ) is available and then remove it in Section V. The vehicle measures a noisy bearing angle m k ( ) from itself to the target, i.e.,

$$m ( k ) = \phi ( k ) + e ( k ),$$

where φ k ( ) = arctan ( ( x k ( ) -x T ) / y k ( ( ) -y T ) ) denotes the true azimuth bearing angle, and e k ( ) denotes the measurement noise.

Fig. 1. The bearing-only target search problem with a mobile vehicle.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-3810e18d-6654-4783-805c-0bb640bb6e8b.jpg)

In this work, we adopt the following discrete-time Dubins vehicle [27]

$$\text{ration} \quad \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \\ \text{ion is} \quad \text{$x(k+1)=x(k)+\frac{\frac{2v_{c}}{\omega(k)}\sin\left(\frac{\omega(k)h}{2}\right)\cos\left(\theta(k)+\frac{\omega(k)h}{2}\right)$,} \\ \text{re the} \quad \text{inicate} \quad \text{$y(k+1)=y(k)+\frac{\frac{2v_{c}}{\omega(k)}\sin\left(\frac{\omega(k)h}{2}\right)\sin\left(\theta(k)+\frac{\omega(k)h}{2}\right)$,} \\ \text{r, the} \quad \text{$\theta(k+1)=\theta(k)+\omega(k)h$} \\ \quad \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \, \cdot \end{cases}$$

$$( 2 )$$

to search the target, where v c is the constant forward velocity, ω k ( ) is the angular velocity to be tuned, θ k ( ) is the orientation, and h is the sampling period. Since the forward velocity v c is finite, the vehicle can only move within a certain region in a sampling period, which is called the feasible region, see Fig. 1. The bearing-only target search problem is to sequentially design the angular velocity ω k ( ) only by the noisy bearing angles, such that the vehicle approaches the unknown target position as fast as possible.

Obviously, if p T were known to the vehicle, the optimal angular velocity is the one that drives the vehicle to the closest waypoint to the target within the feasible region, i.e., p c ( k +1) , which is called the control objective. Whereas this waypoint is not achievable for p T is unknown and needs to be estimated. In light of [10], the optimal next waypoint to optimize the estimation performance for the unknown position p T is illustrated as p e ( k +1) in the feasible region, which is called the estimation objective. However, this waypoint even cannot ensure the vehicle to approach the target.

Based on these observations, the optimal angular velocity ω ∗ ( k ) should simultaneously consider both the estimation and the control objectives. To this end, we formulate the following bi-objective optimization problem

$$\max _ { \omega ( k ) } \left [ \begin{matrix} f _ { e } ( \mathbf p ( k + 1 ) ) \\ f _ { c } ( \mathbf p ( k + 1 ) ) \end{matrix} \right ] \text{subject to } ( 2 ), \quad \ \ ( 3 )$$

where the objective functions f e ( p ( k +1)) and f c ( p ( k +1)) represent the estimation and the control objectives, respectively. Notice that both objectives depend on the unknown position p T and are essential to the target search problem.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-e4ac3796-f955-4364-afbb-46890b6b50f3.jpg)

}{√ } ⌈ ⌈ ˜ √

√

Fig. 2. The controller architecture for bearing-only target search.

## B. The controller architecture

To solve the bearing-only target search problem, we design a controller for the mobile vehicle including a recursive leastsquare estimator and an optimization solver. Fig. 2 shows the proposed controller architecture. Firstly, the online measurement m k ( ) is utilized to estimate the unknown target position by the recursive least-square approach. Then, we solve the biobjective optimization problem (3) in ties with the estimate and obtain the angular velocity ω ∗ ( k ) . The details of the controller design are elaborated in later sections.

## III. SOLVING THE BI-OBJECTIVE OPTIMIZATION PROBLEM

To solve the bi-objective optimization problem (3), we adopt the weighted sum method [22] and combine the two objective functions of (3) into the following composite form

$$f _ { e } ( { \mathbf p } ( k + 1 ) ) + \alpha f _ { c } ( { \mathbf p } ( k + 1 ) ), \text{ \quad \ \ } ( 4 )$$

where the weighting parameter α &gt; 0 plays an important role in balancing the two objectives, but lacks physical meaning.

In this section, we first describe the concrete expressions of the two objective functions f e ( p ( k + 1)) and f c ( p ( k + 1)) , and then solve the optimization problem.

## A. The estimation and the control objective functions

Assume the measurement noise e k ( ) in (1) is white Gaussian noise with zero mean and σ 2 variance, i.e., e k ( ) ∼ N (0 , σ 2 ) . The estimation performance is commonly evaluated by the Cramer-Rao lower bound (CRLB) C ( p T ) [32]. We utilize the most recent n measurements at the current time step k glyph[greaterorequalslant] n , denoted as m n ( k ) = [ m k ( -n +1) , m k ( -n + 2) , . . . , m ( k )] ′ , to measure the CRLB, and we adopt the Doptimal design. In light of Theorem 1 in [10], the estimation objective at the time step k is explicitly expressed as

$$\lambda _ { \ } e ( \mathbf p ( k + 1 ) ) = \det ( C ( \mathbf p _ { T } ) ^ { - 1 } ) = \frac { 1 } { \sigma ^ { 4 } } \sum _ { i < j } \frac { \mathbf n ^ { 2 } ( \phi ( j ) - \phi ( i ) ) } { r ^ { 2 } ( i ) r ^ { 2 } ( j ) }, \quad B. \ T r a n _ { \ } n i z a t i o r. \\ \intertext { w h a r a i n i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z a t i z } \text{s} h a r a s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i n s i n t r a n d i m i z a t i o r }$$

where i, j ∈ { k -n +1 , . . . , k +1 } and r ( ) = · ‖ p ( ) · -p T ‖ is the Euclidean distance from p ( ) · to p T , i.e., the range between the vehicle and the target.

For the control performance of the vehicle approaching the target, we consider to use the range difference within a sampling period, which evaluates the approaching speed of the vehicle to the target. We express the control objective by maximizing the sum of the difference of the squared range in n sampling periods

$$f _ { c } ( \mathbf p ( k + 1 ) ) = \sum _ { i = k - n + 2 } ^ { k + 1 } \left ( r ^ { 2 } ( i - 1 ) - r ^ { 2 } ( i ) \right ) \quad \ \ ( 6 ) \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \ \$$

Fig. 3. The decomposition of the forward velocity and geometrical relationship of two sequential waypoints.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-4ee556ad-acc3-4e04-906c-fe501e13752c.jpg)

under the constraint r i ( ) glyph[lessorequalslant] r i ( -1) , since the vehicle is always required to approach the target. Observe that we employ the difference of the squared range in (6) instead of the absolute range difference for the ease of solving the optimization problem.

For brevity, we only utilize the vehicle's current position and measurements to design the angular velocity, i.e., n = 1 in (5) and (6). Then the composite optimization problem in (3)-(6) is expressed as

$$2$$

$$\intertext { n t o r l e } \underset { \omega ( k ) } { \maximize } \, \frac { \sin ^ { 2 } \left ( \phi ( k + 1 ) - \phi ( k ) \right ) } { \sigma ^ { 4 } r ^ { 2 } ( k ) r ^ { 2 } ( k + 1 ) } + \alpha \left ( r ^ { 2 } ( k ) - r ^ { 2 } ( k + 1 ) \right ) } { \text{$n t o r l e $ } } \text{$r(k+1)$} \leqslant r ( k ) \text{ and } ( 2 ),$$

$$( 7 )$$

where the first constraint ensures the vehicle to gradually approach the target, and the second constraint results from the dynamics of the Dubins vehicle.

Remark 1. Both objective functions (5) and (6) depend on the measurements m n ( k ) , which are utilized to estimate the bearing angles and the ranges. In fact, a small number of the measurements, i.e., m n ( k ) with a small n , are sufficient to evaluate both estimation and control performance at the k -th time step, since it has weak correlation with the measurement at the i -th time step for i glyph[lessmuch] k . In this work, we consider a simple case of n = 1 .

In the rest of this section, we show how to solve the composite optimization problem (7).

B. Transformation and normalization of the composite optimization problem (7)

1) Transformation: Observe that it is very difficult to directly solve the composite optimization problem (7). We transform it into a solvable form by decomposing the forward velocity v c into v r ( k ) and v t ( k ) , which denote the radial and tangential velocities of the vehicle relative to the target at the time step k , respectively, see Fig. 3.

Then, the following relation is clear from Fig. 3 that

$$\text{ing in} \quad v _ { c } \begin{bmatrix} \cos ( \theta ( k ) ) \\ \sin ( \theta ( k ) ) \end{bmatrix} = v _ { r } ( k ) \begin{bmatrix} \sin \phi ( k ) \\ \cos \phi ( k ) \end{bmatrix} + v _ { t } ( k ) \begin{bmatrix} \cos \phi ( k ) \\ - \sin \phi ( k ) \end{bmatrix}. \ ( 8 )$$

Assuming that v r ( k , v ) t ( k ) and φ k ( ) are constant within the sampling period from kh to ( k +1) h , we take Euler integral

with respect to the continuous dynamics on the right side of (8) and obtain that

$$x ( k + 1 ) - x ( k ) & = h \left ( v _ { r } ( k ) \sin ( \phi ( k ) ) + v _ { t } ( k ) \cos ( \phi ( k ) ) \right ), \quad \text{which} \\ y ( k + 1 ) - y ( k ) & = h \left ( v _ { r } ( k ) \cos ( \phi ( k ) ) - v _ { t } ( k ) \sin ( \phi ( k ) ) \right ). \quad \text{weight}.$$

Together with the Dubins vehicle dynamics (2), it holds that

$$\text{together with the Dubins vehicle dynamics (2), it holds that } \text{     control} \\ \frac { 2 v _ { c } } { \omega ( k ) } \sin \left ( \frac { \omega ( k ) h } { 2 } \right ) \cos \left ( \theta ( k ) + \frac { \omega ( k ) h } { 2 } \right ) \text{                                                                                                                                                                                                         the two                                                                                                                                                                                                       
      = & h \left ( v _ { r } ( k ) \sin ( \phi ( k ) ) + v _ { t } ( k ) \cos ( \phi ( k ) ) \right ), \quad ( 9 ) \\ \frac { 2 v _ { c } } { \omega ( k ) } \sin \left ( \frac { \omega ( k ) h } { 2 } \right ) \sin \left ( \theta ( k ) + \frac { \omega ( k ) h } { 2 } \right ) \text{                                                                                                                                                                         
      = & h \left ( v _ { r } ( k ) \cos ( \phi ( k ) ) - v _ { t } ( k ) \sin ( \phi ( k ) ) \right ). \quad \text{                                                                                                                                                                                                    
   Furthermore, dividing the two equalities in ( 9 ) leads to that$$

Furthermore, dividing the two equalities in (9) leads to that

$$\tan \left ( \theta ( k ) + \frac { \omega ( k ) h } { 2 } \right ) = \frac { v _ { r } ( k ) \cos ( \phi ( k ) ) - v _ { t } ( k ) \sin ( \phi ( k ) ) } { v _ { r } ( k ) \sin ( \phi ( k ) ) + v _ { t } ( k ) \cos ( \phi ( k ) ) },$$

which implies the relation between ω k ( ) and v r ( k , v ) t ( k ) , i.e.

$$\omega ( k ) & = \frac { \cdot } { h } \arctan \left ( \frac { v _ { r } ( k ) \cos ( \phi ( k ) ) - v _ { t } ( k ) \sin ( \phi ( k ) ) } { v _ { r } ( k ) \sin ( \phi ( k ) ) + v _ { t } ( k ) \cos ( \phi ( k ) ) } \right ) ^ { \quad \text{Append} } ( 1 0 ) ^ { \quad \text{In coi} } \\ & \quad - \frac { 2 } { h } \theta ( k ). ^ { \quad \text{weightiri} } \text{importal} } \text{move a machine to a } \text{the } \text{the } \text{$i$} \, \text{$which} \, i$$

Moreover, the decomposition of v c shows that v r ( k h ) is exactly the projected distance of p ( k + 1) -p ( k ) onto the direction of p T -p ( k ) , i.e.,

$$v _ { r } ( k ) h = \left \| \frac { ( \mathbf p _ { T } - \mathbf p ( k ) ) ( \mathbf p _ { T } - \mathbf p ( k ) ) ^ { \prime } } { \| \mathbf p _ { T } - \mathbf p ( k ) \| ^ { 2 } } ( \mathbf p ( k + 1 ) - \mathbf p ( k ) ) \right \|, \quad \text{the estir} \\ \text{the is uncle} \\ \text{disense}$$

which is denoted as the bold dashed segment on the direction of p T -p ( k ) in Fig. 3. Similarly, v t ( k h ) can be represented as

$$\stackrel { \dots } { v _ { t } ( k ) h } = \left \| \left ( I _ { 2 } - \frac { ( \mathbf p _ { T } - \mathbf p ( k ) ) ( \mathbf p _ { T } - \mathbf p ( k ) ) ^ { \prime } } { \| \mathbf p _ { T } - \mathbf p ( k ) \| ^ { 2 } } \right ) ( \mathbf p ( k + 1 ) - \mathbf p ( k ) ) \right \|, \text{ different } } ^ { \text{lems with} } \text{ the soluti} \\ \dots \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \colon \end{cases}$$

which is denoted as the bold dashed segment on the direction perpendicular to p T -p ( k ) in Fig. 3. Jointly with (11), (12) and the triangular geometric relationship in Fig. 3, it follows that

$$& \sin ^ { 2 } ( \phi ( k + 1 ) - \phi ( k ) ) = \frac { v _ { t } ^ { 2 } ( k ) } { ( r ( k ) / h - v _ { r } ( k ) ) ^ { 2 } + v _ { t } ^ { 2 } ( k ) }, \quad ( 1 3 ) \quad \begin{array} { c } \text{$be seen$} \\ \text{$target $s\alpha$} \\ r ^ { 2 } ( k ) - r ^ { 2 } ( k + 1 ) = 2 r ( k ) v _ { r } ( k ) h - ( v _ { r } ^ { 2 } ( k ) + v _ { t } ^ { 2 } ( k ) ) h ^ { 2 }. \end{array} \quad \begin{array} { c } \text{$the seen$} \\ \text{$three target$} \\ \end{array} \quad \begin{array} { c } \text{$the target$} \\ \end{array} \quad \begin{array} { c } \text{$the target$} \\ \end{array} \quad \begin{array} { c } \text{$the target$} \\ \end{array} \quad \begin{array} { c } \text{$the target$} \\ \end{array} \quad \begin{array} { c } \text{$the target$} \\ \end{array}$$

Besides, it follows from (8) that

$$v _ { r } ^ { 2 } ( k ) + v _ { t } ^ { 2 } ( k ) = v _ { c } ^ { 2 }. \quad \quad \ ( 1 4 ) \ ^ { \ e s t i n }$$

Substituting (13) and (14) into the composite optimization problem (7) yields that

$$\maximize _ { v _ { r } } & f _ { e } ( v _ { r } ) + \alpha f _ { c } ( v _ { r } ) & \quad \text{and} \, i \\ \text{subject to} \, \frac { v _ { c } ^ { 2 } h } { 2 r } \leqslant v _ { r } \leqslant v _ { c }, & C. \ S t \\ \cdots \cdots \cdots \cdots c & \, \nu \, \cdots \cdots \cdots \cdots c \cdots & \, \nu..$$

where the dependence of v r ( k ) and r k ( ) on k is dropped for simplifying notations and

$$\cdot \cdot \cdot \ & \swarrow \\ f _ { e } ( v _ { r } ) & = \frac { ( v _ { c } ^ { 2 } - v _ { r } ^ { 2 } ) h ^ { 2 } } { \sigma ^ { 4 } r ^ { 2 } ( r ^ { 2 } - 2 r v _ { r } h + v _ { c } ^ { 2 } h ^ { 2 } ) ^ { 2 } }, \\ f _ { c } ( v _ { r } ) & = 2 r v _ { r } h - v _ { c } ^ { 2 } h ^ { 2 }.$$

,

2) Normalization: Observe that f e ( v r ) and f c ( v r ) are used to quantify the objectives with different physical meanings, which results in the lack of physical understanding of the weighting parameter α . Consequently, it is unclear how to provide an appropriate α to balance the estimation and the control objectives. To overcome this difficulty, we normalize the two objective functions in (16) to be dimensionless and with values in [0 , 1] . Then, the composite optimization problem (15) is rewritten in the following normalized form

$$\lim _ { v _ { r } } f ( v _ { r } ) \text{ subject to } \frac { v _ { c } ^ { 2 } h } { 2 r } \leqslant v _ { r } \leqslant v _ { c }, \quad \ \ ( 1 7 )$$

where the normalized objective function

$$f ( v _ { r } ) \ = \ \left ( 1 - \frac { ( v _ { c } ^ { 2 } - v _ { r } ^ { 2 } ) ( r ^ { 2 } - v _ { c } ^ { 2 } h ^ { 2 } ) ^ { 2 } } { v _ { c } ^ { 2 } ( r ^ { 2 } - 2 r v _ { r } h + v _ { c } ^ { 2 } h ^ { 2 } ) ^ { 2 } } \right ) \\ + \beta \, \frac { ( v _ { c } - v _ { r } ) ( v _ { c } ^ { 2 } h ^ { 2 } + r ^ { 2 } ) } { v _ { c } ( v _ { c } h - r ) ^ { 2 } } \quad ( 1 8 ) \\ \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \cdots \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \ \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \cdots \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad\, \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \,
 \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \. \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \quad \, \$$

introduces a new weighting parameter β ∈ [0 , + ∞ ) . See Appendix A for the detailed normalization process.

In contrast with the weighting parameter α in (15), the new weighting parameter β in (18) can directly quantify the relative importance of the control objective to the estimation objective, which is also an advantage of the normalized optimization problem (17) over those in [23]-[26] with simple weighted sums as the objective functions. To be more specific, β = 1 in (18) means the control objective has the same importance as the estimation objective, while the physical meaning of α = 1 is unclear. Therefore, the new weighting parameter β facilitates designers to emphasize the importance between the estimation and the control objectives.

Furthermore, solutions to the normalized optimization problems with different values of β lead the vehicle to move in different trajectories. It is interesting to note that if β glyph[greaterorequalslant] 4 , the solution of (17) results in the heading direction of the vehicle directly toward the estimate of the target position. In this case, the controller design ignores the estimation objective as in [18]-[20].

Remark 2. We set a constant β in this work. In fact, it can be seen as an optimization variable to be adjusted during the target search task, which, however, is beyond of the scope of this work. A heuristic idea is to design β k ( ) ∝ 1 /r k ( ) , since the target position may have already been well estimated when the vehicle is close to the target. Another idea is to measure estimation progress by using the following ratio

$$\frac { f _ { e } ( v _ { r } ( k ) ) ) } { f _ { e } ( v _ { r } ( k - 1 ) ) }$$

and increase β k ( ) if this ratio is close to unity, and decrease it otherwise.

## C. Solution to the normalized optimization problem (17)

Before providing the solution, we introduce an auxiliary variable

$$\rho = v _ { c } h / r$$

to denote the level of completing the target search task. Specifically, ρ &gt; 1 , i.e., v h &gt; r c , means that the vehicle's

traveling distance in a sampling period is strictly larger than its distance to the target. This implies that the vehicle is sufficiently close to the target. Hence, we only focus on the nontrivial case that 0 &lt; ρ glyph[lessorequalslant] 1 .

Lemma 1. The equation d d v r f ( v r ) = 0 has at least one root in the interval [ v , v s c ] , where v s = 2 rv 2 c h/ r ( 2 + v h 2 c 2 ) .

## Proof.

$$P r o o f \colon & \text{See Appendix B} \colon & \text{$\Delta$} \colon & \text{$^{\prime}\mathcal{ }^{\prime}$}$$

Let ρ c = (2 -√ 4 β -β 2 ) / β ( -2) and v z denote the largest root of the equation d d v r f ( v r ) = 0 , v r ∈ [ v , v s c ] by Lemma 1. Then, we provide the solution of (17) in the following theorem.

Theorem 1. The optimal solution v ∗ r of the normalized optimization problem (17) is given as follows.

$$\begin{smallmatrix} ( \mathbf a ) & I f \, \beta \in [ 4, \infty ), \, t h e n \, v _ { r } ^ { * } = v _ { c }. \\ ( \mathbf b ) & I f \, \beta \in [ 2, 4 ), \, t h e n \end{smallmatrix}$$

$$v _ { r } ^ { * } = \begin{cases} \ v _ { c }, & w h e n \ \rho \in ( 0, \rho _ { c } ], \\ \ v _ { z }, & w h e n \ \rho \in ( \rho _ { c }, 1 ]. \end{cases}$$

$$( \mathfrak { c } ) \, \underset { \Phi } { \, I f } \, \beta \in [ 1, 2 ), \, \underset { \Phi } { \, t h e n } \, v _ { r } ^ { * } = v _ { z }.$$

$$( \text{d} ) \ I f \, \beta \in [ 0, 1 ), \, t h e n$$

$$v _ { r } ^ { * } = \begin{cases} \ v _ { z }, & \text{when $f(v_{z}) < \beta$}, \\ \ v _ { s }, & \text{otherwise}. \end{cases}$$

Proof. See Appendix C.

glyph[squaresolid]

Note that v s glyph[lessorequalslant] v z glyph[lessorequalslant] v c . Theorem 1 shows that a larger β results in a larger v ∗ r . As we discussed in Section III-B that a larger β emphasizes more the importance of the control objective, a larger radial velocity v r is derived to approach the target faster. This is also consistent with our intuition. Specifically, if β ∈ [4 , ∞ ) , then v ∗ r directly takes its maximum value v c , which implies that the vehicle almost ignores the estimation objective.

In addition, observe that finding v z is not needed in all cases in Theorem 1. This potentially saves computational cost, compared to the works in [23]-[26], [28].

## IV. THE OPTIMIZATION-BASED CONTROLLER

It should be noted from Theorem 1 and (10) that the vehicle's optimal angular velocity depends on the range r k ( ) and the bearing angle φ k ( ) , both of which are unknown and need to be estimated. Accordingly, we design a recursive least-square estimator, and utilize the certainty-equivalence principle for the controller design.

## A. The recursive least-square estimator

Observing that the unknown r k ( ) and φ k ( ) are both originated from the unknown position p T , we only need to obtain an optimal estimate of p T , which is the solution of the following nonlinear least-square minimization problem

$$\min _ { \mathbf P _ { T } } \, \frac { 1 } { 2 } \sum _ { i = 1 } ^ { k } \left \| m ( i ) - \arctan \left ( \frac { x ( i ) - x _ { T } } { y ( i ) - y _ { T } } \right ) \right \| ^ { 2 }. \quad ( 1 9 ) \quad \text{vehicle} ^ { \prime } \\ \text{consider}$$

Similar to Stansfield estimator in [29], we also assume the measurement noise e k ( ) is sufficiently small to justify the replacement of m i ( ) -arctan(( x i ( ) -x T ) / y i ( ( ) -y T )) with sin ( m i ( ) -arctan(( x i ( ) -x T ) / y i ( ( ) -y T )) ) . Using the relation that

$$\ n e \ r o o t \quad & \sin \left ( m ( i ) - \arctan \left ( \frac { x ( i ) - x _ { T } } { y ( i ) - y _ { T } } \right ) \right ) \\ \Lambda \quad & = \frac { 1 } { r ( i ) } \left ( \sin ( m ( i ) ) ( y ( i ) - y _ { T } ) - \cos ( m ( i ) ) ( x ( i ) - x _ { T } ) \right ),$$

we approximate the optimization problem (19) as

$$& [ v _ { s }, v _ { c } ] \\ & \text{in the} \quad \minimize \ \frac { 1 } { 2 } \sum _ { i = 1 } ^ { k } \| x ( i ) \cos ( m ( i ) ) - y ( i ) ( \sin ( m ( i ) ) \\ & d \ o p t i \text{.} & & - [ \cos ( m ( i ) ), - \sin ( m ( i ) ) ] [ x _ { T }, y _ { T } ] ^ { \prime } \| ^ { 2 },$$

where we ignore the rather weak effect of 1 /r i ( ) in (20).

Since the optimization problem (21) is a standard leastsquare problem, we design a recursive algorithm to solve it, see Algorithm 1. In practice, the initial estimate of the target position ̂ p T (1) in Algorithm 1 is given as a solution of the following equations

$$\tan ( m ( i ) ) = \frac { x ( i ) - x _ { T } } { y ( i ) - y _ { T } }, i = 0, 1. \quad \quad ( 2 2 )$$

## Algorithm 1 Least-square estimate

- 1: At time k = 1 , set ̂ p T (1) as the solution of (22), and Q (1) = I 2 .
- 2: At time k glyph[greaterorequalslant] 2 , the vehicle takes a bearing angle m k ( ) and computes

$$\text{ and computes} \quad & \text{$H=[-cos(m(k))$,sin(m(k))]^{\prime}$,} \\ \text{$K=Q(k-1)H(H^{\primeQ(k-1)H+1)^{-1}$,} \\ e _ { T } = \text{$H^{\primep(k)}$,} \\ \text{$Q(k)=(I_{2}-KH^{\prime})Q(k-1),} \\ \widehat { \mathbf p } _ { T } ( k ) = \widehat { \mathbf p } _ { T } ( k - 1 ) + \text{$K(e_{T}-H^{\prime}\widehat{ \mathbf p }_{T}(k-1))$.}$$

Note that Algorithm 1 is developed under the assumption of small e k ( ) and ignorance of 1 /r k ( ) , and might be biased. Nonetheless, it is verified by the simulation results in [33] to perform well with a small estimation error and to be of practical significance.

## B. The controller for bearing-only target search

Now, the detailed controller providing the angular velocity is shown in Algorithm 2.

## V. CONTROL FOR TARGET SEARCH WITHOUT GPS INFORMATION

Observe that the controller in Algorithm 2 requires the vehicle's position p ( k ) in the global coordinate, which is unavailable in GPS-denied environments. In this section, we consider the target search problem without the vehicle's global position, and utilize the local bearing angles defined in the vehicle's coordinate to solve it.

## Algorithm 2 The angular velocity for the vehicle

- 1: At time k = 0 , the vehicle takes a bearing angle m k ( ) , and randomly selects ω k ( ) .
- 2: At time k glyph[greaterorequalslant] 1 , the vehicle takes a bearing angle m k ( ) and completes the following steps:
- (a) runs Algorithm 1 to obtain

$$\widehat { \mathbf p } _ { T } ( k ) = [ \widehat { x } _ { T } ( k ), \widehat { y } _ { T } ( k ) ] ^ { \prime },$$

- (b) computes both the bearing and range estimates ̂ φ k ( ) = arctan ( ( x k ( ) -̂ x T ( k )) / y k ( ( ) -̂ y T ( k )) ) ,

$$\widehat { r } ( k ) = \| \mathbf p ( k ) - \widehat { \mathbf p } _ { T } ( k ) \|, \\ \intertext { \i \ m e c \ T h o o r a m \ 1 \ t o \ c o m u t a }$$

- (d) computes the angular velocity in (10) by replacing φ k ( ) and v r ( k ) by ̂ φ k ( ) and v ∗ r ( k ) , respectively.
- (c) uses Theorem 1 to compute v ∗ r ( k ) where r k ( ) is replaced by its estimate ̂ r k ( ) ,

Fig. 4. The bearing-only target search problem in the vehicle's coordinate.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-857be009-cda6-4896-bc00-d26572a5473c.jpg)

## A. The target search problem without GPS information

Without the vehicle's global position p ( k ) , we consider the target search problem in the vehicle's coordinate Σ l . As Fig. 4 shows, this coordinate is attached to the vehicle with its positive x-axis coincident with the heading of the vehicle. The target's position in this coordinate is denoted as p l T ( k ) = [ x l T ( k , y ) l T ( k )] ′ at the time step k , which is time-varying due to the moving vehicle's coordinate. The orientation θ k ( ) is still accessible by a compass. Moreover, the vehicle measures a local bearing angle m k l ( ) in the vehicle's coordinate relative to the target, i.e.,

$$m ^ { l } ( k ) = \phi ^ { l } ( k ) + e ^ { l } ( k ),$$

where φ l ( k ) = arctan( y l T ( k /x ) l T ( k )) denotes the true local bearing angle, and e l ( k ) ∼ N (0 , σ 2 ) denotes the local measurement noise.

Then, the bearing-only target search problem becomes to sequentially design the angular velocity ω k ( ) for the vehicle (2) to approach the unknwon target position as fast as possible only by local noisy bearing angles. To solve this problem, we still design the controller depending on the solution to the bi-objective optimization problem (3).

## B. The optimization-based controller without GPS information

Similar to Section IV, a recursive least-square algorithm is designed in the proposed controller to estimate the target position p l T ( k ) . We approximately solve the nonlinear leastsquare optimization problem

$$\underset { \mathbf p _ { T } ^ { l } ( k ) } { \min } \, \frac { 1 } { 2 } \sum _ { i = 1 } ^ { k } \| m ^ { l } ( i ) - \arctan \left ( \frac { x _ { T } ^ { l } ( i ) } { y _ { T } ^ { l } ( i ) } \right ) \| ^ { 2 }$$

by a linear least-square form, i.e.,

$$& \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb{v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \text{$\mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb {v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } /\mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb{ v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } /\text{$\mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } /\mathbb {v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb {v } / \text{$\mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb {v } / \mathbb {v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \text{$\mathbb {v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb{ v } / \text{$\mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb{ v } / \mathbb {v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } /\mathbb { v } / \text{$\mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } /\mathbb { v } / \mathbb {v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb { v } / \mathbb {v } / \mathbb { v } / \extractor$$

The difference from Section IV lies in the time-varying target position in the vehicle's coordinate due to the moving vehicle. Thus, the state equation of the target position p l T ( k ) is to be derived. According to the relation between global coordinate and the vehicle's coordinate, it holds that

$$\mathbf p _ { T } ^ { l } ( k ) = \begin{bmatrix} \cos ( \theta ( k ) ) & \sin ( \theta ( k ) ) \\ - \sin ( \theta ( k ) ) & \cos ( \theta ( k ) ) \end{bmatrix} ( \mathbf p _ { T } - \mathbf p ( k ) ). \quad ( 2 5 )$$

Combining with the vehicle dynamics, we further obtain that

$$\mathbf p _ { T } ^ { l } ( k ) = \mathbf A ( k - 1 ) \mathbf p _ { T } ^ { l } ( k - 1 ) + \mathbf b, \quad \ \ ( 2 6 )$$

where

$$\mathbf A ( k - 1 ) = \begin{bmatrix} 1 & \omega ( k - 1 ) \\ - \omega ( k - 1 ) & 1 \end{bmatrix}, \mathbf b = [ - v _ { c } h, 0 ] ^ { \prime }.$$

Jointly with (24) and (26), the recursive least-square estimator with only local bearing angles is given in Algorithm 3, where the initial estimate of the target position in the vehicle's coordinate ̂ p l T (1) is selected as a solution of the following equations

$$\mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak{ n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cd� \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdot \mathfrak { n } \cdoplus \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfra \cdot \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfring { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n } \mathfrak { n }$$

Algorithm 3 Least-square estimation algorithm with local bearing angles

- 1: At time k = 1 , set ̂ p l T (1) as the solution of (27), and Q (1) = I 2 .
- 2: At time k glyph[greaterorequalslant] 2 , the vehicle takes a local bearing angle m k l ( ) and computes

$$\stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim  \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel | \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel. } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stack { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim } \stackrel { \sim }$$

The angular velocity for the vehicle is also given in Algorithm 2 where we use Algorithm 3 to replace Algorithm 1 and the bearing angle and the range estimates are computed as

$$\mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ {\quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \mu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \ \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \\ \mu _ { \quad } \widehat { r } ( k ) & = \frac { \pi } { 2 } - \mu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \ :: } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ _ {\quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nolimits _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \cdot } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad } \nu _ { \quad }$$

Fig. 5. The trajectory of the vehicle and the estimated positions of the target. The error between the vehicle and the target is around 0.41m.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-b5caecb0-3a50-4ac1-a1cd-e62ae101858a.jpg)

and the rest of Algorithm 2 remain unchanged. Thus, we obtain the controller for the target search problem without GPS information.

## VI. SIMULATIONS

In this section, we illustrate the effectiveness of the proposed optimization-based controllers for the bearing-only target search problem. The initial state of the vehicle is set as p (0) = [0 0] , ′ , θ (0) = 0 , the constant forward velocity v c = 4 m/s, the target position p T = [100 100] , ′ m, and its initial estimate ̂ p T (0) = [40 80] , ′ m. We take the weighting parameter β = 1 unless particular specification, which means the control objective has the same importance as the estimation objective. Moreover, the sampling period is set as h = 0 25 . s.

## A. Simulation results with GPS information

With GPS information, the controller in Algorithm 1 is implemented in a target search task to approach an unknown target position. Fig. 5 shows the vehicle's trajectory and the estimated target positions, where the small square in the upper left corner is the magnification of the results in the target's neighbor. Clearly, the vehicle moves in a spiral curve and approaches the true target gradually, while the estimated target positions converges to its true position. It reveals that the vehicle achieves the estimation and the control objectives simultaneously with the weighting parameter β = 1 .

The cases of β = 0 and β = 2 are depicted in Fig. 6, which shows that the vehicle tends to circle the target for the estimation objective and to directly approach the target for the control objective, respectively. This also explains the rationality of the spiral trajectory in Fig. 5, as β = 1 is between these two cases. The videos of the vehicle's motions under the proposed controller in the cases of β = 0 , β = 1 , and β = 2 are uploaded to https://www.dropbox.com/s/9dtpijap65gup50/ video motion.mp4?dl=0.

## B. Simulation results without GPS information

For the target search problem without the vehicle's global position p ( k ) , we apply the controller in Algorithm 3 with local bearing angles. The results are shown in Fig. 7, where both

Fig. 6. The trajectories of the vehicle and the estimated positions of the target with β = 0 and β = 2 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-3c7e1d4f-06d7-4de9-8184-cab1e453ebaf.jpg)

Fig. 7. The trajectory of the vehicle and the estimated positions of the target with local bearing angles. The final error between the vehicle and the target is around 0.43m.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-2177e165-a4e3-4af3-a2cd-d6c80e12788e.jpg)

of the vehicle's trajectory and the estimated target positions gradually converge to the true target position. It is consistent with our results in Section V. Since the true target position is time-varying in the vehicle's coordinate, the estimated target positions have larger fluctuation than those in Section VI-A.

In addition, we compare the two recursive least-square estimators of Algorithm 1 and 3 by denoting the estimation error as e est ( k ) = ‖ p T ( k ) -̂ p T ( k ) ‖ . We observe that the e est ( k ) of Algorithm 3 converges to zero as that of Algorithm 1, although its value is larger in the initial stage and the convergence rate is slower due to lack of the vehicle's GPS information. It implies that the GPS information is not necessary, but might improve the search performance, which is not studied in previous works [10]-[14].

## C. The role of the weighting parameter β

In this subsection, we investigate the role of the weighting parameter β in balancing the estimation and the control objectives, and provide some guidelines for choosing an appropriate β .

In Fig. 9, the time for the estimation error less than 1 m are t = 51 s, t = 64 s, and t = 145 s for the cases of β = 1 , β = 2 , and β = 3 . From Fig. 10, the time of estimated

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-21fd4384-6faf-43e8-8d6a-51e423bf72b9.jpg)

Fig. 8. The estimation errors of the target position with global and local bearing angles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-bb2c1773-d8e9-4af8-a9e6-dc54630be49c.jpg)

Fig. 9. The estimation errors of the target position with different values of β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-edb300cb-ca99-4a36-ac57-6c90e278fe18.jpg)

Fig. 10. The ranges between the vehicle and the estimated target with different values of β .

Fig. 11. The trajectory of the vehicle and the estimated positions of the target with β = 3 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-8bc426f8-37be-4c37-8244-3b79f4e42470.jpg)

Fig. 12. The time for target search with different values of β .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-c759e668-9137-44cb-88d2-6fc66baa5fbe.jpg)

range in Algorithm 2 less than 1 m are t = 90 s, t = 71 s, and t = 28 s. That is, the larger the value of β , the slower the estimation convergence and the faster the vehicle approaching the estimated target position. In particular, we can observe from Fig. 10 that the estimated range easily tends to zero for a large β , i.e., the vehicle tends to approach the estimated target position. This is consistent with Section III-B that a larger β emphasizes more the importance of the control objective. However, the estimate of the target position may be far from the true target position in the case of a larger β . Figure 11 further shows that the vehicle approaches a 'wrong' target due to the large estimation error. Thus, the weighting parameter β must be chosen carefully to balance the estimation and the control objectives.

To further explore the search performance with different values of β , we set a terminal condition of Algorithm 2 as r k ( ) &lt; v h c , indicating the completeness of the target search task. The times for completing the target search task with different values of β are given in Fig. 12, where the search time initially decreases with a larger β , and then increases as β become much larger. Finally, it keeps constant when β glyph[greaterorequalslant] 4 5 . . This implies that neither too small or too large value of β is preferable in the target search problem.

## VII. CONCLUSION

We have derived an optimization-based controller for the discrete-time Dubins vehicle to address the bearing-only target search problem, where the vehicle was to be steered to the unknown target position as fast as possible and the estimation performance for the unknown target position was considered. In the proposed scheme, we estimated the unknown target position via the recursive estimator with bearing angles in global and the vehicle's coordinate. Using the estimated target position, we formulated and solved a bi-objective optimization problem with the estimation and the control objectives. In particular, the solution was in terms of the predefined importance of the control objective relative to the estimation objective. Then, an optimization-based controller was obtained and drove the vehicle to approach the target as fast as possible while optimizing the estimation performance. Finally, simulations were given to verify the effectiveness of the proposed controller.

## APPENDIX A

NORMALIZATION OF THE OPTIMIZATION PROBLEM (15)

In the appendix, we denote ˙ ( f · ) , f ¨ ( · ) as the derivative and second derivative of the function f ( ) · with respect to the variable in the bracket.

To normalize the objective function in (15), we take the derivative of f e ( v r ) and f c ( v r ) in (16) as

$$\dot { f } _ { e } ( v _ { r } ) = \frac { - 2 v _ { r } h ^ { 2 } ( r ^ { 2 } - 2 r v _ { r } h + v _ { c } ^ { 2 } h ^ { 2 } ) + 4 r h ^ { 3 } ( v _ { c } ^ { 2 } - v _ { r } ^ { 2 } ) } { \sigma ^ { 4 } r ^ { 2 } ( r ^ { 2 } - 2 r v _ { r } h + v _ { c } ^ { 2 } h ^ { 2 } ) ^ { 3 } }, \quad \begin{array} { c } \text{obs} \\ \text{i.e.,} \\ \dot { f } ( v \end{array}.$$

˙ f c ( v r ) = 2 rh.

Let ˙ f e ( v r ) = 0 , and we obtain a stationary point of f e ( v r ) denoted as v s = 2 rv 2 c h/ r ( 2 + v h 2 c 2 ) . Thus, the increasing interval and decreasing interval of f e ( v r ) are [ v h/ r, v 2 c 2 s ) and [ v , v s c ] , respectively, and f c ( v r ) is always increasing in the feasible set [ v h/ r, v 2 c 2 c ] . It follows that the Pareto-optimal set of the composite optimization problem (15) is [ v , v s c ] .

Instituting v r = v s into f e ( v r ) and f c ( v r ) , we obtain the ideal value of f e ( v r ) and the nadir value of f c ( v r ) below

$$f _ { e } ^ { * } = \frac { v _ { c } ^ { 2 } h ^ { 2 } } { r ^ { 2 } ( r ^ { 2 } - v _ { c } ^ { 2 } h ^ { 2 } ) ^ { 2 } }, f _ { c } ^ { n a d } = \frac { 3 r ^ { 2 } v _ { c } ^ { 2 } h ^ { 2 } - v _ { c } ^ { 4 } h ^ { 4 } } { r ^ { 2 } + v _ { c } ^ { 2 } h ^ { 2 } }. \quad ( 2 8 ) \quad \text{implyir}$$

Similarly, taking v r = v c , we obtain the nadir value of f e ( v r ) and the ideal value of f c ( v r ) below

$$f _ { e } ^ { \text{nad} } = 0, f _ { c } ^ { * } = 2 r v _ { c } h - v _ { c } ^ { 2 } h ^ { 2 }. \text{ \quad \ \ } ( 2 9 ) \text{ \ \ } \text{$name}.$$

Jointly with (16), (28) and (29), we can normalize the objective functions as

$$\text{number} \, \infty \, \infty \\ f _ { e } ^$$

Then the normalized objective function is given as

$$f ( v _ { r } ) = f _ { e } ^ { \text{norm} }$$

with a new weighting parameter β ∈ [0 , + ∞ ) .

## APPENDIX B

## PROOF OF LEMMA 1

The derivative of f ( v r ) in (17) is computed as

$$\mathfrak { s$$

Let ˙ ( f v r ) = 0 , and we have three solutions denoted as z , z 1 2 , and z 3 . Since we only focus on the interval v r ∈ [ v , v s c ] , the first derivatives of f ( v r ) at v s and v c are computed as follows

$$\intertext { \text{nce} } _ { \text{oe} } \quad \dot$$

Thus, v s is one of the roots of ˙ ( f v r ) = 0 , which implies that ˙ ( f v r ) = 0 has at least one root in the interval v r ∈ [ v , v s c ] .

## APPENDIX C

## PROOF OF THEOREM 1

Obviously, solutions of the normalized optimization problem (17) exist in the Pareto-optimal set, thus we only need to focus on f ( v r ) for v r ∈ [ v , v s c ] . Based on the proof of Lemma 1, the first derivative of ˙ ( f v c ) is indefinite, the sign of which is to be analyzed in the following.

To this end, we firstly consider the continuity of ˙ ( f v r ) and observe that there exists and only exists a discontinuity point, i.e., v r = ( r 2 + v h 2 c 2 ) / 2 rh , which makes the denominator of ˙ ( f v r ) equal to zero. Furthermore, this discontinuity point is larger than v c , implying that the function ˙ ( f v r ) is continuous in the Pareto set [ v , v s c ] .

Then, we take the second derivative of f ( v r ) in v r ∈ [ v , v s c ] and obtain that

$$\inf _ { \substack { \iota } } \, \text{and} \\ \$$

where Θ = r 2 -v h 2 c 2 and Ξ = r 2 -2 v hr r + v h 2 c 2 . Inserting the bounds of the Pareto-optimal set v , v s c into the second derivative yields that

$$\ddot { f } ( v _ { s } ) > 0, \ddot { f } ( v _ { c } ) > 0,$$

implying that the function ˙ ( f v r ) is increasing at v s and v c .

Now, we focus on the sign of

˙ (

f

v

c

)

.

Since the denominator is positive in

[

v , v

s

c

]

,

we only needs to analyze the sign of the numerator

(2

-

β

)(

v h

2

2

+

r

2

) + 4

rv h

c

.

In light of the fundamental inequality a 2 + b 2 glyph[greaterorequalslant] 2 ab, ∀ a, b, we obtain that β glyph[greaterorequalslant] 4 is a sufficient condition of ˙ ( f v c ) glyph[lessorequalslant] 0 . The upper subfigure in Fig.13 approximately depicts ˙ ( f v r ) in this case, where ˙ ( f v r ) is continuous in [ v , v s c ] , ˙ ( f v s ) = 0 , ˙ ( f v c ) glyph[lessorequalslant] 0 , and (33) holds. Then the function curve of f ( v r ) can be approximated corresponding to that of ˙ ( f v r ) , where f ( v s ) -f ( v c ) = β -1 &gt; 0 and f ( v c ) = 1 . Therefore, the minimum of f ( v r ) is f ( v c ) .

c

Similarly, we have that β glyph[lessorequalslant] 2 is a sufficient condition of ˙ ( f v c ) glyph[greaterorequalslant] 0 . There are two cases depicted in the upper subfigures in Fig. 14 and Fig. 15. In the first case, the other two zero points besides v s exist in ( v , v s c ) , and we denote the largest zero point as v z . Since f ( v s ) -f ( v c ) = β -1 ,

Fig. 13. A possible case of ˙ ( f v r ) and f ( v r ) curves.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-a24a8d12-875d-44ab-98b7-0a51f7643ed3.jpg)

it holds that f ( v z ) is the minimum if β glyph[greaterorequalslant] 1 . For β &lt; 1 , the minimum is the smaller one of f ( v s ) and f ( v z ) . Fig. 14 depicts the function curves of ˙ ( f v r ) and f ( v r ) , where ˙ ( f v s ) = 0 , the other two zero points of ˙ ( f v r ) are in [ v , v s c ] , and the corresponding f ( v r ) are marked with β glyph[greaterorequalslant] 1 and β &lt; 1 , respectively. In the second case, the function ˙ ( f v r ) remains positive in ( v , v s c ) , since there is no zero points. It implies that f ( v r ) always increases and β &lt; 1 , as Fig. 15 shows. Therefore, the minimum of f ( v r ) is f ( v s ) .

For the case of β ∈ (2 4) , , it clearly holds that f ( v s ) -f ( v c ) = β -1 &gt; 0 . Let (2 -β )( v h 2 c 2 + r 2 ) + 4 rv h c = 0 , and we obtain the following relation that

$$\rho = \frac { v _ { c } h } { r } = \frac { 2 - \sqrt { 4 \beta - \beta ^ { 2 } } } { \beta - 2 }.$$

It follows that if ρ ∈ ( 0 (2 , -√ 4 β -β 2 ) / β ( -2) ] , then ˙ ( f v c ) glyph[lessorequalslant] 0 . The function curves are shown in Fig.13. Otherwise, the function ˙ ( f v c ) &gt; 0 , and the function curves are similar to that in Fig. 14 marked by β ≥ 1 .

To summarize, we have following conclusions with different values of β .

- (a) If β ∈ [4 , ∞ ) , then ˙ ( f v c ) glyph[lessorequalslant] 0 . It follows from f ( v s ) -f ( v c ) = β -1 &gt; 0 that the minimum of f ( v r ) is f ( v c ) . The function curves are shown in Fig. 13.
- (b) If β ∈ [2 , 4) , together with f ( v s ) -f ( v c ) = β -1 &gt; 0 , there are two conditions that i) If ρ ∈ ( 0 (2 , -√ 4 β -β 2 ) / β ( -2) ] , then ˙ ( f v c ) glyph[lessorequalslant] 0 . Thus, the minimum of f ( v r ) is f ( v c ) . The function curves are shown in Fig. 13.
- ii) If ρ ∈ ( (2 -√ 4 β -β 2 ) / β ( -2) 1 , ) , then ˙ ( f v c ) &gt; 0 . Thus, the minimum of f ( v r ) is f ( v z ) , where v z is the largest root of ˙ ( f v r ) = 0 , v r ∈ [ v , v s c ] . The function curves are shown in Fig. 14.
- (c) If β ∈ [1 , 2) , then ˙ ( f v c ) &gt; 0 . It follows from f ( v s ) -f ( v c ) = β -1 glyph[greaterorequalslant] 0 that the minimum of f ( v r ) is f ( v z ) . The function curves are shown in Fig. 14.
- (d) If β ∈ [0 , 1) , then ˙ ( f v c ) &gt; 0 . Together with f ( v s ) -f ( v c ) = β -1 &lt; 0 , there are two conditions that i) If there is any root of ˙ ( f v r ) = 0 in the interval ( v , v s c ] ,

Fig. 14. A possible case of ˙ ( f v r ) and f ( v r ) curves.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-758c033b-d27c-479b-9891-135c9017e40c.jpg)

√

√

Fig. 15. A possible case of ˙ ( f v r ) and f ( v r ) curves.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-72480de4-2593-4d83-b8de-8e3da12adec5.jpg)

then the minimum of f ( v r ) is the smaller one of f ( v z ) and f ( v s ) . The function curves are shown in Fig. 14. ii) If there is no root of ˙ ( f v r ) = 0 in the interval ( v , v s c ] , then the minimum of f ( v r ) is f ( v s ) . The function curves are shown in Fig. 15.

The optimal solution v ∗ r of the normalized optimization problem (17) can be selected as in Theorem 1.

## REFERENCES

- [1] L. W. Chen, J. H. Cheng, and Y. C. Tseng, 'Optimal path planning with spatial-temporal mobility modeling for individual-based emergency guiding,' IEEE Transactions on Systems Man &amp; Cybernetics: Systems , vol. 45, no. 12, pp. 1491-1501, 2015.
- [2] R. Raja, A. Dutta, and K. Venkatesh, 'New potential field method for rough terrain path planning using genetic algorithm for a 6-wheel rover,' Robotics and Autonomous Systems , vol. 72, pp. 295-306, 2015.
- [3] G. Best, W. Martens, and R. Fitch, 'Path planning with spatiotemporal optimal stopping for stochastic mission monitoring,' IEEE Transactions on Robotics , vol. 33, no. 3, pp. 629-646, 2017.
- [4] P. Lanillos, E. Besada-Portas, G. Pajares, and J. J. Ruz, 'Minimum time search for lost targets using cross entropy optimization,' in International Conference on Intelligent Robots and Systems . IEEE, 2012, pp. 602609.

- [5] J. Song and S. Gupta, ' ε glyph[star] : An online coverage path planning algorithm,' IEEE Transactions on Robotics , vol. 34, no. 2, pp. 526-533, 2018.
- [6] X. C. Ding, A. R. Rahmani, and M. Egerstedt, 'Multi-UAV convoy protection: An optimal approach to path planning and coordination,' IEEE Transactions on Robotics , vol. 26, no. 2, pp. 256-268, 2010.
- [7] B. Zhang and S. Liu, 'Planning expected-time optimal paths for target search by robot,' in World Congress on Intelligent Control and Automation . IEEE, 2012, pp. 3881-3886.
- [8] M. Otte and N. Correll, 'C-FOREST: Parallel shortest path planning with superlinear speedup,' IEEE Transactions on Robotics , vol. 29, no. 3, pp. 798-806, 2013.
- [9] S. Koziol, R. Wunderlich, J. Hasler, and M. Stilman, 'Single-objective path planning for autonomous robots using reconfigurable analog VLSI,' IEEE Transactions on Systems Man &amp; Cybernetics: Systems , vol. 47, no. 7, pp. 1301-1314, 2017.
- [10] A. N. Bishop, B. Fidan, B. D. Anderson, K. Dogancay, and P. N. Pathirana, 'Optimality analysis of sensor-target geometries in passive localization: Part 1 Bearing-only localization,' in International Conference on Intelligent Sensors, Sensor Networks and Information . IEEE, 2007, pp. 7-12.
- [11] A. Mavrommati, E. Tzorakoleftherakis, I. Abraham, and T. D. Murphey, 'Real-time area coverage and target localization using receding-horizon ergodic exploration,' IEEE Transactions on Robotics , vol. 34, no. 1, pp. 62-80, 2018.
- [12] D. Moreno-Salinas, A. M. Pascoal, and J. Aranda, 'Optimal sensor placement for multiple target positioning with range-only measurements in two-dimensional scenarios,' Sensors , vol. 13, no. 8, p. 10674, 2013.
- [13] K. Dogancay, 'UAV path planning for passive emitter localization,' IEEE Transactions on Aerospace and Electronic systems , vol. 48, no. 2, pp. 1150-1166, 2012.
- [14] W. Meng, L. Xie, and W. Xiao, 'Communication aware optimal sensor motion coordination for source localization,' IEEE Transactions on Instrumentation and Measurement , vol. 65, no. 11, pp. 2505-2514, 2016.
- [15] F. Lorussi, A. Marigo, and A. Bicchi, 'Optimal exploratory paths for a mobile rover,' in International Conference on Robotics and Automation . IEEE, 2001, pp. 2078-2083.
- [16] B. T. Hinson, M. K. Binder, and K. A. Morgansen, 'Path planning to optimize observability in a planar uniform flow field,' in American Control Conference . IEEE, 2013, pp. 1392-1399.
- [17] M. Cognetti, P. Salaris, and P. R. Giordano, 'Optimal active sensing with process and measurement noise,' in International Conference on Robotics and Automation . IEEE, 2018, pp. 2118-2125.
- [18] M. Deghat, I. Shames, B. D. Anderson, and C. Yu, 'Localization and circumnavigation of a slowly moving target using bearing measurements,' IEEE Transactions on Automatic Control , vol. 59, no. 8, pp. 2182-2188, 2014.
- [19] S. G¨ uler, s ¸FidanBar, S. Dasgupta, B. D. Anderson, and I. Shames, 'Adaptive source localization based station keeping of autonomous vehicles,' IEEE Transactions on Automatic Control , vol. 62, no. 7, pp. 3122-3135, 2017.
- [20] Y. Liu, Q. Wang, H. Hu, and Y. He, 'A novel real-time moving target tracking and path planning system for a quadrotor UAV in unknown unstructured outdoor scenes,' IEEE Transactions on Systems Man &amp; Cybernetics: Systems , no. 99, pp. 1-11, 2018.
- [21] A. N. Bishop, B. Fidan, B. D. Anderson, K. Do˘ ganc ¸ay, and P. N. Pathirana, 'Optimality analysis of sensor-target localization geometries,' Automatica , vol. 46, no. 3, pp. 479-492, 2010.
- [22] R. T. Marler and J. S. Arora, 'The weighted sum method for multiobjective optimization: new insights,' Structural and multidisciplinary optimization , vol. 41, no. 6, pp. 853-862, 2010.
- [23] H. Ergezer and K. Leblebicioglu, 'Path planning for UAVs for maximum information collection,' IEEE Transactions on Aerospace &amp; Electronic Systems , vol. 49, no. 5, pp. 502-520, 2013.
- [24] H. Duan and P. Qiao, 'Pigeon-inspired optimization: a new swarm intelligence optimizer for air robot path planning,' International Journal of Intelligent Computing and Cybernetics , vol. 7, no. 1, pp. 24-37, 2014.
- [25] D. Sidoti, G. V. Avvari, M. Mishra, L. Zhang, B. K. Nadella, J. E. Peak, J. A. Hansen, and K. R. Pattipati, 'A multiobjective pathplanning algorithm with time windows for asset routing in a dynamic weather-impacted environment,' IEEE Transactions on Systems Man &amp; Cybernetics: Systems , vol. 47, no. 12, pp. 3256-3271, 2017.
- [26] I. Maurovi´ c, M. Seder, K. Lenac, and I. Petrovi´ c, 'Path planning for active SLAM based on the D* algorithm with negative edge weights,' IEEE Transactions on Systems Man &amp; Cybernetics: Systems , vol. 48, no. 8, pp. 1321-1331, 2018.
- [27] S. J. Liu and P. P. Zhang, 'Discrete-time stochastic source seeking via tuning of forward velocity,' IFAC-PapersOnLine , vol. 48, no. 28, pp. 927-932, 2015.
- [28] S. Candido and S. Hutchinson, 'Minimum uncertainty robot navigation using information-guided pomdp planning,' in International Conference on Robotics and Automation . IEEE, 2011, pp. 6102-6108.
- [29] M. Gavish and A. J. Weiss, 'Performance analysis of bearing-only target location algorithms,' IEEE Transactions on Aerospace and Electronic Systems , vol. 28, no. 3, pp. 817-828, 1992.
- [30] C. Lin, Z. Lin, R. Zheng, G. Yan, and G. Mao, 'Distributed source localization of multi-agent systems with bearing angle measurements,' IEEE Transactions on Automatic Control , vol. 61, no. 4, pp. 1105-1110, 2016.
- [31] R. Zheng, Y. Liu, and D. Sun, 'Enclosing a target by nonholonomic mobile robots with bearing-only measurements,' Automatica , vol. 53, pp. 400-407, 2015.
- [32] H. L. V. Trees, Detection, Estimation, and Modulation Theory: Detection, Estimation, and Linear Modulation Theory . Publishing House of Electronics Industry, 2003.
- [33] K. Do˘ ganc ¸ay, 'On the bias of linear least squares algorithms for passive target localization,' Signal Processing , vol. 84, no. 3, pp. 475-486, 2004.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-e6b28b60-0efd-4692-ae2c-69e1034a40c5.jpg)

Zhuo Li received the B.S. degree from the Department of Automation, Harbin Institute of Technology, Harbin, China, in 2016. She is currently pursuing the Ph.D. degree at the Department of Automation, Tsinghua University, Beijing, China. Her research interests include nonlinear control, distributed optimization, source seeking and learning.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-0b85897c-38b3-4138-b3e4-71b049a436a2.jpg)

Keyou You (SM'17) received the B.S. degree in Statistical Science from Sun Yat-sen University, Guangzhou, China, in 2007 and the Ph.D. degree in Electrical and Electronic Engineering from Nanyang Technological University (NTU), Singapore, in 2012. After briefly working as a Research Fellow at NTU, he joined Tsinghua University in Beijing, China where he is now a tenured Associate Professor in the Department of Automation. He held visiting positions at Politecnico di Torino, The Hong Kong University of Science and Technology, The

University of Melbourne and etc. His current research interests include networked control systems, distributed optimization and learning, and their applications.

- Dr. You received the Guan Zhaozhi award at the 29th Chinese Control Conference in 2010, the CSC-IBM China Faculty Award in 2014 and the ACA (Asian Control Association) Temasek Young Educator Award in 2019. He was selected to the National 1000-Youth Talent Program of China in 2014 and received the National Science Fund for Excellent Young Scholars in 2017.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-6454b41f-6e41-49d5-b100-af414087985e.jpg)

Shiji Song received the Ph.D. degree in the Department of Mathematics from Harbin Institute of Technology in 1996. He is a professor in the Department of Automation, Tsinghua University. His research interests include system modeling, control and optimization, computational intelligence and pattern recognition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00380v1-figure-05159c3d-5add-4605-a566-200da6492bd4.jpg)

Anke Xue received the Ph.D. degree from Zhejiang University, in 1997. He is a professor in Hangzhou Dianzi University. His current research interests include robust control theory and applications.
Zhuo Li Keyou You Shiji Song Anke Xue
Optimization-based Control for Bearing-only Target Search with a Mobile Vehicle
This work aims to design an optimization-based controller for a discrete-time
Dubins vehicle to approach a target with unknown position as fast as possible
by only using bearing measurements. To this end, we propose a bi-objective
optimization problem, which jointly considers the performance of estimating the
unknown target position and controlling the mobile vehicle to a known position,
and then adopt a weighted sum method with normalization to solve it. The
controller is given based on the solution of the optimization problem in ties
with a least-square estimate of the target position. Moreover, the controller
does not need the vehicle's global position information. Finally, simulation
results are included to validate the effectiveness of the proposed controller.