On Thu, Oct 30, 2008 at 11:44 PM, frank wang <[EMAIL PROTECTED]> wrote:
> Hi, Bob, > > The problem is that I want to resample my data with another sampling rate. > the two rates is very close. I use the formula: > > s(t)=sum(a_k*sinc(t-kTs)). > > the new sampling rate is Ts', so I have > s(nTs')=sum(a_k*sinc(nTs'-kTs)). The sum index k is over the (-P, P), > Centered at n. The n is start from zero. THe code is using two for loops and > it is slow. The length of s(nTs) is very long, so it takes quite long time > to do it. > > I think you can use some form of chirpz/Bluestein for this. Just think of the data as the Fourier transform of its spectrum. Chuck
_______________________________________________ Numpy-discussion mailing list Numpy-discussion@scipy.org http://projects.scipy.org/mailman/listinfo/numpy-discussion