[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Next Article in Journal
Time-Aware Service Ranking Prediction in the Internet of Things Environment
Next Article in Special Issue
A Double Rate Localization Algorithm with One Anchor for Multi-Hop Underwater Acoustic Networks
Previous Article in Journal
Horizontal Directional Drilling-Length Detection Technology While Drilling Based on Bi-Electro-Magnetic Sensing
Previous Article in Special Issue
Energy-Efficient Channel Coding Strategy for Underwater Acoustic Networks
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Artificial Measurements-Based Adaptive Filter for Energy-Efficient Target Tracking via Underwater Wireless Sensor Networks

1
State Key Laboratory of Industrial Control Technology, Hangzhou 310027, China
2
College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China
3
School of Marine Science and Technology, Northwestern Polytechnical University, Xi’an 710072, China
*
Author to whom correspondence should be addressed.
This article is an expanded version of an earlier conference paper: Zhang, S.; Chen, H.; Liu, M.; Zhang, Q. Artificial Measurements Based Filter for Energy-Efficient Target Tracking in Underwater Wireless Sensor Networks. In Proceedings of the 2016 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Hong Kong, China, 5–8 August 2016.
Sensors 2017, 17(5), 971; https://doi.org/10.3390/s17050971
Submission received: 21 March 2017 / Revised: 19 April 2017 / Accepted: 23 April 2017 / Published: 27 April 2017
(This article belongs to the Special Issue Advances and Challenges in Underwater Sensor Networks)
Figure 1
<p>Conventional distributed fusion architecture for target tracking.</p> ">
Figure 2
<p>Measurement residual indicator based distributed fusion architecture.</p> ">
Figure 3
<p>Artificial measurement-based distributed fusion architecture.</p> ">
Figure 4
<p>Flow chart of artificial measurements-based adaptive filter.</p> ">
Figure 5
<p>Target tracking performance: <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics> </math> versus <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics> </math>.</p> ">
Figure 6
<p>Performance comparison: <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics> </math> versus <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics> </math>. (<b>a</b>) Target tracking error: <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics> </math> versus <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics> </math>. (<b>b</b>) Energy consumptions: <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>0</mn> </mrow> </semantics> </math> versus <math display="inline"> <semantics> <mrow> <mi>δ</mi> <mo>=</mo> <mn>1</mn> </mrow> </semantics> </math>.</p> ">
Figure 7
<p>Impacts of normalized threshold <math display="inline"> <semantics> <mi>δ</mi> </semantics> </math>. (<b>a</b>) Target tracking error with different <math display="inline"> <semantics> <mi>δ</mi> </semantics> </math>. (<b>b</b>) Energy consumptions with different <math display="inline"> <semantics> <mi>δ</mi> </semantics> </math>.</p> ">
Figure 8
<p>Impacts of pre-given reference value <math display="inline"> <semantics> <msub> <mi mathvariant="normal">Θ</mi> <mi>r</mi> </msub> </semantics> </math>. (<b>a</b>) Target tracking error with different <math display="inline"> <semantics> <msub> <mi mathvariant="normal">Θ</mi> <mi>r</mi> </msub> </semantics> </math>. (<b>b</b>) Energy consumptions with different <math display="inline"> <semantics> <msub> <mi mathvariant="normal">Θ</mi> <mi>r</mi> </msub> </semantics> </math>.</p> ">
Figure 9
<p>Performances of different sensor groups. (<b>a</b>) Target tracking error with different sensor groups. (<b>b</b>) Energy consumptions with different sensor groups.</p> ">
Figure 10
<p>Performances of different number of selected sensors. (<b>a</b>) Target tracking error with different number of selected sensors. (<b>b</b>) Energy consumptions with different number of selected sensors.</p> ">
Figure 11
<p>Number of cases needed to try of different search algorithms.</p> ">
Versions Notes

Abstract

:
We study the problem of energy-efficient target tracking in underwater wireless sensor networks (UWSNs). Since sensors of UWSNs are battery-powered, it is impracticable to replace the batteries when exhausted. This means that the battery life affects the lifetime of the whole network. In order to extend the network lifetime, it is worth reducing the energy consumption on the premise of sufficient tracking accuracy. This paper proposes an energy-efficient filter that implements the tradeoff between communication cost and tracking accuracy. Under the distributed fusion framework, local sensors should not send their weak information to the fusion center if their measurement residuals are smaller than the pre-given threshold. In order to guarantee the target tracking accuracy, artificial measurements are generated to compensate for those unsent real measurements. Then, an adaptive scheme is derived to take full advantages of the artificial measurements-based filter in terms of energy-efficiency. Furthermore, a computationally efficient optimal sensor selection scheme is proposed to improve tracking accuracy on the premise of employing the same number of sensors. Simulation demonstrates that our scheme has superior advantages in the tradeoff between communication cost and tracking accuracy. It saves much energy while loosing little tracking accuracy or improves tracking performance with less additional energy cost.

1. Introduction

More than 70% of the earth’s surface is covered by seas and oceans. Seas and oceans are mysterious and charismatic to human beings because of the huge amount of unexploited resources. Underwater wireless sensor networks (UWSNs) technologies are developing gradually to enhance our abilities to discover resources in aquatic environments [1,2,3,4]. UWSNs are three-dimensional (3D) networks. The communication between underwater sensors relies on acoustic waves. UWSNs have a broad range of applications such as environmental monitoring, undersea exploration, disaster prevention, and distributed tactical surveillance, etc. We study the problem of accurately and energy-efficiently tracking a maneuvering target via UWSNs. UWSNs are the extending of wireless sensor networks (WSNs) which are applied to terrestrial environments [5,6,7]. One of the significant differences [8] between UWSNs and WSNs is the cost. Since underwater sensors need to work in the extreme underwater environment, they are much more expensive than terrestrial sensors. Underwater sensors use acoustic waves while terrestrial sensors use radio frequency waves. The energy consumption for communication between underwater sensors is higher than terrestrial sensors. Moreover, the sensors of UWSNs are battery-powered and it is impracticable to replace batteries when exhausted. This means that the battery life affects the lifetime of the whole network. Compared with the energy cost of sensing and processing, communication cost dominates the whole energy cost according to the energy model shown in [9]. Thus, in this paper, we improve the energy efficiency of target tracking by cutting down less helpful communications between local sensors and the fusion center.
This paper addresses the issue of implementing the tradeoff between the communication rate and target tracking accuracy. Local sensors need to figure out whether to send their information to the fusion center or not. This idea was inspired by some research about remote state estimation under communication constraints [10,11,12,13]. When a local sensor obtains a measurement about the target, it needs to figure out whether the measurement residual is large enough. A large measurement residual means the new measurement has enough value to be sent to the fusion center. If the measurement residual is larger than the threshold, the fusion center receives information from the local sensor and works as usual. If the measurement residual is smaller than the threshold, the fusion center receives nothing from the local sensor and generates an artificial measurement to approximate the unsent one, which makes full use of the information of the unsent measurement. Then, we derive the corresponding artificial measurements-based recursive form the filter. A preliminary version of the present paper appeared as a conference paper in [14]. The current version extends the conference version by providing an adaptive method for determining proper criteria which are used to tell local sensors whether their measurements have enough value to be sent to the fusion center. Moreover, a computationally efficient optimal sensor selection scheme is proposed to improve tracking accuracy on the premise of employing the same number of sensors.
The main contributions of this paper are threefold. First, we derive an artificial measurements-based filter, which has advantages in energy-efficiency. Second, in order to exploit advantages of our filter, we propose an adaptive method for determining proper criteria, which results in our artificial measurements based adaptive filter. Last, an optimal sensor selection scheme is proposed to further improve the energy-efficiency.
The rest of the paper is organized as follows. In Section 2, we discuss the related work in the area of target tracking in UWSNs. In Section 3, we formulate the problem and introduce some propaedeutic. In Section 4, we introduce our artificial measurements-based adaptive filter. In Section 5, we present our simulation results to verify our adaptive filter and discuss its characteristics. Finally, in Section 6, we provide the conclusions.

2. Related Work

Target tracking is a focused application for underwater defense systems. Intended targets to be tracked are unmanned underwater vehicles (UUVs) and submarines. As an emerging research interest, only a few works about target tracking in UWSNs can be found in the literature. In early work, a simple target tracking method utilizing only measurement information for 3D underwater is presented by Isbitiren et al. [15]. Based on the time of arrival of the echoes from the target after transmitting acoustic pulses from the sensors, the ranges of the nodes to the target are determined, and trilateration is used to obtain the location of the target. This method tracks the target only based on current measurements, which is adverse in terms of achieving high target tracking accuracy. In sparse networks, it results in tracking failure if not enough sensors are involved. In order to get better target tracking performance, Wang et al. [16] proposed an algorithm that combines the interacting multiple model (IMM) with the particle filter (PF) to cope with uncertainties in target maneuvers. To realize energy-effective target tracking, Yu et al. [17] provided an algorithm named wake-up/sleep (WuS) increasing energy efficiency of each sensor by using a distributed architecture. At each time step, WuS means waking up sensors that have an opportunity to detect the target and sending those that do not to sleep. However, it wastes energy by employing all candidate sensors without the survival of the fittest. Later, Zhang et al. [18] proposed an adaptive sensor scheduling scheme which saves energy by changing the sampling interval. The sampling interval is variable according to whether the tracking accuracy is satisfactory or not at each time step. The main distinction between this paper and the reference [18] is that they improve energy-efficiency from different dimensions. This paper focus on the tradeoff between communication rate and tracking accuracy at each time step, which is from a spatial dimension. Zhang et al. [18] focused on the tradeoff between sampling interval and tracking accuracy, which is from the temporal dimension. Recently, Zhang et al. [19] studied the effect of sensor topology on the target tracking in UWSNs with quantized measurements. They proposed a sensor selection method which selects the optimal topology by minimizing the posterior Cramer–Rao lower bound (PCRLB). This method improved the target tracking performance under the premise of employing the same number of sensors. However, the computation of PCRLB is complicated. In our work, we use the trace of predicted estimate covariance to select the optimal sensor group, which is more convenient than PCRLB.

3. Problem Formulation

This section formulates the problem of single target tracking via distributed UWSNs. The issues to be covered include system model, distributed fusion architectures and measurement residual-based sensor scheduling. For ease of reference, we list notations that will be used frequently in Table 1.

3.1. System Model

We consider the conventional target motion model, which is defined as
X k = F k X k 1 + w k 1 ,
where X k denotes the target state (positions and velocities) at time k, F k is the state transition matrix at time k and w k 1 is the process noise with zero-mean white Gaussian distribution N ( 0 , Q k 1 ) .
UWSNs consist of N wireless acoustic sensors floating at different seawater layers. The positions of sensors in Cartesian coordinates are denoted by s i = ( x i s , y i s , z i s ) ,   i = 1 , N . Sensors measure the distance to the target by transmitting acoustic pulses (ping) and calculating the time-of-arrival (ToA) of the pings and echoes.
The measurement model of the sensor s i at time k is given by
Z k i = h k i ( X k ) + v k i ,
where h k i ( X k ) is the measurement function, and v k i is the measurement noise with zero-mean white Gaussian distribution N ( 0 , R k i ) . The measurement function is given by
h k i ( X k ) = ( x k x i s ) 2 + ( y k y i s ) 2 + ( z k z i s ) 2 .
where ( x k , y k , z k ) is the location of the target at time k. The corresponding Jacobian matrix H k i , which is a useful approximation technique from the well-known extended Kalman filter (EKF) of the measurement function h k i ( · ) is given by
H k i ( X k ) = ( x k x i s ) / d , 0 , ( y k y i s ) / d , 0 , ( z k z i s ) / d , 0 .
where d = ( x k x i s ) 2 + ( y k y i s ) 2 + ( z k z i s ) 2 is the distance between the target and the sensor i.

3.2. Distributed Fusion Architectures

At the same time, different local sensors have different measurements about the same target. The information comes from different sensors must be fuse together to acquire more accurate estimates of target states. There are two types of fusion architectures, distributed fusion architectures and centralized fusion architectures. Distributed fusion architectures have advantages over centralized architectures in lower communications and processing costs. Therefore, distributed fusion architectures are preferential for application in resource-limited UWSNs. Figure 1 shows the normal structure of the distributed fusion system. Local sensors sample the measurements ( Z k 1 , Z k 2 , , Z k N ) from the target periodically. Then, based on new measurements and past information ( X ^ k 1 ), local sensors obtain local estimates ( X ^ k 1 , X ^ k 2 , , X ^ k N ) and transmit local estimates to the fusion center. Finally, the fusion center collects all local estimates and fuses them together to get the fusion estimate ( X ^ k ). The fusion estimate will be sent back to local sensors to predict future target states.

3.3. Measurement Residual-Based Sensor Scheduling

For the purpose of saving communication costs, we want local sensors to think carefully before sending their local estimates to the fusion center periodically as usual. If some local estimates have low values in updating the target state estimate, we should leave them at local sensors to reduce energy costs. The value of a local estimate can be measured by a measurement residual before we calculate the local estimate as
Z ˜ k i = Z k i Z ^ k | k 1 ,
Z ^ k | k 1 = h k i ( X ^ k | k 1 ) ,
where Z ˜ k i is the measurement residual and Z ^ k | k 1 is the predicted measurement. These make sense since the larger the measurement residual is, the larger the difference between the measurement updated estimate X ^ k and the predicted estimate X ^ k | k 1 = F k X ^ k 1 will be. A small measurement residual means it can change the predicted estimate X ^ k | k 1 only a little, so the fusion center can simply keep the predicted estimate X ^ k | k 1 or do some approximation.
The fusion center should formulate criteria to tell local sensors whether their estimates are needed or not. We adopt the standard proposed in [11] and make some changes in formulations for further convenience. We define an indicator function as
λ k i = 0 ,   | Z ˜ k i | δ E k i T 1 ,   otherwise
where δ is the normalized threshold and the weight E k i is determined by:
E k i T E k i = H k i P k | k 1 H k i T + R k i ,
where H k i is the Jacobian matrix of the measurement function with the predicted estimate X ^ k | k 1 and P k | k 1 is the error covariance of X ^ k | k 1 . Once a local sensor i obtains fresh measurements at time k, it should calculate the corresponding indicator value λ k i as Equation (7). If λ k i = 0 , sensor i needs to do nothing but keep silent to save energy. If λ k i = 1 , sensor i will calculate the local estimate and send it to the fusion center. Figure 2 illustrates how these indicator values work. It should be noticed that the measurement residual based fusion framework includes a feedback path from fusion center to local sensors. However, this feedback path adds to the negligible energy consumption of local sensors because the cost of receiving energy is much smaller than the cost of transmissitting energy according to [9].

4. Artificial Measurement Based Adaptive Filter

4.1. Artificial Measurement Model

Even if the fusion center did not receive local estimates from some local sensors, it obtained useful information that their measurement residuals are smaller than the threshold. For instance, if a fusion center did not get a packet from local sensor i at time k, then λ k i = 0 and
Z ^ k | k 1 δ E k i T Z k i Z ^ k | k 1 + δ E k i T .
Based on the well-known Bayesian formula, Z k i conditioned on λ k i = 0 and X ^ k | k 1 i has the distribution
f ( Z k i | λ k i = 0 , X ^ k | k 1 i ) f ( λ k i = 0 | Z k i , X ^ k | k 1 i ) f ( Z k i | X ^ k | k 1 i ) .
where f ( λ k i = 0 | Z k i , X ^ k | k 1 i ) has the uniform distribution
f ( λ k i = 0 | Z k i , X ^ k | k 1 i ) U ( [ Z ^ k | k 1 i δ E k i T , Z ^ k | k 1 i + δ E k i T ] ) ,
and before we get λ k i ,
E ( Z k i Z ^ k | k 1 i ) = 0 ,
Cov ( Z k i Z ^ k | k 1 i ) = H k i P k | k 1 H k i T + R k i .
f ( Z k i | X ^ k | k 1 i ) has the normal distribution as
f ( Z k i | X ^ k | k 1 i ) N ( Z ^ k | k 1 i , H k i P k | k 1 H k i T + R k i ) .
So f ( Z k i | λ k i = 0 ,   X ^ k | k 1 i ) is simply a truncated normal distribution as
f ( Z k i | λ k i = 0 , X ^ k | k 1 i ) N [ a , b ] ( Z ^ k | k 1 i , H k i P k | k 1 H k i T + R k i ) .
where a = Z ^ k | k 1 i δ E k i T and b = Z ^ k | k 1 i + δ E k i T . We do not want to drop above information about Z k i . Thus, we define an artificial measurement model as
Z ¯ k i = Z k i + u k i .
where Z ¯ k i is the artificial measurement and u k i is the zero-mean measurement noise. This model can be regarded as a measurement model of the real measurement and we want to use the artificial measurement to approximate the real measurement. The Equation (15) can be rewritten as
Z ^ k | k 1 i = Z k i + ξ k i .
and
ξ k i N [ a , b ] ( 0 , H k i P k | k 1 H k i T + R k i ) .
It is obvious that Equation (17) matches our artificial measurement model such as Equation (16) well. If we let Z ^ k | k 1 i stand as the realization of our artificial measurement, the measurement noise u k i has the same distribution as ξ k i . According to the characteristics of the truncated normal distribution, the variance of u k i can be calculated as
Var ( u k i ) = ( H k i P k | k 1 H k i T + R k i ) ( 1 Δ ( δ ) ) .
Δ ( δ ) = 2 / 2 δ e δ 2 / 2 [ 2 Φ ( δ ) 1 ] 1 .
where Φ ( · ) is the standard Φ -function defined by
Φ ( δ ) = δ 1 / 2 π e x 2 / 2 d x .
So our artificial measurement model Equation (16) has a unique realization Z ^ k | k 1 i and a truncated normal noise u k i .

4.2. Artificial Measurement Based Filter

Based on previous discussion about our artificial measurement model, we can derive our artificial measurement-based filter as follows
(1) Predict:
Assume that we already obtain the fusion estimate X ^ k 1 and the corresponding error covariance P k 1 at time k. The predicted state estimate and the corresponding error covariance at time k can be calculated as
X ^ k | k 1 = F k X ^ k 1 ,
P k | k 1 = F k P k 1 F k T + Q k .
(2) Update:
If λ k i = 1 , local sensor i can update the estimate and corresponding error covariance as EKF as follows.
Measurement residual:
Z ˜ k i = Z k i h k i ( X ^ k | k 1 ) ,
covariance of Measurement residual:
S k i = H k i P k | k 1 H k i T + R k i ,
where H k i = H k i ( X ^ k | k 1 i ) refer to Equation (4). Optimal Kalman gain:
K k i = P k | k 1 H k i T S k i 1 ,
Updated state estimate and updated estimate covariance:
X ^ k i = X ^ k | k 1 i + K k i Z ˜ k i ,
P k i = P k | k 1 i K k i H k i P k | k 1 i .
If λ k j = 0 , the fusion center will not receive the local estimate of j, it can update the artificial local estimate and corresponding covariance itself with the help of the artificial measurement model as follows. Measurement residual:
Z ¯ ˜ k j = h k j ( X ^ k | k 1 ) h k j ( X ^ k | k 1 ) = 0 ,
covariance of Measurement residual:
S ¯ k j = ( H k j P k | k 1 H k j T + R k j ) ( 2 Δ ( δ ) ) ,
The above equation can be obtained by
S ¯ k j = Cov ( Z ¯ ˜ k j ) = Cov ( Z ¯ k j Z ^ k | k 1 ) = Cov ( Z k j Z ^ k | k 1 + u k j ) = Cov ( Z k j Z ^ k | k 1 ) + Cov ( u k j ) = H k j P k | k 1 H k j T + R k j + ( H k j P k | k 1 H k j T + R k j ) ( 1 Δ ( δ ) ) = ( H k j P k | k 1 H k j T + R k j ) ( 2 Δ ( δ ) ) ,
Optimal Kalman gain:
K ¯ k j = P k | k 1 H k j T S ¯ k j 1 ,
Updated state estimate and updated estimate covariance:
X ^ k j = X ^ k | k 1 j + K ¯ k j Z ¯ ˜ k j = X ^ k | k 1 j ,
P k j = P k | k 1 j K ¯ k j H k j P k | k 1 j .
Finally, the fusion center will fuse these local estimates and artificial local estimates together to get the fusion estimate using the fusion algorithm proposed in [20].
P k = [ i = 1 N P k i 1 ( N 1 ) P k | k 1 1 ] 1 .
X ^ k = P k [ i = 1 N P k i 1 X ^ k i ( N 1 ) P k | k 1 1 X ^ k | k 1 ] .
Figure 3 shows the structure of our artificial measurements-based distributed fusion filter. The red part tells us that the fusion center makes full use of the information of λ k 1 = 0 and generates the artificial local estimate to compensate for the unsent real local estimate.

4.3. Adaptive δ Determination

The normalized threshold δ plays a key role in our artificial measurements based filter. Indicate function refereing to Equation (7), which is a function of δ , gives local sensors a criterion to decide whether to send local estimates to the fusion center or not. The larger the δ is, the lower probability local sensors have to send local estimates. That means δ affects the communication frequency of local sensors. In other words, δ can determine the energy consumption of local sensors. In addition, the estimate covariance in Equation (34) is a function of δ , which means δ affects not only the energy consumption but also the estimate accuracy of the target state. Both energy cost and target tracking accuracy are important in underwater target tracking issues. In this section, we proposes an adaptive δ determination method to make a better tradeoff between the communication cost and target tracking accuracy.
From Equations (7) and (8), we can obtain the probability distribution function of the indicator λ k i as
p ( λ k i = 0 ) = 2 Φ ( δ ) 1 .
p ( λ k i = 1 ) = 2 2 Φ ( δ ) .
Consequently, the expectation of communication rate of local sensor i at time k is
E ( λ k i ) = 2 2 Φ ( δ ) = Energy ( δ ) .
In this paper, we use the number of packets sent to the fusion center to measure energy cost at local sensors, which is defined as Energy( δ ).
From Equations (28) and (34), we can calculate the expectation of estimate covariance of local sensor i at time k as
E ( P k i ) = p ( λ k i = 0 ) ( P k | k 1 i K ¯ k i H k i P k | k 1 i ) + p ( λ k i = 1 ) ( P k | k 1 i K k i H k i P k | k 1 i ) .
we rewrite above equation as
E ( P k i ) = P k | k 1 i + Error ( δ ) K k i H k i P k | k 1 i ,
and
Error ( δ ) = [ p ( λ k i = 0 ) / ( 2 Δ ( δ ) ) + p ( λ k i = 1 ) ] .
It is clear that the estimate covariance increases with the increase of Error ( δ ). Equations (39) and (41) formulate how δ affects the energy costs of local sensors and target tracking accuracy. For the purpose of selecting proper δ during target tracking missions, we define an objective function as
O ( δ ) = Energy ( δ ) + α k Error ( δ ) ,
where α k is a coefficient to adjust the weight of tracking error. The optimal δ can be determined by
δ o p t = arg min O ( δ ) .
Since we should guarantee target tracking performance first, the coefficient α k should be large when tracking performance is bad. Conversely, α k should be small if tracking accuracy meets our demand. In this paper, we use the trace of estimate covariance to measure tracking performance and α k is determined as
α k = Θ k / Θ r ,
where Θ k is the trace of P k in Equation (35) and Θ r is a pre-given reference value. It is not easy to get an analytical solution of Equation (44). Therefore, we provide an efficient numerical way to get the proper δ . Since δ is equivalent to the standard deviation of standard normal distribution, we select δ from 0 to 3 according to the well-known Pauta criterion. Then we take n uniformly-spaced samples from [0, 3] and form a set as
δ { L 1 , , L n 1 , L n } .
Thus, Equation (44) can be rewritten as
δ o p t = arg min δ { L 1 , , L n 1 , L n } O ( δ ) .
It should be noticed that Energy ( L i ) and Error ( L i ) can be calculated off-line to improve the online computational efficiency.

4.4. Optimal Sensor Group Selection

Assume that filtering results X ^ k and P k are given at time k. Then, target position at time k + 1 can be predicted referring to Equation (1) and the distance d k + 1 i from sensor i to the predicted target position can be calculated. Sensor i will have a chance to track the target and be a candidate sensor at time k + 1 if d k + 1 i is smaller than its sensing range. However, in 3D networks, we need four sensors to locate the position of a target [15], which means that adopting more than four sensors is not worthwhile if we consider their energy consumption. Thus, we select the best four sensors if there are more than four candidate sensors. Certainly, we employ all the candidate sensors when the number is less than or equal to four. We had proposed a posterior Cramer-Rao lower bound (PCRLB) based sensor selection scheme for particle filters in [21]. It calculated the PCRLBs of different sensor groups to evaluate how they contribute to tracking performance. However, in this paper, the fusion estimate covariance is given in Equation (35). That means we can evaluate how sensor groups contribute to tracking performance by P k + 1 . This is more convenient than PCRLBs.
Given that we have a set of N c candidate sensors as C = { c 1 , c 2 , , c N c } at time k + 1 , the sensor selection problem can be formulated by
G o p t = arg min G C Ψ k + 1 ,
Ψ k + 1 = Trace ( P k + 1 ) = Trace { [ i G P k + 1 i 1 ( N 1 ) P k + 1 | k 1 ] 1 } ,
where G = { c g 1 , c g 2 , c g 3 , c g 4 } stands for sensor groups selected from set { c 1 , c 2 , , c N c } . Equation (48) tells us that the sensor group that minimizes the trace of the fusion estimate covariance will be the best one. Since the rank of the traces of sensor groups in Equation (48) will not be changed with δ , P k + 1 i is given by Equation (28) for simplified calculation. The exhaustive search is the most direct way to find the optimal sensor group and there will be N c ! 4 ! ( N c 4 ) ! groups. However, when N c is large, the number of the groups increase rapidly and the exhaustive search has a heavy computational burden. Therefore, we use the generalized Breiman, Friedman, Olshen and Stone (GBFOS) proposed in [22] to find the optimal sensor group. Initially, N = N c and the GBFOS algorithm keeps finding the optimal N 1 elements subset from the last optimal N elements set until N = 4 . Thus, GBFOS needs to try ( N c 4 ) ( N c + 5 ) 2 sensor groups to find the optimal one. Table 2 lists some numerical examples to compare exhaustive search with GBFOS. It is obvious that GBFOS can reduce much more computational burden when N c becomes larger and larger.
We should mention the other search algorithm, called the greedy search, which is the reverse of GBFOS. The greedy algorithm keeps taking one optimal sensor out from candidate sensors until four sensors are taken out. Thus, the greedy search needs to try ( 4 N c 6 ) sensor groups to find the optimal one. It seems better than GBFOS when N c is larger than 8. However, it is infeasible for our sensor selection scheme because there is no optimal sensor if we only choose one. That is according to
Trace ( P k + 1 i ) = Trace ( P k + 1 j ) , i , j C .
The flow chart of our artificial measurements-based adaptive filter is shown in Figure 4. It shows how information flows between the fusion center and local sensors. The dash line means the information flow from local sensors to the fusion center is non-existant when λ = 0 and artificial measurements are introduced to compensate this missing information flow.

5. Simulation and Results

5.1. Simulation Scenario

We employ our artificial measurement-based filter to a target tracking mission for verification. In order to get more realistic performance measures, the target is assumed to move in a 3D underwater environment. The monitored field is 1000 m×1000 m×1000 m and sensors are deployed as a 5 × 5 × 5 uniform grid. All local sensors are identical. Their detection radius and measure covariance are 300 m and 10 m 2 , respectively. The initial state of the target is assumed to be [ 300 , 10 , 300 , 2 , 10 , 2 ] T . From 1 s to 40 s, it moves at constant velocity (CV). From 41 s to 80 s, it makes a coordinate turn (CT) with turn rate 0.052 rad/s. From 81 s to 100 s, it moves at CV. CV and CT can be formulated as
X k = F CV ( X k 1 ) + w k ,
X k = F CT ( X k 1 ) + w k ,
where F CV and F CT are state transition matrixes. w k is the process noise with zero-mean white Gaussian distributions N ( 0 , Q k ) . F CV , F CT and Q k are given by
F CV = 1 T 0 0 0 0 0 1 0 0 0 0 0 0 1 T 0 0 0 0 0 1 0 0 0 0 0 0 1 T 0 0 0 0 0 1 ,
F CT = 1 sin ( w T ) w 0 cos ( w T ) 1 w 0 0 0 cos ( w T ) 0 sin ( w T ) 0 0 0 1 cos ( w T ) w 1 sin ( w T ) w 0 0 0 sin ( w T ) 0 cos ( w T ) 0 0 0 0 0 0 1 T 0 0 0 0 0 1 ,
Q k = q 2 T 3 3 T 2 2 0 0 0 0 T 2 2 T 0 0 0 0 0 0 T 3 3 T 2 2 0 0 0 0 T 2 2 T 0 0 0 0 0 0 T 3 3 T 2 2 0 0 0 0 T 2 2 T .
where q is the intensity of the process noise. For an underwater target, we consider that only on the xoy plane does it move as a CT model and it moves as a CV model in Z-axis direction.

5.2. Performance Verification

Simulation results are averaged over 100 Monte-Carlo runs. We adopt root mean square error (RMSE) to assess the accuracy of target tracking and the number of packets sent from local sensors to indicate the energy consumption.

5.2.1. Performance Comparison

In our simulation, we compare performances between the conventional target tracking scheme referring to δ = 0 and our artificial measurements-based energy-efficient target tracking scheme with δ = 1 to see how our algorithm achieves the goal of energy-efficiency. Figure 5 shows the real trajectory of the target and performances of different tracking schemes. All real measurements mean that the normalized threshold is equal to 0 and the fusion center has all local estimates. Containing artificial measurements means that the normalized threshold is not equal to 0 (in this case δ = 1 ) and the fusion center updates the estimate with the help of artificial measurements. Both target tracking schemes can successfully track the target with high accuracy. The detailed tracking error and communication costs of the two schemes are displayed in Figure 6a,b respectively. They tell us that our target tracking scheme can save much energy (about 80%) while only loosing a small amount of tracking accuracy (about 40%). This is a worthwhile trade in the pursuit of saving energy.

5.2.2. Impacts of δ

We also want to know how the normalized threshold δ affects the target tracking performance and communication costs of our algorithm. So we change the δ from 0 to 2 with the increment a of 0.1 at each step. Figure 7a illustrates the target tracking performance with different δ and we mark some points for further discussion. It is clear that target tracking error increases very slowly with a small δ . However, it increases faster and faster with the increase of δ . Compared with δ = 0 , it increases only 5.6% when δ = 0 . 5 and 41.2% when δ = 1 . The error turns out to be extremely high (more than 300%) when δ = 2 . Correspondingly, the impact of different δ on communication costs is displayed in Figure 7b. Inversely, the communication rate falls rapidly when δ is small and it changes slowly when δ is lager than 1.5. The communication cost decreases by 46.5% when δ = 0 . 5 and 79% when δ = 1 . Through Figure 7, we can find that our algorithm has a nice property with δ , which shows the potential that we can save much energy while only loosing a little tracking accuracy. Take δ = 0 . 5 as a powerful example. We save 46.5% energy while loosing only 5.6% tracking accuracy. Therefore, the artificial measurements-based adaptive filter is proposed to find the proper δ .

5.2.3. Performance of Adaptive Filter

Since the normalized threshold δ has opposite effects on target tracking accuracy and energy consumptions, it is important to find the proper δ to exploit the advantages of our artificial measurements-based filter. The core idea is to reduce energy consumption as much as possible under the premise of sufficient target tracking performance. In this paper, tracking performance is represented by α k in Equation (45). A large α k means tracking error is too large and we should put more effort into decreasing tracking error. In contrast, a small α k means tracking performance meets our demand and we should pay attention to reducing energy cost. Θ r in Equation (45) is a pre-given reference value to represent our demand for target tracking performance. In order to show the superiority of our artificial measurements-based adaptive filter, we change Θ r from 10 to 100 with the increment of 10 at each step and set n in Equation (46) to 30. Figure 8a displays the target tracking performance with different Θ r . It is clear that target tracking error increases slowly with Θ r , which means our adaptive filter can guarantee target tracking performance even with a large Θ r . in contrast, the target tracking error becomes extremely high with a large δ in Figure 7a. Correspondingly, the impact of different Θ r on communication costs is shown in Figure 8b. The communication cost decreases observably as Θ r increases even with a small Θ r , which means our adaptive filter can effectively reduce energy consumption. Marked points in Figure 8 illustrate that our adaptive filter can exploit its advantages in energy saving and guarantee good target tracking performance. For example, compared with δ = 0 , we save 29.1% energy while lost only 2.11% tracking accuracy when Θ r = 10 and we save 74.35% energy while loosing 32.9% tracking accuracy when Θ r = 100 . By setting different Θ r , the artificial measurements-based adaptive filter can achieve varying degrees of pursuits of energy saving.

5.2.4. Performance of Sensor Group Selection

Sensor group selection is a feasible method to improve the energy-efficiency of target tracking in UWSNs. We compare performances of three sensor group selection schemes to support this opinion. Here, we set Θ r = 40 . Target tracking error and energy consumptions of three schemes are displayed in Figure 9a,b, respectively. Selecting the best four sensors means that the sensor group is optimized using Equation (48). In contrast, selecting the worst four sensors means the sensor group is generated by
G w o r s t = arg max G C Ψ k + 1 .
Selecting a random four sensors means the sensor group is generated at random. Obviously, the performance of selecting the worst four sensors is much worse than other schemes. This must be avoided in terms of energy-efficiency and target tracking accuracy. The sensor group optimized by our sensor selection scheme has the best performance in both target tracking and energy saving. Average tracking accuracy and the total number of packets are listed in Table 3. Compared with the random sensor group, selecting the best four sensors can improve tracking accuracy by 19.61% and save 8.65% of energy.
Moreover, in some cases, we want to improve tracking accuracy with less additional energy consumption. This goal can be realized by selecting more sensors and using our artificial measurements-based adaptive filter. In order to have more candidate sensors, the sensors are deployed as a 7 × 7 × 7 uniform grid. The number of selected sensors N s at each step is changed from 4 to 11. We compare performances between artificial measurements-based adaptive filter with Θ r = 40 and conventional all real measurements-based filter. The averaged target tracking error and communication costs of the two schemes are shown in Figure 10a,b respectively. It is clear that selecting more sensors can improve target tracking accuracy effectively with regard to both schemes. However, the energy consumption of an artificial measurements-based adaptive filter increases much more slowly than a conventional all real measurements-based filter. That means our target tracking scheme can improve target tracking performance with less additional energy consumption. Through marked points in Figure 10, the energy cost of our scheme with N s = 10 is 4.25% less than conventional scheme with N s = 4 and the target tracking error of our scheme with N s = 10 is 28.1% lower than conventional scheme with N s = 4 . That means our scheme can have similar energy to a conventional scheme but get much better tracking performance than a conventional scheme.
In a dense sensor network, the exhaustive search needs to try too many cases to find the best sensor group, which deteriorates the real-time performance of our optimal sensor group selection scheme. Hence, we use the GBFOS algorithm to reduce the number of cases and improve the computation efficiency of our optimal sensor group selection scheme. The number of cases needed to try to find the best sensor group of exhaustive search and GBFOS is plotted in Figure 11. Here we set N s = 4 . Compared with the exhaustive search, the GBFOS algorithm reduces the number of cases by about two orders of magnitude. Hence, the GBFOS algorithm can remarkably improve the real-time performance of our optimal sensor group selection scheme.
Overall, this work focus on providing an energy-efficient target tracking algorithm for resource limited UWSNS. Our artificial measurements-based adaptive filter is easy to implement because it has widely applied the Kalman filter structure and low online computation demand.

6. Conclusions

This paper proposes an artificial measurements-based energy-efficient target tracking scheme in UWSNs. The basic idea of our approach is that, under the distributed fusion framework, we abandon low value local measurements to decrease the communication rate from local sensors to the fusion center to save energy. We guarantee the tracking accuracy by generating corresponding artificial measurements in the fusion center to compensate for unsent measurements. Then, we derive an adaptive filter based on these artificial measurements. In addition, we propose an optimal sensor selection scheme to further improve the energy-efficiency. Through simulation results, we can draw the following conclusions. Firstly, the artificial measurements based adaptive filter can save much energy while loosing less tracking accuracy. Secondly, by setting a different pre-given reference value Θ r , this adaptive filter can achieve varying degrees of pursuit of energy saving. Thirdly, our computationally efficient optimal sensor selection algorithm can efficiently improve target tracking performance under the premise of employing the same number of sensors. Finally, with the increase of the number of selected sensors, our artificial measurements-based adaptive filter better utilizes its advantages in energy-efficiency.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grants 61374021, 61531015, U1609204, and 61673345.

Author Contributions

Huayan Chen and Senlin Zhang conceived the main idea; Senlin Zhang and Qunfei Zhang designed the simulation; Huayan Chen and Meiqin Liu performed the simulation; Huayan Chen wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaime, L. Underwater sensor nodes and networks. Sensors 2013, 13, 11782–11796. [Google Scholar]
  2. Climent, S.; Capella, J.V.; Meratnia, N.; Serrano, J.J. Underwater sensor networks: A new energy efficient and robust architecture. Sensors 2012, 12, 704–731. [Google Scholar] [CrossRef] [PubMed]
  3. Cui, J.H.; Kong, J.J.; Zhou, S.L. The challenges of building scalable mobile underwater wireless sensor networks for aquatic applications. IEEE Netw. 2006, 20, 12–18. [Google Scholar]
  4. Climent, S.; Sanchez, A.; Capella, J.V.; Meratnia, N.; Serrano, J.J. Underwater acoustic wireless sensor networks: advances and future trends in physical, MAC and routing layers. Sensors 2014, 14, 795–833. [Google Scholar] [CrossRef] [PubMed]
  5. Souza, E.L.; Nakamura, E.F.; Pazzi, R.W. Target Tracking for Sensor Networks: A Survey. ACM Comput. Surv. (CSUR) 2016, 49, 30. [Google Scholar] [CrossRef]
  6. Han, G.; Jiang, J.; Zhang, C.; Duong, T.Q.; Guizani, M.; Karagiannidis, G.K. A survey on mobile anchor node assisted localization in wireless sensor networks. IEEE Commun. Surv. Tutor. 2016, 18, 2220–2243. [Google Scholar] [CrossRef]
  7. Zhou, W.; Shi, W.R.; Zhang, H.D.; Wang, X.G.; Kang, D.J. Study on distributed target tracking in wireless sensor networks. Chin. J. Sci. Instrum. 2013, 34, 1485–1491. [Google Scholar]
  8. Catipovic, J. Performance limitations in underwater acoustic telemetry. IEEE J. Ocean. Eng. 1990, 15, 205–216. [Google Scholar] [CrossRef]
  9. Bhardwaj, M.; Chandrakasan, A.P. Bounding the lifetime of sensor networks via optimal role assignments. In Proceedings of the IEEE INFOCOM, New York, NY, USA, 23–27 June 2002; pp. 1587–1596. [Google Scholar]
  10. Li, L.; Lemmon, M.; Wang, X. Event-triggered state estimation in vector linear processes. In Proceedings of the American Control Conference, Baltimore, MD, USA, 30 June–2 July 2010; pp. 2138–2143. [Google Scholar]
  11. Wu, J.; Jia, Q.S.; Johansson, K.H.; Shi, L. Event-based sensor data scheduling: Trade-off between communication rate and estimation quality. IEEE Trans. Autom. Control 2013, 58, 1041–1046. [Google Scholar] [CrossRef]
  12. Wu, J.; Li, Y.; Quevedo, D.E.; Lau, V.; Shi, L. Data-driven power control for state estimation: A bayesian inference approach. Automatica 2015, 54, 332–339. [Google Scholar] [CrossRef]
  13. Shi, D.; Chen, T.; Shi, L. On set-valued Kalman filtering and its application to event-based state estimation. IEEE Trans. Autom. Control 2015, 60, 1275–1290. [Google Scholar] [CrossRef]
  14. Zhang, S.; Chen, H.; Liu, M.; Zhang, Q. Artificial measurements based filter for energy-efficient target tracking in underwater wireless sensor networks. In Proceedings of the 2016 IEEE International Conference on Signal Processing, Communications and Computing (ICSPCC), Hong Kong, China, 5–8 August 2016; pp. 1–6. [Google Scholar]
  15. Isbitiiren, G.; Akan, O.B. Three-dimensional underwater target tracking with acoustic sensor networks. IEEE Trans. Veh. Technol. 2011, 60, 3897–3906. [Google Scholar] [CrossRef]
  16. Wang, X.; Xu, M.; Wang, H. Combination of interacting multiple models with the particle filter for three-dimensional target tracking in underwater wireless sensor networks. Math. Probl. Eng. 2012, 2012, 829451. [Google Scholar] [CrossRef]
  17. Yu, C.H.; Lee, J.C.; Choi, J.W.; Park, M.K.; Kang, D.J. Energy efficient distributed interacting multiple model filter in UWSNs. In Proceedings of the International Conference on Control, Automation and Systems, JeJu Island, Korea, 17–21 October 2012; pp. 1093–1098. [Google Scholar]
  18. Zhang, S.; Chen, H.; Liu, M. Adaptive sensor scheduling for target tracking in underwater wireless sensor networks. In Proceedings of the 2014 IEEE International Conference on Mechatronics and Control (ICMC), Jinzhou, China, 3–5 July 2014; pp. 55–60. [Google Scholar]
  19. Zhang, Q.; Liu, M.; Zhang, S. Node Topology Effect on Target Tracking Based on UWSNs Using Quantized Measurements. IEEE Trans. Cybern. 2015, 45, 2323–2335. [Google Scholar] [CrossRef] [PubMed]
  20. Zhu, Y.M.; You, Z.S.; Zhao, J.; Zhang, K.S.; Li, X.R. The optimality for the distributed Kalman filtering fusion with feedback. Automatica 2001, 37, 1489–1493. [Google Scholar] [CrossRef]
  21. Zhang, Q.; Liu, M.; Zhang, S.; Chen, H. Node topology effect on target tracking based on underwater wireless sensor networks. In Proceedings of the 2014 IEEE 17th International Conference on Information Fusion (FUSION), Salamanca, Spain, 7–10 July 2014; pp. 1–8. [Google Scholar]
  22. Riskin, E. Optimal bit allocation via the generalized bfos algorithm. IEEE Trans. Inf. Theory 1991, 37, 400–402. [Google Scholar] [CrossRef]
Figure 1. Conventional distributed fusion architecture for target tracking.
Figure 1. Conventional distributed fusion architecture for target tracking.
Sensors 17 00971 g001
Figure 2. Measurement residual indicator based distributed fusion architecture.
Figure 2. Measurement residual indicator based distributed fusion architecture.
Sensors 17 00971 g002
Figure 3. Artificial measurement-based distributed fusion architecture.
Figure 3. Artificial measurement-based distributed fusion architecture.
Sensors 17 00971 g003
Figure 4. Flow chart of artificial measurements-based adaptive filter.
Figure 4. Flow chart of artificial measurements-based adaptive filter.
Sensors 17 00971 g004
Figure 5. Target tracking performance: δ = 0 versus δ = 1 .
Figure 5. Target tracking performance: δ = 0 versus δ = 1 .
Sensors 17 00971 g005
Figure 6. Performance comparison: δ = 0 versus δ = 1 . (a) Target tracking error: δ = 0 versus δ = 1 . (b) Energy consumptions: δ = 0 versus δ = 1 .
Figure 6. Performance comparison: δ = 0 versus δ = 1 . (a) Target tracking error: δ = 0 versus δ = 1 . (b) Energy consumptions: δ = 0 versus δ = 1 .
Sensors 17 00971 g006
Figure 7. Impacts of normalized threshold δ . (a) Target tracking error with different δ . (b) Energy consumptions with different δ .
Figure 7. Impacts of normalized threshold δ . (a) Target tracking error with different δ . (b) Energy consumptions with different δ .
Sensors 17 00971 g007
Figure 8. Impacts of pre-given reference value Θ r . (a) Target tracking error with different Θ r . (b) Energy consumptions with different Θ r .
Figure 8. Impacts of pre-given reference value Θ r . (a) Target tracking error with different Θ r . (b) Energy consumptions with different Θ r .
Sensors 17 00971 g008
Figure 9. Performances of different sensor groups. (a) Target tracking error with different sensor groups. (b) Energy consumptions with different sensor groups.
Figure 9. Performances of different sensor groups. (a) Target tracking error with different sensor groups. (b) Energy consumptions with different sensor groups.
Sensors 17 00971 g009
Figure 10. Performances of different number of selected sensors. (a) Target tracking error with different number of selected sensors. (b) Energy consumptions with different number of selected sensors.
Figure 10. Performances of different number of selected sensors. (a) Target tracking error with different number of selected sensors. (b) Energy consumptions with different number of selected sensors.
Sensors 17 00971 g010
Figure 11. Number of cases needed to try of different search algorithms.
Figure 11. Number of cases needed to try of different search algorithms.
Sensors 17 00971 g011
Table 1. List of notations.
Table 1. List of notations.
NotationsExplanations
X k Target state at time k
X ^ k Estimate of target state at time k
X ^ k | k 1 Predicted estimate of target state at time k
F k State transition matrix at time k
ω k 1 Process noise at time k 1
Q k Covariance of process noise at time k
Z k i Measurement of sensor i at time k
υ k i , u k i , ξ k i Measurement noise of sensor i at time k
R k i Covariance of measurement noise of sensor i at time k
Z ^ k | k 1 i Predicted measurement at time k
Z ˜ k i Measurement residual of sensor i at time k
Z ¯ k i Artificial measurement of sensor i at time k
h k i ( · ) Measurement function of sensor i at time k
H k i ( · ) Jacobian matrix of sensor i at time k
( x k , y k , z k ) Target location at time k
( x i s , y i s , z i s ) Location of Sensor i
δ Normalized threshold
λ k i Indicator value of sensor i at time k
P ^ k | k Estimate error covariance at time k
P ^ k | k 1 Predicted estimate error covariance at time k
f ( · ) Distribution of random variable
E ( · ) Expectation of random variable
Cov ( · ) Covariance of random variable
p ( · ) Probability of random variable
S k i Covariance of measurement residual of sensor i at time k
S ¯ k i Covariance of measurement residual of sensor i at time k with artificial measurement
K k i Kalman gain of sensor i at time k
K ¯ k i Kalman gain of sensor i at time k with artificial measurement
Θ k Trace of P k
Θ r Pre-given reference value
Table 2. The number of cases needed to try to find the best one.
Table 2. The number of cases needed to try to find the best one.
N c = 6 N c = 8 N c = 10 N c = 15 N c = 20
Exhaustive Search157021013654845
GBFOS112645110200
Table 3. Performances of different sensor groups.
Table 3. Performances of different sensor groups.
Worst Sensor GroupRandom Sensor GroupBest Sensor Group
Target tracking error10.63085.39764.3389
Number of packets292.65210.34192.16

Share and Cite

MDPI and ACS Style

Chen, H.; Zhang, S.; Liu, M.; Zhang, Q. An Artificial Measurements-Based Adaptive Filter for Energy-Efficient Target Tracking via Underwater Wireless Sensor Networks . Sensors 2017, 17, 971. https://doi.org/10.3390/s17050971

AMA Style

Chen H, Zhang S, Liu M, Zhang Q. An Artificial Measurements-Based Adaptive Filter for Energy-Efficient Target Tracking via Underwater Wireless Sensor Networks . Sensors. 2017; 17(5):971. https://doi.org/10.3390/s17050971

Chicago/Turabian Style

Chen, Huayan, Senlin Zhang, Meiqin Liu, and Qunfei Zhang. 2017. "An Artificial Measurements-Based Adaptive Filter for Energy-Efficient Target Tracking via Underwater Wireless Sensor Networks " Sensors 17, no. 5: 971. https://doi.org/10.3390/s17050971

APA Style

Chen, H., Zhang, S., Liu, M., & Zhang, Q. (2017). An Artificial Measurements-Based Adaptive Filter for Energy-Efficient Target Tracking via Underwater Wireless Sensor Networks . Sensors, 17(5), 971. https://doi.org/10.3390/s17050971

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop