Next Article in Journal
Parameter Estimation of a Ground Moving Target Using Image Sharpness Optimization
Previous Article in Journal
Classification of Kiwifruit Grades Based on Fruit Shape Using a Single Camera
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Joint Parameter Estimation for the Two-Wave with Diffuse Power Fading Model

by
Jesus Lopez-Fernandez
*,
Laureano Moreno-Pozas
,
Francisco Javier Lopez-Martinez
and
Eduardo Martos-Naya
Departamento de Ingeniería de Comunicaciones, ETS Ingeniería de Telecomunicación, Universidad de Málaga, Málaga 29071, Spain
*
Author to whom correspondence should be addressed.
Sensors 2016, 16(7), 1014; https://doi.org/10.3390/s16071014
Submission received: 14 April 2016 / Revised: 23 June 2016 / Accepted: 25 June 2016 / Published: 30 June 2016
(This article belongs to the Section Sensor Networks)

Abstract

:
Wireless sensor networks deployed within metallic cavities are known to suffer from a very severe fading, even in strong line-of-sight propagation conditions. This behavior is well-captured by the Two-Wave with Diffuse Power (TWDP) fading distribution, which shows great fit to field measurements in such scenarios. In this paper, we address the joint estimation of the parameters K and Δ that characterize the TWDP fading model, based on the observation of the received signal envelope. We use a moment-based approach to derive closed-form expressions for the estimators of K and Δ, as well as closed-form expressions for their asymptotic variance. Results show that the estimation error is close to the Cramer-Rao lower bound for a wide range of values of the parameters K and Δ. The performance degradation due to a finite number of observations is also analyzed.

1. Introduction

Rician fading model is extensively used to characterize the random fluctuations of the received signal amplitude in line-of-sight (LOS) environments [1]. In such scenarios, the scattering waves arriving at the receiver can be split into the dominant LOS component plus a diffuse component, which accounts for the effect of the non-LOS (NLOS) propagation. The relative strength of the LOS component with respect to the NLOS one is measured by the Rician K-factor, defined as the ratio between the powers of both components.
The estimation of the Rician K-factor is of paramount importance in the context of wireless communications systems, as the proper system operation heavily relies on the quality of the estimation of K [2,3]. As a matter of fact, a flurry of papers have addressed this estimation problem from different perspectives [4,5,6,7,8] in the last years.
However, in some indoor and outdoor LOS environments the Rician distribution falls short to accurately modeling small-scale fading effects. This is of special relevance in the case of wireless sensor networks (WSN) deployed on the inner surface of cavities (e.g., tunnels, plane and helicopter airframes, buses, shipping containers) used to measure data for maintenance, comfort, health and security purposes [9,10,11] or in vehicle-to-vehicle links [12], or in general, machine-to-machine systems like WSN deployed inside and/or around complex structures. In all the aforementioned scenarios, field measurements show that fading in the presence of LOS components may even degenerate into a regime more severe than Rayleigh fading [9,10,11,12].
Among the different distributions in the literature which may be useful to account for this special propagation regime, the Two-Wave with diffuse power (TWDP) fading model is the preferred alternative because of its flexibility [11] and clear physical interpretation. This fading model was first presented in [13] as a generalization of the Rician fading model, by incorporating a second dominant LOS component with uniformly distributed phase. The effect of this new LOS component is captured through the parameter Δ, which measures the relative magnitudes of the two LOS components to one another. This model has succeeded on characterizing a wide variety of fading behaviors, from purely Rician to worse than Rayleigh fading [14,15], and has been considered by many authors [16,17,18,19] in order to evaluate the performance of communication systems operating under this peculiar fading condition.
Despite its relevance, the estimation of the parameters K and Δ for the TWDP fading model has not been addressed in the literature to the best of our knowledge, except for the preliminary results in [20] based on maximum-likelihood estimation. Thence, we aim to fill this gap by studying the problem of jointly estimating K and Δ from the observation of the received signal envelope. Specifically, we design a moment-based estimator for the TWDP fading parameters leveraging the closed-form expressions of the moments recently proposed in [15]. We also derive the asymptotic variance of the estimator, and compare it to the Cramer-Rao lower bound (CRB) via Monte Carlo simulations. The effect of a finite number of observations on the estimator performance is also quantified.
The remainder of this paper is structured as follows: Section 2 introduces the main aspects of the TWDP fading model, which will be of use in the following derivations. The parameter estimation is addressed in Section 3, whereas the CRB and the asymptotic variance of the proposed estimators are presented in Section 4. The effect of having a finite number of observations in the estimation performance is studied in Section 5, and the main conclusions are outlined in Section 6.

