International Journal of Distributed Sensor Networks

Volume 2012 (2012), Article ID 438090, 15 pages

doi:10.1155/2012/438090

Research Article

Cooperative Group Localization Based on Factor Graph for Next-Generation Networks

Xuefei Zhang, Qimei Cui, Yulong Shi, and Xiaofeng Tao

Key Laboratory of Universal Wireless Communications, Ministry of Education, Beijing University of Posts and Telecommunications, Beijing 100876, China

Received 3 May 2012; Revised 11 August 2012; Accepted 25 August 2012

Academic Editor: Yonghui Li

Copyright © 2012 Xuefei Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In ill-conditioned communication environment, multiple target localization is of important practical significance. The cooperative group localization (CGL) model was firstly put forward, which has verified the effectiveness of localization performance gain and simultaneous multiple target localization in ill conditions. However, there exist two inherent difficulties: the strict demand for CGL topology and the high complexity. By the rational use of information to relax restrictions on topology and by dividing the complex problem into some simple local ones, the factor graph (FG) together with the sum-product algorithm is a perfect candidate for the problems above. In order to solve the two problems, we propose the weighted FG-based CGL (WFG-CGL) algorithm which incorporates the optimal weights based on the information reliability. In order to further reduce the complexity, we propose the low-complexity FG-based CGL (LCFG-CGL) algorithm. The Cramer-Rao lower bound (CRLB) of the localization error in CGL is first derived. Theoretical analysis and numerical results indicate that the proposed algorithms not only perform better in relaxing CGL topology requirement, but also enjoy high localization accuracy under low complexity in comparison with the existing CGL algorithm.

1. Introduction

Wireless localization technologies, which are designated to estimate the position of a mobile terminal (MT), have drawn a lot of attention over the past few decades. Federal Communication Commission (FCC) has mandated the cellular network operators to estimate the position of emergency caller with the error of less than 100 meters in the enhanced 911 (E-911) [1]. In recent years, there are increasing demands for commercial applications to adopt the localization information within their system design, such as the navigation system, the health care system, the wireless sensor network (WSN) [2–4], and the intelligent transportation system (ITS) [5, 6]. Moreover, more and more people rely on the instant localization, which may make the communication services more user-friendly. With the emergent interests in the localization-based services (LBSs), the wireless localization with enhanced accuracy becomes necessary for the applications under different circumstances.

The research of the localization techniques is an important part of wireless localization fields [7]. A variety of wireless localization techniques have been investigated, such as basic localization algorithms and their corresponding Cramer-Rao lower bound (CRLB) analysis [8–10]. CRLB provides a useful means for the analysis of the limits of localization accuracy. CRLB results for localization, based on different types of internode measurements, can be found in [2, 3, 11–15]: the localization error is shown to depend on several elements, in particular on the measurement reliability (which may be impaired by multipath, non-line of sight, synchronization errors, etc.) and the topology of the network.

Most of current academic achievements in localization field focus on the enhancement localization techniques in the well condition. The well condition means that mobile terminal (MT) to be localized has the sufficient signal resources. For example, if one MT can be accurately located by the measurement parameters of time of arrival (TOA), time difference of arrival (TDOA), angle of arrival (AOA), or received signal strength (RSS), the MT connects well with at least three base stations (BSs) and line-of-sight (LOS) propagation paths existing between BSs and MT [16]. Nevertheless, the scenarios with sufficient signal resources do not always happen in real circumstances. Sometimes, the parameter measurements may be corrupted by non-line-of-sight (NLOS) error, or MT can not simultaneously connect with three BSs [9, 16]. These ill conditions with insufficient signal resources can make the accurate localization unrealistic.

The cooperative group localization (CGL), which can achieve the multiple target localization simultaneously, is proposed as an approach to solving the above ill-conditioned localization problem in [17, 18]. The existing research indicates that the CGL scheme not only has higher localization accuracy over the traditional methods in ill conditions, but also can effectively enhance the robustness of localization. However, there also exists two inherent difficulties: (1) the strict demand for CGL topology and (2) the high complexity. According to the rigid graph theory [19], the CGL topology must be global rigidity so as to satisfy the unique solvability. Furthermore, CGL has to face the problem of high complexity, especially with the increasing number of cooperative MTs.

