Two-Phase Robust Target Localization in Ocean Sensor Networks Using Received Signal Strength Measurements

Target localization plays a vital role in ocean sensor networks (OSNs), in which accurate position information is not only a critical need of ocean observation but a necessary condition for the implementation of ocean engineering. Compared with other range-based localization technologies in OSNs, the received signal strength (RSS)-based localization technique has attracted widespread attention due to its low cost and synchronization-free nature. However, maintaining relatively good accuracy in an environment as dynamic and complex as the ocean remains challenging. One of the most damaging factors that degrade the localization accuracy is the uncertainty in transmission power. Besides the equipment loss, the uncertain factors in the fickle ocean environment may result in a significant deviation between the standard rated transmission power and the usable transmission power. The difference between the rated and actual transmission power would lead to an extra error when it comes to the localization in OSNs. In this case, a method that can locate the target without needing prior knowledge of the transmission power is proposed. The method relies on a two-phase procedure in which the location information and the transmission power are jointly estimated. First, the original nonconvex localization problem is transformed into an alternating non-negativity-constrained least square framework with the unknown transmission power (UT-ANLS). Under this framework, a two-stage optimization method based on interior point method (IPM) and majorization-minimization tactic (MMT) is proposed to search for the optimal solution. In the first stage, the barrier function method is used to limit the optimization scope to find an approximate solution to the problem. However, it is infeasible to approach the constraint boundary due to its intrinsic error. Then, in the second stage, the original objective is converted into a surrogate function consisting of a convex quadratic and concave term. The solution obtained by IPM is considered the initial guess of MMT to jointly estimate both the location and transmission power in the iteration. In addition, in order to evaluate the performance of IPM-MM, the Cramer Rao lower bound (CRLB) is derived. Numerical simulation results demonstrate that IPM-MM achieves better performance than the others in different scenarios.


