数字信号处理实验指导书.docVIP

班级:11 电信 姓名：牟万明 学号：1141301037 Experiment 2 Spectral analysis using FFT 实验2 用FFT对信号作频谱分析 1．实验目的 学习用FFT对连续信号和时域离散信号进行频谱分析（也称谱分析）的方法，了解可能出现的分析误差及其原因，以便正确应用FFT。  2. 实验原理 用FFT对信号作频谱分析是学习数字信号处理的重要内容。 经常需要进行谱分析的信号是模拟信号和时域离散信号。 对信号进行谱分析的重要问题是频谱分辨率D和分析误差。 频谱分辨率直接和FFT的变换区间N有关，因为FFT能够实现的频率分辨率是2π/N， 因此要求2π/N≤D。 可以根据此式选择FFT的变换区间N。 误差主要来自于用FFT作频谱分析时，得到的是离散谱，而信号（周期信号除外）是连续谱，只有当N较大时，离散谱的包络才能逼近于连续谱，因此N要适当选择大一些。 周期信号的频谱是离散谱，只有用整数倍周期的长度作FFT，得到的离散谱才能代表周期信号的频谱。 如果不知道信号周期，可以尽量选择信号的观察时间长一些。 对模拟信号进行谱分析时，首先要按照采样定理将其变成时域离散信号。如果是模拟周期信号，也应该选取整数倍周期的长度，经过采样后形成周期序列，按照周期序列的谱分析进行。 3． 实验内容及步骤 （1）对序列进行谱分析: 选择FFT的变换区间N为8和16的两种情况进行频谱分析。分别打印其幅频特性曲线，并进行对比、分析和讨论。  （2）对周期序列进行谱分析 选择FFT的变换区间N为8和16的两种情况分别对以上序列进行频谱分析。分别打印其幅频特性曲线。并进行对比、分析和讨论。 （3）对模拟周期信号进行谱分析: x6(t)=cos8πt+cos16πt+cos20πt 选择采样频率Fs=64 Hz，变换区间N=16，32，64的三种情况进行谱分析。分别打印其幅频特性，并进行分析和讨论。 4.实验程序 实验1程序：lab 2.m %实验1：对序列进行频谱分析 x1n=[4,3,2,1,1,2,3,4,]; x2n=[2,2,2,2]; x1k8=fft(x1n,8); x2k8=fft(x2n,8); x1k16=fft(x1n,16); x2k16=fft(x2n,16); N=8; f=2/N*(0:N-1); figure(1); subplot(2,2,1); stem(f,abs(x1k8),.); title(a.8点FFT[x1n]); xlabel(W/pi); ylabel(幅度); N=16; f=2/N*(0:N-1); figure(1); subplot(2,2,2); stem(f,abs(x1k16),.); title(a.16点FFT[x1n]); xlabel(W/pi); ylabel(幅度); N=8; f=2/N*(0:N-1); figure(1); subplot(2,2,3); stem(f,abs(x2k8),.); title(b.8点FFT[x1n]); xlabel(W/pi); ylabel(幅度); N=16; f=2/N*(0:N-1); figure(1); subplot(2,2,4); stem(f,abs(x2k16),.); title(b.16点FFT[x1n]); xlabel(W/pi); ylabel(幅度); %实验2：对周期序列进行频谱分析 N=8; n=0:N-1; x3n=cos(pi*n/4); x4n=cos(pi*n/4)+cos(pi*n/8); x3k8=fft(x3n,8); x4k8=fft(x4n,8); x3k16=fft(x3n,16); x4k16=fft(x4n,16); figure(2); N=8; f=2/N*(0:N-1); subplot(1,2,1); stem(f,abs(x3k8),.); title(c.8点FFT[x3n]); xlabel(W/pi); ylabel(幅度); N=16; f=2/N*(0:N-1); figure(2); subplot(1,2,2); stem(f,abs(x3k16),.); title(c.16点FFT[x3n]); xlabel(W/pi); ylabel(幅度); figure(2); N=8; f=2/N*(0:N-1); subplot(1,2,1); stem(f,abs(x4k8),.); title(d.8点FFT[x4n]); xlabel(W/pi); ylabel(幅度); N=16; f=2/N*(0:N-1); figure(2); subplot(1,2,2)