北京理工大学数字信号处理实验的报告.docVIP

精品 本科实验报告 实验名称： 数字信号处理实验 课程名称： 数字信号处理实验 实验时间： 任课教师： 实验地点： 4—423 实验教师： □ 原理验证 实验类型： □综合设计 学生姓名： □ 自主创新 学号/ 班级： 组 号 ： 学 院： 信息与电子学院 同组搭档： 专 业： 信息工程 成 绩 ： -可编辑 - 精品 实验 1 利用 DFT 分析信号频谱 一、实验目的 1.加深对 DFT 原理的理解。 2.应用 DFT 分析信号频谱。 3.深刻理解利用 DFT 分析信号频谱的原理，分析现实过程现象及解决办法。 二、实验原理 1 、 DFT 和 DTFT 的关系 有限长序列 x(n) 的离散时间傅里叶变换 X (e j ) 在频率区间 (0 2 ) 的 N 个等分点 { x(0), x(1),? x(k ),? x( N 1)} 上的 N 个取样值可以由下式表示： N 1 j 2 X (ej ) | 2 k x(n)e kn X (k) 0 k N 1 (2 1) N k 0 由上式可知，序列 x(n) 的 N 点 DFT X (k) ，实际上就是 x(n) 序列的 DTFT 在 N 个等间隔频率点 { X (0), X (1),?X (k ),? X ( N 1)} 上样本 X (k) 。 2 、利用 DFT 求 DTFT 方法 1 ：由 X (k) 恢复出 X (ej ) 的方法如图 2.1 所示： 图 2.1. 由 N 点 DFT 恢复频谱 DTFT 的流程由图 2.1 所示流程图可知： X (e j ) x(n)e j n 1 X (k)W kn e j n (2 2) N n N n k 0 由式 2-2 可以得到 -可编辑 - 精品 x(ej ) N 2 k) X (k) ( (2 3) k 1 N 其中 (x) 为内插函数 sin(N ) N 1 j ( ) 2 ?e 2 (2 4) N sin( 2) 方法 2 ：然而在实际 MATLAB 计算中，上诉插值公式不见得是最好的方法。由于 DFT 是 DTFT 的取样值，其相邻的两个频率样本点的间距为 2 ,所以如果我们增加数据的长 N 度 N ，使得得到的 DFT 谱线就更加精细， 其包络就越接近 DTFT 的结果，这样可以利用 DFT 来近似计算 DTFT 。如果没有更多的数据，可以通过补零来增加数据长度。 3 、利用 DFT 分析连续时间信号的频谱 采用计算机分析连续时间信号的频谱， 第一步就是把连续时间信号离散化， 这里需要进 行连个操作：一是采样，二是截断。 对于连续非周期信号 x (t ) ，按采样间隔 T 进行采样，截取长度为 M ，那么 M 1 x (nT )e j nT X ( j ) x (t )e j t dt T (2 5) n 0 对 X ( j ) 进行 N 点的频率采样，得到 M 1 j 2 kn X ( j ) | 2 T x ( nT )e N TX M (k) (2 6) k n 0 NT 因此，可以将利用 DFT 分析连续非周期信号频谱的步骤归纳如下： （ 1）确定时域采样间隔 T，得到离散序列 x(n) ； （ 2）确定截取长度 M ，得到 M 点离散序列 xM (n) x(n)w(n)，这里的 w(n) 为窗 函数。 （ 3）确定频域采样点数 N ，要求 N M 。 （ 4）利用 FFT 计算离散序列的 N 点 DFT ，得到 X M (k) 。 -可编辑 - 精品 （ 5）根据式（ 2-6 ）由 X M (k) 计算 X ( j ) 采样点的近似值。 采用上诉方法计算的频谱，需要注意如下三点问题： 1 ）频谱混叠。如果不满足采样定理的条件，频谱会很出现混叠误差。对于频谱无限宽的信号，应考虑覆盖大部分主要频率的范围。 （ 2）栅栏效应和频谱分辨率。使用 DFT 计算频谱，得到的结果只是 N 个频谱样本值， 样本值之间的频谱是未知的，就像通过一个栅栏观察频谱，称为“栅栏效应” 。频谱分辨率 与记录长度成正比，提高频谱分辨率，就要增加记录时间。 3 ）频谱泄露。对于信号截断会把窗函数的频谱会引入到信号频谱中，造成频谱泄露。解决这问题的主要办法是采用旁瓣小的窗函数，频谱泄露和窗函数均会引起误差。 因此，要合理选取采样间隔和截取长度，必要时还需考虑适当的窗。 对于连续周期信号， 我们在采用计算机进行计算时， 也总是要进行截断， 序列总是有限长的，仍然可以采用上诉方法近似计算。 4 、可能用到 MATLAB 函数与代码 实验中的 DFT 运算可以采用 MATLAB 中提供的 FFT 来实现。 DTFT 可以利用 M