用FFT对信号作频谱分析Matlab程序.docVIP

对以下序列进行FFT分析 x1(n)=R4(n) x2(n)= x3(n)= x1n=[ones(1,4)]; %产生R4(n)序列向量 X1k8=fft(x1n,8); %计算x1n的8点DFT X1k16=fft(x1n,16); %计算x1n的16点DFT %以下绘制幅频特性曲线 N=8; f=2/N*(0:N-1); （不懂） figure(1); subplot(1,2,1);stem(f,abs(X1k8),r,.); %绘制8点DFT的幅频特性图，abs求得Fourier变换后的振幅 title((1a) 8点DFT[x_1(n)]);xlabel(ω/π);ylabel(幅度); N=16; f=2/N*(0:N-1); subplot(1,2,2);stem(f,abs(X1k16),.); %绘制8点DFT的幅频特性图 title((1b) 16点DFT[x_1(n)]);xlabel(ω/π);ylabel(幅度); %x2n 和 x3n M=8;xa=1:(M/2); xb=(M/2):-1:1; %从M/2到1每次递减1 x2n=[xa,xb]; %产生长度为8的三角波序列x2(n) x3n=[xb,xa]; X2k8=fft(