fft_lisanfuliyebianhuan.m

来自「fft离散傅立叶变化」· M 代码 · 共 22 行

%对连续的单一频率周期信号 按采样频率fs=8fa 采样，截取长度N分别选N =20和N =16，观察其DFT结果的幅度谱。 
%解 :此时离散序列x(n)=sin(2*pi*nfa/fs)=sin(2pi*n/8)，即k=8。用MATLAB计算并作图，函数fft用于计算离散傅里叶变换DFT，程序如下
    k=8; 
	n1=[0:1:19]; 
	xa1=sin(2*pi*n1/k); 
	subplot(2,2,1) 
	plot(n1,xa1) 
	xlabel('t/T');ylabel('x(n)'); 
	xk1=fft(xa1);xk1=abs(xk1); 
	subplot(2,2,2) 
	stem(n1,xk1) 
	xlabel('k');ylabel('X(k)'); 
	n2=[0:1:15]; 
	xa2=sin(2*pi*n2/k); 
	subplot(2,2,3) 
	plot(n2,xa2) 
	xlabel('t/T');ylabel('x(n)'); 
	xk2=fft(xa2);xk2=abs(xk2); 
                        subplot(2,2,4) 
	stem(n2,xk2) 
	xlabel('k');ylabel('X(k)');