MATLAB中fft

Posted 2020-08-15 JokerMan

   一般的fft需要将点数补成2的整数次幂，MATLAB中有fft函数，输入N点序列则输出也是N点序列，其中N不一定为2的整数次幂。所以会疑惑MATLAB做的fft和N点序列对应的dft是否一样。经过验证，MATLAB中对N点序列做fft的结果与N点dft的结果是完全一样的，没有对N点序列进行补0后再做操作。

  下面是MATLAB代码和结果。

n = 0: 49;
A = 444.128;
a = 50 * sqrt(2.0) * pi;
w0 = 50 * sqrt(2.0) * pi;
T = 0.001;
x = A * exp( -a * n *T) .* sin(w0 * n * T);
k = 0 : 49;
W = 2 * pi / 50;
X = x *(exp(-j * W)).^ (n\' * k);
magX = abs(X);%X的幅度谱
close all;
subplot(311);stem(x);title(\'理想信号采样1\');set(gca,\'box\',\'off\');%采样完成
subplot(312);stem(magX);title(\'理想采样信号的幅频特性1(dft)\');set(gca,\'box\',\'off\');
subplot(313);stem(abs(fft(x)));title(\'理想采样信号的幅频特性1(fft)\');set(gca,\'box\',\'off\');

结果显示：