Next Article in Journal
Limitations of Foot-Worn Sensors for Assessing Running Power
Next Article in Special Issue
Efficient Channel Feedback Scheme for Multi-User MIMO Hybrid Beamforming Systems
Previous Article in Journal
A Framework for Maternal Physical Activities and Health Monitoring Using Wearable Sensors
Previous Article in Special Issue
Uplink Sparse Channel Estimation for Hybrid Millimeter Wave Massive MIMO Systems by UTAMP-SBL
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Asymptotic Efficiency of the M2M4 Signal-to-Noise Estimator for Deterministic Complex Sinusoids

Department of Engineering, University of Campania “L. Vanvitelli”, 81031 Aversa, CE, Italy
Sensors 2021, 21(15), 4950; https://doi.org/10.3390/s21154950
Submission received: 4 June 2021 / Revised: 9 July 2021 / Accepted: 16 July 2021 / Published: 21 July 2021
(This article belongs to the Special Issue Massive MIMO and mm-Wave Communications)

Abstract

:
The moment-based M 2 M 4 signal-to-noise (SNR) estimator was proposed for a complex sinusoidal signal with a deterministic but unknown phase corrupted by additive Gaussian noise by Sekhar and Sreenivas. The authors studied its performances only through numerical examples and concluded that the proposed estimator is asymptotically efficient and exhibits finite sample super-efficiency for some combinations of signal and noise power. In this paper, we derive the analytical asymptotic performances of the proposed M 2 M 4 SNR estimator, and we show that, contrary to what it has been concluded by Sekhar and Sreenivas, the proposed estimator is neither (asymptotically) efficient nor super-efficient. We also show that when dealing with deterministic signals, the covariance matrix needed to derive asymptotic performances must be explicitly derived as its known general form for random signals cannot be extended to deterministic signals. Numerical examples are provided whose results confirm the analytical findings.

1. Introduction

In [1], the authors address the problem of estimating the signal-to-noise ratio (SNR) in the case of a complex sinusoidal signal with a deterministic but unknown phase corrupted by additive Gaussian noise. A general signal model that, for example, represents the complex envelope of frequency modulated signals with deterministic modulating signals. The problem is found in many fields of practical interest such as in digital communication and radar systems (see, for example, [1] and references therein). The authors propose an estimator based on the sample-moment estimators of the second and fourth orders ( M 2 M 4 SNR estimator) and study its performances through numerical examples. The results let them conclude that the proposed SNR estimator is asymptotically efficient and presents finite sample super-efficiency in some cases that depend on the values of the signal and noise power.
In this paper, we show that the proposed SNR estimator is neither (asymptotically) efficient nor super-efficient, in any case. In fact, we study the performances of the estimator through asymptotic analytical analysis besides numerical examples and show that, though the estimator is not efficient, the performances are acceptable in many cases of interest and not too distant from the Cramer–Rao lower bounds (CRLB).
The M 2 M 4 SNR estimator belongs to a wider class of moment-based SNR estimators that have been proposed mainly in the context of modulated telecommunication signals. As reminded in [2], a similar estimator was first proposed by Benedict and Soong in real AWGN channels [3]. Subsequently, in [4], Matzner derived a similar expression and, in [5], together with Englberger, obtained the same result using a different approach. In [6], the estimation of SNR for non-constant modulus constellations over a frequency-flat fading channel is considered, and a family of estimators that use higher-order moments is derived, and their performances studied. More specifically, the authors show that the M 1 M 2 and M 2 M 4 estimators contained in the family achieve the corresponding CRLBs for constant modulus constellation, while for multilevel constellations, performances degrade as the SNR increases. A different family of moment-based SNR estimators that makes use of the second, fourth and sixth moments is proposed in [7], with the purpose to improve the performances in the case of non-constant modulus constellations for increasing SNRs. Subsequently, in [8], the authors propose an extended family of moment-based SNR estimators that makes use of higher-order moments to further improve the performances with non-constant modulus constellations. In [1], a different signal model is considered, as the authors assume a complex exponential with a deterministic but unknown phase rather then randomly modulated signals. In practice, the authors consider a signal model very similar to the M-PSK vector model for the AWGN channel, with the only difference being the deterministic nature of the unknown phase. Under this assumption, the derivation in [1] is quite cumbersome. We show in this paper that a more general method can be used instead, as suggested by Kay [9]. The result is an M 2 M 4 SNR estimator that has the same form of the similar estimator obtained for M-PSK constellations.
The M 2 M 4 SNR estimator is simple to implement, has low computational complexity and is blind as it does not require prior knowledge of signal or noise. However, the performance of this estimator must be evaluated by numerical examples with Monte Carlo simulations and/or asymptotic analytical analyses. In fact, exact analysis is difficult as mathematics quickly becomes intractable as it involves the transformation of random variables to obtain the estimator pdf. In [6,7,8], the performances of the proposed moment-based estimators are studied in terms of asymptotic variance and CRLBs. However, care is needed when analytical asymptotic performances are derived in the deterministic case, since, in general, it is not possible to extend the results obtained for the random signals to deterministic signals, and in this paper, we extend to the complex case what has been already shown in the case of real deterministic a sinusoid in [10]. Since the derivation of asymptotic variances depends on a covariance matrix that, in turn, depends on the signal model, we show that, in general, the covariance matrix must be derived explicitly in the case of deterministic signals. Nevertheless, it turns out that, in the specific case of complex sinusoids with deterministic phase, the two covariance matrices coincide.
The paper is organized as follows: in Section 2, it is shown the derivation of even-order moments and estimators for signal and noise power as well as the SNR; in Section 3, asymptotic variances for all the estimators under investigation are derived; details of derivations are presented in Appendix A and Appendix B; we provide some numerical examples in Section 4; finally, in Section 5, we draw the conclusions.