In [20], factor graph (FG) theory is proposed, which can simplify a complex problem into multiple simple sub problems. The optimal or the near-optimal solution can be obtained by solving the sub problems, respectively. According to its feasibility and simplicity, it has been introduced into cellular or sensor network based on measurement parameters TOA, TDOA, AOA, and RSS [21–25] and reaches the good tradeoff between the localization accuracy and complexity. In [26], Wymeersch et al. analyzes the feasibility of FG theory applied into the cooperative localization.

In this paper, we combine FG theory with CGL. The FG-CGL graph can adapt its constitutions and architecture based on the different measurement information, considered as self-organizing. First, a weighted FG-based CGL algorithm (WFG-CGL) is proposed to effectively solve the above problems of strict demand for CGL topology and high complexity in the ill condition. It enhances the accuracy by adopting the different information reliability of BSs and MTs into its formulation. Then, a low-complexity FG-based CGL (LCFG-CGL) algorithm is proposed to further diminish the complexity. Moreover, CRLB analysis is extended to CGL scenarios. Numerical simulation results verify that compared to the existing CGL method, the proposed algorithms can relax CGL topology requirement and provide the high localization accuracy under low complexity. And the important factors for the FG-based CGL algorithms are analyzed based on the simulation results.

The paper is organized as follows. In Section 2, CGL model is introduced. In Section 3, the principles of factor graph and sum-product algorithm are presented. Section 4 proposes WFG-CGL algorithm and LCFG-CGL algorithm. Section 5 completes CRLB analysis of CGL. Section 6 evaluates the performance and the complexity of the proposed algorithms by numerical simulation. Finally, the conclusion is included in Section 7.

2. Mathematical Model of CGL

For the description convenience, take the scenario shown in Figure 1 for example, which presents the principle of CGL in the following sections. In CGL, the conception of the terminal group (TG) is firstly put forward. TG refers to a group of MTs supporting peer-to-peer communication mode in the next-generation communication networks. To represent the most general situation, all MTs are supposed to be the localization unknown nodes, and the relative distances between some certain MT pairs are known. For the next-generation communication networks, the advantages of the cooperation between MTs can be applied widely. As shown in Figure 1, each MT in TG connects with less than three BSs, which can be included in the ill conditions mentioned above.

Figure 1: CGL topology diagram.

The CGL graph can be divided into three sub-graphs: the first one is BS subgraph , where represents the BSs, and represents the existing edges between any two vertices of ; the second one is TG subgraph , where includes all the MTs in TG, and represents the distances determined by direct communication among the MTs; the third is cooperative sub-graph , where includes the BSs or the cooperative MTs in TG, and represents the edges existing between and . The operator of the sets denotes the number of the elements in a set.

Given the CGL graph, the aim of CGL is to simultaneously solve the positions of all the MTs in TG. However, not all topologies of CGL have unique solvability. The unique solvability means that the positions of MTs in a TG can be located without ambiguity. It can be investigated by the rigid graph theory [19].

3. Factor Graph and Sum-Product Algorithm

3.1. Factor Graph

A factor graph [20] is a bipartite graph that illustrates how a complicated global function with many variables is simplified into the product of several simple local functions. Each local function is a function with few variables. A factor graph has the variable node for each variable, the function node for each local function, and the edge connecting variable node to factor node if and only if the variable is an argument of the function. For example, is a real-valued function with five variables, which is divided into the product of five local functions and , as is shown in Figure 2, as follows:

Figure 2: An example of factor graph.

3.2. Sum-Product Algorithm

Sum-product algorithm [20] is one of the important algorithms in graphical model, which provides the convenient solution to marginal probability. The algorithm works by the passing messages along the edges between the nodes. There are two types of messages.

With two rules above, information is passed among neighboring nodes. Each variable node modifies its own value according to the received message. After a few times, the information converges. The variable can be estimated according to the product of overall message from its neighboring nodes as follows:

4. FG-Based CGL Algorithms

In the following analysis, the positions of several MTs are estimated simultaneously based on TOA measurements. In order to reduce the complexity of the two-dimensional (2D) problem, the problem is divided into two one-dimensional (1D) problems [21], -coordinate group and -coordinate group, as illustrated in Figure 3. In Figure 3, and represent the coordinates for the th MT; and represent the relative distances between the th MT and the th BS; and represent the relative distances between the th MT and the th BS after the weight adaptation; and represent relative distances between the th MT and the th cooperative MT; and represent relative distances between the th MT and the th cooperative MT after the weight adaptation; and are the measurement distance from the th BS to the th MT, and the distance from the th cooperative MT to the th MT; and are produced by the Gaussian distribution with the means and , and with their corresponding variances. The meanings of function nodes are described below. The messages are processed and passed between variable nodes and function nodes.