2. TWDP Fading Model

In the TWDP fading model, the complex baseband signal V r at the receiver side can be expressed as ([13] Equation (4))
V r = V 1 exp ( j ϕ 1 ) + V 2 exp ( j ϕ 2 ) + V d
We observe two LOS components with uniformly distributed phases ϕ 1 , ϕ 2 U ( 0 , 2 π ) (The symbol ∼ reads as statistically distributed as.) and constant amplitudes V 1 and V 2 , plus a diffuse component V d = X + j Y regarded as a complex Gaussian random variable, being X , Y N ( 0 , σ 2 ) and independent. The TWDP model is fully characterized by the parameters K, Δ and Ω defined as
K = V 1 2 + V 2 2 2 σ 2 , Δ = 2 V 1 V 2 V 1 2 + V 2 2 , Ω = V 1 2 + V 2 2 + 2 σ 2
The parameter K [ 0 , ) has an equivalent interpretation as in the Rician case, whereas the parameter Δ [ 0 , 1 ] indicates whether the two LOS components have equal amplitude ( Δ = 1 ) or not ( Δ < 1 ), degenerating in the Rician fading model for Δ = 0 . The pdf of r = | V r | is given by ([13] Equation (32))
f ( r ) = r 2 π σ 2 e - r 2 2 σ 2 - K 0 2 π e K Δ cos ( θ ) I 0 r σ c θ d θ
where I 0 ( · ) is the modified Bessel function of the first kind and order zero, and c θ = 2 K [ 1 + Δ cos ( θ ) ] . Finally, the parameter Ω corresponds to the average received signal power and hence is merely a scale factor for the model. Although we focus our study exclusively in the estimation of the K and Δ parameters, the parameter Ω is also unknown and related to K and Δ and must therefore be estimated as well.
In ([15] Equation (30)), the moments of the squared signal envelope γ = r 2 were calculated as a finite-range integral involving the Kummer confluent hypergeometric function 1 F 1 ( · , · ; · ) . Using the connection between the Kummer function and the Laguerre polynomials L k ( · ) [21] [Equation 8.972], we have
E ( γ k ) = k ! Ω k ( 1 + K ) k 2 π 0 2 π L k - K ( 1 + Δ cos θ ) d θ
Noting that the Laguerre polynomials can be expressed as [21] [Equation 8.970]
L k ( - x ) = m = 0 k k m x m m !
the moments can be computed in closed-form by direct integration of Equation (4) as
E ( γ k ) = k ! Ω k ( 1 + K ) k 2 π m = 0 k k m K m m ! 0 2 π ( 1 + Δ cos θ ) m d θ
by using a simple change of variables t = cos θ . Besides, the first moments of γ coincide with the first even moments of r 2 , and can be computed from Equation (6) as
E ( γ ) = μ 2 = Ω
E ( γ 2 ) = μ 4 = Ω 2 ( 4 + K ( 8 + K ( 2 + Δ 2 ) ) ) 2 ( 1 + K ) 2
E ( γ 3 ) = μ 6 = Ω 3 ( 12 + K ( 3 + K ) ( 12 + K ( 2 + 3 Δ 2 ) ) ) 2 ( 1 + K ) 3
where μ n is the n-th order moment of the received signal envelope, i.e., μ n = E [ r n ] . The second, fourth and sixth-order moments μ 2 , μ 4 , μ 6 are expressed in closed-form in terms of the parameters K and Δ and Ω. We will use these expressions to design a moment-based joint estimator.

3. Moment-Based Estimation of K and Δ