Introduction
In recent years, wireless sensor networks (WSNs) have been widely used in many fields [1][2][3], such as data collection, environmental monitoring, target localization and tracking, and traffic control. One of the critical applications is the use of ocean sensor networks (OSNs) to locate ocean targets [4]. In order to achieve effective and accurate positioning, some sensor nodes that are sensitive to location indicators are usually deployed on the ocean surface or underwater, including anchor nodes and unknown nodes. It deserves indicators are usually deployed on the ocean surface or underwater, including anchor nodes and unknown nodes. It deserves to be noted that the location of some anchor nodes is fixed and pre-known, and some are mobile and equipped with the localization hardware that can provide location information. The goal is to make full use of the anchor node information to locate the target node. In this paper, buoy nodes on the ocean surface, the anchor nodes underwater with the known position, and unknown target nodes are considered in OSNs where they can communicate using the acoustic signal [Error! Reference source not found.,Error! Reference source not found.]. The whole system model for OSNs is depicted in Figure 1.  It is well known that, on account of the dynamics and complicated communication channels in the ocean environment, underwater localization is more complicated and challenging than ground positioning [Error! Reference source not found.]. Different from terrestrial radio-frequency wireless communication, the use of acoustic signal seems to be a better way in OSNs [Error! Reference source not found.]. Unfortunately, it is difficult to estimate a target with the acoustic signal due to various adverse factors, including the dynamics of the environment, multipath and absorption loss of the signal, and noise interference [7,9,10]. In this context, how to locate a target in such a complex ocean environment is a challenge.
In terms of localization methods in OSNs, it mainly includes cooperative localization (CL) scheme and non-cooperative localization (NL) scheme. The unknown target nodes can communicate with anchor nodes and other unknown nodes within a certain communication range for the former, whereas the latter cannot. Frankly speaking, the CL scheme consumes more energy and computational time, together with a higher complexity compared to NL scheme [Error! Reference source not found.]. In this case, we concentrate on the NL scheme in OSNs. Among the techniques in terms of localization in NL scheme, received signal strength (RSS)-based method has gained the upsurge of interest by researchers with its low complexity, easy implementation, cost-effective, and synchronization-free, compared with time of arrival (TOA), time difference of arrival (TDOA), angle of arrival (AOA)-based techniques, and their variants [12][13][14][15]. It is well known that, on account of the dynamics and complicated communication channels in the ocean environment, underwater localization is more complicated and challenging than ground positioning [7]. Different from terrestrial radio-frequency wireless communication, the use of acoustic signal seems to be a better way in OSNs [8]. Unfortunately, it is difficult to estimate a target with the acoustic signal due to various adverse factors, including the dynamics of the environment, multipath and absorption loss of the signal, and noise interference [7,9,10]. In this context, how to locate a target in such a complex ocean environment is a challenge.
In terms of localization methods in OSNs, it mainly includes cooperative localization (CL) scheme and non-cooperative localization (NL) scheme. The unknown target nodes can communicate with anchor nodes and other unknown nodes within a certain communication range for the former, whereas the latter cannot. Frankly speaking, the CL scheme consumes more energy and computational time, together with a higher complexity compared to NL scheme [11]. In this case, we concentrate on the NL scheme in OSNs. Among the techniques in terms of localization in NL scheme, received signal strength (RSS)-based method has gained the upsurge of interest by researchers with its low complexity, easy implementation, cost-effective, and synchronization-free, compared with time of arrival (TOA), time difference of arrival (TDOA), angle of arrival (AOA)-based techniques, and their variants [12][13][14][15].
In what concerns the RSS-based technique, various approaches have been developed to solve the localization problem for terrestrial wireless sensor networks. In [16], the authors proposed a novel RSS-based indoor localization method using belief function theory, in which they modeled RSS by non-Gaussian probability density functions and put forward decision-making modes to improve the localization accuracy. In [17], the authors formulated the localization problem into a nonlinear weighted least squares (WLS) problem, where the unscented transformation and the bisection method were adopted to estimate the transmit power and the target location jointly. In [18], a suboptimal estimator has been derived from dealing with the localization problem. Unfortunately, the suboptimal estimator seems to perform poorly under the visible bias and variance in a noisy environment. In this context, considerable efforts have been made to improve performance [19][20][21][22][23]. The authors in [20] have employed a hybrid model that combined the distance and angle measurements to reduce the estimation error. Nevertheless, the hybrid information extraction may increase the network complexity, albeit the localization accuracy has been improved. Chan et al. in [22] have designed a WLS estimator to solve the 3-D AoA/RSS Difference (RSSD) non-cooperative localization problem in a relatively small-scale WSN with low noise power. In addition, to multiple-target localization in 3-D WSNs, a linear quadratic estimator was presented in [23], which integrated RSS and AOA measurements and utilized spherical representation of the considered problem to acquire a closed solution.
However, different from the terrestrial WSNs, it remains challenging to locate targets in OSNs using RSS-based methods due to the dual factors of path loss and absorption loss for the underwater acoustic signals. The dual factors may negatively impact localization in OSNs and limit the capability of the RSS-based methods. In this case, substantial efforts have been made to improve the localization performance using RSS-based methods in OSNs [24][25][26][27][28]. In [24], the maximum likelihood (ML) problem regarding the RSSbased underwater target localization was discussed, and a series of methods were derived using the l 1 norm, but the accuracy is not satisfactory. To further improve accuracy, the authors in [25] transformed the original localization problem into a generalized trust-region subproblem (GTRS) and developed a novel weighted least squares (NWLS) with a known transmit power scheme (NWLS-K) and an unknown transmit power scheme (NWLS-U). Moreover, the authors in [26] presented a different way to solve the localization problem with an SD/SOCP scheme at the tremendous iteration cost to achieve high location accuracy. In [28], the authors have designed a robust, non-cooperative localization algorithm (RNLA); although it can reduce the estimated error in outliers, the signal loss resulting from the absorption was not considered. Overall, some methods like NWLS or SD/SCOP could improve the localization accuracy in a way, the performance might not be suitable for the dynamic scenario. Some methods like RNLA or LS may provide a fast and feasible solution, but the localization accuracy seems unsatisfactory. In this case, we would like to develop an exclusive method for the dynamic ocean environment with a relatively high localization accuracy for target localization.
In the paper, a two-phase robust target localization method in OSNs under unknown transmission power is proposed. First, the original localization problem is described by taking the maximum likelihood criterion; then, it is transformed into an unknown transmittingpower alternating non-negativity-constrained least square problem (UT-ANLS) by applying some approximations. Furthermore, a method based on the interior point method (IPM) is developed to acquire the solution. However, the accuracy of the solution obtained by IPM may be far away from the accurate one. In this case, a majorization-minimization tactic (MMT) is proposed to refine the transmission power and target location simultaneously in the iteration in which the solution of IPM is as the initial guess.
The main contributions of the paper can be summarized as follows: (1) The considered localization problem is converted into a UT-ANLS framework by exploiting certain approximations. (2) A two-phase optimization method, i.e., IPM-MM, is described to solve the problem, wherein the first is the IPM, and the second is the MMT. The rest of the paper is arranged as follows. The problem is formulated in Section 2. In Section 3, the proposed method, IPM-MM, is illustrated. Numerical simulation results in different scenarios are shown and discussed in Section 4. In the last section, we make a summary of the paper.

