实验三：用FFT对信号作频谱分析.docxVIP

实验三:用FFT对信号作频谱分析实验报告 -、 实验目的与要求 学习用FFT对连续信号和时域离散信号进行谱分析的方法， 了解可能出现的 分析误差及其原因，以便正确应用 FFT。 二、 实验原理 用FFT对信号作频分析是学习数字信号处理的重要内容， 经常需要进行分析 的信号是模拟信号的时域离散信号。 对信号进行谱分析的重要问题是频谱分辨率 D和分析误差。频谱分辨率直接和 FFT的变换区问N有关，因为FFT能够实现的 频率分辨率是2兀/N,因此要求2兀/N小丁等丁 D。可以根据此式选择FFT的变 换区问N。误差主要来自丁用FFT作频谱分析时，得到的是离散谱，而信号（周 期信号除外）是连续谱，只有当 N较大时，离散谱的包络才能逼近连续谱，因此 N要适当选择大一些。 三、 实验步骤及内容（含结果分析） （1）对以下序列进行FFT分析： xi(n)=R4(n) 「n+1 0 n 3 X2（n）= \ 8-n 4 n 7 0 其它n 4-n 0 n 3 X3（n）= n-3 4n 7 0 其它n 选择FFT的变换区间N为8和16两种情况进行频谱分析，分别打印出幅频特性曲线, 并进行讨论、分析与比较。 【实验结果如下】： IOPa）场点 IO Pa）场点D叩咐 201 15 ? W w 哩 CJ/TT 印）日点DFT[电㈣辨 ■ ■ 15 10 5 0 0.5 1 1.5 2 实验结果图形与理论分析相符 对以下周期序列进行谱分析: X4（n）=cos[（兀 /4）*n] X5（n）= cos[（兀 /4）*n]+ cos[（兀 /8）*n] 选择FFT的变换区间N为8和16两种情况进行频谱分析，分别打印出幅频特性曲 线，并进行讨论、分析与比较。 tto(4a) 3 tto (4a) 3点 DFT(n)] 0 ?——?—*—?—*—I— TOC \o 1-5 \h \z 0 0.6 1 1.5 2 cjAt 真司B点DFTK(n)] G,— ■—~ D C 0.5 1 1.5 2 w/tt （3）对模拟周期信号进行频谱分析: x 6（n）= cos（8 兀 t）+ cos（16 兀 t）+ cos（20 兀 t） 选择采样频率Fs=64H乙FFT的变换区间N为16、32、64三种情况进行频谱分析, 分别打印出幅频特性曲线，并进行讨论、分析与比较。 【实验结果如下】： 四、 【附录】（实验中代码） 产生R4（n）序列向量 产生R4（n）序列向量 计算x1n的8点DFT 计算x1n的16点 DFTX1k8=fft(x1n,8); % DFT X1k16=fft(x1n,16); % %以下绘制幅频特性曲线 N=8; f=2/N*(0:N-1); 绘制8点 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; subplot(1,2,1);stem(f,abs(X1k8),.); % title((1a) 8 点 DFT[x_1(n)]);xlabel( N=16; 绘制8 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; subplot(1,2,2);stem(f,abs(X1k16),.); % title((1a) 16 点 DFT[x_1(n)]);xlabel( %x2n 和 x3n M=8;xa=1:(M/2); xb=(M/2):-1:1; x2n=[xa,xb]; % 产生长度为8的三角波序列x2(n) x3n=[xb,xa]; X2k8=fft(x2n,8); X2k16=fft(x2n,16); X3k8=fft(x3n,8); X3k16=fft(x3n,16); figure(2); N=8; 绘制8点DFT 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; 绘制8点DFT的幅频特性图 ①/ 兀）;ylabel（ 幅度）; f=2/N*(0:N-1); subplot(2,2,1);stem(f,abs(X2k8),.); % title((2a) 8 点 DFT[x_2(n)]);xlabel( subplot(2,2,3);stem(f,abs(X3k8),.); % title((3a) 8 点 DFT[x_3(n)]);xlabel( N=16; f=