Dividing Equation (8) by μ 2 2 and Equation (9) by μ 2 3 , we obtain
μ 4 μ 2 2 = 4 + K ( 8 + K ( 2 + Δ 2 ) ) 2 ( 1 + K ) 2
μ 6 μ 2 3 = 12 + K ( 3 + K ) ( 12 + K ( 2 + 3 Δ 2 ) ) 2 ( 1 + K ) 3
Using the sample moments μ ^ n = 1 N l = 1 N r l n instead of the ensemble averages μ n , Equations (10) and (11) can be solved for K and Δ resulting in moment-based estimators (A similar procedure is used in [6] for designing moment-based estimators of the Rician K factor where a function to denote the ratio of different moments is introduced and denoted as f n , m . The terms μ 4 μ 2 2 and μ 6 μ 2 3 on the left hand side of Equations (10) and (11) would correspond to f 4 , 2 and f 6 , 2 respectively.). In this procedure notice that the third parameter Ω is also being estimated when using μ ^ 2 instead of μ 2 as indicated in Equation (7). Solving Equation (10) for Δ and substituting it in Equation (11) results in the cubic polynomial equation on K ^
K ^ 3 + a 1 K ^ 2 + a 2 K ^ + a 3 = 0
with coefficients
a 1 = 6 μ ^ 6 - 30 μ ^ 4 μ ^ 2 + 24 μ ^ 2 3 2 μ ^ 6 - 6 μ ^ 4 μ ^ 2 + 4 μ ^ 2 3 a 2 = 6 μ ^ 6 - 42 μ ^ 4 μ ^ 2 + 48 μ ^ 2 3 2 μ ^ 6 - 6 μ ^ 4 μ ^ 2 + 4 μ ^ 2 3 a 3 = 2 μ ^ 6 - 18 μ ^ 4 μ ^ 2 + 24 μ ^ 2 3 2 μ ^ 6 - 6 μ ^ 4 μ ^ 2 + 4 μ ^ 2 3
The closed-form solution of a cubic polynomial results in the estimator
K ^ = p + p 2 + q 3 1 2 1 3 + p - p 2 + q 3 1 2 1 3 - a 1 3
with
p = 1 54 ( 9 a 1 a 2 - 27 a 3 - 2 a 1 3 ) , q = 1 9 ( 3 a 2 - a 1 2 )
Note that, while Equation (12) has three possible solutions, only one of them yields a valid estimation of K (i.e., real and positive); this is the solution stated in Equation (14). Plugging the estimate of K in Equation (10) yields a quadratic polynomial in Δ, for which the following estimator is obtained:
Δ ^ = 2 μ ^ 4 μ ^ 2 2 ( 1 + K ^ ) 2 - 2 K ^ 2 - 8 K ^ - 4 K ^
where the positive-valued solution has been selected, as Δ 0 by definition.

4. Asymptotic Variance and Cramer-Rao Bound

