2020年新版实验一利用DFT分析信号频谱.docxVIP

精选资料，欢迎下载 精选资料，欢迎下载 实验一利用DFT分析信号频谱 一、 实验目的 加深对DFT原理的理解。 应用DFT分析信号的频谱。 深刻理解利用DFT分析信号频谱的原理，分析实现过程中出现的现象及解决方法。 二、 实验设备与环境 计算机、MATLAB^件环境。 三、 实验基础理论 DFT与DTFT的关系 DFT实际上是 DFT实际上是 DTFT在单位圆上以 的抽样，数学公式表示为: N-1 N-1 _j 空 k X(k) = X(z)| 耳八 x(n)e N z” N n=0 (2 — 1) 2、利用 DFT求DTFT 方法一：利用下列公式： 2rk X(ej )二、X(k)( ) k=0 N k= 0,1,..N - 1 (2 — 2) Sn(N， Sn(N，/2) Nsin(，/2) .N A e 2为内插函数 方法二：实际在MATLAB十算中，上述插值运算不见得是最好的办法。 由于DFT是DTFT的取 样值，其相邻两个频率样本点的间距为 —,所以如果我们增加数据的长度 N,使得到的 N DFT谱线就更加精细，其包络就越接近 DTFT的结果，这样就可以利用 DFT计算DTFT如果 没有更多的数据，可以通过补零来增加数据长度。 3、利用DFT分析连续时间函数 利用DFT分析连续时间函数是，主要有两个处理：①抽样，②截断 对连续时间信号xa (t) 一时间T进行抽样，截取长度为 M则 址 ML (2— 3)(2—4)Xa(N)「-xa(t)e4dt 二「xa (2— 3) (2—4) n=0 再进行频域抽样可得 Xa(r1)| Xa(r1)| 送，T Xa(nT)e N =TXM(k) NT n=0 因此，利用DFT分析连续时间信号的步骤如下： (1 )、确定时间间隔，抽样得到离散时间序列 x(n). 、选择合适的窗函数和合适长度 M得到M点离散序列xM (n) = x(n)w(n). (3 )、确定频域采样点数 N,要求NA M。 、利用FFT计算N点DFT,得到XM (k)。 、根据式(2 — 4)计算Xa(j0)的近似值。 利用上述方法分析连续连续时间时，应该注意以下问题： (1 )、频谱泄露 (2 )、频谱混叠 、栅栏效应和频谱分辨率 四、实验内容 1、已知x(n)={2,-1,1,1} ，完成如下要求: 、计算他的DTFT并画出卜n , n ]区间的波形。 、计算4点DFT并把结果显示在(1)所画的图形中。 、对x(n)补零，计算64点DFT并显示结果。 、 根据实验结果，分析是否可以由 DFT计算DTFT如果可以，如何实现 (2)实验代码如下： x=[2 -1 1 1]; n=0:3; w=0:0.01*pi:pi*2; X1=x*exp(-j* n*w); X2=fft(x) subplot(211); plot(w,abs(X1)); hold on; stem( n*pi/2,abs(X2),filled); axis tight; subplot(212); plot(w,a ngle(X1)); hold on; stem( n*pi/2,a ngle(X2),filled); axis tight; MATLAB图形如下： (3 )实验代码如下： N=0:63; x=[-2 -1 1 1 zeros(1,60)]; Y=fft(x); subplot(211); stem(abs(Y),filled); subplot(212); stem(a ngle(Y),filled); MATLA图像如下： 570 5 70 精选资料，欢迎下载 精选资料，欢迎下载 432io2030405060io 4 3 2 i o 20 30 40 50 60 io 答：可以由DFT计算DTFT由实验结果波形看出，序列补零后，长度越长， DFT点数越多, DFT当做其DTFT其DFT越逼近其 DFT当做其DTFT 2、考察序列 x(n)=cos(0.48 n n)+cos(0.52 n n) 0=n=10时，用DFT估计x(n)的频谱；将x(n)补零加长到长度为 100点序列用DFT 估计x(n)的频谱，要求画出相应波形。 0=n=100时，用DFT估计x(n)的频谱。并画出波形。 根据实验结果，分析怎样提高频谱分辨率 (1)实验程序代码如下： n=0:10; k=0:10; x=cos(0.48*pi* n)+cos(0.52*pi* n); Y=fft(x); subplot(211); stem(k,abs(Y),filled); subplot(212); stem(k,a ngle(Y),filled); MATLA波形如下： 精选资料，欢迎下载 精选资料，欢迎下载 将x(n)补零至100点