Resample using bandlimited interpolation.
signal to be resampled.
not need to be integers since this routine does not use a polyphase
rate change algorithm, but instead uses bandlimited interpolation,
wherein the continuous time signal is estimated by summing the sinc
functions of the nearest neighbouring points up to distance
Note that resample computes all samples up to but not including time n+1.
If you are increasing the sample rate, this means that it will
generate samples beyond the end of the time range of the original
signal. That is why
xf must go all the way to 10.95 in the example below.
Nowadays, the signal version in Matlab and Octave contain more modern code for resample that has not been ported to the signal R package (yet).
The resampled signal, an array of length
ceiling(length(x) * p / q).
Original Octave version by Paul Kienzle email@example.com. Conversion to R by Tom Short.
J. O. Smith and P. Gossett (1984). A flexible sampling-rate conversion method. In ICASSP-84, Volume II, pp. 19.4.1-19.4.2. New York: IEEE Press.
Octave Forge https://octave.sourceforge.io/
1 2 3 4 5 6 7 8 9 10 11 12 13
xf <- seq(0, 10.95, by=0.05) yf <- sin(2*pi*xf/5) xp <- 0:10 yp <- sin(2*pi*xp/5) r <- resample(yp, xp, xf) title("confirm that the resampled function matches the original") plot(xf, yf, type = "l", col = "blue") lines(xf, r[1:length(xf)], col = "red") points(xp,yp, pch = 19, col = "blue") legend("bottomleft", c("Original", "Resample", "Data"), col = c("blue", "red", "blue"), pch = c(NA, NA, 19), lty = c(1, 1, NA), bty = "n")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.