Figure 3: FG in CGL (take one of the MT as an example).

4.1. The Functions of All Nodes in FG-Based CGL

(1) Variable and (in the th iteration)

The message from variable node to function node or is a Gaussian probability density function (PDF) of and can be expressed as

where , , , . is the number of whole MTs to be located simultaneously. and represent the number of BSs and other cooperative MTs connecting to the th MT, respectively. and are Gaussian PDF of . Note that the product of any Gaussian pdf is also Gaussian and can be derived as [20]

where , , and . Similarly, variable node can be calculated.

(2) Variable Nodes , , , , , , , , , and

Variable node and can transmit the receiving message directly as follows:

Similarly, variable node , , , , , , , and all play the role as a message transfer between the corresponding two function nodes.

(3) Function Nodes , , , and

The function of these nodes is to convert the relative position information into the absolute position information, and vice versa

where and represent the position of the th BS connecting to the th MT.

where and are the variance of the Gaussian soft information from and . Similarly, the message from , , and can be calculated.

(4) Function Nodes and

The function of these nodes is to merge the separate information from -coordinate group and the -coordinate group and then compare it with real measurement data as follows:

updates itself according to the messages from and ; updates itself according to the messages from and .

where and are the measurement information between BSs and MTs, and between MTs. Similarly, and can be calculated.

(5) Function Nodes and

The function of these nodes is to perform the weighting process based on information reliability as follows:

where is the weighting coefficient and (without weighing) is for the information from BS. Similarly, , , , , , and can be calculated.

(6) Function Nodes and

The measurement information and enter FG through the function node and as follows:

After iterations, information in variable nodes and converges. and can be estimated.

4.2. The Characteristics of FG-Based CGL

