Signal Processing Stack Exchange is a question and answer site for practitioners of the art and science of signal, image and video processing. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

Looking for a solution to the following problem:

A signal $x(t)$ is band limited to $B$ Hz, and sampled above the Nyquist rate, with corresponding $f_s = 1/T$. If the sampled signal is given by

$$ y(n) = x(nT), $$

describe how to compute the phase-shifted sampled signal

$$ z(n) = x(nT + T/2), $$

from the samples $y(n)$ without first reconstructing $x(t)$ from $y(n)$ and resampling.

Have not been able to find one on the web.

share|improve this question

Because the sampling frequency is above the Nyquist rate, the original signal can be written in terms of its samples $y(n)$:


By setting $t=nT+T/2$ you get from (1)


and since

$$\sin[\pi(n-m)+\pi/2]=(-1)^{n-m}$$ you finally get


So for computing $z(n)$ you do not need to reconstruct $x(t)$ but you need to compute a sum over all values $y(m)$. In practice a few samples $y(m)$ around the value $n$ will be sufficient because the weighting factors in the sum decrease with increasing distance $|n-m|$.

share|improve this answer

Yet another possibility is to zero-pad the signal sample vector, FFT it, rotate the phase of each FFT result bin linearly with the bin index, and IFFT a time shifted result.

share|improve this answer

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

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