北理工数字信号处理实验报告.docVIP

专业资料 word完美格式 本科实验报告 实验名称： 数字信号处理实验 课程名称： 数字信号处理实验 实验时间： 周二8~10节 任课教师： 范哲意 实验地点： 4—423 实验教师： 范哲意 实验类型： □ 原理验证 □ 综合设计 □ 自主创新 学生姓名： 任科飞 学号/班级： 1120121159组 号： 学 院： 信息与电子 同组搭档： 专 业： 信息工程 成 绩：  实验1 利用DFT分析信号频谱 一、实验目的 1.加深对DFT原理的理解。 2.应用DFT分析信号频谱。 3.深刻理解利用DFT分析信号频谱的原理，分析现实过程现象及解决办法。 二、实验原理 1、DFT和DTFT的关系 有限长序列的离散时间傅里叶变换在频率区间的N个等分点上的N个取样值可以由下式表示： 由上式可知，序列的N点DFT，实际上就是序列的DTFT在N个等间隔频率点上样本。 2、利用DFT求DTFT 方法1：由恢复出的方法如图2.1所示： 图 2.1.由 N点DFT恢复频谱DTFT的流程 由图2.1所示流程图可知： 由式2-2可以得到 其中为内插函数 方法2：然而在实际MATLAB计算中，上诉插值公式不见得是最好的方法。由于DFT是DTFT的取样值，其相邻的两个频率样本点的间距为,所以如果我们增加数据的长度N，使得得到的DFT谱线就更加精细，其包络就越接近DTFT的结果，这样可以利用DFT来近似计算DTFT。如果没有更多的数据，可以通过补零来增加数据长度。 3、利用DFT分析连续时间信号的频谱 采用计算机分析连续时间信号的频谱，第一步就是把连续时间信号离散化，这里需要进行连个操作：一是采样，二是截断。 对于连续非周期信号，按采样间隔T进行采样，截取长度为M，那么 对进行N点的频率采样，得到 因此，可以将利用DFT分析连续非周期信号频谱的步骤归纳如下： （1）确定时域采样间隔T，得到离散序列； （2）确定截取长度M，得到M点离散序列，这里的为窗函数。 （3）确定频域采样点数N，要求。 （4）利用FFT计算离散序列的N点DFT，得到。 （5）根据式（2-6）由计算采样点的近似值。 采用上诉方法计算的频谱，需要注意如下三点问题： （1）频谱混叠。如果不满足采样定理的条件，频谱会很出现混叠误差。对于频谱无限宽的信号，应考虑覆盖大部分主要频率的范围。 （2）栅栏效应和频谱分辨率。使用DFT计算频谱，得到的结果只是N个频谱样本值，样本值之间的频谱是未知的，就像通过一个栅栏观察频谱，称为“栅栏效应”。频谱分辨率与记录长度成正比，提高频谱分辨率，就要增加记录时间。 （3）频谱泄露。对于信号截断会把窗函数的频谱会引入到信号频谱中，造成频谱泄露。解决这问题的主要办法是采用旁瓣小的窗函数，频谱泄露和窗函数均会引起误差。 因此，要合理选取采样间隔和截取长度，必要时还需考虑适当的窗。 对于连续周期信号，我们在采用计算机进行计算时，也总是要进行截断，序列总是有限长的，仍然可以采用上诉方法近似计算。 4、可能用到MATLAB函数与代码 实验中的DFT运算可以采用MATLAB中提供的FFT来实现。 DTFT可以利用MATLAB矩阵运算的方法进行计算。 三、实验内容 1. ，完成如下要求： （1）计算其DTFT，并画出区间的波形。 （2）计算4 点DFT，并把结果显示在（1）所画的图形中。 （3）对补零，计算64 点DFT，并显示结果。 （4）是否可以由DFT 计算DTFT，如果可以，请编程实现。 2. 考察序列 （1）时，用DFT 估计的频谱；将补零加长到长度为100点序列用DFT 估计的频谱。要求画出相应波形。 （2）时，用DFT 估计x(n)的频谱，并画出波形。 3. 已知信号 ，其中，，。从的表达式可以看出，它包含三个频率的正弦波，但是，从其时域波形来看，似乎是一个正弦信号，利用DFT做频谱分析，确定适合的参数，使得到的频谱的频率分辨率符合需要。 4.利用DFT近似分析连续时间信号xt 四、实验代码及实验结果 第一题 实验代码 x = [2, -1, 1, 1]; n = 0 : 3; n1 = 0 : 63; z = zeros(1,60); x1 = [x z]; y = fft(x); y1 = fft(x1); w = -pi : 0.01 * pi : pi; X = x * exp(-j * n * w); subplot(211); plot(w, abs(X)); xlab