From Figure 3, it can be concluded that the passing messages from BSs or cooperative MTs only merge in variable nodes and , while the passing message in other variable nodes is just related with the information in their local branches. The estimation for (or ) is , where . represents the mean of the information in the branch (no matter the BS branch or the MT branch) connecting to the th MT, and, represents the variance of the information in the branch connecting to the th MT. Moreover, it is supposed that the variances of the information from BSs are the same (denoted as , and the variances of the information from cooperative MTs are the same as well (denoted as . Therefore, can be written as , where denotes the number of BSs connecting to the th MT, and denotes the number of cooperative MTs connecting to the th MT. The estimated localization error is

where is the real value.

According to the form of (26), the localization error does not decrease linearly with the growing number of BSs or the number of cooperative MTs, and it is affected by the value of . When the number of BSs and MTs is definite, the localization error can reduce with the increasing value of . Moreover, the influence of and on localization accuracy weakens with the increasing value of . Therefore, the strict demand for topology is relieved which embodies the robustness of the algorithms.

4.3. The Formulation of WFG-CGL Algorithm

Obviously, the information provided by the cooperative MT is less reliable than that of BS. Generally, the variances information represents the reliability. The information reliability lowers with the larger variance. If the variance information is multiplied by a weighted coefficient based on its reliability before performing the position estimation, the higher localization accuracy can be obtained. Considering the influence from the weighted coefficient, a weights optimization method is proposed, which can calculate the optimal weights for multiple MTs, respectively.

Denote the localization error as and

where and are the functions of ; and are the real position for the th MT.

Take one element of the matrix in (27) as an example to calculate the corresponding weight.

Denote and represent the mean and variance of the Gaussian variable , , and represent the mean and variance of the Gaussian variable ; and are the constants.

Therefore, and follow the noncentral Chi-squared distribution. , where is the freedom degrees of the th MT in -coordinate and . Due to the fact that is independent of as follows:

where , , , and are the functions of .

The optimal weight adaptation achieves minimum localization error for multiple MTs, respectively,

where and are based on (15) and (30).

where only and are the functions of . Numerical methods can be utilized to obtain estimation.

According to the above analysis, the procedures of WFG-CGL algorithm are described as follows.

Step 1 (initialization). All variable nodes need deterministic initial values, except and . Initial values can be made at random or by some specific algorithms, not detailed here. Set .

Step 2 (entrance of measurement information in the th iteration). The measurement information and from TOA measurement enter FG through function nodes and , according to (24).

Step 3 (uplink calculation in the th iteration). Function node updates and according to (14), and (15). Similarly, function node updates and . is calculated in (35), and function nodes and update , , , and according to (18) and other similar formulae. Function nodes , , , and update the messages for and according to (12) and other similar formulae.

Step 4 (downlink calculation in the th iteration). Variable nodes transmit message to and according to (3). Similarly, transmits message to and . Function node updates according to (11). Similarly, function nodes , , and update , , and . Function node updates and according to (22) and (23). Similarly, function node updates and .

Step 5. In the th iteration, Variable nodes and are updated and estimate the position of the th MT according to (25).

Step 6. Targeting for the th MT who has not updated the position in the th iteration, we make , return to Step 2. Else, when all of the cooperative MTs in a TG have completed the position updating in the th iteration, we make .

Step 7. If , return to Step 2. Else, and , the algorithm stops.

4.4. The Formulation of LCFG-CGL Algorithm

In WFG-CGL algorithm, the solution to the weights increases the calculation complexity. In order to further diminish the complexity, the weights for all MTs in TG are not differentiated, leading to the fact that all weights are considered as 1. Therefore, we propose LCFG-CGL algorithm to solve the problem, and the procedures of LCFG-CGL algorithm are described as follows.

Step 1 (initialization). All variable nodes need deterministic initial values, except and . Initial values can be made at random or by some specific algorithms, not detailed here. Set .

Step 2 (entrance of measurement information in the th iteration). The measurement information and from TOA measurement enter FG through function nodes and , according to (24).

Step 3 (uplink calculation in the th iteration). Function node updates and according to (14) and (15). Similarly, function node updates and . All are set as 1, and function nodes , and work as transfers to , , , and . Function nodes , , , and update the messages for and according to (12) and other similar formulae.

Step 4 4 (downlink calculation in the th iteration). Variable nodes transmits message to and according to (3). Similarly, transmits message to and . Function node updates according to (11). Similarly, function nodes , , and update , , and ; Function node updates and according to (22)and (23). Similarly, function node updates and .

Step 5. In the th iteration, variable nodes and are updated and estimate the position of the th MT according to (25).

Step 6. Targeting for the th MT who has not updated the position in the th iteration, we make , return to Step 2. Else, when all of the cooperative MTs in a TG have completed the position updating in the th iteration, we make .

Step 7. If , return to Step 2. Else, and , the algorithm stops.

5. CRLB Analysis

The CRLB is commonly used as a performance benchmark of an estimator because it gives the lowest possible variance [27]. Although the CRLB analysis is a classical localization metric, the CRLB analysis targeting CGL has not been performed in detail. Here, we make the CRLB analysis of CGL firstly.

Considering a vector of all nodes parameters in CGL, nodes represent the blindfolded MTs, and nodes represent the BSs. The unknown MT parameter vector is , where for . Note that the BS parameter vector is known. Nodes and make pair-wise observations with density . Let =: node makes pair-wise observations with node . By convention, a node cannot make pair-wise observation with itself, so that . By symmetry, if , then .

It is assumed by reciprocity that . Thus, it is sufficient to consider only the lower triangle of the observation matrix when formulating the joint likelihood function. We assume that are statistically independent for . This assumption can be somewhat oversimplified but necessary for analysis. The log of the joint conditional PDF is

The CRLB on the covariance matrix of any unbiased estimator is cov, where the Fisher information matrix (FIM) is defined as

where the diagonal elements for of can be simplified to a single sum over since there are term in (38) that depend on . The off-diagonal elements can be further reduced: when for ; , there is at most one summand in (38) that is a function of both and .

where is an indicator function: 1 if or 0 otherwise.

For a 2D system, , where . The relative localization problem corresponds to the estimation of blindfolded device coordinates , where , , and the known BSs coordinates are . In the cellular network, , ; is the measured TOA between BSs and MTs. , ; ; is the measured TOA between MTs. Assuming that both and are independent Gaussian distributed, which is denoted as

where ; .

where ; ; .

where is the speed of light propagation, and denote the variance of TOA measurements between BSs and MTs, and between MTs, respectively.

For simplification, the FIM will have a similar form to (39) if partitioned into blocks

where is given by (39) using only the parameter vector , and is given by (39) using only . The off-diagonal blocks and are similarly derived. According to this cooperative structure, the elements of the submatrics of (47) are derived as follows:

And the FIM is for

for ,

where is 1 when node and node can make pair-wise measurements, else is zero. is 0 when or , else is 1. is 1 when , else is 0.

The formula derivation of FIM has been completed in (47)~(52), (53) and (54). Then, the trace of the covariance for the th blind MT satisfies

6. Simulation Results

In this section, numerical simulations include two parts. Part I is the research on important factors for the proposed algorithms. Part II is performed to show the effectiveness of the proposed WFG-CGL and LCFG-CGL algorithms under two typical topologies for localization research, including global rigidity which has the unique solvability shown in Figure 4(a), and rigidity which has no unique solvability and converges in the finite positions shown in Figure 4(b). The topologies with three MTs are taken as examples due to the fact that the simple and specific topologies are easy to follow. The comparison of localization accuracy and complexity between the proposed algorithms and the existing Taylor-series-based CGL algorithm [18] is performed. Moreover, CRLB analysis is conducted.

(a) CGL topology(global rigidity)

(b) CGL topology (rigidity)

Figure 4: CGL topology.

Any topology scenario is constructed by the random birth of MTs in the efficient range of BSs. The average error is evaluated over 10,000 independent trials. While in the different scenarios with the same number of MTs, the positions of MTs and BSs are the same. First, we consider the edges existing between BSs and MTs or between MTs LOS paths. The distance measurements derived from TOA between MTs and BSs are , the distance measurements between MTs are , where denotes the real distance between BS and MT, denotes the real distance between MTs, and and denote the corresponding measurement error and follow the Gaussian distribution with zero mean value. Moreover, in urban environment, we consider NLOS propagation case. The NLOS error was modeled as the exponential distribution random variable in the literatures [28, 29]. The probability density function (PDF) of NLOS error follows

where denotes the mean of . Referring to the literature [29], is set as 100 m.

The basic parameter configuration is listed in Table 1. The localization accuracy is analyzed by means of root mean square error (RMSE) between the real MTs’ position and the estimated location , that is,

Table 1: Parameters configuration.

The number of BSs ()
3

The number of MTs in a TG ()
3

Cell radius (m)
3000

TG radius (m)
750

6.1. The Important Factors for FG-Based CGL

(1) The Iteration Times

for the FG-based localization, the iteration times is predefined and has strong influence on the localization accuracy. Apparently, more iteration times can make the localization more precise. However, excessive iterations will lead to the increasing complexity. Therefore, the reasonable value of should be investigated.

The intuitive results can be obtained from Figures 5 and 6 which shows the localization error in terms of the iteration times from to with the topologies of Figures 4(a) and 4(b) in LOS and NLOS scenarios. It can be observed from Figures 5 and 6 that no matter LOS or NLOS scenario, the performance has no significant improvement with the growth of after in both the rigidity and the global rigidity. Considering the high accuracy and low complexity simultaneously, is adopted in the following simulations.

(a) Iteration times

(b) Iteration times

Figure 5: The influence of iteration times in LOS.

(a) Iteration times

(b) Iteration times

Figure 6: The influence of iteration times in NLOS.

(2) The Influence of and

Figures 7 and 8 indicate the localization error of LOS and NLOS paths between MTs and BSs, respectively. As shown in Figure 7(a) and Figure 8(a) where  m² in LOS and NLOS scenarios, and Figure 7(b) and Figure 8(b), where  m² in LOS and NLOS scenarios, they indicate that the standard deviation between MT and BS has stronger effects on the localization accuracy than the standard deviation between MTs for the proposed algorithms in both LOS and NLOS scenarios. The main reason is attributive to the fact that BS can provide the information of higher accuracy, compared to the cooperative MT. However, the information in NLOS scenario is more inaccurate than that in LOS scenario due to the fact the influence of in NLOS is slighter than that in LOS. This can be illustrated in Table 2.

Table 2: Statistics analysis for difference .

:

Global rigidity
Rigidity

LCFG
LOS
(0.134, 0.143)
(0.200, 0.162)

NLOS
(0.131, 0.002)
(0.198, 0.002)

WFG
LOS
(0.219, 0.057)
(0.175, 0.101)

NLOS
(0.159, 0.002)
(0.167, 0.001)

(a)

(b)

Figure 7: RMSE of difference in LOS.

(a)

(b)

Figure 8: RMSE of difference in NLOS.

Denote , where represents the RMSE in , and represents the RMSE in . and represent the mean and variance of .

The statistic data can be obtained from Table 2 which shows and in terms of from  m to  m for the two proposed algorithms in different scenarios. No matter the LCFG-CGL algorithm or the WFG-CGL algorithm, no matter the global rigidity topology or the rigidity topology, and no matter the LOS or NLOS, the values of and are all small which verifies the little performance difference between  m² and  m². Thereinto, we just take  m² as an example in the following simulations.

(3) The Influence of

According to the analysis in Section 4.2 and the simulation results in Figure 9, the localization error decreases with the growing number of the cooperative BSs, while the decline rate of the localization error becomes smaller with the growing number of the cooperative BSs. It is assumed that the variance of the measurement error between MTs is definite. Figure 9 illustrates that the value of reduces with the rising of , and the influence of on localization accuracy weakens with the increasing value of , no matter in LOS and NLOS scenarios. The simulation results consist with the analysis in Section 4.3, which can verify the robustness of our algorithms for the topology of CGL graph. It also clarifies the localization accuracy cannot be reduced through increasing the amount of measurement information (the number of BSs) when the measurement information is inaccurate. Therefore, the localization accuracy can be improved, only through optimizing the algorithm. Fortunately, our proposed algorithms can achieve this.

(a)

(b)

Figure 9: RMSE of difference .

6.2. The Localization Accuracy Analysis

The simulation results are shown in Figures 10 and 11, where the average localization errors versus of the range measurement errors are compared in LOS and NLOS scenarios. The comparison is made among the Taylor-series-based CGL algorithm, the LCFG-CGL algorithm, the WFG-CGL algorithm and CRLB in the same scenario, where the Taylor-series-based CGL algorithm is analyzed in [18]. As expected, no matter the global rigidity or the rigidity, (1) the average localization errors increase when the range measurement error increases, (2) the proposed algorithms are superior to the traditional Taylor series CGL algorithm, (3) the WFG-CGL algorithm has higher localization accuracy than LCFG-CGL algorithm primarily due to its exploitation of the information reliability as the weights, (4) in NLOS scenario, LCFG-CGL algorithm and WFG-CGL algorithm perform far well than the traditional Taylor-series CGL algorithm, and (5) especially, the WFG-CGL algorithm is approximate to CRLB results in both global rigidity and rigidity regarding the LOS scenario.

(a)

(b)

Figure 10: RMSE of different topologies in LOS.

(a)

(b)

Figure 11: RMSE of different topologies in NLOS.

Denote , where represents the RMSE in rigidity topology, and represents the RMSE in global rigidity topology. and represent the mean and variance of .

Remarkably, Taylor-series-based CGL algorithm performs better in global rigidity topology than in rigidity topology. It is obvious that the global rigidity topology can provide more localization information. However, the two proposed algorithms in rigidity have the similar localization accuracy in global rigidity, as shown in Figures 10 and 11. Table 3 derived from Figures 10 and 11 which shows the and in terms of from  m to  m in different scenarios and different algorithms. It denotes that the two proposed algorithms in rigidity have the similar localization accuracy in global rigidity, which demonstrates the strong robustness of the proposed algorithms for different CGL topologies.

Table 3: Statistics analysis for different topologies.

Taylor series
LOS
0.41
0.0021

NLOS
0.409
0.002

LCFG
LOS
0.12
0.0002

NLOS
0.11
0

WFG
LOS
0.09
0.0002

NLOS
0.10
0

CRLB
LOS
0.08
0

NLOS
0.002
0

The localization accuracy depends on how to utilize the available information. Generally, the variances can represent the practical environment effectively. The proposed algorithms make full use of the variances in nearly each node of FG, which is considered as a great help for performance upgrade in rigidity topology, even with less information than global rigidity. Furthermore, the power of the proposed algorithms lies in the ability to distributively process the problem and to provide the optimum or near-optimum solutions, which makes it a good approach to solving distributed problem. On the other hand, CGL is a distributed problem, that is, measurement parameters are only locally related. Therefore, the excellent combination between FG and CGL can relax the CGL topology requirement under ensuring the localization accuracy to some extent.

6.3. Complexity Analysis

This section analyzes the computation complexity of the proposed algorithms and Taylor-series-based algorithm, respectively.

(1) Complexity Analysis of Taylor-Series-Based CGL Algorithm

In Taylor-series-based CGL algorithm [18], the most complex step is the solution to pseudoinverse of large scale matrix. In each iteration, the computation complexity is , where is the number of edges in CGL topology. Moreover, the iteration number is not definite, because it is determined by the tolerable localization error and the distance measurement error. Especially when the TG includes more cooperative MTs, the complexity of CGL-based Taylorseries localization algorithm increases sharply.

(2) Complexity Analysis of WFG-CGL Algorithm and LCFG-CGL Algorithm

In the two proposed FG-based CGL algorithms, the amount of messages passed between nodes is limited. And the iteration times is pre-defined, commonly . Because after , the performance has no significant improvement with the growth of . The performance of the proposed algorithms approaches high accuracy after 7 iterations, which shows that the convergence speeds of the two proposed algorithms are reasonably high.

Although the solution of the weights in WFG-CGL algorithm results in the increasing calculation complexity to some extent, they can be calculated distributively, which means that it just solves multiple equations, respectively, instead of the solution to pseudo-inverse of large scale matrix. The computation complexity of weights solution is , which is less than the Taylor-series-based CGL algorithm.

On the other hand, the convergence probabilities of the two proposed algorithms are shown in Figure 12. The proposed algorithms have very high convergence probabilities in both LOS and NLOS scenarios. Even when the standard derivation is big or the topology is not a unique solvability, the convergence probabilities of the two proposed algorithms are still acceptable. Those all lead to the low computational complexity and transmission loads while ensuring the localization accuracy.

From the above analysis and simulation results, it can be concluded that (1) both of the proposed algorithms can relax the CGL topology requirement and provide high localization accuracy. (2) WFG-CGL algorithm possesses higher accuracy than LCFG-CGL algorithm while LCFG-CGL algorithm possesses lower complexity than WFG-CGL algorithm. Although the complexity of WFG-CGL algorithm increases to some extent, it still performs with lower complexity than the traditional Taylor-series-based CGL algorithm.

(a)

(b)

Figure 12: Convergence probability of FG-based CGL.

7. Conclusion

In this paper, the WFG-CGL algorithm and the LCFG-CGL algorithm are proposed in ill conditions. Especially for WFG-CGL algorithm, the weights based on the information reliability are incorporated into the algorithm in order to obtain higher localization accuracy. The corresponding analysis and simulations show that the proposed algorithms can relax the demand for CGL topology and provide high localization accuracy under low complexity in comparison with the existing CGL algorithm. Moreover, the important factors and CRLB for FG-based localization are analyzed.

Acknowledgments

This work is supported by the National Nature Science Foundation of China Project (Grant 61001119) and Fundamental Research Funds for the Central Universities (G470712).

References

F. C. Commission, Revision of the Commissions Rules to Insure Compatibility with Enhanced 911 Emergency Calling Systems, 1996

N. Patwari, J. N. Ash, S. Kyperountas, A. O. Hero III, R. L. Moses, and N. S. Correal, “Locating the nodes: cooperative localization in wireless sensor networks,” IEEE Signal Processing Magazine, vol. 22, no. 4, pp. 54–69, 2005.

S. Gezici, Z. Tian, G. B. Giannakis et al., “Localization via ultra-wideband radios: a look at positioning aspects of future sensor networks,” IEEE Signal Processing Magazine, vol. 22, no. 4, pp. 70–84, 2005.

S. Hara, D. Zhao, K. Yanagihara et al., “Propagation characteristics of IEEE 802.15.4 radio signal and their application for location estimation,” in Proceedings of the IEEE 61st Vehicular Technology Conference (VTC '05), pp. 97–101, June 2005.

S. Feng and C. L. Law, “Assisted GPS and its impact on navigation in intelligent transportation systems,” in Proceedings of the IEEE International Conference on Intelligent Transportation Systems, pp. 926–931, June 2002.

P. Farradyne, Vehicle Infrastructure Integration (VII) Architecture and Functional Requirements, draft, version 1.0, 2005

M. Shafl, W. Aamer, and K. Butterworth, “Next generation wireless networks,” in Proceedings of the The IEEE International Symposium on Circuits and Systems (ISCAS '01), pp. 721–729, Sydney, Australia, May 2001.

Y. Zhao, “Standardization of mobile phone positioning for 3G systems,” IEEE Communications Magazine, vol. 40, no. 7, pp. 108–116, 2002.

A. H. Sayed, A. Tarighat, and N. Khajehnouri, “Network-based wireless location: challenges faced in developing techniques for accurate wireless location information,” IEEE Signal Processing Magazine, vol. 22, no. 4, pp. 24–40, 2005.

J. Chaffee and J. Abel, “GDOP and the Cramer-Rao bound,” in Proceedings of the IEEE Position Location and Navigation Symposium, pp. 663–668, Las Vegas, Nev, USA, April 1994.

E. G. Larsson, “Cramér-rao bound analysis of distributed positioning in sensor networks,” IEEE Signal Processing Letters, vol. 11, no. 3, pp. 334–337, 2004.

Y. Shen, H. Wymeersch, and M. Z. Win, “Fundamental limits of wideband cooperative localization via fisher information,” in Proceedings of the IEEE Wireless Communications and Networking Conference (WCNC '07), pp. 3951–3955, March 2007.

Y. Qi, H. Kobayashi, and H. Suda, “On time-of-arrival positioning in a multipath environment,” IEEE Transactions on Vehicular Technology, vol. 55, no. 5, pp. 1516–1526, 2006.

D. B. Jourdan, D. Dardari, and M. Z. Win, “Position error bound and localization accuracy outage in dense cluttered environments,” in Proceedings of the IEEE International Conference on Ultra-Wideband (ICUWB '06), pp. 519–524, September 2006.

J. N. Ash and R. L. Moses, “On the relative and absolute positioning errors in self-localization systems,” IEEE Transactions on Signal Processing, vol. 56, no. 11, pp. 5668–5679, 2008.

Y. Zhao, “Standardization of mobile phone positioning for 3G systems,” IEEE Communications Magazine, vol. 40, no. 7, pp. 108–116, 2002.

Y. Zhang, Q. Cui, and X. Tao, “Cooperative group localization for 4G wireless networks,” in Proceedings of the IEEE 70th Vehicular Technology Conference Fall (VTC '09), Anchorage, Alaska, September 2009.

Q. Cui and X. Zhang, “CGL-based modified taylor-series localization algorithms for future wireless networks,” in Proceedings of the IEEE International Conference on Communication Technology (ICCT '11), Jinan, China, September 2011.

J. Aspnes, T. Eren, D. K. Goldenberg et al., “A theory of network localization,” IEEE Transactions on Mobile Computing, vol. 5, no. 12, pp. 1663–1678, 2006.

F. R. Kschischang, B. J. Frey, and H. A. Loeliger, “Factor graphs and the sum-product algorithm,” IEEE Transactions on Information Theory, vol. 47, no. 2, pp. 498–519, 2001.

J. C. Chen, Y. C. Wang, C. S. Maa, and J. T. Chen, “Network-side mobile position location using factor graphs,” IEEE Transactions on Wireless Communications, vol. 5, no. 10, pp. 2696–2704, 2006.

C. Mensing and S. Plass, “TDoA positioning based on factor graphs,” in Proceedings of the IEEE 17th International Symposium on Personal, Indoor and Mobile Radio Communications (PIMRC '06), pp. 1–5, September 2006.

B. Omidali and S. A. A. B. Shirazi, “Performance improvement of AOA positioning using a two-step plan based on factor graphs and the Gauss-Newton method,” in Proceedings of the 14th International CSI Computer Conference (CSICC '09), pp. 305–309, October 2009.

C. T. Huang, C. H. Wu, Y. N. Lee, and J. T. Chen, “A novel indoor RSS-based position location algorithm using factor graphs,” IEEE Transactions on Wireless Communications, vol. 8, no. 6, pp. 3050–3058, 2009.

M. Üney and M. Çetin, “Target localization in acoustic sensor networks using factor graphs,” in Proceedings of the IEEE 16th Signal Processing, Communication and Applications Conference (SIU '08), pp. 1–4, April 2008.

H. Wymeersch, J. Lien, and M. Z. Win, “Cooperative localization in wireless networks,” Proceedings of the IEEE, vol. 97, no. 2, pp. 427–450, 2009.

N. Patwari, A. O. Hero, M. Perkins, N. S. Correal, and R. J. O'Dea, “Relative location estimation in wireless sensor networks,” IEEE Transactions on Signal Processing, vol. 51, no. 8, pp. 2137–2148, 2003.

K. Yu, I. Sharp, and Y. J. Guo, Non-Line-of-Sight Identification, Ground-Based Wireless Positioning, John Wiley & Sons, Chichester, UK, 2009.

P. C. Chen, “A non-line-of-sight error mitigation algorithm in location estimation,” in Proceedings of the IEEE Wireless Communications and Networking Conference, pp. 316–320, New Orleans, La, USA, September 1999.