How to resample a signal by a fraction

Hello,

I have the time-domain received signal v(t) and I want to re-sample it at as v(t/(1+a)), where a is not an integer and it is in the order of ~10^-3. How can I do that in MATLAB, since the resample function requires that both P and Q be integer?

Thanks

Answer 1

[Matt Kindig]

You could use interpolation instead. Something like this:

vnew = interp1( t, v, t/(1+a), 'linear');

You could use other interpolation approaches ('spline','cubic', etc.) as well, as appropriate.

[Saed]

This seems to work just fine.

Thanks

Answer 2

[David Barry]

I suggest you reciprocate your a variable and round it to the nearest integer. You then just need to make sure you get the arguments the correct way in the resample function.

y = resample(x,p,q) resamples the sequence in vector x at p/q times the original sampling rate

[Saed]

This is the problem, I do not want to round it to the nearest integer, since it always be 1. I want the resampling factor 1/(1+a) as is.

[David Barry]

Apologies, I missed your 1+a bit. In that case you need to change both P and Q in the resample function by multiplying up until you get the accuracy you need.

Example: Suppose I want to re-sample by 1/(1+a) where a was 1e-3.

1+1e-3 = 1.0010 so I might want to scale up by factor of 1000 and use values of 1000 and 1001 for P and Q.

[Saed]

What do you mean by scale up by 1000? Like

if true
  resample(x,1000*(1+a),1001);
end
?

[David Barry]

I mean having P and Q values of 1 and 1.0010 is the same as having P and Q values of 1000 and 1001.

[Kir]

P and Q are integers

Answer 3

[Malcolm Lidierth]

Is this any help?

% SincResample returns the data convolved with a set of time-shifted windowed sinc functions, one for each of the samples [1..size(x,1)] in the input signal.

http://www.mathworks.co.uk/matlabcentral/fileexchange/12268-sinc-resample