[SciPy-user] Chirp Z transform
Stéfan van der Walt
Thu Jan 22 09:05:09 CST 2009
2009/1/22 Nadav Horesh <email@example.com>:
> Chirp Z transform is a generalization of the Fourier transform.
> Attached here a module for chirp z transform written by Paul Kienzle and I. We tried to follow scipy's coding-style directions. Is it possible (and how) to make it a part of the scipy project?
Thanks for working on this; I, for one, would like to see it in SciPy.
Recently I referred you to another implementation at
Your version is much more complete, but the following struck me as
data = np.random.random(10000)
a = czt.czt(data, w=np.exp(-2*1j*np.pi/float(len(data)))
b = chirpz_s.chirpz(data, 1, np.exp(-2*1j*np.pi/float(len(data))), len(data))
target = np.fft.fft(data)
err_a = np.sum(np.abs(a - target))
err_b = np.sum(np.abs(b - target))
In : err_a / err_b
The only reason I mention this is because you speak about the
inaccuracy in the docstring. The errors are, on average, in the
vicinity of 1e-10 vs. 5e-11 respectively, so I'm probably on a wild
More information about the SciPy-user