In this post, we continue our study of the symmetries of CSP parameters. The second-order parameters–spectral correlation and cyclic correlation–are covered in detail in the companion post, including the symmetries for ‘auto’ and ‘cross’ versions of those parameters.

Here we tackle the generalizations of cyclic correlation: cyclic temporal moments and cumulants. We’ll deal with the generalization of the spectral correlation function, the cyclic polyspectra, in a subsequent post. It is reasonable to me to focus first on the higher-order temporal parameters, because I consider the temporal parameters to be much more useful in practice than the spectral parameters.

This topic is somewhat harder and more abstract than the second-order topic, but perhaps there are bigger payoffs in algorithm development for exploiting symmetries in higher-order parameters than in second-order parameters because the parameters are multidimensional. So it could be worthwhile to sally forth.

As we did with the second-order parameters, we’ll want to determine the symmetries, if any, relating to the lag (delay) vector in the temporal parameters and the cycle frequency.

The th-order temporal moment and cumulant functions are functions of the variables in the set

where is the th optionally applied conjugation, and by convention we say that of these conjugates are used. This notation is often sufficient for our purposes, but not in this post. Let’s define a *binary indicator vector*

and corresponding *conjugation vector* , where

so that our set of variables can be expressed as

and therefore the specific relationship between the th delay and the conjugation status of is made clear through . Key to our understanding of the symmetries of the cyclic moments and cumulants is the notion of permuting the indicator vector. Let denote a permutation of the index set . For example, for , there are permutations, and these are easily enumerated or obtained in MATLAB through the function perms.m, as shown here:

We want to consider permutations of the indicator function , so we’ll use the notation

to indicate that is a permutation of specified by the reordering implied by . For example, suppose

and

Then

that is, no reordering is implied by this . For

we have

Also, we’ll use a bar to denote an indicator vector that is the complement of another indicator vector

So when , , and when , we have .

### Symmetry in the Delay Variables

If we consider the th-order cyclic temporal moment function,

it seems intuitive that if we permute the delays such that each ends up in a with the same conjugation status that it originally had, then the cyclic moment will be unchanged. For example, consider and , which implies

If we use the permutation , then we obtain

and we can rearrange the factors inside the integral to yield exactly (12). So in this case, the cyclic moment is equal for the two delay vectors and .

Taking a step back, if we consider the set of random variables defined by

and we permute the delays so that the set is the same before and after the permutation, then the time-varying temporal moment functions corresponding to the two sets are identical. But we can also consider independently permuting the indicator vector as well. Permute the delay vector by the permutation and the indicator vector by the permutation . Define the set of variables that is implied by these two permutations:

If , then the temporal moment functions for the two sets are equal, and if the temporal moment functions are equal, so are the cyclic moment functions. * This same argument holds for the temporal cumulant and temporal cyclic cumulant functions*.

We can devise a test for whether or not a particular pair of permutations and result in and therefore in equal cyclic moments and cumulants. If after applying the two permutations, the resulting set of variables in can be reordered so that it is identical to an ordering of the variables in , then we have equality, and therefore have found a symmetry.

Define a matrix by

and a corresponding post-permutation matrix by

If there is a permutation of the columns of such that

then and the moments and cumulants are equal.

I’ve estimated the fourth-order cyclic moments and cumulants for the rectangular-pulse BPSK signal with bit rate , carrier offset , and unit power for all indicator vectors of length and all integer-valued delay vectors on the three-dimensional hypercube . I can use those stored sets of moments and cumulants to verify the symmetries we’ve looked at here, at least for .

In particular, I implemented the / test to find the symmetries, then verified by looking up the corresponding estimates in the stored moment or cumulant files. The results of this verification process do not lend themselves to pretty plots (but do continue reading). I get things like this for the cyclic cumulants:

and similar results for the cyclic moments.

### Symmetry in Cycle Frequency

The symmetries developed in the previous section relate to a single cycle frequency value and consider all the different delay vectors that produce identical values of the temporal parameters.

Here we wonder how the cyclic moment for cycle frequency relates to the cyclic moment for , and similarly wonder about the cyclic cumulants and negated cycle frequencies.

Looking back at the variable set

if we permute the indicator vector such that we obtain the complement

then each variable in that is conjugated becomes unconjugated, and each variable that is unconjugated becomes conjugated, and the delay vector is unaffected. So, in this sense, we would have

From elementary considerations, such as the definition of the th-order moment function, we know that the temporal moment and temporal cumulant for such a are the conjugated functions for . So we expect a cycle-frequency symmetry of the form

Again using the stored sets of cyclic cumulant estimates for the rectangular-pulse BPSK signal, I plotted various cases corresponding to (22) and arranged them in a movie. These all confirm the basic symmetry: negated cycle frequency, complementary indicator function (complementary conjugation configuration), conjugation, and identical delay vector. The individual frames of the movie can also be used to verify correct operation of your cyclic-cumulant estimator. Remember that my cyclic cumulant estimator has been validated by comparing its output to the theoretical cyclic cumulants for those signals for which we know them, such as rectangular-pulse BPSK and several different kinds of QAM and PSK signal types.

### Visualizing the Symmetries

What follows are several attempts at visualizing the cyclic moments and cumulants. There are a lot of dimensions to try to include: order , indicator vector , harmonic number , and delay vector .

First let’s look at a couple cyclic moments and cumulants as surfaces, holding the fourth delay at zero, using two other delays as the x- and y-axes variables, and incrementing the final delay to produce each surface. The movie file names indicate whether the movie shows cyclic moments (‘ctmf’) or cyclic cumulants (‘ctcf’).

Here are the cyclic moment and cumulants for order , , and :

Note how the cyclic moment does not appear to be contained in the range of delays , whereas the cyclic cumulant does. The cyclic moment has the property that in certain directions in the four dimensional delay domain, it does not ever die out (this is intimately related to the notion of impure sine waves).

Here are a couple more in that style:

Here are attempts to simultaneously visualize multiple cyclic-cumulant and cyclic-moment magnitudes and phases:

Finally, here are some movies that show the cyclic cumulant magnitudes and phases for several complementary indicator vectors:

Please don’t hesitate to leave a Comment with corrections, suggestions, or questions. Especially corrections.