[SciPy-user] Chirp Z transform
Stéfan van der Walt
stefan at sun.ac.za
Thu Jan 22 10:05:09 EST 2009
Hi Nadav
2009/1/22 Nadav Horesh <nadavh at visionsense.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
http://www.mail-archive.com/numpy-discussion@scipy.org/msg01812.html
Your version is much more complete, but the following struck me as
slightly strange:
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 [152]: err_a / err_b
Out[152]: 1.6138562461610748
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
goose chase.
Regards
Stéfan
More information about the SciPy-User
mailing list