Problem Formulation
Consider a three-dimensional (3-D) OSN with plenty of buoy nodes deployed on the ocean surface including M anchor nodes and an unknown node. Without loss of generality, the ith anchor can be denoted as a i = [a i1 , a i2 , a i3 ] T , where i = 1, 2, . . . , M. The unknown target node can be represented as u = [u 1 , u 2 , u 3 ] T , where T stands for the transpose operation. We suppose that the RSS information propagates between the ith anchor node and target node following the shadowing model, which could be modeled as [24,26,29]: where P si is the signal strength received through the ith anchor node from the target node, d 0 (d 0 ≤ u − a i ) is the reference distance and it is usually set as 1 m, P 0 denotes the transmit power of the target when the distance is d 0 , α is the path loss exponent under different environments, · is 2 norm, and n i is the noise with zero mean and variance σ 2 i . α f is the absorption coefficient which can be obtained from Thorp's formula using a frequency function, denoted as [30]: (2) Figure 2 shows the relationship of α f with the variable frequency f. According to [31,32], most acoustic systems operate in the 0-30 kHz frequency band. In addition, it deserves to be noted that a standard terrestrial localization using RSS measurements can be obtained while α f equals to 0. ploiting certain approximations.
(2) A two-phase optimization method, i.e., IPM-MM, is described to solve the problem, wherein the first is the IPM, and the second is the MMT. (3) A closed-form expression of the Cramer-Rao lower bound (CRLB) of the measurement model is derived under the unknown transmission power.
The rest of the paper is arranged as follows. The problem is formulated in Section 2. In Section 3, the proposed method, IPM-MM, is illustrated. Numerical simulation results in different scenarios are shown and discussed in Section 4. In the last section, we make a summary of the paper.

Problem Formulation
Consider a three-dimensional (3-D) OSN with plenty of buoy nodes deployed on the ocean surface including M anchor nodes and an unknown node. Without loss of generality, the i th anchor can be denoted as The unknown target node can be represented as , where T stands for the transpose operation. We suppose that the RSS information propagates between the i th anchor node and target node following the shadowing model, which could be modeled as [24,26,29]: where si P is the signal strength received through the i th anchor node from the target node, 4100 (2) Figure 2 shows the relationship of αf with the variable frequency f. According to [Error! Reference source not found.,Error! Reference source not found.], most acoustic systems operate in the 0-30 kHz frequency band. In addition, it deserves to be noted that a standard terrestrial localization using RSS measurements can be obtained while αf equals to 0.  Let z i = α f u − a i , then Equation (1) can be written as (3). Given the observation matrix P = [P s1 , P s2 , . . . , P s M ] T , and n i follows a Gaussian distribution with zero-mean and variance σ 2 i , the probability density function (PDF) can be expressed as [28]: Maximizing the joint PDF, the maximum likelihood estimator can be derived as: The above problem is difficult to solve directly due to its highly non-convexity, which motivates us to simplify it using an approximation method to find the sub-optimal solution.

