1、实 验 报 告实验名称： 利用 DFT 分析连续信号频谱课程名称_数字信号处理 _院 系 部:电气与电子工程 专业班级：信息 1002学生姓名：王萌 学 号: 11012000219同 组 人: 实验台号:指导教师：范杰清 成 绩： 实验日期:华北电力大学实验六 利用 DFT 分析模拟信号频谱一、实验目的应用离散傅里叶变换 DFT 分析模拟信号 x(t)的频谱，深刻理解利用 DFT 分析模拟信号频谱的原理、分析过程中出现的现象及解决方法。二、实验原理连续周期信号相对于离散周期信号，连续非周期信号相对于离散非周期信号，都可以通过时域抽样定理建立相互关系。因此，在离散信号的 DFT 分析方法基础上
2、，增加时域抽样的步骤，就可以实现连续信号的 DFT 分析。利用 DFT 计算连续周期信号 )(tx的频谱 分析步骤为：(1) 确定周期信号的基本周期 T0；(2) 计算一个周期内的抽样点数 N。若周期信号的最高次谐频为 p 次谐波 pw0 ，则频谱中有 2p+1 根谱线；若周期信号的频谱无限宽，则认为集中信号 90%以上(或根据工程允许而定)能量的前( p+1)次谐波为近似的频谱范围，其余谐波忽略不计。取 N=2p+1；(3) 对连续周期信号以抽样间隔 T= T0 /N 进行抽样，得到 xk ；(4) 利用 FFT 函数对 xk作 N 点 FFT 运算，得到 Xm；(5) 最后求得连续周期信号
3、的频谱为 X(nw0)=Xm/N。利用 DFT 计算连续非周期信号 x(t) 的频谱 分析步骤为：(1)根据时域抽样定理，确定时域抽样间隔 T，得到离散序列 xk；(2) 确定信号截短的长度 M 及窗函数的类型，得到有限长 M 点离散序列 xMk=xkwk；(3) 确定频域抽样点数 N，要求 N=M；(4) 利用 FFT 函数进行 N 点 FFT 计算得到 N 点的 Xm；(5) 由 Xm可得连续信号频谱 X(jw)样点的近似值 )j(2mTXN三、实验内容1. 利用 FFT 分析信号 )(e)(2tutx的频谱。(1) 确定 DFT 计算的各参数（抽样间隔，截短长度，频谱分辨率等）；(2)
4、比较理论值与计算值，分析误差原因，提出改善误差的措施。解答：选取 fm=25Hz 为近似的最高频率，则抽样间隔 T= )2/(1mf=0.02s6pT进行分析，则截短点数为 N= p/300采用矩形窗，确定频域抽样点数为 512 点。频谱分辨率=25HZ/512= 0.048fsam=50;Tp=6; N=512; T=1/fsam;t=0:T:Tp; x=exp(-2*t);X=T*fft(x,N);subplot(2,1,1);plot(t,x);xlabel(t);title(时域波形);w=(-N/2:N/2-1)*(2*pi/N)*fsam;y=1./(j*w+2);subplot(
5、2,1,2);plot(w,abs(fftshift(X),w,abs(y),r-.);title(幅度谱 );xlabel(w);axis(-10,10,0,1.4)2. 分析周期信号的 )18sin(2)10cos()( tttx 频谱时，如果分析长度不为整周期，利用 fft 函数计算并绘出其频谱，总结对周期信号进行频谱分析时，如何选取信号的分析长度T0=1; N=36; T=T0/N; % 周期 T0、FFT 的点数 N、抽样间隔 Tt=0:T:T0;x=cos(10*pi*t)+2*sin(18*pi*t); %周期信号 Xm=fft(x,N)/N; %利用 FFT 计算其频谱f=(-
6、(N-1)/2:(N-1)/2)/N/T; %若 N 为偶数 f=1/T/N*(-N/2:(N/2-1);subplot(2,1,1);stem(f,abs(fftshift(Xm); %画出幅度谱xlabel(f (Hz);ylabel(magnitude);title(幅度谱 N=36); T0=1; N=90; T=T0/N; % 周期 T0、FFT 的点数 N、抽样间隔 Tt=0:T:T0;x=cos(10*pi*t)+2*sin(18*pi*t); %周期信号 Xm=fft(x,N)/N; %利用 FFT 计算其频谱f=(-(N-1)/2:(N-1)/2)/N/T; %若 N 为偶数
7、 f=1/T/N*(-N/2:(N/2-1);subplot(2,1,2);stem(f,abs(fftshift(Xm); %画出幅度谱xlabel(f (Hz);ylabel(magnitude);title(幅度谱 N=90);3. 假设一实际测得的一段信号的长度为 0.4s，其表达式为x(t)=cos(2pf1t)+0.75 cos(2pf2t)其中 f1=100Hz， f2=110Hz。当利用 FFT 近似分析该信号的频谱时，需要对信号进行时域抽样。试确定一合适抽样频率 fsam，利用DFT 分析信号 x(t)的频谱。若在信号截短时使用 Hamming 窗，由实验确定能够分辨最小谱峰
8、间隔 Df 和信号长度 Tp 的关系。若采用不同参数的 Kaiser 窗，重新确定能够分辨最小谱峰间隔 Df 和信号长度 Tp 的关系。解答：%对信号进行频谱分析fsam=440;Tp=0.4; N=55; T=1/fsam;t=0:T:Tp; f1=100;f2=110; x=cos(2*pi*f1*t)+sin(2*pi*f2*t); %周期信号 Xm=fft(x,N)/N; %利用 FFT 计算其频谱f=(-(N-1)/2:(N-1)/2)/N/T; %若 N 为偶数 f=1/T/N*(-N/2:(N/2-1);subplot(2,1,1);stem(f,abs(fftshift(Xm)
9、; %画出幅度谱xlabel(f (Hz);ylabel(magnitude);title(幅度谱 N=440); %使用 hamming 对信号进行频谱分析fsam=440;Tp=0.4; N=55;T=1/fsam;t=0:T:Tp; N=Tp/T+1;f1=100;f2=110; y=cos(2*pi*f1*t)+0.75*sin(2*pi*f2*t); %周期信号 %选择非矩形窗 hamming 窗分析w=0.54-0.46*cos(2*pi*k/(N-1);x=y.*w;Xm=fft(x,N)/N; %利用 FFT 计算其频谱f=(-(N-1)/2:(N-1)/2)/N/T; %若
10、N 为偶数 f=1/T/N*(-N/2:(N/2-1);subplot(2,1,2);stem(f,abs(fftshift(Xm); %画出幅度谱xlabel(f (Hz);ylabel(magnitude);title(幅度谱怎加 hamming 窗后分析 N=?); 4. 产生一个淹没在噪声中的信号 x(t)，例如由 50Hz 和 120Hz的正弦信号以及一个零均值的随机噪声叠加而成。确定抽样间隔和信号截短长度，分析信号的频谱，指出 50Hz 和 120Hz 的正弦成分对应的谱峰位置，详细写出检测信号的步骤和原理。解答：%对信号进行频谱分析fsam=480;Tp=0.4; N=55; T
11、=1/fsam;t=0:T:Tp; f1=50;f2=120; x=cos(2*pi*f1*t)+0.75*sin(2*pi*f2*t); %周期信号 Xm=fft(x,N)/N; %利用 FFT 计算其频谱f=(-(N-1)/2:(N-1)/2)/N/T; %若 N 为偶数 f=1/T/N*(-N/2:(N/2-1);stem(f,abs(fftshift(Xm); %画出幅度谱xlabel(f (Hz);ylabel(magnitude);title(幅度谱 N=55);四、思考题1. 既然可直接由 Fourier 变换的定义计算连续信号的傅里叶变换，为何利用 DFT 分析连续信号的频谱？
12、答：根据定义是可以根据傅里叶变换的定义直接计算连续信号的福利叶变换，但是定义区间是无限长，这在计算上是不可实施的，无论是人工计算还是通过计算机进行计算。而 DFT 是有限长的序列的傅里叶变换，在计算机上容易实现。再者，在数字信号处理中，希望能够利用数字方法直接计算常见的四种信号的频谱函数，这是需要的时域信号为有限长，其频谱也为有限项。因此常常利用 DFT 对序列进行频谱分析。2. 若信号持续时间无限，且无解析表达式，如何利用 DFT 分析其频谱？答：因为只要是有限长的离散序列都可以通过 DFT 对其进行频谱分析。故对于持续时间无限的信号，因首先得到能表征信号特征的离散序列，在对该离散序列进行
13、DFT 变换分析即可。这就涉及到模拟信号的数字化过程，具体方法是：（1）采样，根据信号的采样定理对该持续时间无限的序列进行采样。 （2）量化，将采样得到的序列进行量化得到原模拟信号对应的离散序列，让后采用 DFT 对该序列进行频谱分析即可。3. 在利用 DFT 分析连续信号频谱时，会出现哪些误差？如何克服或改善这些误差？答：在利用 DFT 分析连续信号的频谱时，会出现的误差可以分为三类。分别是混叠现象，泄露现象，栅栏现象。克服或改善的方法分别是：（1） 混叠现象：对于带限连续信号，只要提高抽样频率使之满足时域抽样定理；对于非带限信号，可以根据实际情况对其进行低通滤波，使之成为带限信号。工程中的
14、信号一般都不是带限信号，连续信号在抽样前通常都进过一个低通滤波器（即抗混叠滤波器）进行低通滤波，以减少混叠误差，提高频谱分析精度。（2） 泄漏现象：在选择矩形窗口的长度时，适当增加窗的长度，可以提高频谱分辨率，但是不能减小旁瓣引起的频谱泄露，因此可以选择旁瓣幅度很小甚至为零的非矩形窗对信号进行加窗处理，就可以降低频谱泄露。（3） 栅栏现象：改善栅栏现象最常用的方法是在离散序列之后补零，得到一个比原有序列更长的序列，这样就可以增加频谱图中的很多细节，降低栅栏现象。4. 在利用 DFT 分析连续信号频谱时，如何选择窗函数？答：在用 DFT 分析连续信号频谱时，选择窗函数一般首选矩形窗，因为对信号进
15、行加窗处理的目的是去截断信号，故一般情况下选择矩形窗就可以了。但是在对频谱分析精度要求高的情况下，就要合理选择非矩形窗，选择旁瓣幅值小甚至为零的非矩形窗以满足要求，提高频谱分析精度。5. 讨论序列后补零对频谱分析结果的影响。答：在序列后补零直接的影响就是增加了序列的长度。但是却提高了频谱分析的精度。因为序列补零后，序列长度增加了，由于抽样频率没有改变，因此频谱图中谱线之间的间隔变小了，从而显示出了更多的细节，提高了频谱分析精度。6. 窗函数对频谱分辨率有何影响？如何提高频谱分辨率？答：窗函数对频谱分析精度的影响主要是有窗函数的主瓣宽度和旁瓣幅度影响的。具体而言是：主瓣宽度越窄，精确度越高；旁瓣幅度越小，精确度越高。因此要想提高频谱分析精度，就必须选择合适的窗函数，应该使所选窗函数的长度尽量长，这样就可以减小窗函数主瓣的宽度；此外尽量选则旁瓣幅度较小甚至为零的窗函数，这样也可以减小频谱泄露，提高频谱分子精度。
