Saya ingin menerapkan Fast Cosine Transform. Saya membaca di wikipedia , bahwa ada versi cepat dari DCT yang juga dihitung untuk FFT. Saya mencoba membaca makalah yang dikutip Makhoul *, untuk implementasi FTPACK dan FFTW yang juga digunakan dalam Scipy , tetapi saya tidak dapat mengekstrak algoritma yang sebenarnya. Inilah yang saya miliki sejauh ini:
Kode FFT:
def fft(x):
if x.size ==1:
return x
N = x.size
x0 = my_fft(x[0:N:2])
x1 = my_fft(x[0+1:N:2])
k = numpy.arange(N/2)
e = numpy.exp(-2j*numpy.pi*k/N)
l = x0 + x1 * e
r = x0 - x1 * e
return numpy.hstack([l,r])
Kode DCT:
def dct(x):
k = 0
N = x.size
xk = numpy.zeros(N)
for k in range(N):
for n in range(N):
xn = x[n]
xk[k] += xn*numpy.cos(numpy.pi/N*(n+1/2.0)*k)
return xk
Uji coba FCT:
def my_fct(x):
if x.size ==1:
return x
N = x.size
x0 = my_fct(x[0:N:2]) # have to be set to zero?
x1 = my_fct(x[0+1:N:2])
k = numpy.arange(N/2)
n = # ???
c = numpy.cos(numpy.pi/N*(n+1/2.0)*k)
l = x0 #???
r = x0 #???
return numpy.hstack([l,r])
* J. Makhoul, "Transformasi kosinus cepat dalam satu dan dua dimensi," IEEE Trans. Acoust. Pidato Sig. Proc 28 (1), 27-34 (1980).