UT-ANLS Framework
First, we conduct a transformation from (1) as [20,33,34]: , then when the noise is relatively small Equation (6) can be rewritten as: Then, the estimation problem in Equation (5) is converted to: We square Equation (8) and expand each term, then: where κ = u 2 .
Let y = u T , 10 P 0 5α , κ T be the estimation vector, the problem in (9) is finally transformed into an equivalent ANLS problem under UT in (10): where: The problem can be further changed to minimize a quadratic objective function under the positive constraint.

PHASE1: Interior Point Method (IPM)
IPM [35] is numerically an effective solution for convex quadratic programs, which can obtain a feasible possible solution for the UT-ANLS framework. A finite number of iterations are conducted to achieve an ideal value. In this section, we introduce the IPM to solve the target localization problem based on the UT-ANLS framework.
According to IPM, we first redescribe the formula (10) by one logarithmic barrier function to transform the original problem into an unconstrained one: where µ is a penalty factor (µ > 0) and the feasible region is denoted by otherwise, it tends to be infinite, which results in φ(y) being infinite. We get the estimated value through multiple iterations, and the barrier parameter of the k th iteration can be expressed as µ k . The above formula can be rewritten as: The first-order and second-order partial derivative can be denoted as ∇φ y (y, . . , y M ), e = (1, 1, . . . , 1) T (the number is M), then the Newton-Raphson method is adopted to search the KKT solution. Let , the KKT condition can be: Then, the problem can be solved by: where W k = ∇φ 2 yy (y k , µ k ), d y k and d z k are the Newton step. Once an initial point that satisfies the above system is found, a line search could be carried out along d y k to find a new iteration value, following y k+1 = y k + α k · d y k , where α k is the step length. Similarly, the update of z satisfies z k+1 = z k + α k · d z k . When the tolerance of all conditions holds, the final solution y * can be obtained. The whole process of IPM is shown in Algorithm 1.

Algorithm 1
The process of IPM Although the interior point method has a relatively low memory consumption, the accuracy of its solution may be unsatisfactory because the barrier function prevents the iteration from approaching the boundary of the inequality constraint.

PHASE2: Majorization-Minimization Tactic (MMT)
For IPM, the performance is good only when it satisfies the given constraint, once falling on the edge of restriction or out the range, it may directly tend to be an approximate local optimum solution or no solution. To acquire the best estimation, an iteration strategy, i.e., MMT, is used to optimize the solution further using the initial point computed by IPM. The whole process can be seen in Figure 3. Although the interior point method has a relatively low memory consumption, the accuracy of its solution may be unsatisfactory because the barrier function prevents the iteration from approaching the boundary of the inequality constraint.