2. Moment-Based Estimators

Let us consider the sequence of K samples
y k = S x k + w k ,
where x k = e j ϕ k , and ϕ k [ 0 , 2 π ) is a deterministic unknown sequence, S is a real positive scalar, and w k is an additive complex Gaussian random noise with zero mean and variance equal to N, i.e., w k CN ( 0 , N ) . The general model (1) includes, as special case, ϕ k = 2 π ν 0 k , i.e., a pure complex sinusoid with normalized frequency ν 0 ; other choices for ϕ k are possible as shown, for example, in [1]. The base assumption made here is that ϕ k is deterministic in contrast, for example, to the case of phase modulated signals where the phase is random [2].
The signal-to-noise ratio (SNR), defined as ρ S / N , is unknown and represents the parameter we wish to estimate when the instantaneous phase ϕ k is also unknown. In [1], an estimator for ρ based on the method of moments, built as a ratio of the estimators of S and N and that employs the second and fourth moment of y k ( M 2 M 4 estimator has been proposed, and its performances presented through numerical examples only. In this paper, we show that the same M 2 M 4 estimator can be derived in an alternative, more general way, and we present asymptotic analytical performances, along with new numerical examples, that contradict some of the results obtained in [1].
The method of moments is a well-known general statistical method to derive estimators as a function of high-order sample-moments [9]. The key idea is to express the parameter to be estimated as a function of the moments of y k of different orders, i.e., μ m = E y k m , where m is an integer, and use the natural estimators of the moments μ ^ m ( 1 / K ) k = 0 K 1 y k m in place of the true moments μ m to obtain an estimate of the parameter from the observed samples. Such estimators typically present good statistical performances that can be derived through either numerical simulations or asymptotic analytical methods. In general, the asymptotic efficiency of an estimator derived using this method cannot be guaranteed; however, in many cases, the estimator turns out to be consistent.
The method relies on the possibility of expressing the moments μ m as a function of the parameters to be estimated, i.e., S and N in this case. However, when in the observed signal model, one term is deterministic, as in the signal model (Equation (1)), the derivation of the moments might be cumbersome, as shown, for example, in [1]. For these cases, Kay [9] suggests to assume a random nature for the “deterministic” unknown phase. In this paper, we follow such an approach to obtain the same M 2 M 4 estimator derived in [1] by assuming that ϕ k is an r.v. independent of w k . As consequence, the general closed form expression for the the even-order moments of y k , μ 2 m , obtained in [7] for general QAM signals, applies also for the signal model (Equation (1)), i.e.,
μ 2 m E [ | y k | 2 m ] = n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E | x k | 2 n S n N ( m n ) .
Note that Equation (2), obtained under the random phase assumption, is equivalent to the expression
μ 2 m = n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 1 K k = 0 K 1 | x k | 2 n S n N ( m n ) ,
under the assumption that K . Since with x k = e j ϕ k , we have E [ | x k | 2 n ] = 1 regardless of the distribution of ϕ k , Equation (2) simplifies to
μ 2 m = n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 S n N ( m n ) .
The above formula gives, as special cases, the second and fourth moments needed to form the M 2 M 4 estimator,
μ 2 = S + N
μ 4 = S 2 + 4 S N + 2 N 2 ,
from which we obtain both S and N as functions of μ 2 and μ 4 as follows
S = 2 μ 2 2 μ 4
N = μ 2 S .
Then, by using the sample moments
μ ^ 2 = 1 K k = 0 K 1 y k 2
μ ^ 4 = 1 K k = 0 K 1 y k 4
in place of the true moments in Equations (7) and (8), we obtain the following estimators for S and N
S ^ = 2 μ ^ 2 2 μ ^ 4
N ^ = μ ^ 2 S ^ ,
where we apply the modulus to ensure that the estimates are always real and positive. Finally, the above estimators are then used to form the M 2 M 4 SNR estimator
ρ ^ = S ^ N ^ .
Not surprisingly, the estimator takes on the same form as the M 2 M 4 estimator proposed for phase modulated signals [2]. However, it is important to note that when it comes to deriving asymptotic performances, as we show in the following section, we cannot assume a random nature for the deterministic signals, and results obtained under such an assumption cannot in principle be readily extended to the signal model (1).

3. Asymptotic Performances

Analytical performances in terms of squared biases and variances of moment-based estimators can only be derived asymptotically for sufficiently large number of observed samples, i.e., K . A method based on a first-order approximation analysis is suggested, for example, in [9], with the purpose of overcoming the mathematical intractability of the method that uses the transformation of random variables to obtain the pdf of the estimator. Alternatively, performances of the moment-based estimators may be obtained by estimating squared biases and variances through Monte Carlo simulations. We adopt both approaches, and we will show that, as expected, the analytical results confirm what we obtain through numerical simulations.
Let us start by defining the sample moment estimator over K observations as
T m y = 1 K k = 0 K 1 y k m
and g T is the function that maps the sample moments to the estimate, and T is the vector of the sample moments. More specifically, we define for the estimator S ^ , the function g T = 2 T 2 2 T 4 , with T = T 2 , T 4 , as we use the second and fourth moments; the estimator N ^ is defined by the function h T = T 2 g T ; the estimator ρ ^ is defined by the function f T = g T / h T .
In general, from the first-order Taylor expansion of g ( T ) about the point T = E T = μ , we know that the mean of the estimator g T is equal to g μ , and by simple substitution using Equations (5) and (6), it can be shown that the estimators are asymptotically unbiased. The same applies to the other estimators defined by h ( T ) and f T , respectively.
We now turn our attention to variances. The first-order approximation of the variance of the estimator S ^ is given by
var S ^ g T T = μ T C T g T T = μ ,
where C T is the covariance matrix of T that by definition is equal to
C T = E T μ T μ T .
In general, for random signals x k in Equation (1), it is already known that each element of the covariance matrix is given by
C T i j = 1 K μ 2 ( i + j ) μ 2 i μ 2 j ,
where μ 2 m is the true even-order moment. However, in the case of generic deterministic signals x k , Equation (17) cannot be readily applied. Nevertheless, for the specific signal under consideration, i.e., the complex exponential signal with a deterministic unknown phase, it turns out that the elements of the covariance matrix can still be expressed by Equation (17), which is valid for any K, even for a small K. Details of derivation are shown in Appendix A.
In Equation (15), we also need the partial derivatives of g T , i.e.,
g T 2 = 2 T 2 g T 2 , T 4
and
g T 4 = 1 2 1 g T 2 , T 4 .
Similarly, we know that the first-order approximation of the variance of the estimator N ^ is given by
var N ^ h T T = μ T C T h T T = μ
where
h T 2 = 1 g T 2
h T 4 = g T 4 .
Finally, we have that the variance of the estimator ρ ^ is well approximated by
var ρ ^ f T T = μ T C T f T T = μ ,
where
f T 2 = g T 2 1 h g h T 2 1 h 2
f T 4 = g T 4 1 h g h T 4 1 h 2 .
Based on the same approach, it is possible to derive asymptotic performances for the sample moment estimators of order m. Let g 2 m ( T 2 m ) = T 2 m , the function that maps the sample moments to the sample moment estimate, and from Equation (15), we have
var μ ^ 2 m g 2 m T 2 m T 2 m = μ 2 m T C T 2 m g 2 m T 2 m T 2 m = μ 2 m ,
where
g 2 m T 2 m T 2 m = μ 2 m = T 2 m T 2 m T 2 m = μ 2 m = 1
and the covariance matrix is
C T 2 m = 1 K μ 4 m μ 2 m 2 .
Putting it all together
var μ ^ 2 m 1 K μ 4 m μ 2 m 2 .
For convenience, we conclude this section by providing the Cramer–Rao lower bounds (CRLB) for all the estimators under analysis because they provide a significant benchmark for the performances of all corresponding estimators.
CRLBs for S ^ and N ^ are immediately given by the Fisher information matrix as derived in [1]. Note that in the derivation of the CRLB, the phases ϕ k should be regarded as nuisance parameters. Thus, the true Fisher Information Matrix (FIM) has size ( K + 2 ) × ( K + 2 ) . However, since ϕ k ’s are uncoupled with S and N, the FIM turns out to be diagonal, and the final result does not change. The CRLBs are as follows:
I ( S , N ) = K 2 S N 0 0 K N 2
and we have
var ( S ^ ) 2 S N K
var ( N ^ ) N 2 K .
For ρ ^ expressed in dB, i.e., ρ ^ d B , we have
var ( ρ ^ d B ) 1 K 10 ln 10 2 2 N S + 1 dB 2 .
CRLBs for μ ^ 2 m , with m = 1 , 2 , , are derived in Appendix B and can be written as
var ( μ ^ 2 m ) μ 2 m S 2 2 S N K + μ 2 m N 2 N 2 K ,
where
S μ 2 m ( S , N ) = n = 0 m ( m ! ) 2 n ( m n ) ! ( n ! ) 2 S n 1 N ( m n )
and
N μ 2 m ( S , N ) = n = 0 m ( m ! ) 2 ( m n ) ( m n ) ! ( n ! ) 2 S n N ( m n 1 ) .

4. Numerical Examples

In this section, we present some examples that show the performances of the proposed moment-based estimators for S and N and signal-to-noise ratio ρ in terms of squared bias and variance. We report estimated squared biases and variances obtained through Monte Carlo simulations, and we show that numerical results confirm the asymptotic performances derived in Section 3. We also compare them to the corresponding CRLBs, and we make some comments on the property of efficiency of the estimators.
For all estimators under study, we have estimated the squared bias and variance by averaging over L = 2 13 = 8192 estimates on a varying number of observations K. The mean of each estimator is estimated by means of a standard sample mean estimator, e.g., for S ^
E [ S ^ ] μ ^ S = 1 L l = 1 L S ^ l ,
where S ^ l is a single estimate of S with K observations. The squared bias is then obtained by using the estimated mean in place of the true mean, e.g., b 2 S ^ = E [ S ^ ] S 2 . The variance of estimators is estimated by means of the standard sample variance estimator, e.g., for S ^ .
var ( S ^ ) σ ^ S 2 = 1 L 1 l = 1 L ( S ^ l μ ^ S ) 2 .
The estimation of squared bias and variance of estimators has been performed at different SNRs that, without loss of generality, are obtained by varying the power of the complex sinusoid S while maintaining the noise power fixed at N = 4 . The instantaneous phase of the complex exponential has been chosen to be ϕ k = 2 π ν 0 k , where k is the discrete-time index and ν 0 = 0.1234 . Note that the results do not depend on the specific form of the instantaneous phase, and other choices are possible, such as the phase used in [1].
In the following, we report results not only for the SNR estimator ρ ^ but also for estimators of the second and fourth moments and for estimator of the signal power S ^ and the noise power N ^ .

4.1. Performances of μ ^ 2

The sample moment estimator μ ^ 2 is, under first order approximation, unbiased. On the other hand, our focus is on the variance of μ ^ 2 , which has been estimated for several SNRs obtained by varying S with N = 4 . We report the results in Figure 1, where the corresponding CRLB and asymptotic performances are shown as well. For all the SNRs under consideration, we see that we have complete superposition of the three curves even for a small number of observations K and that the sample moment estimator is, then, always consistent and efficient. Though performances are always optimal, the variance tends to increase as the SNR increases.

4.2. Performances of μ ^ 4

The sample moment estimator μ ^ 4 is, under first order approximation, unbiased. The performances in terms of variance of μ ^ 4 are shown in Figure 2, where estimated variance, CRLB and asymptotic variance are shown at different SNRs obtained by varying S with N = 4 . In all cases under consideration, the estimated variance and the asymptotic variance are superimposed even for a small numbers of observation K, confirming, also in this case, the validity and usefulness of the asymptotic analysis carried out above. While in all cases the sample estimator μ ^ 4 results are consistent, we see that at low SNRs, the estimator can be regarded as nearly efficient because the variance of μ ^ 4 is not superimposed on the CRLB, though it is very close (as shown, for example, in Figure 2a–f). Note that the variance of the estimator increases to significantly large values as the SNR becomes larger by increasing the power S.

4.3. Performances of S ^

In Figure 3, we report the results of the estimation of the squared bias of S ^ at different SNRs obtained by setting N = 4 and varying S. We see that in all reported cases, the squared bias is sufficiently small even for a small number of observed samples K, and then we find that S ^ is practically (asymptotically) unbiased.
Estimates of the variance of S ^ are shown in Figure 4 together with the corresponding asymptotic variances and CRLBs. First, the results confirm the derivation of the asymptotic performances obtained in Section 3. Though valid only for K , it is interesting to note that the asymptotic results are achieved also for a small number of samples in most of the cases under consideration. In all cases, with K = 2 18 , the variance of S ^ is well approximated by the asymptotic variance that shows that the estimator is consistent. It is not always efficient, especially at a relatively small SNR. Consider, for example, Figure 4a, which shows the variance for S = 1 and N = 4 , corresponding to the SNR equal to −6.0206 dB. The asymptotic variance, which is achieved by the estimator with at least K = 2 12 samples, cannot achieve the corresponding CRLB. This is also true in Figure 4b,c. On the other hand, starting from Figure 4d, we see that the asymptotic variance is superimposed to the corresponding CRLB, and then the estimator becomes efficient. For a sufficiently large SNR, the asymptotic variance is able to predict the effective variance of the estimator S ^ even for a very small number of samples K.

4.4. Performances of N ^

In Figure 5, the estimated squared bias of the estimator N ^ at different SNRs with varying S and N = 4 is shown. As for S ^ , we find that, in practice, we can consider N ^ as an unbiased estimator, as expected by the first-order approximation analysis.
Under the same setup, we have obtained the estimation of the variance of N ^ , and we show the results in Figure 6 together with the asymptotic variances and the corresponding CRLBs obtained in Section 3. Note that we maintain the noise power fixed while varying the signal power S, and therefore, we have the same CRLB on all the figures.
In all cases, the estimated variances achieve the corresponding asymptotic variances, confirming the analytical results obtained in Section 3 and showing that the estimator is consistent.
In general, the asymptotic variance describes the current variance of N ^ very well, even for a small number of samples, with the exception of some cases at low SNRs, as shown in Figure 6a–c. Furthermore, the variance never achieves the corresponding CRLB, and then the estimator is not efficient.

4.5. Performances of ρ ^

The estimated squared biases and variances of the SNR estimator ρ ^ have been obtained for a number of different SNRs with Monte Carlo simulations run with N = 4 and varying the signal power S, as conducted previously with S ^ and N ^ . The results are reported in Figure 7 and Figure 8.
From the plots in Figure 7, we have confirmation that the estimator is practically unbiased.
In Figure 8, we report the estimated variances along with the corresponding asymptotic variances and CRLBs as derived in Section 3. First, the results confirm once again the validity of the analysis carried out in Section 3. In all cases, the estimated variances achieve the corresponding asymptotic variances, showing that the estimator is consistent. Note that the asymptotic variance is achieved with a small number of samples, with the exception of small SNRs, as we have already seen for the estimator of noise power N ^ . Clearly, the limits of the noise power estimator are also reflected in the performances of the SNR estimator as the SNR estimator is obtained as the ratio of S ^ and N ^ . For the same reason, the estimator ρ ^ cannot achieve the performances given by the CRLB, and then it is not efficient, though the difference between the variance of ρ ^ and the corresponding CRLB is very small, especially at high SNRs.
The results obtained in this paper do not confirm what had been concluded in [1], where the authors stated that the same estimator was not only efficient, but in some cases, it was also super-efficient, i.e., with the variance of the estimator below the CRLB. The authors in [1] admit that the conclusion is based only on numerical results obtained in the case of S = 100 and N = 4 that corresponds to ρ = 13.9794 dB and reported in Figure 4 in [1]. Specifically, in Figure 8e, it is shown that we have obtained in this paper the variance of ρ ^ both numerically and mathematically with the same parameters used in Figure 4 in [1]. We have proven that the asymptotic value, which represents with very good approximation the variance of ρ ^ even for a small number of samples, is always above the CRLB, and therefore, even in this specific case, the estimator is neither efficient or super-efficient.

5. Conclusions

The M 2 M 4 SNR estimator for unknown deterministic complex phase signals proposed in [1] represents a useful tool for blind estimation of SNR. However, its performances, while acceptable in many applications of interest, are not those stated in [1]. In fact, through detailed analytical and numerical analysis of the performances presented in this paper, we could not confirm the finite sample super-efficiency or (asymptotically) efficiency of the M 2 M 4 SNR estimator claimed in [1], which was based on numerical examples only.
In this paper, we have also derived analytical asymptotic performances for all the intermediate estimators that are used to form the M 2 M 4 SNR estimator, namely even-order sample moments, signal and noise power. The results show that though the second and fourth sample moment estimators are efficient, as well as the estimator of the signal power for sufficiently high SNRs, the noise power estimator is not efficient. Consequently, the SNR is not efficient though its performances are close to the corresponding CRLB. In contrast to the case of signal and noise power, the performances of the SNR estimator do not degrade as the true SNR increases.
The derivation of analytical asymptotic performances revealed that, in general, it is not possible to readily extend the results obtained for random signals to the case of deterministic signals. More specifically, the covariance matrix required to compute the asymptotic variances needs to be explicitly derived. Nevertheless, in the specific case considered in this paper, i.e., complex sinusoid with deterministic phase, it turns out that the covariance matrix has the same form as that in the random case.

Funding

This research received no external funding.

Conflicts of Interest

The author declares no conflict of interest.

Appendix A. Derivation of Covariance Matrix

In this section, we provide details of the derivation of Equation (17). Let us start from the definition of the covariance matrix, where each element is given by
C T i j = cov T 2 i , T 2 j
E 1 K k = 0 K 1 | y k | 2 i 1 K p = 0 K 1 | y p | 2 j E 1 K k = 0 K 1 | y k | 2 i E 1 K p = 0 K 1 | y p | 2 j .
Now, consider the first term of the right-hand side. We get
E 1 K k = 0 K 1 | y k | 2 i 1 K p = 0 K 1 | y p | 2 j = 1 K 2 k = 0 K 1 p = 0 K 1 E y k 2 i y p 2 j
= 1 K 2 k = 0 K 1 E y k 2 i + j + 1 K 2 k = 0 K 1 p = 0 , p k K 1 E y k 2 i y p 2 j
= 1 K 2 k = 0 K 1 E y k 2 i + j + 1 K 2 k = 0 K 1 p = 0 , p k K 1 E y k 2 i E y p 2 j
Similarly, we obtain that the second term on the right-hand side of (A2) is
E 1 K k = 0 K 1 | y k | 2 i E 1 K p = 0 K 1 | y p | 2 j = 1 K 2 k = 0 K 1 E | y k | 2 i E | y k | 2 j + 1 K 2 k = 0 K 1 p = 0 , p k K 1 E | y k | 2 i E | y p | 2 j .
Putting it all together, we have that each element of the covariance matrix is
C T i j = 1 K 2 k = 0 K 1 E y k 2 i + j 1 K 2 k = 0 K 1 E | y k | 2 i E | y k | 2 j .
The first term can be written explicitly as
1 K E 1 K k = 0 K 1 | y k | 2 i + j = p = 0 i + j ( ( i + j ) ! ) 2 ( i + j p ) ! ( p ! ) 2 1 K k = 0 K 1 | x k | 2 p S p N ( i + j p ) ,
while the second term becomes
1 K 2 k = 0 K 1 E | y k | 2 i E | y k | 2 j = = 1 K k = 0 K 1 p = 0 i ( i ! ) 2 ( i p ) ! ( p ! ) 2 | x k | 2 p S p N ( i p ) q = 0 j ( j ! ) 2 ( j q ) ! ( q ! ) 2 | x k | 2 q S q N ( j q ) = = 1 K k = 0 K 1 p = 0 i q = 0 j ( i ! ) 2 ( i p ) ! ( p ! ) 2 ( j ! ) 2 ( j q ) ! ( q ! ) 2 | x k | 2 p | x k | 2 q S p S q N ( i p ) N ( j q ) = = p = 0 i q = 0 j ( i ! ) 2 ( i p ) ! ( p ! ) 2 ( j ! ) 2 ( j q ) ! ( q ! ) 2 1 K k = 0 K 1 | x k | 2 p | x k | 2 q S p S q N ( i p ) N ( j q ) .
For x k = e j ϕ k , we have an exact expression valid even for small K (we do not need to assume K ), and we obtain
1 K k = 0 K 1 E | y k | 2 ( i + j ) = p = 0 i + j ( ( i + j ) ! ) 2 ( i + j p ) ! ( p ! ) 2 S p N ( i + j p ) = μ 2 ( i + j ) .
and
1 K k = 0 K 1 E | y k | 2 i E | y k | 2 j = p = 0 i q = 0 j ( i ! ) 2 ( i p ) ! ( p ! ) 2 ( j ! ) 2 ( j q ) ! ( q ! ) 2 S p S q N ( i p ) N ( j q )
= μ 2 i μ 2 j .
In summary, we have that, for signals in the form of complex exponential with deterministic phase sequence, the elements of the covariance matrix are given by
C T i j = 1 K μ 2 ( i + j ) μ 2 i μ 2 j .
Note that we obtain the same result as in the general random signal case and that Equation (A12) holds for all K. This is not true in all cases, as it represents a specific result due to the specific form of the deterministic signal x k . In general, the elements of the covariance matrix in the case of deterministic signal is not given by Equation (A12) because clearly
1 K k = 0 K 1 | x k | 2 p | x k | 2 q 1 K k = 0 K 1 | x k | 2 p 1 K k = 0 K 1 | x k | 2 q .
It is still possible to write the covariance matrix in the general case, but its expression requires derivation of the terms
1 K k = 0 K 1 | x k | 2 p
and
1 K k = 0 K 1 | x k | 2 p | x k | 2 q
that depends obviously on signal x k and typically is only possible under the assumption that K . For an example, see [10], where the covariance matrix for deterministic real sinusoids has been derived.

Appendix B. Derivation of CRLBs for μ ^ 2 m

The CRLBs for μ ^ 2 m can be derived from the expression of even-order moments given above through the derivatives of μ 2 m with respect to S and N. We have for random signals
S μ 2 m ( S , N ) = S n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S n N ( m n )
= n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S S n N ( m n )
= n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] n S n 1 N ( m n )
= n = 0 m ( m ! ) 2 n ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S n 1 N ( m n )
Analogously, we have
N μ 2 m ( S , N ) = N n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S n N ( m n )
= n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S n N N ( m n )
= n = 0 m ( m ! ) 2 ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S n m n N ( m n 1 )
= n = 0 m ( m ! ) 2 ( m n ) ( m n ) ! ( n ! ) 2 E [ | x k | 2 n ] S n N ( m n 1 )
The CRLB is obtained as
var ( μ ^ 2 m ) μ 2 m S μ 2 m N I 1 S , N μ 2 m S μ 2 m N T
that can be written as
var ( μ ^ 2 m ) μ 2 m S 2 2 S N K + μ 2 m N 2 N 2 K .

