Sign in

username:

password:



Not a member?

Search compdsp



Search tips

comp.dsp by Keywords

Adaptive Filter | ADPCM | ADSP | ADSP-2181 | Aliasing | AMR | Anti-Aliasing | ARMA | Autocorrelation | AutoCovariance | Beamforming | Bessel | Blackfin | Butterworth | C6713 | CCS | Chebyshev | CIC Filter | Circular Convolution | Code Composer Studio | Comb Filter | Compression | Convolution | Cross Correlation | DCT | Decimation | Deconvolution | Demodulation | DM642 | DSP Boards | DSP/BIOS | DTMF | Echo Cancellation | Equalization | Equalizer | ETSI | EZLITE (Ez-kit Lite) | FFT | FFTW | FIR Filter | Fixed Point | FSK | G.711 | G.723 | G.729 | Gaussian Noise | Goertzel | GPIO | Hilbert Transform | IFFT | IIR Filter | Interpolation | Invariance | JTAG | Kalman | Laplace Transform | Levinson | LPC | McBSP | MIPS | Modulation | MPEG | Multirate | Notch Filter | Nyquist | OFDM | Oversampling | Pink Noise | Pitch | PLL | Polyphase | QAM | QDMA | Quantization | Quantizer | Radar | Random Noise | Reed Solomon | Remez | Resampling | RTDX | Sampling | Sharc | TI C6711 | Undersampling | Viterbi | Wavelets | White Noise | Wiener Filter | Windowing | XDS510PP | Z Transform


Discussion Groups

Free Online Books

See Also

Embedded SystemsFPGAElectronics

Discussion Groups | Comp.DSP | Fractional decimation

There are 3 messages in this thread.

You are currently looking at messages 0 to 3.


Fractional decimation - tharris00 - 2009-03-01 17:41:00

Hello,

I've been studying polyphase decomposition, decimation, interpolation and
fractional decimation. I think I have a handle on what I'm doing, but
there's one thing I don't understand...

Do you have to have greater than an L*Mth order filter in order to use
polyphase decomposition to put the compressor and expander in their most
efficient places (first and last, respectively)?

For example: Let's say I want to reduce my sample rate to Fs'= 0.8333Fs. L
= 5 and M = 6, so they're relatively prime, which is good. 

But then I do Type II decomposition and get R0(z^5), R1(z^5)...R4(z^5).
What if all I needed was a 3rd order FIR filter in the first place? Then
wouldn't R4 =0? 

Whats more, when I go to do the TypeI decompsition, I'm going to get R00,
01, etc, but won't all but the first row of the R filter coefficient matrix
equal zero?

So, either 
(A)it's okay to have zero-terms in the polyphase elements, or 
(B)the minimum order of the original filter, H(z) has to be L*M (per the
example, N=30).

Which is correct?

Thank you for your kind consideration of my question.




______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.

Re: Fractional decimation - John - 2009-03-02 04:25:00



On Mar 1, 5:41=A0pm, "tharris00" <ted.harris.h...@gmail.com> wrote:
> Hello,
>
> I've been studying polyphase decomposition, decimation, interpolation and
> fractional decimation. I think I have a handle on what I'm doing, but
> there's one thing I don't understand...
>
> Do you have to have greater than an L*Mth order filter in order to use
> polyphase decomposition to put the compressor and expander in their most
> efficient places (first and last, respectively)?
>
> For example: Let's say I want to reduce my sample rate to Fs'=3D 0.8333Fs=
. L
> =3D 5 and M =3D 6, so they're relatively prime, which is good.
>
> But then I do Type II decomposition and get R0(z^5), R1(z^5)...R4(z^5).
> What if all I needed was a 3rd order FIR filter in the first place? Then
> wouldn't R4 =3D0?
>
> Whats more, when I go to do the TypeI decompsition, I'm going to get R00,
> 01, etc, but won't all but the first row of the R filter coefficient matr=
ix
> equal zero?
>
> So, either
> (A)it's okay to have zero-terms in the polyphase elements, or
> (B)the minimum order of the original filter, H(z) has to be L*M (per the
> example, N=3D30).
>
> Which is correct?
>
> Thank you for your kind consideration of my question.

I find the polyphase terminology in books to be a huge distraction
from the very simple and obvious concept of avoiding unnecessary
calculations. If you start with a concatenated interpolator and
decimator and elimatate wasted work, you'll end up at the right place.
The only trick is the book-keeping.

John
______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.

Re: Fractional decimation - 2009-03-02 04:47:00

Not 100% sure (would need to think for that :-)). But here is what I
(don't) think.
Since polyphase structure is exactly (theortically and practically)
equivalent to the normal interpolation/filtering/decimation structure
without polyphase, having some terms equal to 0 should be OK if the
filter is very small order.
But, maybe I would really think over it more seriously sometime.

Regards
Piyush


On Mar 2, 3:41=A0am, "tharris00" <ted.harris.h...@gmail.com> wrote:
> Hello,
>
> I've been studying polyphase decomposition, decimation, interpolation and
> fractional decimation. I think I have a handle on what I'm doing, but
> there's one thing I don't understand...
>
> Do you have to have greater than an L*Mth order filter in order to use
> polyphase decomposition to put the compressor and expander in their most
> efficient places (first and last, respectively)?
>
> For example: Let's say I want to reduce my sample rate to Fs'=3D 0.8333Fs=
. L
> =3D 5 and M =3D 6, so they're relatively prime, which is good.
>
> But then I do Type II decomposition and get R0(z^5), R1(z^5)...R4(z^5).
> What if all I needed was a 3rd order FIR filter in the first place? Then
> wouldn't R4 =3D0?
>
> Whats more, when I go to do the TypeI decompsition, I'm going to get R00,
> 01, etc, but won't all but the first row of the R filter coefficient matr=
ix
> equal zero?
>
> So, either
> (A)it's okay to have zero-terms in the polyphase elements, or
> (B)the minimum order of the original filter, H(z) has to be L*M (per the
> example, N=3D30).
>
> Which is correct?
>
> Thank you for your kind consideration of my question.

______________________________
New DSP Code Snippets Section now Live.   Learn more about the reward program for contributors here.