北理个dsp实验报告.docxVIP

? Is信息与电子学院 ? Is 信息与电子学院 導豪北京理工大4* Se oBeijing Institute of Technology Se o 本科实验报告 实验名称：基于MATLAB数字信号处理实验 课程名称： 数字信号处理 实验时间： 任课教师： 实验地点： 4-423 实验教师： 范哲意 实验类型： ■原理验证 □综合设计 □自主创新 学生姓名： 捋卿： 组 号： 学 院： 信息与电子学院 同组搭档： 专 业： 电子科学勺技术 成 绩： 实验1利用DFT分析信号频谱 -实验目的 加深对DFT原理的理解。 84应用DFT分析信号的频谱。 深刻理解利用DFT分析信号频谱的原理，分析实现过程中出现的现象及解决方法。 二、 实验设备与环境 计算机、MATLAB软件环境。 三、 实验基础理论 DFT与DTFT的关系 有限长序列x(n)(O n N - 1)的离散时间傅里叶变换X (e加)在频率区间(0 co 2兀)的N个等间隔分布的点ke = 2 7T k/N(0 kN- 1)上的个取样值可以由下式表示 N-1 X(e j0))L=2Kk/N = x(n)e~^kn = x(/c) 0/cN-l (2-1) k=0 由上式可知，序列x(n)的N点DFT的X(k),实阪上就是x (n)序列的DTFT在N个等间隔频率点 kco = 2 兀 k/N(0 kN- 1)上样本 X (k)。 利用DFT求DTFT 方法1:由X(k)恢复出X @加)的方法如图2. 1所示： )图2. ) 图2. 1由N点DFT恢复频谱DTFT的流程 由图2. 1所7JX流程可知: OO (2-2)X(ei?)=》x(n) e~;wn (2-2) n=-8 由式2-2可以得到 其中(x)为内插函数X(/c)(2-3)k=l(2 — 4) 其中(x)为内插函数 X(/c) (2-3) k=l (2 — 4) 方法2:然而在实际MATLAB计算中，上述插值运算不见得是最好的办法。由于DFT是DTFT 的取样值，其相邻两个频率样本点的间距为2tt/N ,所以如果我们增加数据的长度N,使得 到的MT谱线就更加精细，其包络就越接近DTFT的结果，这样就可以利用MT来近似计算DTFT。 如果没有更多的数据，可以通过补零来增加数据长度。 利用DFT分析连续时间信号的频谱 采用计算机分析连续时间信号的频谱，第一步就是把连续时间信号离散化，这里需要逬 行两个操作：一是采样，二是截断。 对于连续时间非周期信号X°(t),按采样间隔T逬行采样，截取长度为九那么 Xa Xa(jQ)= f xa(t)e~jQtdt = T 丿一8 (2 — 5)》xa(nT)e_7Qnr (2 — 5) n=0 对Xa(jQ)逬行N点频域采样，得到 X(eQM-l=7 X(eQ M-l =7》Xa(n^n n=0 =TXm (Zc) (2-6) 因此，可以将利用DFT分析连续非周期信号频谱的步骤归纳如下： 确定时域采样间隔T,得到离散序列x(n); 确定截取长度M,得到M点离散序列XM(n) = x(n)w(n),这里\v(n)为窗函数。 确定频域采样点数N,要求N M o 利用FFT计算离散序列的N点DFT,得到XM(k)o 根据式(2-6)由X“(k)计算Xa(jQ)采样点的近似值。 采用上述方法计算X a(t)的频谱，需要注意如下三个冋题： 频谱混叠。如果不满足采样定理的条件，频谱会出现混叠误差。对于频谱无限宽 的信号，应考虑覆盖大部分主要频率分量的范围。 栅栏效应和频谱分辨率。使用DFT计算频谱，得到的结果只是N个频谱样本值， 样本值之间的频谱是未知的，像通过一个栅栏观察频谱，称为“栅栏效应” O频谱分辨率与 记录长度成反比，要提高频谱分表率，就要增加记录时间。 频谱泄露。对信号截断会把窗函数的频谱引入信号频谱，照成频谱泄露。解决这 个问题的主要办法是采用旁瓣小的窗函数，频谱泄露和窗函数均会引起误差。 因此，要合理选取采样间隔和截取长度，必要时还需考虑加适当的窗。 对于连续时间周期信号，我们在采用计算机逬行计算时，也总是要逬行截断，序列总是 有限长的，仍然可以采用上述方法近似计算。 4.可能用到的MATLAB函数与代码 实验中DFT运算可采用MATLAB中提供的函数fft来实现。 €-加 e_jn2 Q■■ €-加 e_jn2 Q ■ ■ ■ e_jnN Q (2-7) X (e；Q)=》x[n]e_yQn = …,x[71n]] n=n〔 实验内容 1.已知x(n)二{2, T, 1, 1},完成如下要求: 计算其DTFT,并画出［-it f7i］区间的波形。 计算4点DFT,并把结果显示在(1)所画的图形中。 对x(n)补零，计算64点DFT,并显示结果。 根据实验结果