To assess the performance of the proposed estimators, we will first calculate the asymptotic variance. Then, we will compare it to the CRB bound, which determines the minimum achievable variance of any unbiased estimator.
It is well-known that for any moment-based estimator, as the number of independent and identically distributed (i.i.d.) observations N increases, the estimator bias tends asymptotically to zero. Similarly, in this situation the estimator variance tends to vary proportionally to 1 / N (it is a N -consistent estimator). In our problem, the proposed estimators are a function of three sample moments, i.e., K ^ = h 1 ( μ ^ 2 , μ ^ 4 , μ ^ 6 ) and Δ ^ = h 2 ( μ ^ 2 , μ ^ 4 , μ ^ 6 ) where the function h 1 ( · ) is obtained by substituting Equations (13) and (15) in Equation (14) and the function h 2 ( · ) is obtained by substituting Equation (14) in Equation (16). The asymptotic variance for the estimation of K is given by [22] [Equation 9.16]:
AsV K = g K Σ g K T
where g K is the vector of derivatives of K with respect to each sample moment, evaluated at the corresponding statistical moment value, i.e.,
g K = K ^ μ ^ 2 , K ^ μ ^ 4 , K ^ μ ^ 6 μ ^ 2 = μ 2 , μ ^ 4 = μ 4 , μ ^ 6 = μ 6
and Σ is the 3 × 3 covariance matrix of the ensemble averages whose entries are defined as [ Σ ] i j = cov ( μ ^ 2 i , μ ^ 2 j ) , i , j = 1 , 2 , 3 . It can be easily shown that
[ Σ ] i j = 1 N ( μ 2 i + 2 j - μ 2 i μ 2 j ) , i , j = 1 , 2 , 3
Evaluation of Equation (19) requires the use of all the even order moments up to order 12. The first three moments are shown in Equations (7)–(9). Closed-form expressions for μ 8 , μ 10 and μ 12 can be easily obtained by evaluating Equation (6), and are shown in Appendix A. Replacing K by Δ in Equations (17) and (18) yields the corresponding expression for A s V Δ .
The Cramer-Rao lower bound is directly related to the Fisher Information Matrix (FIM), which will be denoted by I ( θ ) with entries [ I ( θ ) ] i j where θ = [ K , Δ , Ω ] is the parameter vector to be estimated. For i.i.d. observations the FIM matrix is a 3 × 3 matrix with
[ I ( θ ) ] i j = N E ln f ( r ) θ i · ln f ( r ) θ j i , j = 1 , 2 , 3
where N is the number of observations and f ( r ) is the TWDP pdf in Equation (3).
The CRBs for the estimation of K and Δ will be denoted here as CRB K and CRB Δ respectively. According to [22] , these can be expressed as the elements (1,1) and (2,2) of the inverse FIM matrix, respectively, i.e., CRB K = [ I ( θ ) - 1 ] 11 and CRB Δ = [ I ( θ ) - 1 ] 22 . See that the CRBs decrease with the number of observations as 1 / N . In order to compute the FIM entries we must first obtain the derivatives in Equation (20), and then perform numerical integration. As stated before, although we are interested in the estimation of K and Δ, the parameter Ω must likewise be estimated and accounted for in the CRB calculation. Whether Ω is known or not indeed affects the CRB, but the specific value of Ω is irrelevant because the CRB is scale invariant.
Before presenting the results for CRB and AsV, we would like to remark an important aspect relative to normalization. In the context of moment-based estimators it is common (see Figure 2 in [6]) to normalize the AsV and the CRB by N because the (asymptotic) behavior of the proposed estimator when N is high can be seen in one single figure and compared to the CRB. Therefore in the sequel, results for both the CRB and the AsV will be normalized with respect to N. A further normalization with respect to the estimated parameter true value will be carried out in order to have bounds on relative errors (instead of absolute errors) which give better information of the goodness of the estimation. Finally, it is usual in the literature to represent the bound on the standard deviation of the error instead of on the variance. Hence, the square root of CRB and AsV will be considered. For the sake of compactness, in the text we will refer to the square root of the CRB for the estimation of K, normalized to both K and N as sqrt-normalized CRB K which stands for CRB K N K 2 and the same holds for the sqrt-normalized CRB Δ , AsV K and AsV Δ .
In Figure 1, the sqrt-normalized CRB K and the sqrt-normalized AsV K are plotted as a function of K. Different plots correspond to different values of the unknown parameter Δ. Focusing first on the CRB, see that the sqrt-normalized CRB K grows as Δ decreases. This means that the error bound in estimating the ratio of the LOS components power to the diffuse components power is higher when there is only one LOS component. The estimation can improve as the two LOS components become of similar magnitude. See also that the error bound decreases drastically as K increases from 0 to a moderate value (around K = 2 for Δ < 0 . 5 ); this is, as the power of the diffuse components diminishes. Finally, for larger K the error converges to a constant value as K goes to infinity. The value of the sqrt-normalized AsV K is remarkably close to the sqrt-normalized CRB K for the entire considered range of parameters. This means that the proposed estimator of the K parameter is almost asymptotically efficient.
Proceeding in an analogous manner, in Figure 2, we have represented the sqrt-normalized CRB Δ and the sqrt-normalized AsV Δ as a function of Δ. Different plots correspond to different values of the unknown parameter K. With regard to the sqrt-normalized CRB Δ , we notice once again that the estimation improves when K grows large. We also observe that for a given K, the sqrt-normalized CRB Δ increases as Δ 0 ; this implies that it is difficult to determine the relative amplitudes of the two LOS components as one of them vanishes. The converse can be argued when both tend to be equal (i.e., Δ 1 ). We also see that the values of the sqrt-normalized AsV Δ are relatively close to the sqrt-normalized CRB Δ for the range of parameters analyzed, but not as close as when estimating K.

