n161.m

来自「FIR滤波器设计」· M 代码 · 共 21 行

%用所编写的FFT程序
%计算其16点离散傅立叶变换， 并根据DFT的对称性，由     求出                               和                   

%源程序：
x=x2+x3;
xf=kqsfft(x);
xff= abs (xf);
x2k=real(xf);
x3k=imag(xf);
x2kf= abs (x2k);
x3kf= abs (x3k);
subplot(3,1,1);
stem(n,xff)
xlabel('n');ylabel('幅度');title('x=x2+x3谱分析的幅度谱');
subplot(3,1,2);
stem(n,x2kf)
xlabel('n');ylabel('幅度');title('x2k谱分析的幅度谱');
subplot(3,1,3);
stem(n,x3kf)
xlabel('n');ylabel('幅度');title('x3k谱分析的幅度谱');