实验3用FFT对信号作频谱分析(预习报告).doc

实验3 用FFT对信号作频谱分析 1．实验程序及运行结果 实验内容1：对非周期序列进行谱分析 对以下序列进行谱分析 　　　　　　 选择FFT的变换区间N为8和16 两种情况进行频谱分析。分别打印其幅频特性曲线， 并进行对比、分析和讨论。 源程序：shzxhchlshiyan3_1 %用FFT对非周期序列进行谱分析 clear all,close all,clc,clf; 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];%产生长度为8的序列x3(n) X1k8=fft(x1n,8); %计算x1n的8点DFT X1k16=fft(x1n,16); %计算x1n的16点DFT X2k8=fft(x2n,8); %计算x2n的8点DFT X2k16=fft(x2n,16); %计算x2n的16点DFT X3k8=fft(x3n,8); %计算x3n的8点DFT X3k16=fft(x3n,16); %计算x3n的16点DFT Xk=fft(x1n,1024); %1024点FFT[x(n)], 用于近似序列x(n)的TF k=0:1023;wk=k/1024; %0-2π周期内的归一化频率 subplot(3,1,1);plot(wk,abs(Xk));%画FT[x1(n)]a axis([0,1,0,1.1*max(Xk)]);set(gca,Xtick,[0,0.25,0.5,0.75,1]); title((a) 原序列的频谱FT[x_1(n)]); ylabel(|X(e^j^\omega)|);%xlabel(f=\omega/2\pi); nX1k8=0:7;nX1k16=0:15; subplot(3,1,2);stem(nX1k8,abs(X1k8),.); %绘制8点DFT的幅频特性图 axis([0,8,0,1.1*max(X1k8)]); title((b) 8点DFT[x_1(n)]);ylabel(幅度);%xlabel(n); subplot(3,1,3); stem(nX1k16,abs(X1k16),.); %绘制16点DFT的幅频特性图 axis([0,16,0,1.1*max(X1k16)]); title((c) 16点DFT[x_1(n)]);ylabel(幅度);%xlabel(n); figure(2) Xk=fft(x2n,1024); %1024点FFT[x(n)], 用于近似序列x(n)的TF k=0:1023;wk=k/1024; %0-2π周期内的归一化频率 subplot(3,1,1);plot(wk,abs(Xk));%画FT[x2(n)] axis([0,1,0,1.1*max(Xk)]);set(gca,Xtick,[0,0.25,0.5,0.75,1]); title((a) 原序列的频谱FT[x_2(n)]); ylabel(|X(e^j^\omega)|);%xlabel(f=\omega/2\pi); nX2k8=0:7;nX2k16=0:15; subplot(3,1,2);stem(nX2k8,abs(X2k8),.); %绘制8点DFT的幅频特性图 axis([0,8,0,1.1*max(X2k8)]); title((b) 8点DFT[x_2(n)]);ylabel(幅度);%xlabel(n); subplot(3,1,3); stem(nX2k16,abs(X2k16),.); %绘制16点DFT的幅频特性图 axis([0,16,0,1.1*max(X2k16)]); title((c) 16点DFT[x_2(n)]);ylabel(幅度);%xlabel(n); figure(3) Xk=fft(x3n,1024); %1024点FFT[x(n)], 用于近似序列x(n)的TF k=0:1023;wk=k/1024; %0-2π周期内的归一化频率 subplot(3,1,1);plot(wk,abs(Xk));%画FT[x3(n)] axis([0,1,0,1.1*max(Xk)]);set(gca,Xtick,[0,0.25,0.5,0.75,1]); title((a) 原序列的频谱FT[x_3(n)]); ylabel(|X(e^j^\omega)|);%xlabel(f=\omega/2\pi); nX3k8=0:7;nX3k16