数字信号处理-实验二-FFT频谱分析.docxVIP

实验三：用FFT对信号作频谱分析 10.3.1实验指导 实验目的 学习用FFT对连续信号和时域离散信号进行谱分析的方法，了解可能出现的分析 误差及其原因，以便正确应用 FFT。 实验原理 用FFT对信号作频谱分析是学习数字信号处理的重要内容。经常需要进行谱分析的信号是 模拟信号和时域离散信号。对信号进行谱分析的重要问题是频谱分辨率 D和分析误差。频谱分 辨率直接和 FFT的变换区间N有关，因为FFT能够实现的频率分辨率是 2 /N，因此要求 2 /N D。可以根据此式选择 FFT的变换区间N。误差主要来自于用 FFT作频谱分析时， 得到的是离散谱，而信号(周期信号除外)是连续谱，只有当 N较大时离散谱的包络才能逼近 于连续谱，因此 N要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作 FFT，得到的离散谱才能代表周期 信号的频谱。如果不知道信号周期，可以尽量选择信号的观察时间长一些。 对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期 信号，也应该选取整数倍周期的长度， 经过采样后形成周期序列，按照周期序列的谱分析进行。 3?实验步骤及内容 (1)对以下序列进行谱分析。 X1 (n) RHn) n 1, 0 n 3 X2 (n) 8 n, 4 n 7 0 ,其它n 4 n, 0 n 3 X3( n) n 3, 4 n 7 0, 其它n 选择FFT的变换区间N 为8和16两种情况进行频谱分析。 分别打印其幅频特性曲线。 并 进行对比、分析和讨论。 (2)对以下周期序列进行谱分析。 x4(n) cos—n 4 4 x5(n) cos( n/4) cos( n/8) 选择FFT的变换区间N为8和16两种情况分别对以上序列进行频谱分析。 分别打印其幅 频特性曲线。并进行对比、分析和讨论。 (3)对模拟周期信号进行谱分析 x6(t) cos8 t cos16 t cos20 t 选择 采样频率 Fs 64Hz ，变换区间 N=16,32,64 三种情况进行谱分析。分别打印其幅频特 性，并进行分析和讨论。 4．思考题 ( 1)对于周期序列，如果周期不知道，如何用 FFT 进行谱分析？ ( 2)如何选择 FFT 的变换区间？(包括非周期信号和周期信号) (3)当N=8时，X2(n)和X3(n)的幅频特性会相同吗？为什么？ N=16呢？ 5．实验报告要求 ( 1 )完成各个实验任务和要求。附上程序清单和有关曲线。 ( 2 )简要回答思考题。 实验内容( 1) X1n=[ones(1,4)]; X1k8=fft(X1n,8); X1k16=fft(X1n,16); N=8; f=2/N*(0:N-1); subplot(1,2,1);stem(f,abs(X1k8), .); title( (la) 8 点 DFT[X_1(n)] ); Xlabel( 频谱特性 );ylabel( 幅度 ); N=16; f=2/N*(0:N-1); subplot(1,2,2);stem(f,abs(X1k16), .); title( (la) 16 点 DFT[X_1(n)] ); Xlabel( 频谱特性 );ylabel( 幅度 ); 孟一 8 舫BFTXS-)一 聖H旗 度 ts fh二 6^bFTT^n )1 Muooxaul =M、2)bM、2T1^ xznHxa-xs xsnHxb-xar x2k8Aff(x2n-8= X2k16u3(x2n-16= x3k8u3(x3n-8)- X3k16uax3n-16)- figure? NG 7r2、N* (0 乏「1 = subp-of(2-2-1 )kem(f-abs(x2k8)- m_e( -(2a)8lurDFTl2(n)」-「 x_abe一 ( -n )y_abe一 (?M) su bp-of(2-2 - 3)kem(f-abs(x3k8)- m_e( -(3a)8lurDFTl3(n)」-「 x_abe一 ( -n )y_abe一 (?M) NU1Q 7r2、N* (0 乏「1 = subplot(2,2,2);stem(f,abs(X2k16), .); title( (2a) 16 点 DFT[x_2(n)]); xlabel(频谱特性);ylabel(幅度); subplot(2,2,4);stem(f,abs(X3k16), .); title( (3a) 16 点 DFT[x_3(n)]); xlabel(频谱特性);ylabel(幅度); itii-e-(2a) 8f^ DFT[x2[n|]201 1510* -f-5 itii -e- (2a) 8f^ DFT[x2[n|] 201 15 10 * -f- 5 0 *—