5. Effect of a Finite Number of Observations

The performance of the proposed estimators for finite N is studied by resorting to Monte-Carlo simulations. For every fixed pair of K and Δ, we have generated 500 sets of N i.i.d. realizations of the TWDP random variable. As an initial sanity check, Figure 3 shows the estimated K using Equation (14) vs. the true K value for a fixed value of Δ. The dashed line on the graphic results from averaging the estimate over the trials while the solid lined corresponds to the unit slope line that serves as a reference. From this type of plot we can extract qualitative results about the behavior of the estimators. As can be seen in Figure 3, the estimator K ^ shows a reduced bias for the range of K values shown in the figure. See also that the dispersion of the estimated values grows with growing K. Accordingly, Figure 4 shows the estimated Δ using Equation (16) vs. the true Δ value for a fixed value of K. In this case the estimator shows a growing bias and dispersion of values as Δ decreases.
In order to perform a quantitative analysis of the estimators performance we have computed the sample mean squared errors MSE K and MSE Δ where MSE K = 1 500 i = 1 500 ϵ 2 with ϵ = K ^ - K and where MSE Δ is defined in an analogous way. The double normalization used with the CRB and AsV is also applied to the MSE, and hence term sqrt-normalized MSE K denotes the square root of the MSE in the estimation of K normalized to both K and N and corresponds to MSE K N K 2 . The same holds for sqrt-normalized MSE Δ .
In Figure 5, the sqrt-normalized MSE K using Equation (14) is plotted as a function of K for Δ = 0 . 5 and different values of N. Accordingly, in Figure 6, the sqrt-normalized MSE Δ using Equation (16) is plotted as a function of Δ for K = 3 and different values of N. In both figures the corresponding sqrt-normalized CRB and AsV have been included as a reference.
See that for a very high number of observations ( N = 10 6 ) the sqrt-normalized MSE is close to the sqrt-normalized AsV in both cases, as expected. As N decreases, the sqrt-normalized MSE falls below the sqrt-normalized CRB for values of K 4 in Figure 5 and for values of Δ 0 . 6 in Figure 6. The reason for this is that the bias of the proposed estimators grows when decreasing N. It is well-known that the CRB sets a limit on the variance of unbiased estimators. However, in case the estimators are biased the CRB is not practical and the MSE (which accounts for both variance and bias) may take a lower value. This behavior is also observed in the moment-based Rice parameter estimators proposed in the literature for values of the Rician parameter K close to zero and finite N (e.g., see Figure 5 in [6]). Although a CRB and a MSE bound for biased estimators can be also resorted [23], they are impractical since they require an a priori choice of the bias gradient. Although not explicitly included in Figure 5 and Figure 6 a similar behavior is observed for other values of K and Δ. Due to the double normalization applied to the sqrt-normalized MSE, care must be taken when interpreting the relative values of the different sqrt-normalized MSE curves shown in Figure 5 and Figure 6. For instance, in Figure 6, the sqrt-normalized MSE curve corresponding to N = 5 × 10 3 takes a higher value than that of N = 10 3 in most of the range where Δ varies. De-normalizing the sqrt-normalized MSE with respect to N reveals that the sqrt-normalized MSE with N = 5 × 10 3 is lower than that with N = 10 3 for the whole range of Δ; i.e., the estimation error is in fact lower for higher number of observations, as expected.
This fact is explicitly shown in Figure 7 and Figure 8 where in this case both the CRB and MSE are not normalized by N. See how both the CRB and the MSE decrease as N grows from 500 to 10 4 for all the values of K and Δ in Figure 7 and Figure 8, respectively.