PHASE2: Majorization-Minimization Tactic (MMT)
For IPM, the performance is good only when it satisfies the given constraint, once falling on the edge of restriction or out the range, it may directly tend to be an approximate local optimum solution or no solution. To acquire the best estimation, an iteration strategy, i.e., MMT, is used to optimize the solution further using the initial point computed by IPM. The whole process can be seen in Figure 3. According to [Error! Reference source not found.], we redescribe the problem (6) to construct a surrogate function:    According to [36], we redescribe the problem (6) to construct a surrogate function: where g(u) is obviously a convex function, and h(u) is concave. Geometrically, the concave function h(u) can have an upper bound by using a linear function. Then, Q(u, P 0 ) is made up of g(u) and the linear term. To acquire the linear term, first order Taylor expansion [37] is used for h(u) at u (d) where the iteration number is d.
Therefore, we denote Q u u (d) , P 0 (d) by: Then, P 0 (d) is fixed as P 0 using IPM, the u (d+1) can be updated as: We can obtain the concrete u (d+1) as: Then we fix u to derive P 0 (d+1) . The update rule for P 0 (d+1) is: Enough attention should be paid to two conditions for its convergence: (a) The function Q(u, P 0 ) is bounded below; (b) the value of Q(u, P 0 ) is non-increasing, i.e., which implies it has a convergence limit Q (l) . To prove the two conditions (a) and (b), we describe the process in detail as follows.
To (a): According to the formula (16), one has Q(u, P 0 ) ≥ −∑ M i=1 d 2 0 10 2γ i ,which shows that it has a lower bound.
To (b): Q(u, P 0 ) consists of the concave function h(u) and convex function g(u). We use Θ = [u, P 0 ] T , then Q(u, P 0 ) and h(u, P 0 ) can be denoted by Q(Θ) and h(Θ). To find the upper bound of Q(Θ), we take the difference between formula (16) and formula (18) at Θ (d) , then we can have: Because: We find that (23) is an affine function plus a convex term, which shows that

Cramer-Rao Lower Bound (CRLB) Analysis
CRLB [38] is a covariance matrix that represents a lower bound of any unbiased linear estimators. It is usually used as the benchmark for the estimators and can be defined as the trace of the inverse of the Fisher information matrix (FIM): where ∑ is diag(σ 2 1 , . . . , σ 2 M ) = diag(σ 2 ): where: ϑ = − 10α ln 10 . Let: Then, the FIM can be denoted as: Suppose that ŷ i − y i = er, the root mean square error (RMSE) can be related to CRLB by:

Numerical Simulation Results
In this section, several numerical experiments are conducted in different scenarios to evaluate the performance of IPM-MM. To verify the effectiveness of our proposed method, we compare it with LS [19], SRWLS [20], ASM [33], RNLA [28], and CRLB in MATLAB under 3-Dimension. Due to the ocean environment dynamics, we assume that the target and anchor nodes' positions are not fixed at a cube area with side length Sel in each Monte Carlo trial. Unless otherwise stated, the maximum number of iterations I m is set to be 1000, the reference distance d 0 is 1 m, and Monte Carlo (MC) runs are 1000. The root mean square error (RMSE) is adopted to calibrate localization accuracy, which is defined as: wherex andP 0 are the estimates of the location x and power P 0 in the i th Monte Carlo trial, MC is the total number of Monte Carlo trials. Before presenting the simulation experiment, some necessary instructions must be made. The scenarios that RMSE versus varying path loss exponent α, SNR, anchor nodes M, absorption coefficient α f , side length Sel and cumulative distribution function (CDF) of different algorithms are performed meeting the conditions of SNR = 10 · log 10 x − a i / Mσ 2 , whose details can be found in [20,33]. It is worth noting that the environmental noise in the ocean mainly includes the noise generated by factors such as tides and sea static pressure, seismic disturbance, boats and sea waves, and they usually obey a Gaussian distribution. Therefore, we use Gaussian noise to simulate the ocean noise. However, the influence of noise on localization cannot be seen intuitively from the previous simulation experiments. Hence, we supplement the experimental scenario in Section 4.7, wherein the relationship between parameter SNR and noise α 2 does not hold.

Scenario with Varying Path Loss Exponent α
The RMSE versus varying path loss exponent α is illustrated in Figure 4. The specific parameters are set in Table 1. It is worth mentioning that some unlisted parameters remain the same as before. Interestingly, RNLA seems more sensitive than others. It can be observed that almost all methods have acquired the relatively acceptable robustness with increasing α except RNLA, in which IPM-MM outperforms the others for most cases. It may be less satisfactory for the performance of the estimate of x and P 0 when α = 3.5, however, with the increase of α, our method shows a better performance than the rest. For the single IPM, without MMT, it is less ideal compared with IPM-MM according to the simulation results, its performance is just better than RLNA when α ∈ [3.5, 6] shown in Figure 4a, and similarly in Figure 4b, it may not perform well among all methods, but when α = 4, its performance of estimating P 0 is at a moderate level.

