文档介绍:数字信号处理 实验三:用FFT对信号作频谱分析
%
% 用FFT对信号作频谱分析
clear all;close all
%实验内容(1)===================================================
x1n=[ones(1,4)];    %产生序列向量x1(n)=R4(n)
M=8;xa=1:(M/2);  xb=(M/2):-1:1; x2n=[xa,xb];    %产生长度为8的三角波序列x2(n)
x3n=[xb,xa];
X1k8=fft(x1n,8);        %计算x1n的8点DFT
X1k16=fft(x1n,16);      %计算x1n的16点DFT
X2k8=fft(x2n,8);        %计算x1n的8点DFT
X2k16=fft(x2n,16);        %计算x1n的16点DFT
X3k8=fft(x3n,8);        %计算x1n的8点DFT
X3k16=fft(x3n,16);        %计算x1n的16点DFT
%以下绘制幅频特性曲线
subplot(2,2,1);mstem(X1k8); %绘制8点DFT的幅频特性图
title('(1a) 8点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X1k8))])
subplot(2,2,3);mstem(X1k16); %绘制16点DFT的幅频特性图
title('(1b)16点DFT[x_1(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X1k16))])
figure(2)
subplot(2,2,1);mstem(X2k8); %绘制8点DFT的幅频特性图
title('(2a) 8点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X2k8))])
subplot(2,2,2);mstem(X2k16); %绘制16点DFT的幅频特性图
title('(2b)16点DFT[x_2(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X2k16))])
subplot(2,2,3);mstem(X3k8); %绘制8点DFT的幅频特性图
title('(3a) 8点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X3k8))])
subplot(2,2,4);mstem(X3k16); %绘制16点DFT的幅频特性图
title('(3b)16点DFT[x_3(n)]');xlabel('ω/π');ylabel('幅度');
axis([0,2,0,*max(abs(X3k16))])
 
%实验内容(2) 周期序列谱分析==================================
N=