6. Conclusions

In this work, the joint estimation of the two parameters of the TWDP fading model has been addressed. A simple but accurate closed-form moment-based estimator for both parameters was used for the estimation. Results show that the proposed estimators operate relatively close to the CRB for a wide range of parameter values. However, the estimator of K is comparatively better than the estimator of Δ, as the AsV is closer to the CRB. The results here presented can set the underpinnings for alternative solutions for the parameter estimation of the TWDP fading model.

Acknowledgments

This work has been funded by the Consejería de Economía, Innovación, Ciencia y Empleo of the Junta de Andalucía, the Spanish Government and the European Regional Development Fund (projects P2011-TIC-7109, P2011-TIC-8238, TEC2013-44442-P and TEC2014-57901-R).

Author Contributions

The estimation method was conceived by Francisco Javier Lopez-Martinez and Eduardo Martos-Naya. The derivation of the analytical expressions and the computer simulations were carried out by Jesus Lopez-Fernandez. Results were supervised by Laureano Moreno-Pozas. Figures were elaborated by Laureano Moreno-Pozas. The paper was written by Jesus Lopez-Fernandez and Francisco Javier Lopez-Martinez. Revision and proofreading of the final manuscript were conveyed by all authors.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A: Expressions for the Higher Order Moments μ8, μ10 and μ12

Computation of the covariance matrix Σ in Equation (19) requires the use of the higher order moments μ 8 , μ 10 and μ 12 . By direct integration of Equation (6), the following expressions for the higher order moments of the signal envelope can be obtained as
μ 8 = Ω 4 8 ( 1 + K ) 4 ( 192 + K ( 8 ( 4 + K ) ( 24 + K ( 12 + K ) ) + 24 K ( 2 + K ) ( 6 + K ) Δ 2 + 3 K 3 Δ 4 ) )
μ 10 = Ω 5 8 ( 1 + K ) 5 ( 960 + K ( 8 ( 20 + K ( 10 + K ) ) ( 30 + K ( 15 + K ) ) + 40 K ( 60 + K ( 60 + K ( 15 + K ) ) ) Δ 2 + 15 K 3 ( 5 + K ) Δ 4 ) )
μ 12 = Ω 6 16 ( 1 + K ) 6 ( 16 ( 720 + K ( 6 + K ) ( 720 + K ( 780 + K ( 270 + K ( 30 + K ) ) ) ) ) + 120 K 2 ( 360 + K ( 480 + K ( 180 + K ( 24 + K ) ) ) ) Δ 2 + 90 K 4 ( 30 + K ( 12 + K ) ) Δ 4 + 5 K 6 Δ 6 )

