# Tag Info

## Hot answers tagged multirate

9

It's quite benefical to view the mathematical expression of the process as well. In discrete-time interpolation, a mere $N$ times repetition of a single sample (pixel) into the new enlarged block is an interpolation with zero-order hold kernel. This is mostly applied for simple fast but crude image zooming. In 1D, the kernel is $$h_0[n] = \begin{cases}{ 1 ~... 5 For example a delay of 0.75 sampling periods would have an impulse response like this (red squares sampled from the blue delayed sinc): Time (horizontal axis) unit is sampling periods. That kind of a filter is not causal as there are non-zero samples at negative times. None of the samples is zero, no matter how far left or right you go. That's what they ... 5 Assuming you have completed the basic course in DSP with a book like Oppenheim's Discrete-Time Signal Processing book then the next step could be anyone of the followings: Multiresolution Signal Decomposition: Transforms, Subbands,Wavelets, A.Akansu Multirate Digital Signal Processing, Rabiner Multirate Systems and Filter Banks, P.P.Vaidyanathan ... 4 The problem with downsampling is that it can be lossy -- since you're reducing the sampling rate, you can introduce aliasing. So, you can reverse the order whenever downsampling does not result in aliasing. For example, say your discrete-time signal x[n] contains energy in frequencies up to f_N/3, where f_N is the Nyquist frequency. Then, downsampling ... 3 There are two ways (I see right now) to model this system in terms of other known systems. First that comes to my mind is a zero-order sample and hold system. Although those are usually described in a continuous time domain, it might be interesting to look at as well. Another model is a combinition of a usual upsampler followed by an FIR with all N taps ... 3 For one practical example, I'll point to GNURadio's frequency translating FIR filter block: https://github.com/gnuradio/gnuradio/blob/master/gr-filter/lib/freq_xlating_fir_filter_XXX_impl.cc.t It's a channelizer & decimator block that gains an efficiency by spinning a user provided FIR LPF up to be a complex BPF where the channel is and filtering first ... 3 This is a problem from Multirate Systems and Filter Banks by P. P. Vaidyanathan. In Section 4.6.2.E, there is a discussion on what he calls Euclidean Complementary Functions. Specifically, he discusses one of the results of Euclid's theorem: ... if H_0(z) and H_1(z) are relatively prime, there exists polynomials F_0(z) and F_1(z) such that$$H_0(...

3

Answer: You will see residual images of $X(f)$ at multiples $f_s$, $2f_s$ and $3f_s$, and distorted image of $X(f)$ at non-zero multiples of $4f_s$, when sampling in the manner you explained. Depending on value $e$, the size of residual will change. I have explained how in detail below. Ideally, sampling at $4f_s$ would have completely cancelled those ...

2

If you are running a causal/real-time process, note that any upsampling/downsampling sub-process involves anti-alias filtering that will add a delay (a linear delay in the case of a symmetric FIR filter). So, if you want to re-combine resampled and non-resampled processing chains, you might want to add delays to the non-resampled processing paths to match ...

2

You are on the right track. You can use decimation and interpolation (or resampling, which is a combination of both decimation and interpolation) to change the sample rate. For instance, if you start with 48 kHz samples you can decimate by 6 to get down to 8 kHz, do some processing, then interpolate by 2 to get to 16 kHz, do some processing, then ...

2

