I'm trying to fully understand the spectral effects in a polyphase channelizer / polyphase analysis filter bank. I don't really have a signal processing background, so please excuse some incorrect terminology ;-).
So let's consider a very simple case with two output channels. I'll be referring to Polyphase Channelizer Demystified [Lecture Notes] in this question. Let's start with this block diagram:
and this input spectrum:
I want to understand what spectra at A1
, A2
, B1
and B2
look like. So obviously as the input signal is sampled we have a repeating spectrum somewhat depicted in the graph above (showing one repeated half of the grey channel on each side) the paper explains downsampling by factor 2 will scale the signals by factor 1/2 and linearly add one shifted version.
At A1 we have such a simple down-sampling operation, so we should get this:
the signal is downsampled by factor 2, so we have a new sample frequency. Still, the linearly combined signal is plausible, as the spectrum repeats according to the new output sample rate.
Equation (6) and (7) in the paper describe the spectrum of a downsampled signal with one unit delay. This should be the signal B1 and it should look like this:
However, this is not plausible, as the signal spectrum does not repeat properly according to the new sample frequency. What did I miss here?
The next question is: What is the effect of the downsampled filter phases? The filter-downsampling will lead to a 'repeated spectrum' filter transfer function as well, correct? So is this correct for A2
?:
And what's the effect of the offset/shift of the filter coefficients in p1
caused by the polyphase partition?
Now if we only add up the Nyquist band of the A1
and B1
pictures we end up with the orange spectrum representing channel 0. If we do the same but with a phase shift A1
and -B1
gives the grey spectrum, channel 1. So the end result is plausible, but the intermediate signal at B1 is not.
EDIT: To add some equations to this topic: According to R. E. Crochiere and L. R. Rabiner. Multirate Digital Signal Processing, pp. 84–88
Filter partitions: $p_\rho(n) = h(nM + \rho)$
Equation for polyphase branches: $x_\rho(n) = x(nM - \rho)$