1. Introduction
Small-scale fading severely degrades performance of wireless communication systems [
1]. Thus, in order to design highly reliable and efficient transceivers for communication in deep fading conditions, it is of profound significance to have an accurate and tractable fading model [
2,
3]. Traditionally, the signal affected by small-scale fading in nonline-of-sight (NLOS) environments has been modeled as a sum of many diffuse components and described by Rayleigh distribution. Withal, the Rician fading model has been used for the description of received signal variations in better-than-Rayleigh channel conditions, since, except for many diffuse components, it also assumes the presence of one dominant specular component [
4]. However, in some multipath sparse channels, traditional fading models fall short of accurately describing small-scale variations of the received complex envelope [
4], which makes them nonviable for accurate modeling of all propagation conditions [
5].
The aforementioned is especially pronounced in mmWave 5G communication networks equipped with directional antennas or arrays [
6] as well as in wireless sensor networks deployed in cavity environments [
3], where measurements have shown that the signal in multipath sparse propagation conditions may, in some cases, experience fading with characteristics worse than those in Rayleigh channels [
4].
In such scenarios, the TWDP fading model—which assumes the presence of two specular components instead of just one considered by Rician model—appears to be appropriate for modeling small-scale fading effects [
7]. In this regard, it is shown that TWDP distribution describes small-scale fading more accurately than conventional distributions in static sensor networks with their nodes placed within cavity environments, such as aircrafts and buses [
3] and abroad a large transport helicopter [
8]. It is also shown that its PDF shape significantly coincides with empirical results obtained for train-to-infrastructure wireless communication [
9], vehicle-to-vehicle 60 GHz urban communication [
10], and mmWave communication in indoor environments [
11,
12].
Thereat, the TWDP model itself represents generalization of Rayleigh and Rician fading models and can be used for modeling both better-than-Rayleigh and worse-than-Rayleigh fading conditions. It is based on the assumption that the waves arriving at the receiver can be observed as a sum of two specular components with constant magnitudes
and
and uniformly distributed phases, plus many diffuse components treated as a complex zero-mean Gaussian random process with average power
. As such, the model is conventionally parameterized by
K and
, defined as [
13]:
where parameter
K (
) characterizes the ratio of the average power of specular components to the average power of remaining diffuse components (like the Rician parameter
) and parameter
characterizes the relation between magnitudes of specular components. Consequently, the average received signal power
is equal to
.
However, as it is elaborated in [
14], the definition of parameter
is not in accordance with the model’s underlying physical mechanisms. Namely, according to the model’s definition, specular components have constant magnitudes and are propagating in a linear medium. Consequently,
can be nothing but the linear combination of
, wherefore the function which characterizes the ratio between
and
has to be linear [
14]. However,
-based parameterization introduces a nonlinear relation between magnitudes of specular components, since
. This hinders accurate observation of the impact of the ratio between
and
on a system’s performance metrics and causes anomalies within the expressions obtained by integration or derivation with respect to parameter
[
14].
Considering that previously mentioned,
-based parameterization was recently proposed in [
14], by introducing parameters
K and
:
where
, for
obviously ensures linear dependence between
and
. On the other side, the definition of parameter
K remains unchanged, but the definition expression of normalized
with respect to the ratio between
and
, e.g.,
(where
is the Rician parameter), is affected by the choice of second parameter (
vs.
). Obviously, parameter
completely changes its character since
, while
does not, since
[
14].
Despite that aforementioned,
-based parameterization has up to now been considered only in [
14], by elaborating its benefits on ASEP observation accuracy. However, no other benefits of
-based parameterization or the anomalies caused by the nonphysical definition of a
parameter within the expressions which involve integration or derivation with respect to
have been presented yet. Accordingly, in this paper, one such anomaly is identified within
-parameterized TWDP moment-based estimators and overcome by introducing the appropriate estimators for
-parameterized TWDP fading model. In that sense, this paper focuses on performance analysis of moment-based parameter estimators, for both
- and
-based TWDP parameterizations.
Considering that aforementioned, in
Section 2, an overview of the results related to estimation of TWDP parameters is foremost presented, elaborating in detail anomalies caused by conventional
-based parameterization and indicating the absence of those related to a
-parameterized TWDP model. In
Section 3, a closed-form moment-based estimator is derived for
-based parameterization, and their behavior is examined for various values of proposed parameters. In
Section 4, performance analysis of proposed estimators is performed. First, corresponding asymptotic variances are derived and graphically presented and then the limits of estimation problem are explored and determined in terms of Cramer–Rao bound. Finally,
- and
-based estimators are compared in
Section 5, in order to gain insight into the differences in their behavior and their estimation accuracy. Conclusions are outlined in
Section 6.
2. Related Works
The estimation of parameters that characterize fading channel is of practical importance in a variety of wireless scenarios. It includes not only delay insensitive channel characterization and link budget calculations, but also online adaptive coding/modulation for which the estimation of parameters must be both accurate and prompt [
15]. In that sense, inaccurate estimation of channel parameters leads to a non-optimal channel capacity utilization, which might significantly affect highly mobile communication systems such as vehicular-to-vehicular and vehicular-to-infrastructure. Accordingly, different approaches in estimation of TWDP parameters are proposed to make a deal between their computational complexity and estimation accuracy. These approaches usually assume that the value of parameter
can be directly estimated from the data set as a second moment [
10], so the estimation problem has become focused on determination of values of
K and
.
Among the investigated approaches used to estimate parameters of
-based parameterization, the distribution fitting approach is used for measurements performed in air-craft and buses at 2.4 GHz [
3], while the maximum likelihood procedure (ML) is used for measurements performed at 60 GHz in the indoor environment [
11] and in vehicular-to-vehicular propagation scenarios [
10]. However, it is shown that both approaches are computationally very complex and inappropriate for online applications. Accordingly, the moment-based approach is considered in [
4,
7,
12], as a compromise between estimator’s complexity and its accuracy. Thereat, in [
7], estimators are derived only as conditional expressions which can not be used for practical estimations in which both parameters are unknown. To overcome the issue, the exact joint estimators of TWDP parameters are derived in [
4] as computationally simple expressions. However, for certain combinations of parameters
K and
, the estimator of parameter
,
, derived in [
4], provides physically unsubstantiated results, which can be demonstrated by applying ([
4], Equation (
16)) on Monte Carlo simulated samples.
Accordingly, for each combination of
K and
, where
and
, 500 realizations of the TWDP process with
i.i.d. samples are generated. The estimate of parameter
from
j-th realization,
, is obtained by ([
4], Equation (
16)) and, for each combination of
K and
, sample mean value
is calculated as
. The results are presented in
Figure 1.
From
Figure 1, it can be observed that, for large values of
(i.e., for
1),
estimates may exceed 1 regardless of the value of parameter
K. However, according to definition of
:
the parameter has to be lower or equal to one (
) for any
. Accordingly, estimates of
greater than 1 are physically unsubstantiated and can not be used to gain any insight into the relation between
and
.
Presented results related to sample means, when
is in the vicinity of 1, are also in contrast with underlying physical mechanisms. Namely, as magnitudes of two dominant waves become approximately the same (
), i.e., when
, it gets harder to discern the two rays [
16]. However,
Figure 1 shows that error in estimation of
(expressed by the difference between
and
and the dispersion of estimated values) in the vicinity of 1 is smaller than those obtained for moderate values of
. Accordingly, estimated values of
greater than 1, except for being useless for estimation of channel conditions, also lead to the false accurate sample means when
, obtained by averaging potentially accurate values of
smaller than 1 and unsubstantiated estimated values greater than 1.
On the other side, in the region of small and moderate values of , the proposed estimator provides meaningful results. However, its usage leads to huge estimation error of parameter , which is especially pronounced for small and moderate values of K.
Accordingly, due to the nonphysical definition of parameter
, although being analytically correct,
([
4], Equation (
16)) provides irrelevant results for certain combinations of parameters
K and
, and huge errors in estimation of
for some others. Therefore, it is very desirable to derive estimators for physically justified
-based parameterization and to investigate their behavior in different channel conditions.
3. Moment-Based Estimators for the -Parameterized TWDP Model
In order to derive moment-based estimators for the
-parameterized TWDP model, the expression for
n-th moments of a signal envelope
r given by ([
4], Equation (
6)) is expressed in terms of
, as:
and for even values of
n, obtained as a simple closed form expression:
where
is the expectation operator.
Obviously, the
n-th moment of signal’s envelope
depends on three unknown parameters:
K,
, and
. Consequently, estimators of these parameters can be constructed from at least three different moments. Thereby, it is shown that estimation accuracy is the largest when the lower-order moments are used [
4]. Accordingly, moment-based estimators of the
-parameterized TWDP model should be generated using second-, fourth-, and sixth-order moments, since only even moments can be obtained from (
8). Thereat, in order to further reduce the estimation complexity, the impact of parameter
on
K and
could be canceled out by properly defining ratios between the fourth- and the second-order as well as the sixth- and the second-order moments of an envelope [
4]. It finally leads to the system of two equations:
which could be solved for
and
if sample moments
are used instead of the ensemble averages
. In this regard, after (
9) is inserted in (
10), (
10) can be expressed as: is expressed as:
i.e., as:
Then, (
12) is inserted in (
9), transforming it into the following polynomial:
where
a and
b are defined as:
By analyzing the discriminant of the polynomial (
13), it can be shown that, for any values of
K and
, the considered polynomial has one real root and one pair of complex conjugate ones. Thus, following the [
17], one can find the real root of polynomial (
13) as:
where
gives the real part of the complex number, and
Z is defined as:
where
Considering the aforementioned, moment-based estimators of parameters
K and
can be obtained as simple, closed form expressions:
where parameters
a,
b,
Z,
p, and
q can be easily calculated by inserting the second-, the fourth-, and the sixth-order sample moment of an envelope into the Equations (
14), (
15), (
17), (
18) and (
19).
4. AsV and CRB
To further assess the performance of proposed estimators, corresponding asymptotic variances for
K and
,
and
, are calculated from (
20) and (
21), as [
4]:
where:
and CvM is a covariance matrix with elements:
In addition, Cramer–Rao lower bounds are numerically calculated as [
4]:
where the elements of Fisher Information Matrix
,
, are determined as:
and where
N is the number of observations,
is the closed-form TWDP envelope PDF given by [
14] [Equation (
7)], and
Thereby, following the approach presented in [
4,
7], instead of using CRB and AsV, the square root of CRB and AsV normalized to
N and the true value of estimating parameter are used in estimator performance assessment. Thus, to assess the estimation error, the sqrt-normalized CRB and AsV of
,
and
, are plotted in
Figure 4, and sqrt-normalized CRB and AsV of
,
and
, are plotted in
Figure 5.
From
Figure 4, it can be observed that the estimation error of
K increases with the decrease of parameter
K. Thereby, when the power of specular components
is small with respect to the power of diffuse components
(i.e., when the TWDP channel becomes Rayleigh-like), the error in estimation of
K is very large. However, as the value of parameter
K increases, i.e., as
overrides
, the estimation of parameter
K becomes very accurate.
Figure 4 also shows that the error in estimation of
K grows with the reduction of
, indicating that it becomes harder to accurately estimate
K as the specular component with the magnitude
becomes less significant with respect to one with the magnitude
.
From
Figure 4, it can also be observed that the values of sqrt-normalized
are remarkably close to the sqrt-normalized
for the entire considered range of parameters
K and
, indicating almost asymptotic efficiency of the proposed estimator of a parameter
K.
Figure 5 shows that the estimation error of parameter
behaves similarly as the estimation error of
, with respect to
K and
. Hence, the estimation of
deteriorates with the reduction of
K, i.e., as the power of diffuse components becomes more significant with respect to
. The estimation error of
is large for small values of
, indicating that it is hard to estimate the values of
when
is insignificant with respect to
. For moderate values of
,
given by (
21) starts to provide pretty accurate results, especially for large values of
K. However, as
approaches one, estimation of
becomes more inaccurate. In these conditions, the magnitudes of specular components,
and
, become similar. Considering the aforementioned and the fact that the phase difference between these components is uniform, the probability of destructive superposition of specular components becomes more likely, making their overall magnitude often insignificant. Thus, as
approaches one, it gets harder to accurately determine the value of
, especially when the power of diffuse components is large with respect to
.
When it comes to the observation of the proximity of the proposed
estimator to its CRB, it can be concluded from
Figure 5 that the values of the sqrt-normalized
are remarkably close to the sqrt-normalized
for
in the entire range of
, making the proposed estimator asymptotically efficient for the considered values of
K.
Accordingly, except for providing estimation errors significantly close to the corresponding CRBs, moment-based estimators (
20) and (
21) provide accurate estimates obtained from a relatively small number of samples, which can be clearly observed from
Figure 4 and
Figure 5. For example, if we assume that the sufficient accuracy of estimation process is
, it can be achieved with
samples for
and
(obtained by multiplying estimation errors
and
by
).
If necessary, the estimation accuracy in the determined region can be increased, or the region itself can be further expanded, by involving more samples within the estimation process (e.g., by employing
samples, relative estimation error in the considered region of parameters
K and
could be reduced to
, or the estimation error of
could be achieved for the wider range of
K and
, i.e.,
and
). In this way, the procedures used to create
Figure 4 and
Figure 5 can be used to determine the number of samples needed to obtain desired estimation accuracy within the desired range of parameters
K and
.
5. Comparison of Moment-Based Estimators for the - and -Parameterized TWDP Model
In order to observe qualitative differences between
- and
-based parameterization and to gain more precise insight into the relation between estimation errors for considered parameterizations,
and
are normalized to the same parameter ratio
and presented in
Figure 6. This enables us to compare absolute values of AsV for
and
and to discover differences in their estimation errors for considered ratios between
and
.
Figure 6 shows that, for
(which corresponds to
), estimation error of parameter
is up to two times smaller than the error obtained for
parameter estimation. On the other side, for
(i.e.,
) and
, there is no significant difference in accuracy of
and
. Finally, for
, error in estimation of
starts to increase with the increment of
, thus being in line with the model’s physical mechanisms. On the contrary, in the considered region,
-based parameterization provides false accurate results, obtained by also considering values of
greater than one in calculation of
.
Accordingly, in the region of , i.e., , estimation accuracy is significantly improved by using -based parameterization instead of -based, while in the region of , i.e., , -based parameterization prevents the occurrence of nonphysical solutions obtained by estimating parameter .
Except for benefits of
-based parameterization observed with respect to
, it also enables reduction of the estimation error of a parameter
K, for a much wider set of values of a parameter, which reflects the relation between
and
. Namely, based on the expression of parameter
K given in terms of
and the results presented in [
14] [Figure 2], it can be concluded that
is almost constant for the entire range of small and medium values of
, implying that values of
make almost no impact on the value of parameter
K. This causes quite pronounced errors in estimation of
K for the entire range of small and medium values of
(i.e.,
), which can be clearly observed from [
4] [
Figure 1]. On the contrary,
Figure 4 shows that, when
K is expressed in terms of
, no such anomaly can be observed. In these circumstances, errors in estimation of
K are huge only for small values of
(i.e.,
).