It is often useful to know how a signal processing operation affects the probabilistic parameters of a random signal. For example, if I know the power spectral density (PSD) of some signal , and I filter it using a linear time-invariant transformation with impulse response function , producing the output , then what is the PSD of ? This input-output relationship is well known and quite useful. The relationship is

In (1), the function is the transfer function of the filter, which is the Fourier transform of the impulse-response function .

Because the mathematical models of real-world communication signals can be constructed by subjecting idealized textbook signals to various signal-processing operations, such as filtering, it is of interest to us here at the CSP Blog to know how the spectral correlation function of the output of a signal processor is related to the spectral correlation function for the input. Similarly, we’d like to know such input-output relationships for the cyclic cumulants and the cyclic polyspectra.

Another benefit of knowing these CSP input-output relationships is that they tend to build insight into the meaning of the probabilistic parameters. For example, in the PSD input-output relationship (1), we already know that the transfer function at scales the input frequency component at by the complex number . So it makes sense that the PSD at is scaled by the squared magnitude of . If the filter transfer function is zero at , then the density of averaged power at should vanish too.

So, let’s look at this kind of relationship for CSP parameters. All of these results can be found, usually with more mathematical detail, in My Papers [6, 13].

### Addition of Signals

First let’s consider an old friend, the sum of signals. We’ve encountered this operation before, when we discussed and illustrated signal selectivity, the key property of cyclostationary probabilistic parameters like the spectral correlation function. Suppose we have the sum of statistically independent signals with zero mean values

When we look at the second-order moment of ,

we’ll see a lot of cross terms such as which are zero because of the assumptions of independence and zero mean value. The expectation picks out only the ‘auto’ terms such as . We end up with the sum of autocorrelations

It follows that the cyclic autocorrelation functions are additive too, for the sum of independent signals. Since the spectral correlation function is just the Fourier transform of the cyclic autocorrelation, and the Fourier transform is linear, then the spectral correlation function for the sum of independent signals is the sum of their spectral correlation functions:

These second-order relations are actually special cases of the more general relations involving th-order cumulants. It is known that the cumulant function is cumulative for the sum of independent variables, so that we immediately have

from which it follows that the cyclic cumulants are cumulative too,

And from this result we know that the reduced-dimension cyclic cumulant and the cyclic polyspectra are also additive. These relations form the basis of the signal-selectivity property that is so useful in cyclostationary signal processing.

### Linear Time-Invariant Transformations (Filters)

Consider a linear time-invariant system with impulse-response function and transfer function , where and are a Fourier-transform pair. Such systems are usually referred to as simply *filters*. It is straightforward (but tedious) to use the convolution integral representation of the filter’s input-output characteristic to determine the input-output relations for the cyclic cumulants and cyclic polyspectra. First, let’s define the input and output. The input is and the output is , which are related in the time domain by

or

In the frequency domain, using the convolution theorem, we have the input-output relation

assuming, with good reason, that exists.

The cyclic cumulant input-output relation for is given by

and

### Delay (Special Case of Filtering)

A delay (or advance) can be represented by a linear time-invariant system with an impulsive impulse response function,

Using this in (12) yields the effect on the cyclic cumulant of a delay (or advance, depending on the sign of ):

The lag-shifted cyclic cumulant can be shown to be the unshifted cyclic cumulant multiplied by a phase factor that depends on the delay ,

So a delay induces a phase shift of the cyclic cumulant, but does not affect its magnitude. And the cyclic cumulant is still centered at whatever its center point was in -dimensional lag () space that it had before the delay operation was used.

### Product Modulation (Multiplication by Another Signal)

Next, consider the multiplication of two statistically independent signals,

We’ll look at some examples below. It is easy to show what happens with product modulation to the th-order cyclic moment functions. This is due to the assumption of statistical independent between and . We have

which follows from the properties of the expectation operator and statistical independence. We can substitute the expression relating the temporal moment function to the cyclic temporal moment functions for each of and to obtain

Now, we know we can extract the cyclic moment from the temporal moment by Fourier-series analysis,

,

where the angle brackets denote infinite time averaging. This leads to a formula that shows the mixing of the cyclic features for the two signals,

where must equal a cycle frequency for . So the sum is over all pairs of cycle frequencies for and that sum to .

