复件 zuoye2.m

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

%实现X2(n)和X3(n)的DFT
x2n=[1,2,3,4,4,3,2,1];
x3n=[4,3,2,1,1,2,3,4];
M=[3,4];
N=2.^M; 
for i=1:2
    x2k=fft(x2n,N(i));
    subplot(4,2,i);
    stem(0:N(i)-1,abs(x2k),'.');
    grid on
    xlabel('k');
    ylabel('幅度');
    title(['|X2(k)|  N=' num2str(N(i))]);
    subplot(4,2,2+i)
    stem(0:N(i)-1,angle(x2k),'.')
    grid on
    xlabel('k');
    ylabel('相位');
    title(['X2(k)  N=' num2str(N(i))]);
end

for i=1:2
    x3k=fft(x3n,N(i));
    subplot(4,2,4+i);
    stem(0:N(i)-1,abs(x3k),'.');
    grid on
    xlabel('k');
    ylabel('幅度');
    title(['|X3(k)|  N=' num2str(N(i))]);
    subplot(4,2,6+i)
    stem(0:N(i)-1,angle(x3k),'.')
    grid on
    xlabel('k');
    ylabel('相位');
    title(['X3(k)  N=' num2str(N(i))]);
end