matlab中fft定点运算,DSP中的定点FFT运算注意事项(转载)

引用：http://publishblog.blogchina.com/blog/tb.b?diaryID=3811777

在DSP运算中，经常需要把输入时域信号在频域进行处理之后，再还原为时域信号，这样就需要进行FFT和IFFT运算：

x(n) -> FFT -> X(f) -> 频域处理 -> Y(f) -> IFFT -> y(n)

而一般的DSP芯片只支持整数运算，也就是说只能进行定点小数计算。

N点FFT计算出0… N-1，N个复数： 0，A，N/2，A*，A为(N/2-1)个复数，A*为A的共轭复数。

FFT的公式为：

N

X(k) =       sum  x(n)*exp(-j*2*pi*(k-1)*(n-1)/N)、1 < =  k < =  N.

n = 1

IFFT的公式为：

N

x(n) = (1/N) sum  X(k)*exp( j*2*pi*(k-1)*(n-1)/N)、1 < =  n < =  N.

k = 1

假设我们对ADC转换器转换的数字信号进行FFT运算，若输入数据为16bit的短整型数，我们可以把它看作Q15的从-1到1之间的小数。根据FFT的公式我们可以知道，FFT变换之后的结果将超出这个范围。

例如在matlab中输入fft(sin([1:8]*0.5))，可以看到结果：

2.8597，-0.8019 - 3.0216i，0.4312 - 0.8301i