References

  1. Simon, M.K.; Alouini, M.-S. Digital Communication over Fading Channels; Wiley-IEEE Press: New York, NY, USA, 2005. [Google Scholar]
  2. Catreux, S.; Erceg, V.; Gesbert, D.; Heath, R.W. Adaptive modulation and MIMO coding for broadband wireless data networks. IEEE Commun. Mag. 2002, 40, 108–115. [Google Scholar] [CrossRef]
  3. Jayaweera, S.K.; Poor, H.V. On the capacity of multiple-antenna systems in Rician fading. IEEE Trans. Wirel. Commun. 2005, 4, 1102–1111. [Google Scholar] [CrossRef]
  4. Greenstein, L.; Michelson, D.; Erceg, V. Moment-method estimation of the Ricean K-factor. IEEE Commun. Lett. 1999, 3, 175–176. [Google Scholar] [CrossRef]
  5. Abdi, A.; Tepedelenlioglu, C.; Kaveh, M.; Giannakis, G. On the estimation of the K parameter for the Rice fading distribution. IEEE Commun. Lett. 2001, 5, 92–94. [Google Scholar] [CrossRef]
  6. Tepedelenlioglu, C.; Abdi, A.; Giannakis, G. The Ricean K factor: Estimation and performance analysis. IEEE Trans. Wirel. Commun. 2003, 2, 799–810. [Google Scholar] [CrossRef]
  7. Bhattacharjea, R.; Durgin, G.D.; Anderson, C.R. Estimation of Rician K-Factors from Block-Averaged Channel Measurements. IEEE Trans. Wirel. Commun. 2012, 11, 4231–4236. [Google Scholar] [CrossRef]
  8. Medawar, S.; Handel, P.; Zetterberg, P. Approximate Maximum Likelihood Estimation of Rician K-Factor and Investigation of Urban Wireless Measurements. IEEE Trans. Wirel. Commun. 2013, 12, 2545–2555. [Google Scholar] [CrossRef]
  9. Frolik, J.; Sipal, V.; Edwards, D.J. Leveraging depolarization to increase the link reliability for wireless sensors operating in hyper-rayleigh environments. IEEE Sensors J. 2014, 14, 2442–2446. [Google Scholar] [CrossRef]
  10. Frolik, J. A case for considering hyper-Rayleigh fading channels. IEEE Trans. Wirel. Commun. 2007, 6, 1235–1239. [Google Scholar] [CrossRef]
  11. Ketcham, R.; Frolik, J.; Zivanovic, B.; Melais, S.; Weller, T. Compact and simple diversity methods for mitigating severe fading. In Proceedings of the IEEE Annual Wireless and Microwave Technology Conference, 2006, WAMICON ’06, Clearwater Beach, FL, USA, 4–5 December 2006; pp. 1–5.
  12. Matolak, D.; Frolik, J. Worse-than-Rayleigh fading: Experimental results and theoretical models. IEEE Commun. Mag. 2011, 49, 140–146. [Google Scholar] [CrossRef]
  13. Durgin, G.; Rappaport, T.; de Wolf, D.A. New analytical models and probability density functions for fading in wireless communications. IEEE Trans. Commun. 2002, 50, 1005–1015. [Google Scholar] [CrossRef]
  14. Wang, L.; Yang, N.; Elkashlan, M.; Yeoh, P.L.; Yuan, J. Physical Layer Security of Maximal Ratio Combining in Two-Wave With Diffuse Power Fading Channels. IEEE Trans. Inf. Forensics Sec. 2014, 9, 247–258. [Google Scholar] [CrossRef]
  15. Rao, M.; Lopez-Martinez, F.J.; Alouini, M.-S.; Goldsmith, A. MGF Approach to the Analysis of Generalized Two-Ray Fading Models. IEEE Trans. Wirel. Commun. 2015, 14, 2548–2561. [Google Scholar]
  16. Saberali, S.; Beaulieu, N. New expressions for TWDP fading statistics. IEEE Wirel. Commun. Lett. 2013, 2, 643–646. [Google Scholar] [CrossRef]
  17. Lu, Y.; Yang, N. Symbol error rate of decode-and-forward relaying in two-wave with diffuse power fading channels. IEEE Trans. Wirel. Commun. 2012, 11, 3412–3417. [Google Scholar] [CrossRef]
  18. Lu, Y.; Wang, X.; Yang, N. Outage probability of cooperative relay networks in two-wave with diffuse power fading channels. IEEE Trans. Commun. 2012, 60, 42–47. [Google Scholar] [CrossRef]
  19. Dixit, D.; Sahu, P. Performance of QAM signaling over TWDP fading channels. IEEE Trans. Wirel. Commun. 2013, 12, 1794–1799. [Google Scholar] [CrossRef]
  20. Henderson, A.H.; Durkin, C.J.; Durgin, G.D. Finding the right small-scale fading distribution for a measured indoor 2.4 GHz channel. In Proceedings of the 2008 IEEE Antennas and Propagation Society International Symposium, San Diego, CA, USA, 5–11 July 2008; pp. 1–4.
  21. Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series and Products, 7th ed.; Academic Press Inc.: New York, NY, USA, 2007. [Google Scholar]
  22. Kay, S.M. Fundamentals of Statistical Signal Processing. (Volume I), Estimation Theory; Prentice Hall: Upper Saddle River, NJ, USA, 1993. [Google Scholar]
  23. Eldar, Y. Mse bounds with affine bias dominating the cramer-rao bound. IEEE Trans. Signal Proc. 2008, 56, 3824–3836. [Google Scholar] [CrossRef]
