%计算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;是什么意思? |