References

  1. Sekhar, S.C.; Sreenivas, T.V. Signal-to-noise ratio estimation using higher-order moments. Signal Process. 2006, 86, 716–732. [Google Scholar] [CrossRef]
  2. Pauluzzi, D.R.; Beaulieu, N.C. A comparison of SNR estimation techniques for the AWGN channel. IEEE Trans. Commun. 2000, 48, 1681–1691. [Google Scholar] [CrossRef]
  3. Benedict, T.R.; Soong, T.T. The Joint Estimation of Signal and Noise from the Sum Envelope. IEEE Trans. Inf. Theory 1967, 13, 447–454. [Google Scholar] [CrossRef]
  4. Matzner, R. An SNR estimation algorithm for complex baseband signals using higher order statistics. Facta Univ. 1993, 6, 41–52. [Google Scholar]
  5. Matzner, R.; Englberger, F. An SNR estimation algorithm using fourth-order moments. In Proceedings of the 1994 IEEE International Symposium on Information Theory, Trondheim, Norway, 27 June–1 July 1994; p. 119. [Google Scholar]
  6. Ping, G.; Tepedelenlioglu, C. SNR estimation for nonconstant modulus constellations. IEEE Trans. Signal Process. 2005, 53, 865–870. [Google Scholar] [CrossRef]
  7. López-Valcarce, R.; Mosquera, C. Sixth-order statistics-based non-data-aided SNR estimation. IEEE Commun. Lett. 2007, 11, 351–353. [Google Scholar] [CrossRef]
  8. Álvarez-Díaz, M.; López-Valcarce, R.; Mosquera, C. SNR estimation for multilevel constellations using higher-order moments. IEEE Trans. Signal Process. 2010, 58, 1515–1526. [Google Scholar] [CrossRef]
  9. Kay, S. Fundamentals Of Statistical Signal Processing; Prentice Hall: Hoboken, NJ, USA, 2001; Volume 1. [Google Scholar]
  10. Romano, G. Asymptotic Performances of a Signal-To-Noise Ratio Moment-Based Estimator for Real Sinusoids in Additive Noise. Appl. Sci. 2019, 9, 5556. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Asymptotic and estimated variances of the moment estimator μ ^ 2 versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 1. Asymptotic and estimated variances of the moment estimator μ ^ 2 versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g001
