zuoye1.m

序列的dft实现代码对序列快速分析,代码很人性化的

%实现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(2,3,i);
stem(0:N(i)-1,abs(x1k),'.');
grid on
xlabel('k');
ylabel('幅度');
title(['X1(n)  N=' num2str(N(i))]);
subplot(2,3,3+i)
stem(0:N(i)-1,angle(x1k),'.')
grid on
xlabel('k');
ylabel('相位');
title(['X1(n)  N=' num2str(N(i))]);
end