I'm implementing a upsample function in Matlab but it's not perfect right now,for reasons I'm not sure. Here is my code:

U=5; %upsampling rate
N = U*length(x);

% Init
x_up = zeros(1,N);

% Zero-pad
for n=1:N
    if mod(n,U)==0
       ind = ind+1;

%% Design a LP filter
tol = (1/U)/2;
[n,fo,mo,w] = firpmord([1/(2*U)-tol 1/(2*U)+tol],[1 0],[0.01,0.01]); 
 b = firpm(n,fo,mo,w);

x_up = [zeros(1,n),x_up]; %To avoid group-delay
x_up = conv(b,x_up,'full');
% Remove leading and tailing zeros
x_up = x_up(find(x_up,1,'first'):find(x_up,1,'last'));    
x_up = x_up(1:N);

And these are the results: enter image description here

enter image description here

enter image description here

enter image description here

From theory perspective, everything looks like fine to me but how can I make my upsampler perfect?

  • $\begingroup$ Is there a reason you are using firpmord() instead of simply using a $sinc()$ filter? $\endgroup$ – havakok Jan 21 '20 at 13:40
  • $\begingroup$ @havakok not actually, i'll try with sinc. $\endgroup$ – kubicwerke Jan 21 '20 at 15:02
  • $\begingroup$ I usually use MATLAB's designfilt tool. I think you'd find it pretty helpful for this $\endgroup$ – Engineer Jan 21 '20 at 18:42

View your frequency response after your low pass filter on a dB scale to better show the limitations of your filter.

Use a multiband filter with the least squares algorithm for an optimized rejection filter for zero-fill interpolation. This will concentrate the rejection to be specifically where the images are that need to be removed.

Given your original signal extends to the band edge, you will have significant challenges interpolating this with zero insert and filtering since the filter complexity is driven by the width of your transition band. I give details below on such an interpolator filter design - notice specifically what defines the signal you want to pass and the signal you want to reject. To simplify that and make this approach feasible, consider first low pass filtering your waveform to reduce the high frequency components while still being acceptable for what you want to have in the time domain. Then you can interpolate using the zero-insert and multi-band filter design approach. The low pass step is not necessary (as the interpolate filter can be designed to do the same thing) but helpful as you will be able to observe your signal prior to further processing to access if the time domain is still acceptable.

To the extent you can pass your signal within the pass band without any distortion, and completely eliminate the images, you will achieve perfect interpolation (this is not possible without an infinitely long filter and infinite delay so you can only approximate this-- given enough taps you can do very well).

This first slide shows how images are created through zero-insert. The digital spectrum is already periodic, zero insert just extends the sampling rate without changing the periodicity or original shape of the spectrum-- so when we interpolate by 4 for example as done here, the images that were originally around $F_s$, $2F_s$ become part of our spectrum in the first Nyquist zone at $F_s/4$ and $F_s/2$. These must be filtered out to complete the interpolation.


An optimized interpolation filter need only reject the images while not distort the passband signal of interest. A least squares (firls in MATLAB/Octave) FIR filter is an excellent choice and support multiband filters where we can concentrate the rejection where it is needed most, providing greater rejection for the same number of taps.



Your Answer

By clicking “Post Your Answer”, you agree to our terms of service, privacy policy and cookie policy

Not the answer you're looking for? Browse other questions tagged or ask your own question.