As you progress through the various stages of learning CSP (intimidation, frustration, elucidation, puzzlement, and finally smooth operation), the symmetries of the various functions come up over and over again. Exploiting symmetries can result in lower computational costs, quicker debugging, and easier mathematical development.
What exactly do we mean by ‘symmetries of parameters?’ I’m talking primarily about the evenness or oddness of the time-domain functions in the delay and cycle frequency variables and of the frequency-domain functions in the spectral frequency and cycle frequency variables. Or a generalized version of evenness/oddness, such as , where and are closely related functions. We have to consider the non-conjugate and conjugate functions separately, and we’ll also consider both the auto and cross versions of the parameters. We’ll look at higher-order cyclic moments and cumulants in a future post.
You can use this post as a resource for mathematical development because I present the symmetry equations. But also each symmetry result is illustrated using estimated parameters via the frequency smoothing method (FSM) of spectral correlation function estimation. The time-domain parameters are obtained from the inverse transforms of the FSM parameters. So you can also use this post as an extension of the second-order verification guide to ensure that your estimator works for a wide variety of input parameters.
The Non-Conjugate Cyclic Autocorrelation Function
Let’s start our symmetry odyssey with arguably the simplest CSP function: the non-conjugate cyclic autocorrelation function. Recall that this function is equal to the conventional autocorrelation function when the cycle frequency is zero. The non-conjugate cyclic autocorrelation is defined by
where the notation denotes the infinite time average:
The symmetry in the lag parameter is determined by examining for negative ,
A little algebra yields
Note that (4) also reveals the symmetry in the cycle frequency parameter, since it can be rearranged as
Let’s illustrate our sequence of symmetry results using our old friend, the rectangular-pulse BPSK signal. It has ten samples per bit, or Hz, and a carrier frequency offset of Hz. The carrier frequency offset renders the signal complex valued, and we also add complex-valued white Gaussian noise. The power of the signal is unity and the power of the noise is . The true parameters for this signal are known: frequency domain and time domain.
First, let’s take a look at the autocorrelation, which yields the symmetry relation (4) with ,
This implies that the real part and the magnitude are even functions of and the imaginary part is an odd function of . Figure 1 shows the estimated autocorrelation, confirming these even/odd predictions:
Moving on to the cycle frequency of , we obtain the estimates shown in Figure 2. The result in (4) implies that the real part of the CAF for at should be equal to the real part of the CAF for at (verify with ). It also implies that the imaginary part for and should be the negative of that for and (verify with ). So the equation and the estimates are consistent.
The Conjugate Cyclic Autocorrelation Function
A similar analysis applies to the conjugate cyclic autocorrelation function. Let’s go through it. The function is defined by
Symmetry in follows easily
because we don’t have any conjugations to deal with. This simple result implies that the real and imaginary parts (and the magnitude) are even functions of . This symmetry is validated using estimates in Figure 3.
The symmetry relation (7) tells us nothing about the symmetry in cycle frequency . So let’s look at ,
But (10) is just the conjugated conjugate CAF for the conjugated signal ,
The symmetry relation (12) is validated by estimating the conjugate CAF of for and the conjugate CAF of for and plotting the results. Both functions should be even functions of , and they should be conjugates of each other. This is shown in Figure 4.
The Non-Conjugate Cyclic Cross Correlation Function
Let’s stay with the time-domain parameters and introduce cross functions, starting with the non-conjugate cyclic cross correlation. The basic definition is
Consider the function ,
The relation (17) also contains the symmetry relation for cycle frequency and the symmetry in ‘function order,’ that is, the symmetry relating to switching the order of and in the integrand. The symmetry relation for the cross cyclic CAF is illustrated for in Figure 5, which shows the cross correlation symmetry because . Here is the same rectangular-pulse BPSK signal as we’ve used in the previous illustrations, and is the same signal shifted by samples and subject to a different carrier phase and independent Gaussian noise. We therefore expect the cross correlation to peak at in the XY cross function and at in the YX version.
A more complicated example is shown in Figure 6, which features the cycle frequencies of .
The Conjugate Cyclic Cross Correlation Function
A similar analysis yields the symmetry relations for the conjugate cyclic cross correlation function,
The symmetry in is validated in Figure 7, whereas the symmetry in cycle frequency is validated in Figure 8.
The Non-Conjugate Spectral Correlation Function
The symmetry derivations and illustrations for the cyclic correlation functions were fun, yes, but let’s move on to the important function: spectral correlation. First, our oldest and dearest friend, the non-conjugate spectral correlation function.
The SCF is the temporal correlation between the time-series obtained from two narrowband spectral components of a signal. The narrowband components can be obtained using a simple sliding Fourier transform,
Their idealized correlation is
where we first average over all time, and then let the bandwidth of the narrowband components approach zero. As we go along looking at the symmetry properties of the non-conjugate, conjugate, cross, and conjugate cross SCFs, we’ll need an expression for the narrowband component of the conjugate of ,
Let’s look at the spectral correlation for the conjugated signal,
So the non-conjugate SCF doesn’t possess any symmetry in on its own, but if we include the conjugated-signal SCF, we obtain a sort of symmetry.
Symmetry in cycle frequency is a bit easier,
So the information in the non-conjugate SCFs for negative cycle frequencies is redundant with that for positive cycle frequencies. The two forms of symmetry for the non-conjugate SCF are illustrated in Figures 9 and 10.
The Conjugate Spectral Correlation Function
Using an analysis approach similar to that for the non-conjugate SCF, we obtain the symmetry relations for the conjugate SCF:
The conjugate SCF is even in , and has no cycle-frequency symmetry itself, but possesses a symmetry relation with the conjugate SCF for the conjugated signal. Figures 12-14 illustrate and validate the symmetries.
The Non-Conjugate Cross Spectral Correlation Function
The non-conjugate cross SCF is defined by
Using the same kinds of analysis as above, we obtain a symmetry relation for frequency of
and for cycle frequency of
Equation (36) also provides the symmetry in function order. So there is no symmetry in for the function itself, but there is between the function and the non-conjugate cross SCF for conjugated and . Similarly, there is no symmetry in for the function itself, but there is between (34) and the conjugated function-reversed version of (34) with negated cycle frequency.
The symmetries are illustrated in Figures 15 and 16.
The Conjugate Cross Spectral Correlation Function
If you can bear it, there is one more case. The conjugate cross spectral correlation function is defined by
The symmetry relations are actually relatively easy to derive here. For frequency and function order we have
and for cycle frequency we have
Once again, no simple symmetry involving just the function itself (that is, (37)); all symmetry relations depend on reversing the function order and/or involve the conjugate cross spectral correlation function for conjugated inputs. The relations (38) and (39) are illustrated and validated in Figures 17 and 18.
Symmetries for the Spectral Coherence Function
I can hear you saying, “Yeah, great, but what about the coherence function?” And to be honest, that makes me tired, as this has been a long post. But I’m here for you. The non-conjugate spectral coherence function, greatly useful in blind CSP, is defined by
The denominator of (40) is symmetric in : it is the same when is replaced by . But it is not symmetric in . Recall that the numerator is not symmetric in , but obeys the cycle-frequency symmetry
So we have for the non-conjugate coherence
For the conjugate coherence,
both numerator and denominator are symmetric in , but not in . So we have the symmetry
The cross coherences don’t seem to possess much symmetry. I’ll leave those as an exercise for the reader.
I think the main point for CSP practitioners is that the non-conjugate spectral correlation for is redundant with respect to that for . If you’ve estimated the non-conjugate SCF for , you’ve got all the non-conjugate information you need. On the other hand, you have to estimate the conjugate spectral correlation (or coherence) for the entire range of cycle frequency to ensure you find all the valid cycle frequencies.
Comments, corrections, suggestions, and compliments welcome below!