数字信号处理实验报告三__用FFT对信号作频谱分析.docVIP

实验三用FFT对信号作频谱分析 姓名： 班级： 学号： 一、 实验目的 学习用FFT对连续倍号和时域离散倌号进行谱分析的方法，了解可能出现的分析误差 及其原W，以便正确应用FFT。 二、 实验原理与方法 用FR对倌号作频谱分析是学习数字倌号处理的秉要P、j容。经常需要进行谱分析的信 号是模拟信号和时域离散信号。对信号进行谱分析的重耍问题是频谱分辨率D和分析误差。 频谱分辨率直接和FFT的变换区闽N柯关，因为FFT能够实现的频率分辨率是2;z7；V ,因 此要求27T/7VSZ）。可以根裾此式选择FFT的变换区间N。误差主要來自于川FFT作频谱 分析吋，得到的足离散谱，而信号（周期信号除外）足连续谱，只侖当N较人吋离散谱的包 络才能逼近于连续谱，因此N要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的讼度作FFT,得到的离散谱冰能代表周 期信兮的频谱。如果不知道信兮周期，可以尽fi选择信兮的观察时间长一些。 对模拟信号进行谱分析时，首先要按照采样定理将W变成时域离散信号。如果是模拟 周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分 析进行。 三、 实验内容及步骤 （1） 对以下序列进行谱分析。 打 +1， 0 打 3 x2（n） = S-n, 4nl 0，其他n 4-n, 0 n 3 X3（7?） = H - 3, 4 Z2 7 0，其他n 选择FFT的变换区叫N为8和16两种情况进行频谱分析。分别打印K幅频特性曲线。并 进行对比、分析和讨论。 （2） 对以下周期序列进行谱分析。 X4(Z?) = COS^Zt x5 () = cos(7rn / 4) + cos(;rn / 8) 选择FFT的变换区间N为8和16两种情况分别对以上序列进行频谱分析。分別打印 其幅频特性曲线。丼进行对比、分析和讨论。 对模拟周期信号进行谱分析 x6(t) = cos 8/n + cos 167rt + cos 2Q7rt 选择采样频率T7、=64Hz，变换区间N=16,32,64三种情况进行谱分析。分别打印艽幅频特 性，并进行分析和讨论。 、实验结果 实验源程序 %用FFT对信号作频谱分析 clear all;close all %实脸内游(1)============= %产生长度为8的三角波序列x2(n)xln=[ones(l, 4)]; %产生长度为8的三角波序列x2(n) M=8;xa=l:(M/2); xb=(M/2)1:1; x2n=[xa, xb]; x3n=[xb, xa]; %计算xln的 %计算xln的8点DFT %计算xln的16点DFT %计算xln的8点DFT ?算 xln 的 16 点 DPT %计算xln的8点DFT %计算xln的16点DFT Xlkl6=fft(xln, 16); X2k8=fft(x2n, 8); X2kl6=fft(x2n, 16): X3k8=fft(x3n，8): X3kl6=fft(x3n, 16): %以下绘制幅频特性曲线 subplot (3, 2,1):mstem(Xlk8); %给制 8 点 DFT 的幅频特性图 xlabcl ({’ w/n ’ ；’(la) 8 点 DFT[x_l (n)]’}):ylabelC 幅度’): axis([0, 2, 0,1. 2*max(abs(Xlk8))]) subplot (3, 2,2) ;mstem(Xlkl6); %绘制 16 点 DFT 的幅频特性图 xlabel ({’ ⑴/ it ’ ；’(lb) 16 点 DFT[x_l (n)]’}):ylabelC 幅度’); axis([0, 2, 0, 1. 2*inax(abs(X1k16))]) subplot (3, 2, 3):mstem(X2k8): %绘制 8 点 DFT 的描频特性图 xlabel (f w/n’ ：’(2a) 8 点 DFT[x_2(n)]}) :ylabel (’幅度’ k axis([0, 2, 0,1. 2*max (abs (X2k8))]) subplot (3, 2,4) ;mstem(X2kl6); %绘制 16 点 DFT 的幅频特性图 xlabel({’ u/ it ’ ：’(2b) 16 点 DET[x_2(n)]’}):ylabelC 幅度): axis ([0, 2, 0, 1. 2*raix (abs (X2kl6))]) subplot (3, 2, 5) ;mstem(X3k8); %绘制 8 点 DFT 的幅频特性图 xlabel ({’ w/ ji ’ ；’(3a) 8 点 DFT[X-3(n)]’ } ;ylabel (’ 幅度’); ax is ([0,