Figure 2. Asymptotic and estimated variances of the moment estimator μ ^ 4 versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 2. Asymptotic and estimated variances of the moment estimator μ ^ 4 versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g002
Figure 3. The estimated squared bias b 2 ( S ) ^ of signal power estimator S ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 3. The estimated squared bias b 2 ( S ) ^ of signal power estimator S ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g003
Figure 4. Asymptotic and estimated variances of the signal power estimator S ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 4. Asymptotic and estimated variances of the signal power estimator S ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g004aSensors 21 04950 g004b
Figure 5. The estimated squared bias b 2 ( N ) ^ of the noise power estimator N ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 5. The estimated squared bias b 2 ( N ) ^ of the noise power estimator N ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g005aSensors 21 04950 g005b
Figure 6. Asymptotic and estimated variances of the noise power estimator N ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 6. Asymptotic and estimated variances of the noise power estimator N ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g006
Figure 7. The estimated squared bias b 2 ( ρ ) ^ of the SNR estimator ρ ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 7. The estimated squared bias b 2 ( ρ ) ^ of the SNR estimator ρ ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g007
Figure 8. Asymptotic and estimated variances of the SNR estimator ρ ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Figure 8. Asymptotic and estimated variances of the SNR estimator ρ ^ versus the number of observations K. Each sub figure (al) shows results for a given nominal value of SNR obtained by changing the signal power S and with noise power N = 4 .
Sensors 21 04950 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Romano, G. On Asymptotic Efficiency of the M2M4 Signal-to-Noise Estimator for Deterministic Complex Sinusoids. Sensors 2021, 21, 4950. https://doi.org/10.3390/s21154950

AMA Style

Romano G. On Asymptotic Efficiency of the M2M4 Signal-to-Noise Estimator for Deterministic Complex Sinusoids. Sensors. 2021; 21(15):4950. https://doi.org/10.3390/s21154950

Chicago/Turabian Style

Romano, Gianmarco. 2021. "On Asymptotic Efficiency of the M2M4 Signal-to-Noise Estimator for Deterministic Complex Sinusoids" Sensors 21, no. 15: 4950. https://doi.org/10.3390/s21154950

APA Style

Romano, G. (2021). On Asymptotic Efficiency of the M2M4 Signal-to-Noise Estimator for Deterministic Complex Sinusoids. Sensors, 21(15), 4950. https://doi.org/10.3390/s21154950

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