Scenario with Varying Path Loss Exponent α
The RMSE versus varying path loss exponent α is illustrated in Figure 4. The specific parameters are set in Table 1. It is worth mentioning that some unlisted parameters remain the same as before. Interestingly, RNLA seems more sensitive than others. It can be observed that almost all methods have acquired the relatively acceptable robustness with increasing α except RNLA, in which IPM-MM outperforms the others for most cases. It may be less satisfactory for the performance of the estimate of x and P0 when α = 3.5, however, with the increase of α, our method shows a better performance than the rest. For the single IPM, without MMT, it is less ideal compared with IPM-MM according to the simulation results, its performance is just better than RLNA when [ ]

3.5, 6
α ∈ shown in 4a, and similarly in 4b, it may not perform well among all methods, but when α = 4, its performance of estimating P0 is at a moderate level.      Figure 5 depicts the RMSE of different methods with growing SNR. It is noteworthy that unless otherwise stated, all fixed parameters are the same as before. Other parameter settings are shown in Table 2. Obviously, the performance of all methods is improved when the size of SNR increases. Compared with LS, SRWLS, ASM, RNLA, our IPM-MM is closer to CRLB. Especially when SNR is 54 dB, our method performs better than others for both estimating x (shown in Figure 5a) and P 0 (shown in Figure 5b). In addition, the performance gap between IPM and IPM-MM seems a little bit big, IPM-MM performs better than IPM, which demonstrates the combination is effective. Although when SNR = 51 dB, it may not evident in Figure 5b, yet it is apparent in Figure 5a. Overall, other schemes are not as good as our method in the estimation of x, and are not as close to CRLB as our scheme for estimating P 0 .

Scenario with Varying SNR
M (The anchor number) 7 Sel (Side length) 100 m Figure 5 depicts the RMSE of different methods with growing SNR . It is noteworthy that unless otherwise stated, all fixed parameters are the same as before. Other parameter settings are shown in Table 2. Obviously, the performance of all methods is improved when the size of SNR increases. Compared with LS, SRWLS, ASM, RNLA, our IPM-MM is closer to CRLB. Especially when SNR is 54 dB, our method performs better than others for both estimating x (shown in 5a) and P0 (shown in 5b). In addition, the performance gap between IPM and IPM-MM seems a little bit big, IPM-MM performs better than IPM, which demonstrates the combination is effective. Although when SNR = 51 dB, it may not evident in 5b, yet it is apparent in 5a. Overall, other schemes are not as good as our method in the estimation of x , and are not as close to CRLB as our scheme for estimating P0.

Scenario with Varying Anchor Nodes M
In this scenario, the simulation of considered methods is carried out with varying anchor nodes M. Except for some previously fixed parameters, other parameters are indicated in Table 3. The RMSE results are shown as follows. It is deserved to be pointed out that the more location information available with the increase of anchor nodes M. Therefore, when anchor nodes increase from 6 to 12, the performance is improved.  Figure 6, when the anchor's number is on the rise, among all the methods, our proposed IPM-MM obtains a relatively satisfactory performance, in which the results are closer to the CRLB compared to others. Especially when the anchor number M = 8 or 12, its outperformance is remarkable. IPM may seem less robust than others under current conditions, but adding MMT surpasses the other methods and gets a perfect improvement.

Scenario with Varying Anchor Nodes M
In this scenario, the simulation of considered methods is carried out with varying anchor nodes M . Except for some previously fixed parameters, other parameters are indicated in Table 3. The RMSE results are shown as follows. It is deserved to be pointed out that the more location information available with the increase of anchor nodes M. Therefore, when anchor nodes increase from 6 to 12, the performance is improved.
From Figure 6, when the anchor's number is on the rise, among all the methods, our proposed IPM-MM obtains a relatively satisfactory performance, in which the results are closer to the CRLB compared to others. Especially when the anchor number M = 8 or 12, its outperformance is remarkable. IPM may seem less robust than others under current conditions, but adding MMT surpasses the other methods and gets a perfect improvement.