The ideal delay has a frequency response of: $$H(e^{j\omega}) = e^{-j\omega D}$$ this has impulse response $$h(n) = \mbox{sinc}(n - D) = \frac{\sin(\pi(n-D)}{\pi(n-D)}.$$ For $D$ an integer, this just becomes: $$h(n) = \delta(n-D)$$ where $\delta$ is the Kronecker delta. That means there is no resampling, so I am not sure where you get that from (it's ...

2

The Sampling operation (both upsample and downsample) depends on two very critical conditions: 1- The existance and applicability of ideal frequency selective filters 2- The operated signal being strictly bandlimited For most practical systems these two critical conditions are only approximately met. Hence the computational results by those practical ...

2

I think you are referring to the symbol shown at the left of the structure: the arm that moves between different phases? See slide 4 of this PDF. The point the Polyphase representation is trying to get across is that we give different samples of the input to different Phases of the Polyphase filter. Polyphase Filtering is a Digital Signal Processing ...

2

In addition to what hotpaw2 said, it's important to understand that sample rate reduction is the corner stone of actual applicability in many systems. Software Defined Radio frontends nowadays produce hundreds of complex-valued Megasamples per second – if you don't reduce that amount of data, no commonly available PC can deal with that. Also, many devices ...

2

I will use Wikipedia notations - Kalman Filter. In most models the state transition model matrix $F$ depends on the interval parameter $T$. The same goes for the Process Noise Covarinace Matrix $Q$. For instance, for the constant velocity model: $$F = \begin{bmatrix} 1 & T \\ 0 & 1 \end{bmatrix}, \; Q = \begin{bmatrix} \frac{ {T}^4 }{4} & ... 2 You would upsample by 441 and downsample by 640 since 11025 \cdot 640 = 16000 \cdot 441 That would be a rather expensive operation, so in practice you would either fugde it a bit: 31/45 would get you to 11022Hz and 11/16 to 11000 Hz or you would use a suitable poly-phase irrational conversion algorithm 2 [Good question, that made me rethink of stuff I believed natural. I shall incorporate them in future lectures] Downsampling the highpass (and the lowpass) provide you with a critically sampled filterbank, in other words, not redundant. Critical sampling is a price to pay in some applications like compression, and in the context of finding an orthogonal ... 2 As you know you have to use an ideal low-pass filter, and zero order hold is an approximation which needs more compensation. Ideal LPF has a flat response with sudden cut-off and zero stop band but your moving average window has a sinc frequency response. So the deep valleys belong to the zeros of sinc function and the remaining which has high attenuation ... 2 The decimator by integer M can be shown to be the following block:$$ x[n] \longrightarrow \boxed{ \downarrow M } \longrightarrow y[n] = x[Mn] $$Assuming that the input x[n] is WSS it has the ACS as$$ r_x[k] = E\{ x[n] x^*[n+k] \} $$Then the output autocorrelation can be defined as:$$ r_y[n,n+k] = E\{ y[n] y^*[n+k] \} = E\{ x[Mn] x^*[Mn+Mk] \} =...

2

Linear and more generally polynomials are pretty common methods for interpolation or extrapolation, easy to implement, and simple in a causal setting (prediction). Parabolic interpolation (in the Fourier domain) or Savitzky-Golay filtering are practically useful examples, especially for peak-like signals. Polynomials however can be poor at modeling ...

1

With every polyphase filter bank I have worked with, the first block in the analysis phase is an IFFT, and the block in the synthesis phase is a DFT. These operations essentially cancel one other, so it should be fairly intuitive. Think of the DFTs acting as sinusoids modulators in this case rather than an operation to convert to the frequency domain. By ...

1

just a WSS signal $x[n]$ passed through an $M$-fold decimator (a multirate signal processing block where $M=2$, used in fractional rate conversion), no filter follows the decimator. So $y[n] = x[Mn]$ WSS means that the autocorrelation function only depends on the distance between the points (and that the mean and variance are time-independent); i.e. the ...

1

@Royi's and user28715's answers are correct. So just add this answer to theirs. If you let the system's output matrix $H$ go to zero, then two things happen. First, you are modeling a time increment where you know nothing about what actually happened. Second, the Kalman gain ($K$ in most representations of the Kalman filter) goes to zero, and the entire ...

1

Put band pass filters (5 MHz wide) … Re-sample the chunks to 4 MHz. Don't do that! If you reduce the sampling rate to 4 MS/s, you need to filter to 4 MHz bandwidth anyway. So you could instead just use 4 MHz wide filters and get rid of the resampling filter.

1

If the only reason for the buffer is so your algorithm has something to process, then the simplest way to implement this is to keep your current algorithm but use a shorter buffer. If the larger buffer was present for other reasons (i.e. data arrives in large chunks), then keep that larger buffer, but feed it into a shorter buffer that is given to your ...

1

25.6 MHz is larger than the largest possible sampling rate of these two USRP models (25MHz). Therefore, this is impossible to implement directly; in any way, you'll need to generate your 25.6 MS/s signal first, and then resample to 25 MS/s. That resampler will be very CPU-consuming, probably more than your whole OFDM signal processing. The underruns happen ...

1

I will try to explain from my understanding. For interpolation, the whole process can be viewed as two parts: up-samping with padding zeros, and the low-pass filtering. The zero padding will increase the length and thus the computational costs. For the low-pass filtering, there are lots of zeros, thus these multiplications could be omitted. For decimation,...

1

So the central question here is what the restrictions for the filters of this class of filterbanks (i.e. $N=2$-filter bank with the HPF generated from a lowpass/highpass transform of the LPF) is. Let us for that first introduce a couple of helping definitions: \begin{align} \mathbf{h}(z) &= \begin{bmatrix}h_0(z)\\h_1(z)\end{bmatrix}&\text{...

1

"Multirate" has multiple acceptations. The main ideas is that data is composed of different components that each different dynamics or rates (your two sines for instance). The standard approach is "he who can do more can do less": one adapts to the fastest rate. The multirate approach is a divide and conquer process: somehow, you can use different (and ...

Only top voted, non community-wiki answers of a minimum length are eligible