文档介绍：羈实验名称螂专业、年级莀学号肀姓名肄用FFT对信号做频谱分析蒄腿膀蒅以下内容由实验指导教师填写羂实验工程完成情况膂实验工程成绩芀指导教师袆时间蚄羁荿芇2018年11月20日肂备注:按照要求独立完成实验内容。蚀实验五用FFT对信号做频谱分析葿一、实验目的蒄学****用FFT对连续信号和时域离散信号进行谱分析的方法,了解可能出现的分析误差及其原因,以便正确应用FFT。GK36S7UQMvb5E2RGbCAP袄二、实验原理葿用FFT对信号作频分析是学****数字信号处理的重要内容,经常需要进行分析的信号是模拟信号的时域离散信号。对信号进行谱分析的重要问题是频谱分辨率D和分析误差。频谱分辨率直接和FFT的变换区间N有关,因为FFT能够实现的频率分辨率是2π/N,因此要求2π/N小于等于D。可以根据此式选择FFT的变换区间N。误差主要来自于用FFT作频谱分析时,得到的是离散谱,而信号<周期信号除外)是连续谱,只有当N较大时,离散谱的包络才能逼近连续谱,因此N要适当选择大一些。GK36S7UQMvp1EanqFDPw蕿三、实验内容<包括代码与产生的图形及分析讨论):节蒂选择FFT的变换区间N为8和16两种情况进行频谱分析。分别打印其幅频特性曲线,并进行对比、分析和讨论。GK36S7UQMvDXDiTa9E3d蕿functionmstem(Xk>芆%mstem(Xk>绘制频域采样序列向量Xk的幅频特性图羄M=length(Xk>。芁k=0:M-1。wk=2*k/M。%产生M点DFT对应的采样点频率<关于pi归一化值)虿stem(wk,abs(Xk>,'.'>。boxon。%绘制M点DFT的幅频特性图蚇xlabel('w/\pi'>。ylabel('幅度'>。蒁axis([0,2,0,*max(abs(Xk>>]>。聿x1n=[ones(1,4>]。%产生序列向量x1(n>=R4(n>蝿X1k8=fft(x1n,8>。%计算x1n的8点DFT肇X1k16=fft(x1n,16>。%计算x1n的16点DFT膃%以下绘制幅频特性曲线肂subplot(3,2,1>。mstem(X1k8>。%绘制8点DFT的幅频特性图衿title('(1a>8点DFT[x_1(n>]'>。xlabel('ω/π'>。ylabel('幅度'>。GK36S7UQMvRTCrpUDGiT膄axis([0,2,0,*max(abs(X1k8>>]>袅subplot(3,2,2>。mstem(X1k16>。%绘制16点DFT的幅频特性图袁title('(1b>16点DFT[x_1(n>]'>。xlabel('ω/π'>。ylabel('幅度'>。GK36S7UQMv5PCzVD7HxA罿axis([0,2,0,*max(abs(X1k16>>]>芅x2n=[12344321zeros(1,50>]。%产生序列向量x1(n>=R4(n>GK36S7UQMvjLBHrnAILg羃X2k8=fft(x2n,8>。%计算x1n的8点DFT芀X2k16=fft(x2n,16>。%计算x1n的16点DFT蚈%以下绘制幅频特性曲线蚆subplot(3,2,3>。mstem(X2k8>。%绘制8点DFT的幅频特性图螅title('(2a>8点DFT[x_2(n>]'>。xlabel('ω/π'>。ylabel('幅度'>。GK36S7UQMvxHAQX74J0X聿axis([0,2,0,*max(abs(X2k8>>]>螈subplot(3,2,4>。mstem(X2k16>。%绘制16点DFT的幅频特性图肇title('(2b>16点DFT[x_2(n>]'>。xlabel('ω/π'>。ylabel('幅度'>。GK36S7UQMvLDAYtRyKfE膂axis([0,2,0,*max(abs(X2k16>>]>肂x3n=[43211234zeros(1,50>]。%产生序列向量x1(n>=R4(n>GK36S7UQMvZzz6ZB2Ltk袈X3k8=fft(x3n,8>。%计算x1n的8点DFT膃X3k16=fft(x3n,16>。%计算x1n的16点DFT袄%以下绘制幅频特性曲线袀subplot(3,2,5>。mstem(X3k8>。%绘制8点DFT的幅频特性图羈title('(3a>8点DFT[x_3(n>]'>。xlabel('ω/π'>。ylabel('幅度'>。GK36S7UQMvdvzfvkwMI1薄axis([0,2,0,*max(abs(X3k8>>]>莂subplot(3,2,6>。mstem(X3k16>。%绘制16点DFT的幅频特性图虿title('(3b>16点DFT[x_3(n>]'>。xlabel('ω/π'