Scenario with Varying the Absorption Coefficient αf
Absorption bias plays an important role in Equation (3). It is necessary to discuss the relationship between the RMSE and absorption coefficient αf, which is displayed in Figure  7. Apart from some previously fixed parameters, other parameters are displayed in Table  4. Obviously, RNLA seems more sensitive to αf than others, and αf has a great influence on each scheme. Nevertheless, we observe that our solution acquires the relatively acceptable robustness among all the schemes and is still optimal to target location estimation. Figure  7b, in which IPM-MM outperforms the others for all cases. For the single IPM, without

Scenario with Varying the Absorption Coefficient α f
Absorption bias plays an important role in Equation (3). It is necessary to discuss the relationship between the RMSE and absorption coefficient α f , which is displayed in Figure 7. Apart from some previously fixed parameters, other parameters are displayed in Table 4. Obviously, RNLA seems more sensitive to α f than others, and α f has a great influence on each scheme. Nevertheless, we observe that our solution acquires the relatively acceptable robustness among all the schemes and is still optimal to target location estimation. Figure 7b, in which IPM-MM outperforms the others for all cases. For the single IPM, without MMT, it is not ideal compared with IPM-MM according to the simulation results; its performance is just at a moderate level among all methods. With the increase of α f , our improvements on the basis of IPM could be more visible. each scheme. Nevertheless, we observe that our solution acquires the relatively acceptable robustness among all the schemes and is still optimal to target location estimation. Figure  7b, in which IPM-MM outperforms the others for all cases. For the single IPM, without MMT, it is not ideal compared with IPM-MM according to the simulation results; its performance is just at a moderate level among all methods. With the increase of αf, our improvements on the basis of IPM could be more visible.

Scenario with Varying Side Length Sel
Due to the fact that the absorption bias is determined by the absorption coefficient f α and the distance between each anchor and the target node, it is essential to carry out the simulation experiment. The related parameters are shown in Table 5. According to Figure 8, when the side length ranges from 70 m to 120 m, except for a sensitive method RNLA, the RMSE of other algorithms is increasing because the adverse impact of the absorption bias would rise synchronously. In light of the results, IPM-MM has a relatively low error among all the considered methods. It can be found that when Sel equals 110 m, its performance is similar to SRWLS, but is better than other methods with side range is less than 110 m, and interestingly, when the range is up to 120 m, it recovers its good performance. It is worth noting that when Sel = 80 m, the IPM-MM performance of estimating 0 P is slightly close to IPM. However, IPM-MM acquires great improvements on the basis of IPM under other conditions. For ensuring location accuracy and good power estimation, our IPM-MM is effective and recommendable to some extent.

Scenario with Varying Side Length Sel
Due to the fact that the absorption bias is determined by the absorption coefficient α f and the distance between each anchor and the target node, it is essential to carry out the simulation experiment. The related parameters are shown in Table 5. According to Figure 8, when the side length ranges from 70 m to 120 m, except for a sensitive method RNLA, the RMSE of other algorithms is increasing because the adverse impact of the absorption bias would rise synchronously. In light of the results, IPM-MM has a relatively low error among all the considered methods. It can be found that when Sel equals 110 m, its performance is similar to SRWLS, but is better than other methods with side range is less than 110 m, and interestingly, when the range is up to 120 m, it recovers its good performance. It is worth noting that when Sel = 80 m, the IPM-MM performance of estimating P 0 is slightly close to IPM. However, IPM-MM acquires great improvements on the basis of IPM under other conditions. For ensuring location accuracy and good power estimation, our IPM-MM is effective and recommendable to some extent.

