实验三：用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 4 n7 0 其它n 选择FFT的变换区间N为8和16两种情况进行频谱分析，分别打印出幅频特性曲线, 并进行讨论、分析与比较。 实验结果如下 I—ILJ I—I 文件］插入j ZLM(皇画有口0菁削■文伽端事I ?i值入工具臭面］邕匚C藉珈用 21。 21。N * fe弋-、凹倒瘁贰， d te /、凹倒变就， 分析：图(1a)和(1b)说明X1 (n) =R4n的8点DFT和16点DFT分别是XI (n)的频 谱函数的8点和16点采样；因 X3 (n) =X2 ((n-3)) 8R8(n).故X3 (n)与X2 (n)的8点DFT 的模相等，如图(2a) (2b)所示。但当N=16时，X3(n)与X2(n)不满足循环移位关系， 故图(2b) (3b)的模不同 屋］Figure 1 — □ X 文件CD 漏港? 直善Q0 插入(D THCD 卖宜曲窗口出帮助闿 □ d』点|加气公跛国厦x 哈| 口目| ■ 口 20 2 0 6 4 2 0 (5a) B^DFT[xfi(n)] (5b)倡京DFTpcJn)] 分析： X4(n)=cos( n/4)的周期为8,故N=8和N=16匀是其周期的整数倍，得到正确的单一频率正弦波的 频谱，仅在0.25 处有1根单一谱线，如图(4a)和图(4b)所示。X5(n)= cos( n/4)+ cos( n/ 8)的周期为16,故N=8不是其周期的整数倍，得到的频谱不正确如图(5a)月机N=16^l其一个周期,得到 正确的频谱，0.25 ji和0.125 ji有2根单一谱线，如图(5b)所示。 ■?j Figure 4 — □ X 文件(B遍辑旧查卷。插入(D IMCD夏面?霞口㈣箱助(B 公fe L %的佥买北▼耳口匡1| ■也 so n 咨到 1Q qL」」」 0 0.5 1 1.5 Ld/n 分析：X6 (t)有3个频率成分。f1=4Hz,f2=8Hz,f3=10H z,故其周期为0.5s。采样频率 Fs=64Hz,f1=Bf2=6.4f3 变换区间 N=64时，观察区间 TP=16T=0.24s，不是X6 (t)的整数倍周 期，故得周期不正确，如图(6a)所示。变换区间 N=32,64时，观察区间Tp=0.5s,1s,时x6(t)得 整数倍周期，所得频率正确。如图( 6b) (6c).图中3根谱线正好分别位于 4,8,10Hz处。 四、【附录】(实验中代码) 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),.); % 绘制 8 点 DFT的幅频特性图 title((1a) 8 点 DFT[x_1(n)]);xlabel( 3 / 兀);ylabel( 幅度); N=16; f=2/N*(0:N-1); 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]; % x3n=[xb,xa]; X2k8=fft(x2n,8); X2k16=fft(x2n,16); X3k8=fft(x3n,8); X3k16=fft(x3n,16); figure(2); N=8; f=2/N*(0:N-1); subplot(2,2,1