Tuesday, April 25, 2017

filter bank - Understanding multirate filtering fundamentals


I am having trouble grasping some of the fundamental concepts of multirate filtering. I see from various sources that the basic building blocks of a multirate filter are the dyadic analysis and synthesis blocks.




  • Question 1:


    The analysis block structure looks like the following, where the wide band signal is split into lowpass and highpass bands, each with a cutoff of FS/4 (Nyquist/2). Each band is then decimated by a factor of 2.


    enter image description here


    How can you accurately represent the signal in the high frequency band when it contains frequency information above the Nyquist limit of the new decimated sample rate?





  • Question 2:


    The analysis block structure looks like the following, where the sub-band signal is interpolated, re-filtered and then summed.


    enter image description here


    What is the purpose of the second filtering?





Answer



I will answer question 2 first, and hopefully that will help explain what is going on with question 1.


When you sample a baseband signal there are implicit aliases of the baseband signal at all integer multiples of the sampling frequency, as shown in the picture below. Assymetric baseband signal w/ aliases The solid image is the original baseband signal, and the aliases are represented by the dashed images. I chose an assymetric (i.e. complex) signal to help demonstrate the inversion that happens at odd multiples of the sampling frequency.


You might ask, "Do the aliases really exist?" It's a bit of a philosophical question. Yes, in a mathematical sense they do exist, because all of the aliases (including the baseband signal) are indistinguishable from each other.



When you upsample by inserting zeros in between the original samples, you are effectively increasing the sampling rate by the upsample rate. So if you upsample by a factor of two (putting one zero in between each sample), you increase your sampling rate and the Nyquist rate by a factor of 2, resulting in the picture below. Upsampled assymetric baseband signal


As you can see, one of the implicit aliases in the earlier image has now become explicit. If you FFT the samples it will show up. A non-rigorous proof that the DFT transform does not fundamentally change is given below.


Now that you have the two explicit aliases, if you just want the baseband alias then you have to low-pass filter to get rid of the other alias. Sometimes, though, people use the other aliases to do their modulating for them. In that case you would high-pass filter to get rid of the baseband signal. I hope that answers question 2.


Question 1 is basically the inverse of question 2. Suppose that you are already in the situation shown in the second picture. There are two ways to get the baseband signal that you want. The first way is to low-pass filter (thereby getting rid of the higher alias) and then decimating by a factor of two. That gets you to picture #1.


The second way is to high-pass filter (getting rid of the baseband alias) and then decimating by a factor of two. The reason that this works is that you are intentionally aliasing the signal into the baseband, thus, once again, getting you to picture #1.


Why would you want to do it that way? Because in most situations the signals will not be the same, so you can choose which signal you want, or do them both separately.


If you are studying multi-rate processing, I highly recommend getting "Multirate Signal Processing for Communication Systems" by Frederic Harris. He does a really good job of explaining the theory without neglecting the math, and giving a lot of practical advice too.


EDIT: Intentionally sampling a signal at less than the Nyquist rate is called undersampling. The following is my attempt at explaining mathematically why the FFT does not change when you upsample. "x[n]" is the original set of samples, "u" is the upsampling factor, and "x'[n]" is the upsampled set of samples.


$$ \begin{eqnarray*} X[k] &=& \sum_{n=0}^{N-1} x[n]e^{-i2\pi kn/N}\\ X'[k] &=& \sum_{n=0}^{uN-1} x'[n]e^{-i2\pi kn/uN} , \{ &x&'[n] = 0, n \neq mu, m \in (0..N-1)\\ &x&'[n] = x[n/u], n = mu\\ &=& \sum_{n=0}^{N-1} x'[un]e^{-i2\pi kun/uN}\\ &=& \sum_{n=0}^{N-1} x[n]e^{-i2\pi kn/N}\\ &=&X[k] \end{eqnarray*} $$


Apologies for the ugly formatting. I am a LaTex noob.



EDIT 2: I should have pointed out that the DFT's of x[n] and x'[n] are not truly identical. The sample rate is higher, which as I explained in the earlier part of the answer, causes aliases to be "exposed". I was trying to point out in my non-mathematician way that the DFT's are, aside from the sample rate, the same.


No comments:

Post a Comment

periodic trends - Comparing radii in lithium, beryllium, magnesium, aluminium and sodium ions

Apparently the of last four, $\ce{Mg^2+}$ is closest in radius to $\ce{Li+}$. Is this true, and if so, why would a whole larger shell ($\ce{...