Figure 1. Sqrt-normalized CRB K and AsV K as a function of K, for different values of Δ.
Figure 1. Sqrt-normalized CRB K and AsV K as a function of K, for different values of Δ.
Sensors 16 01014 g001
Figure 2. Sqrt-normalized CRB Δ and AsV Δ as a function of Δ, for different values of K.
Figure 2. Sqrt-normalized CRB Δ and AsV Δ as a function of Δ, for different values of K.
Sensors 16 01014 g002
Figure 3. Samples of K ^ using Equation (14) vs. true K for Δ = 0 . 5 .
Figure 3. Samples of K ^ using Equation (14) vs. true K for Δ = 0 . 5 .
Sensors 16 01014 g003
Figure 4. Samples of Δ ^ Equation (16) vs. true Δ for K = 10 .
Figure 4. Samples of Δ ^ Equation (16) vs. true Δ for K = 10 .
Sensors 16 01014 g004
Figure 5. Sqrt-normalized MSE K for the estimation of K using Equation (14) for Δ = 0 . 5 and different number of observations N. Also shown are the corresponding sqrt-normalized AsV K and sqrt-normalized CRB K .
Figure 5. Sqrt-normalized MSE K for the estimation of K using Equation (14) for Δ = 0 . 5 and different number of observations N. Also shown are the corresponding sqrt-normalized AsV K and sqrt-normalized CRB K .
Sensors 16 01014 g005
Figure 6. Sqrt-normalized MSE Δ for the estimation of Δ using Equation (16) for K = 3 and different number of observations N. Also shown are the corresponding sqrt-normalized AsV Δ and sqrt-normalized CRB Δ .
Figure 6. Sqrt-normalized MSE Δ for the estimation of Δ using Equation (16) for K = 3 and different number of observations N. Also shown are the corresponding sqrt-normalized AsV Δ and sqrt-normalized CRB Δ .
Sensors 16 01014 g006
Figure 7. Square root of CRB K and MSE K normalized by K for the estimation of K using Equation (14) for Δ = 0 . 5 and different number of observations N.
Figure 7. Square root of CRB K and MSE K normalized by K for the estimation of K using Equation (14) for Δ = 0 . 5 and different number of observations N.
Sensors 16 01014 g007
Figure 8. Square root of CRB Δ and MSE Δ normalized by Δ for the estimation of Δ using Equation (16) for K = 3 and different number of observations N.
Figure 8. Square root of CRB Δ and MSE Δ normalized by Δ for the estimation of Δ using Equation (16) for K = 3 and different number of observations N.
Sensors 16 01014 g008

Share and Cite

MDPI and ACS Style

Lopez-Fernandez, J.; Moreno-Pozas, L.; Lopez-Martinez, F.J.; Martos-Naya, E. Joint Parameter Estimation for the Two-Wave with Diffuse Power Fading Model. Sensors 2016, 16, 1014. https://doi.org/10.3390/s16071014

AMA Style

Lopez-Fernandez J, Moreno-Pozas L, Lopez-Martinez FJ, Martos-Naya E. Joint Parameter Estimation for the Two-Wave with Diffuse Power Fading Model. Sensors. 2016; 16(7):1014. https://doi.org/10.3390/s16071014

Chicago/Turabian Style

Lopez-Fernandez, Jesus, Laureano Moreno-Pozas, Francisco Javier Lopez-Martinez, and Eduardo Martos-Naya. 2016. "Joint Parameter Estimation for the Two-Wave with Diffuse Power Fading Model" Sensors 16, no. 7: 1014. https://doi.org/10.3390/s16071014

APA Style

Lopez-Fernandez, J., Moreno-Pozas, L., Lopez-Martinez, F. J., & Martos-Naya, E. (2016). Joint Parameter Estimation for the Two-Wave with Diffuse Power Fading Model. Sensors, 16(7), 1014. https://doi.org/10.3390/s16071014

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