Scenario with Varying Side Length Sel
Due to the fact that the absorption bias is determined by the absorption coefficient f α and the distance between each anchor and the target node, it is essential to carry out the simulation experiment. The related parameters are shown in Table 5. According to Figure 8, when the side length ranges from 70 m to 120 m, except for a sensitive method RNLA, the RMSE of other algorithms is increasing because the adverse impact of the absorption bias would rise synchronously. In light of the results, IPM-MM has a relatively low error among all the considered methods. It can be found that when Sel equals 110 m, its performance is similar to SRWLS, but is better than other methods with side range is less than 110 m, and interestingly, when the range is up to 120 m, it recovers its good performance. It is worth noting that when Sel = 80 m, the IPM-MM performance of estimating 0 P is slightly close to IPM. However, IPM-MM acquires great improvements on the basis of IPM under other conditions. For ensuring location accuracy and good power estimation, our IPM-MM is effective and recommendable to some extent.   Figure 9 indicates the cumulative distribution function (CDF) of various methods versus x x − when f α = 0.04 dB/m. The other parameter settings are shown in Table 6.

Cumulative Distribution Function (CDF)
It is obvious that our proposed method obtains the best performance. When CDF reaches almost 90%, the error x x − of our IPM-MM is 5.20 m, whereas the errors of SRWLS,  4.6. Cumulative Distribution Function (CDF) Figure 9 indicates the cumulative distribution function (CDF) of various methods versus x − x when α f = 0.04 dB/m. The other parameter settings are shown in Table 6. It is obvious that our proposed method obtains the best performance. When CDF reaches almost 90%, the error x − x of our IPM-MM is 5.20 m, whereas the errors of SRWLS, ASM, LS, IPM, RNLA are 7.05, 8.70, 9.25, 10.94 and 12.01 m, respectively. According to the numerical results, we find that IPM combined with the majorization-minimization strategy brings better improvement than the single IPM and realize 52.5% performance improvement.  Figure 9 indicates the cumulative distribution function (CDF) of various methods versus x x − when f α = 0.04 dB/m. The other parameter settings are shown in Table 6.

Cumulative Distribution Function (CDF)
It is obvious that our proposed method obtains the best performance. When CDF reaches almost 90%, the error x x − of our IPM-MM is 5.20 m, whereas the errors of SRWLS, ASM, LS, IPM, RNLA are 7.05, 8.70, 9.25, 10.94 and 12.01 m, respectively. According to the numerical results, we find that IPM combined with the majorization-minimization strategy brings better improvement than the single IPM and realize 52.5% performance improvement.     Figure 10 depicts the RMSE of various methods with growing σ 2 . The noise parameter range is set from 1 to 10. The relevant parameters are displayed in Table 7. Intuitively, with the increase of σ 2 , the overall performance shows a downward trend. Compared with LS, SRWLS, ASM, RNLA, our IPM-MM has a lower RMSE especially estimating x. Although in Figure 10b, our proposed method has no obvious advantages compared with SRWLS and RNLA, our scheme is more accurate for the position estimation. In addition, IPM-MM always presents a better performance for the estimate of x and P 0 than IPM. Besides, we find that the performance gap between SRWLS and RNLA is small; it occurs because both two converted the location problem into a GTRS. Generally, to ensure good performance in a noisy environment, our method provides the best choice of all.

Conclusions
In this paper, a novel localization method called IPM-MM is proposed in OSNs using RSS measurements, which is applicable for the unknown transmission power case. The objective function of the original localization problem is transformed into the UT-ANLS framework, and the IPM together with the MMT is exploited to find an optimum solution. The simulation results prove that our proposed IPM-MM method can be superior to other state-of-the-art approaches in terms of localization accuracy in a relatively high dynamics environment.
Author Contributions: Y.Z. proposed the main ideas, wrote the paper, designed the whole framework, and carried out the simulations. H.W. provided guidance for the work, reviewed the paper, and acquired funding. X.M. provided guidance for the work, reviewed the paper, and collaborated in the discussion on the proposed system framework and location scheme. J.X., W.W., Q.Z., L.L. assisted in testing the code and checked the paper. All authors have read and agreed to the published version of the manuscript.