zuoye1.asv

来自「序列的dft实现代码对序列快速分析,代码很人性化的」· ASV 代码 · 共 24 行

clear all
clc
close all
%实现R4(n)的DFT
x1n=[1 1 1 1];
M=[3,4,5];
N=2.^M;
for i=1:3
x1k=fft(x1n,N(i));
subplot(3,2,i);
stem(0:N(i)-1,abs(x1k),'.');
grid on
xlabel('k');
ylabel('幅度');
title(['X1(n)  N=' num2str(N(i))]);

subplot(3,2,2+i)
stem(0:N(i)-1,angle(x1k),'.')
grid on
xlabel('k');
ylabel('相位');
title(['X1(n)  N=' num2str(N(i))]);

end