# [SciPy-User] [SciPy-user] sinc interpolation

bdb112 boyd.blackwell@anu.edu...
Thu Jun 7 04:05:18 CDT 2012

```Corrected typo and improved code:

typo:  should be     yy = resample(y,100)

Improved code:
from scipy.signal import resample
from numpy import linspace, sin, pi
from pylab import plot, show

# replace 4*pi to see effect of non-periodic function (e.g. at either end)
x = linspace(0,4*pi,10,endpoint=False)
y = sin(x)
# if you supply t, you get the interpolated t back as well.
(yy,xx) = resample(y, 100, t = x)

plot(x,y,'ro', xx, yy)
show()

Travis Oliphant-5 wrote:
>
>
>
> On Nov 20, 2009, at 2:26 PM, David Trem wrote:
>
>> Hello,
>>
>> Is sinc interpolation available in Scipy ?
>
> Yes, use scipy.signal.resample  which uses a Fourier method to
> downsample or upsample a signal:
>
> from scipy.signal import resample
> from numpy import r_, sin
> from pylab import plot
>
> x = r_[0:10]
> y = sin(x)
> yy = resample(x, 100)
>
> # This is a bit tricky to get the x-samples right
> xx = r_[0:10:101j][:-1]
>
> plot(x,y,'ro', xx, yy)
>
>
> -Travis
>
> _______________________________________________
> SciPy-User mailing list
> SciPy-User@scipy.org
> http://mail.scipy.org/mailman/listinfo/scipy-user
>
>

--
View this message in context: http://old.nabble.com/sinc-interpolation-tp26449760p33974580.html
Sent from the Scipy-User mailing list archive at Nabble.com.

```