实验五 利用DFT分析模拟信号频谱概要1.docxVIP

数字信号处理及实验实验报告实验题目利用DFT分析模拟信号频谱姓名MYT组别班级光电14学号144320200206【实验目的】运用离散傅里叶变换DFT分析模拟信号x(t)的频谱，深刻理解利用DFT分析模拟信号的频谱的原理、分析过程中出现的现象及解决的方法。【实验原理】连续周期信号相对于离散周期信号，连续的非周期信号相对于离散非周期信号，都可以通过时域抽样定理建立相互关系。因此，在离散信号DFT分析方法的基础上，增加时域抽样的步骤，就可以实现连续信号的DFT分析。1.利用DFT分析连续周期信号的频谱?周期为的连续时间周期信号x(t)的频谱函数定义为式中：是信号的周期；?称为信号的基频（基波）；?称为信号的谐频。?连续周期信号的频谱是非周期离散谱，谱线间隔为。相比离散周期信号的DFT分析方法，连续周期信号的DFT分析方法增加了时域抽样的环节。如果不满足抽样定理的约束条件，将会出现混叠误差。连续周期信号的分析步骤为：?（1）确定周期信号的基本周期。?（2）计算一个周期内的抽样点数N。若周期信号的最高次谐频为p次谐波。则频谱中有（2p+1）根谱线；若周期信号的频谱无限宽，则认为集中信号90%以上（或根据工程允许而定）能量的前（p+1）次谐波为近似的频谱范围，其余谐波忽略不计。取。?（3）对连续周期信号以抽样间隔T进行抽样，。?（4）利用FFT函数对x[k]作N点FFT运算，得到X[m]。?（5）最后求得连续周期信号的频谱为=。?因为当对连续周期信号按间隔T进行均匀抽样，每周期抽取N点时，则有t=KT, =NT, , 代入式（1.5.1）可得若能够按照满足抽样定理的抽样间隔抽样，并选取整周期为信号分析长度，则利用DFT计算得到的离散频谱值等于原连续周期信号离散频谱的准确值。2.利用DFT计算连续非周期信号的频谱??连续时间非周期信号x(t)的频谱函数是连续谱，定义为?相比离散非周期信号的DFT分析方法，连续非周期信号的DFT分析方法增加了时域抽样的环节。如果不满足抽样定理的约束条件，会出现混叠误差。如果信号在时域加窗截短过程中，窗口宽度（截断长度）或窗口类型不合适，则会产生较大的频率泄露而影响频谱分辨率。因此，合理地确定抽样间隔T和相应的截断长度是决定DFT能否正确地分析信号频谱的关键。?连续非周期信号的分析步骤为：?（1）根据时域抽样定理，确定时域抽样间隔T，得到离散序列x[k]。?（2）确定信号截断的长度M及窗函数的类型，得到有限长M点离散序列。?（3）确定频域抽样点数N，要求。?（4）利用FFT函数进行N点FFT计算得到N点的X[m]。?（5）由X[m]可得连续信号的频谱样点的近似值。?因为信号按T进行均匀抽样，截断长度M，则有t=KT, , 代入式（1.5.3）可得?对进行N点频域抽样，可得【实验结果与数据处理】1.利用FFT分析信号的频谱。??（1）确定DFT计算的各参数（抽样间隔T，时域截断长度，频谱分辨率△等）。??（2）比较理论值与计算值，分析误差原因，提出改善误差的措施。clc,clear,closeallfsam=50;Tp=6;T=1/fsam;N=512;t=0:T:Tp;x=exp(-2*t);X=T*fft(x,N); subplot(2,1,1); plot(t,x); w=(-N/2:N/2-1)*(2*pi/N)*fsam; y=1./(j*w+2); subplot(2,1,2); plot(w,abs(fftshift(X)),w,abs(y)); axis([-10 10 1 1.5]);由图可见，理论频谱与由DFT近似计算的频谱之间存在一定的误差，由于信号不是限带信号，在时域抽样时产生混叠，可以降低抽样频率，以减少DFT的计算量。2.假设一实际测得的一段信号的长度为0.4s，其表达式为式中：=100Hz，=110Hz。当利用FFT近似分析该信号的频谱时，需要对信号进行时域抽样。试确定一合适的抽样频率，利用DFT分析信号x(t)的频谱。若在信号截断时使用Hamming窗，由实验确定能够分辨最小谱峰间隔△f和信号长度的关系。若采用不同参数的kaiser窗，重新确定能够分辨最小谱峰间隔△f和信号长度的关系。clc,clear,closeallfsam=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;f=(-(N-1)/2:(N-1)/2)/N/T;figure(1)stem(f,abs(fftshift(Xm)));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