compare.m

来自「OFDM系统」· M 代码 · 共 31 行

%这一部分是对于我自己编的FFT函数和MATLAB中自带的函数的计算精度的比较
%利用的只是简单的余弦函数cos(pi*n/4)

clear
n = 0:15;
x=cos(pi*n/4);
y = FFT_self(x,16);
z = fft(x);
if(z == y)
    disp('true');
else
    disp('false');
end

figure(1)
stem(n,imag(z));
hold on;
h = stem(n,imag(y));
set(h,'MarkerFaceColor','red');
title('虚部比较');
ylabel('虚部');
xlabel('n');

figure(2)
stem(n,real(z));
hold on;
h = stem(n,real(y));
set(h,'MarkerFaceColor','red');
title('实部比较');
ylabel('实部');
xlabel('n');