That’s about as far as I know how to take product modulation without adding some further property to either or . The cyclic cumulant doesn’t have a simple formula for the general case. However, if one of the signals, say is non-random, then we can indeed find a formula for the cyclic cumulants of .

First, notice that for a non-random signal (in our case this means a constant, periodic signal, or polyperiodic signal), the moment function is equal to the lag product itself. That is, consider the th-order lag product

We have

So our general temporal moment formula becomes, for any order ,

When we find the expression for the temporal cumulant for by combining lower-order moments, each lower-order moment will have the form (24) for some order , so that we end up with the cumulant for multiplied by the lag product (moment) for ,

or

so that the cyclic cumulant is given by the mixture of the cyclic cumulants for and the cyclic moments (Fourier components of the lag product in this case) of ,

The conclusion is that when a cyclostationary signal is multiplied by periodic (or polyperiodic) function, the resulting signal has cycle frequencies that are shifted versions of the cycle frequencies for the original signal. For example, if has a cycle frequency of , then will have that cycle frequency only if has a cycle frequency of . That is, only if the th-order lag product for has a non-zero average value. Even if it does have a non-zero average value, that value might be small. Thus, product modulation has serious implications for the output cycle frequencies in terms of those of the input.

A useful special case of product modulation is had by setting the non-random signal equal to a complex constant , where is complex, is real and non-negative, and is real. In this case

For example, if we knew the cyclic cumulants for some signal , then we can easily find the cumulants for $s(t)e^{i\phi_0}$, where could represent an unknown carrier phase.

### Time Gating (Special Case of Product Modulation)

By time gating I mean the signal of interest appears only in regularly spaced time windows or gates. This may happen if a communication system allows the signal to transmit only during certain time slots, or if a receiver follows a schedule for signal detection that allows reception only during brief well-separated regularly spaced intervals. In either case, we might model the signal as the ungated (undisturbed) signal multiplied by a periodic binary-valued gating or windowing function. Mathematically, we have the model

where is the cyclostationary signal of interest and is the binary periodic function defined by

Here the function is a unit-height rectangle, centered at the origin, and with width . So is a pulse train, provided $T_1 > T_2$. Because is periodic with period , so are its lag products. Therefore, the lag products have cycle frequencies (Fourier frequencies) equal to .

We used this kind of gating function in My Papers [34], where we were trying to detect the presence of a cyclostationary signal, but were allowed to receive RF data only during short periodically occurring time windows.

Here is a numerical example illustrating the time-gating effect. We consider a textbook BPSK signal with square-root raised-cosine pulses having roll-off of , a symbol rate of , and a carrier offset frequency of . The gating function is a binary rectangular pulse train, with samples and samples. The blindly estimated (using the SSCA) cycle frequencies and their associated spectral correlation magnitudes for the BPSK signal and its gated version are shown here:

So for the ungated signal, we see the prominent BPSK cycle frequency at , as expected. For the gated signal, we see as well, but also several other cycle frequencies offset from by multiples of . Similarly, the power-related non-conjugate cycle frequency of is scaled and shifted to various harmonics of . You can even see a sinc-function like shape to the variation of the gated-signal spectral correlation magnitudes, as one might expect since the gating signal involves rectangles.

### Frequency Translation (Special Case of Product Modulation)

An important special case of product modulation is frequency translation, where a signal is multiplied by a sine wave in order to shift its support in frequency. For example, a baseband signal is translated to an intermediate frequency by multiplying it by a carrier sine wave.

Mathematically, we have

where can be positive or negative. So in this case and the lag product is

or

The first bracketed term is the only one that depends on time, and is equivalent to

because there are conjugations. The second bracketed term is a phase factor that depends on the frequency and the various delays . So, the multiplication of the signal by a sine wave ends up shifting all of ‘s cycle frequencies by .

Note that when , and the cycle frequencies for are the same as those for . In particular, when , the non-conjugate cycle frequencies for are unchanged from those for . Frequency shifting does not affect non-conjugate cycle frequencies.

Another important signal processing operation is periodic sampling. My colleague Antonio Napolitano has examined that topic in great detail. So if you can’t wait for a future post on it here at the CSP Blog, go check out his work! (Start in The Literature.)