z = x+y*i ;Z = fft(z);X = 0.5*(Z+conj(Z(1+mod(0:-1:1-N,N))));Y = 0.5*(Z-conj(Z(1+mod(0:-1:1-N,N))))/i;end一度人们认为，用离散哈特利转换（Discrete Hartley Transform）来处理纯实数的DFT会更有效率，但接着人们发现，对于同样点数的纯实数DFT，经过设计的FFT，可以比DHT省下更多的运算。Bruun演算法是第一个试着从减少实数输入的DFT运算量的演算法，但此法