实数、复数FFT

实数FFT、复数FFT：


FFT是DFT的快速算法，从时域转换到实数FFT和复数FFT中的实数和复数是相对与输入数据而言，从上面两张图中可以看出：

N点的复数FFT需要时域的N点实数和N点虚数（一般设置为0），经FFT转换后，频域得到N点的实数和N点的虚数（因为前后两部分是共轭对称，所以有效数据只有N/2+1个实数和N/2+1个虚数）。

N点的实数FFT只需要时域N点的实数，经FFT转换后，得到频域N/2+1个实数和N/2+1个虚数。

另外N点的复数FFT可以实现2N点的实数FFT：

clear all;
close all;
clc;
% one N-points complex FFT to calculate two N-points real sequences's FFT
N = 1024;
% construct nTap taps FIR filter

h1 = 2*rand(1,N)-1;
h2 = 2*rand(1,N)-1;
H1 = fft(h1);
h4 = zeros(1, 1024);
H4 = fft(h1 + j.*h4);
H2 = fft(h2);
h = h1+j.*h2;
H3 = fft(h);
n = 0:1:N-1;
H3_conj_R = conj(H3(mod(-n,N)+1));
H11 = (H3 + H3_conj_R)/2;
H21 = (H3 - H3_conj_R)/(2*j);

error1 = H11 - H1;
error2 = H21 - H2;
figure(1);
plot(error1);
figure(2);
plot(error2);

实数FFT速度快的原因：
先计算N/2点的复数FFT，再利用FFT变换后频谱的对称性重塑数据从而获得半个FFT频谱即可。