| %计算0≤n≤128时,x(n)的DFT N=128;                  %数据长度
 n=0:N-1;
 fs=20;                   %采样率
 xn=sin(3.1*2*pi*n/fs)+cos(3*2*pi*n/fs);
 Xk=dft(xn);            %x(n)的DFT
 figure;
 subplot(2,1,1),plot(n,xn),title('x(n) (0=<n<128)');
 m=(0:N/2-1)*fs/N;
 subplot(2,1,2),plot(m,abs(Xk(1:N/2)));xlabel('Hz');
 
 
 上面的程序中,为什么要用plot(m,abs(Xk(1:N/2)))而不用plot(n,abs(Xk));程序中的m=(0:N/2-1)*fs/N;是什么意思?
 |