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 (
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
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
and
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
SNR estimator that has the same form of the similar estimator obtained for
M-PSK constellations.
The
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
where
, and
is a deterministic unknown sequence,
S is a real positive scalar, and
is an additive complex Gaussian random noise with zero mean and variance equal to
N, i.e.,
. The general model (
1) includes, as special case,
, i.e., a pure complex sinusoid with normalized frequency
; other choices for
are possible as shown, for example, in [
1]. The base assumption made here is that
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
, is unknown and represents the parameter we wish to estimate when the instantaneous phase
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
(
estimator has been proposed, and its performances presented through numerical examples only. In this paper, we show that the same
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
of different orders, i.e.,
, where
m is an integer, and use the natural estimators of the moments
in place of the true moments
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
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
estimator derived in [
1] by assuming that
is an r.v. independent of
. As consequence, the general closed form expression for the the even-order moments of
,
, obtained in [
7] for general QAM signals, applies also for the signal model (Equation (
1)), i.e.,
Note that Equation (
2), obtained under the random phase assumption, is equivalent to the expression
under the assumption that
. Since with
, we have
regardless of the distribution of
, Equation (
2) simplifies to
The above formula gives, as special cases, the second and fourth moments needed to form the
estimator,
from which we obtain both
S and
N as functions of
and
as follows
Then, by using the sample moments
in place of the true moments in Equations (
7) and (8), we obtain the following estimators for
S and
N
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
SNR estimator
Not surprisingly, the estimator takes on the same form as the
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.,
. 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
and
is the function that maps the sample moments to the estimate, and
is the vector of the sample moments. More specifically, we define for the estimator
, the function
, with
, as we use the second and fourth moments; the estimator
is defined by the function
; the estimator
is defined by the function
.
In general, from the first-order Taylor expansion of
about the point
, we know that the mean of the estimator
is equal to
, 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
and
, respectively.
We now turn our attention to variances. The first-order approximation of the variance of the estimator
is given by
where
is the covariance matrix of
that by definition is equal to
In general, for random signals
in Equation (
1), it is already known that each element of the covariance matrix is given by
where
is the true even-order moment. However, in the case of generic deterministic signals
, 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
, i.e.,
and
Similarly, we know that the first-order approximation of the variance of the estimator
is given by
where
Finally, we have that the variance of the estimator
is well approximated by
where
Based on the same approach, it is possible to derive asymptotic performances for the sample moment estimators of order
m. Let
, the function that maps the sample moments to the sample moment estimate, and from Equation (
15), we have
where
and the covariance matrix is
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
and
are immediately given by the Fisher information matrix as derived in [
1]. Note that in the derivation of the CRLB, the phases
should be regarded as nuisance parameters. Thus, the true Fisher Information Matrix (FIM) has size
. However, since
’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:
and we have
For
expressed in dB, i.e.,
, we have
CRLBs for
, with
, are derived in
Appendix B and can be written as
where
and
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
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
where
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.,
. The variance of estimators is estimated by means of the standard sample variance estimator, e.g., for
.
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
. The instantaneous phase of the complex exponential has been chosen to be
, where
k is the discrete-time index and
. 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 and the noise power .
4.1. Performances of
The sample moment estimator
is, under first order approximation, unbiased. On the other hand, our focus is on the variance of
, which has been estimated for several SNRs obtained by varying
S with
. 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
The sample moment estimator
is, under first order approximation, unbiased. The performances in terms of variance of
are shown in
Figure 2, where estimated variance, CRLB and asymptotic variance are shown at different SNRs obtained by varying
S with
. 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
results are consistent, we see that at low SNRs, the estimator can be regarded as nearly efficient because the variance of
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
In
Figure 3, we report the results of the estimation of the squared bias of
at different SNRs obtained by setting
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
is practically (asymptotically) unbiased.
Estimates of the variance of
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
, 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
, the variance of
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
and
, corresponding to the SNR equal to −6.0206 dB. The asymptotic variance, which is achieved by the estimator with at least
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
even for a very small number of samples
K.
4.4. Performances of
In
Figure 5, the estimated squared bias of the estimator
at different SNRs with varying
S and
is shown. As for
, we find that, in practice, we can consider
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
, 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
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
and varying the signal power
S, as conducted previously with
and
. 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
. 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
and
. 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
and
that corresponds to
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
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
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 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.