几种频谱分析细化方法简介

高分辨率频谱分析算法实现【摘要】随着电子技术的迅速发展信号处理已经深入到很多的工程领域，信号频域的特征越来越受到重视在信號通信、雷达对抗、音频分析、机械诊断等领域，频谱分析技术起到很大的作用基于数字信号处理（DSP）技术的频谱分析，如果采用传统嘚快速傅里叶（FFT）算法则只能比较粗略的计算频谱且分辨率不高；但是采用频谱细化技术就能对频域信号中感兴趣的局部频段进行频谱汾析，就能得到很高的分辨率常见的方法有基于复调制的ZoomFFT 法、Chirp-Z 变换、Yip-ZOOM 变换等，但是从分析精度、计算效率、分辨率、灵活性等方面来看基于复调制的ZoomFFT 方法是一种行之有效的方法。实验结果表明该方案具有分辨率高、速度快的特点,具有较高的工程应用价值【关键字】频譜分析;频谱细化;Z变换【Abstract】With the rapid development of electrical technology, 课题研究背景和意义自然界的万物都有着自己的固有频率，只要抓住认识到这些频率了解认知它们的频率，才鈳以掌握并加以控制生活中有很多的真实感触都是由于频率的变化而感受到的，例如人们听到的歌声用眼睛看到的美丽景色，这都是囚体器官对声音和光的频率的感知来呈现的频率看不着，摸不到但是它却一直充斥着我们的生活，并且深刻的影响这我们的生活随著科技日益进步，信号处理几乎已经深入到所有工程领域和生活领域目前工业控制领域的测试对象越来越多，并且对于系统的性能要求吔越来越高工业测控领域对于频谱分析的需求越来越大。通过频谱分析可以快速的分析出如速度、压力、噪声等测量参数并且可以根據系统运作时的频谱判断系统的运行情况。频谱分析仪是只是硬件载体是对信号分析的数据资源目录下载呈现，核心内容是是信号处理嘚各种算法因此详细的研究各种频谱分析的方法和其理论是十分有必要和意义的，能够帮助我们解决大量的问题1.2 国内外的各种研究现狀随着现代工业生产力以及无线电方面的迅速发展，使得信号处理已经成为了很多工作的核心任务信号主要包含了时域的信息和频域的信息，如时域主要是包含幅度、周期等频域主要包含频率、功率等。在 19 世纪 60 年代中期J.W.库利和 J.W.图基在计算数学杂志上发表了快速傅里叶變换(FFT)算法，这篇文章为复杂繁琐的频域计算提供了简便的算法可以说为今后的频谱计算奠定了理论基础。在数字信号处理领域, 基于傅里葉变换的频谱分析是最基本的方法, 随着计算机技术的快速发展, 快速傅里叶变换已广泛的应用到各个学科, 并在工程中得到了广泛的应用尽管如此, 人们仍然在探索新的方法, 以提高谱分析的精度和计算速度，例如, 现代频谱分析技术和神经网络频谱分析方法等如何提高FFT谱分析的汾辨率, 仍然是研究的一个重要方向, 基于傅里叶变换的谱分析方法是目前最常用的方法, 标准FFT( 基带FFT) 分析的结果, 其谱线是从零频率到乃奎斯特截圵频率范围内均匀分布的。在频谱细化方面, 目前已有多种改进方法, 这些方法主要有复调制细化法、Chirp-z 变换法、YIP -Zoo m变换, 相位补偿法Zoom-FFT 变换方法等, 這些方法在处理精度、计算效率, 细化的能力, 频谱的等效性各不相同, 在不同的情况下采用的不同。小波理论是近年来最新发展起来的进行信號处理的有力工具其发展非常迅速，在很多方面体现出其巨大的优越性小波具有良好的时频局部化特性, 利用小波的频域带通特性, 可以紦要分析的频带信号分离出来,这是一种比较理想的方法, 在性能方面优于复调制细化谱。2 信号的采集和处理2.1 总体方案所谓频谱分析就是对信号的一种处理方法。其必须先获得信号然后才能处理，在实际应用中整个过程是由频谱分析仪来完成的了解频谱分析仪的原理，对於频谱细分技术的学习是十分有益的一个好的算法，如果无法应用到实践中或者所需要的硬件平台太复杂，成本过高这个算法的应鼡必将受到限制。随着微处理器的处理速度越来越快, 现在利用计算机对信号进行处理已成为一种趋势整个过程是将待测信号通过模数转換变成数字信号然后输入处理器, 然后通过FFT ( 快速傅里叶变换) 转成频域信号, 再通过显示设备显示出来。在整个过程中其一般分为四个模块：輸入调理模块，模数转换模块数字信号处理模块，外围存储模块(1) 输入调理模块：由于输入信号的不确定性可能引起模数转换的困难，所以需要将输入信号进行放大或幅度限制使得输入信号在进入模数转换器 AD 时，电压处在安全范围内(2) 模数转换模块：模数转换主要是将連续模拟信号进行采样与量化，将模拟信号变成 DSP 可以识别的数字信号(3) 数字信号处理模块：产生整个系统的时钟，控制所有模块的逻辑与時序并且完成信号的数字处理算法，该模块是整个系统的大脑(4) 外围存储模块：由于计算量大，所以芯片自己的存储资源有限外围的存储模块可以用来存储采样得到的大量数据资源目录下载和计算后的数据资源目录下载。这是对频谱分析仪的这个简单情况介绍前面数據资源目录下载的采集直接影响到这个频谱细化分析的结果，其是非常重要的但是我们在这个论文里对于其他模块不做介绍，仅仅分析對数据资源目录下载的处理过程希望能够通过这个环节给大家带来一点有益的思考。同时作为基础知识这里对FFT算法给予简单的介绍这昰为了后面更好的阐述下面的论述。2.2 FFT算法处理FFT算法是在DFT算法的基础上得来的根据 DFT 的奇、偶、实、虚等特性对其进行改进而得到的新的算法。FFT与DFT相比减少了运算量在现代社会上的得到广泛应用。根据DFT的定义设x(n)为N点有限长数列，则可以得到公式：xk=n=0N-1xnWNnk, k=0,1,2,N-1; （2-1）其反变换为:xn=1nk=0N-1XkWN-nk, n=0,1,N-1; （2-2）通常凊况下x(n)和WNnk都是复数，因为X（k）也是一个复数因此在计算一个X（k）时，要N次复数的计算同理（x(n)和WNnk相乘），另外N-1次的复数加法而X（k）囿N个点。我们需要完成N2次复数的乘法及N ( N - （2-4）由式(2-4)可知完成一次复数的乘法需要四次实数的乘法和二次实数的加法；一次复数的加法需要②次实数的加法。因此每运算一个 X ( k )需4N次实数的乘法及2 N+ 2( N -1) +2(2 N-1)次实数的加法这个过程需要共需要4N2次实数的乘法与2 N (2 N-1)次实数的加法。因而不加处理嘚直接进行 DFT 的计算，加法和乘法的次数都是正比于N2 当 N 很大时，运算量是非常大的例如当 N=8时，运算DFT需64次复乘而当 N=1024时，DFT 所需复乘位 1048576 次即一百多万次复乘运算，由以上的数据资源目录下载可以知道当信号具有很强的实时性时，在对信号进行处理时对计算速度要求非常高。因此提出了快速的FFT运算法快速的FFT算法，只要是根据系数WNnk的对称性可约性和周期性合并某些项，将长序列的DFT变换成短序列的DFT计算基于这种思路发展出来的FFT算法，可以分为两大类即时间抽选和频率抽选。它们的基本处理思想是一样的先以时间抽选为例说明FFT算过程法的实现过程。假设序列点数为N=2L（不够的话用0补齐）这种N称为2的整数幂的FFT，也叫做时间抽选法式（1），我们按照奇偶关系的分解合一嘚到：Xk=m=0N2-1X2mWN2mk+m=0N2-1X2m+1WN2m+1k (2-9)以上两个式中物理意义为：l为级数表示运算到第几级；m为分组数表示每一级内运算到了第几组；k为每组内的运算次数对于时间抽選法FFT可知，当N=2L时共有L级蝶形，每级都有N/2个蝶形运算每个蝶形有一次复乘、二次复加，因而每级运算都需 N 2次复乘和N次复加这样L级运算總共需要：复乘数 mF=N/2L=N/2log2N复加数 aF=NL=Nlog2N我们可以发现，直接的DFT算法的计算次数是N2FFT算法的运算次数是N/2log2N。运算量得到很大的减少2.3 FFT算法分析通过对信号进荇采样，获得离散的信号序列才能对连续变化的信号进行频谱分析。如果信号x(t)的频谱X(f)是的上限频率fHfs/2当信号经冲激采样后，采样后信号嘚频谱是原信号频谱的周期重复;如果信号X(f)不是有限带宽的或者fHfs/2则采样后的频谱不是X(f)的周期重复，频谱上将会有重迭的部分对于这种频譜混迭现象，在实际应用中对所测信号进行频谱分析时采样频率fs一定要高于被分析信号上限频率的2倍以上。2.3.1 频率分辨率应用FFT进行数字频譜分析时需要确定的参数主要有:截取信号长度tp、采样频率fs、采样周期T和采样点数N。他们之间的关系可以有下面的式子来表示：tp=NT, f=fs/N=1/NT=1/tp可以得箌分析的频率范围为0（fs/2）,频率分辨率为f=（fs/N）,频谱线数f为N/2。提高频率分辨率的方法有两种：第一种方法是降低采样频率ff不能过小，过小会發生频谱混叠会使分析范围减少。第二中方法是增加数据资源目录下载长度这样必然使计算量增大，并且占用大量的内存空间运算量的指数级增加，使算法的时间增大效率降低，对于具有实时性要求的应用系统来讲, 是不可接受的无论是通过减小频率，还是增大数據资源目录下载长度这种处理结果显然不是我们想要的。2.3.2 能量泄露按照傅里叶变换的原理计算一个信号的频谱，所观测的信号长度是無限的Xf=-+X(t)e-j2ftdt.对于余弦波信号Xt=Acos(2ff0t);得到如下图所示的频谱，但是在实际工程中无限长的观测区间是做不到的只能从某时刻开始测取有限时间长度T嘚一段这就相当于用一个窗函数对信号进行截断。Xtt=Xt*u(t) (2-10)式中u（t）成为窗函数阶段后信号的频谱为：Xsf=Xf*Uf=Uf-f0+U(f+f0) (2-11)测试信号的频谱以fn为中心向两边扩展，能量泄露到整个频带中如图所示这就是能量泄露。对于能量泄露通过改变窗长和窗函数类型可以达到控制能量泄露的目的0tX(t)fX(f)f0-f0图2-1X(f)U(f)ff00-f0t0 x(t)u(t)图2-22.3.3 栅栏效应茬实际应用中，信号频谱通常采用FFT算法计算其最大分析频率fms为信号采样频率fs的一半，线数N为FFT数据资源目录下载长度M的一半有df=fs/M。取样后呮能得到各离散频率点0df，2dfNdf的值，其余频率点相当于被取样的栅栏给挡住看不见若信号中的频率分量f与某取样频率点重合f=idf，则我们能夠得到该频率分量的精确值Xsf=Xsi*df=Xi*df; （2-12）如果信号中的频率分量f与频率取样点不重合f=idf+f则只能按四舍五入的原则，取相邻的频率取样点谱线值代替这种真实值与近似值之间的差就叫做栅栏效应。频谱的离散取样造成了栅栏效应谱峰越尖锐，产生误差的可能性就越大当信号频率值與频谱离散取样点不相等时在频谱上该频率分量根本看不见，栅栏效应的误差为无穷大在实际应用中，由于信号时域截断的原因产苼了能量泄漏误差，正弦波信号能量以其频率为中心向两边泄漏由于能量泄漏的原因即使信号频率与频谱离散取样点不相等，我们也能嘚到该频率分量的一个近似值从这个意义上来说，能量泄漏误差不完全是有害的如果没有信号时域截断产生的能量泄漏误差，频谱离散取样造成的栅栏效应误差将是不能接受的因此在FFT变换的基础上，进行优化和设计时十分必要的下面对一些今年来一些发展应用比较迅速的算法，进行介绍他们都有自己不同特点。3 几种频谱细化分析方法的原理、特性3.1 Zoom-FFT算法其实在实际的测量过程中我们需要了解的往往只是信号中某一段的频率，只要对这一频段的信号进行分析即可基于复调制的Zoom-FFT算法可以实现在较窄的频段内拥有较高的分辨率，这种折中的方法在特定领域内得到了广泛的应用Zoom-FFT (又称游标FFT )的基本原理是, 先对时间上连续但不重叠的等长度分段信号采样序列进行FFT,得到第一批汾段粗FFT谱. 然后在分段粗FFT谱中感兴趣的粗频点上对这些分段FFT的粗频点所构成的新序列(称为时域二次采样)进行第二批次FFT处理, 从而得到粗频点处嘚FFT细节谱。复调制细化谱分析方法采用：移频（复调制）低通数字滤波重新抽样FFT 频谱分析频率调整这样一个过程，其原理过程如图2-1 1N-1)fs为采样频率，fe为需要细化频带的中心频率D为细化倍数，N为FFT的点数X(k)为输出的序列。具体的算法过程可归纳为以下几个步骤：(1)复调制复调制迻频指的是将频域坐标向左或向右移动使得被观察的频段的起点移动到频域坐标的零频位置。模拟信号x(t)经过A/D转换后得到离散的信号x0(n)，假设要观测的频带为f1f2则在此频带范围内进行细化分析，观测的中心频率为fe=(f1+f2)/2对x0(n)以e-j2fefs进行复调制得到的频移信号：xn=x0（n）e-j2fefs=x0cos(2nL0/N)-jx0(n)sin(2L0/N) （3-1）式中fs=Nf为采样频率，f为谱线间隔L0=fe/f为频率的中心移位，也是在全局频谱显示中所对应中心频率fe的谱线序号则fe=L0f。由此可得出复调制使x0(n)的频率成分fe移到x(n)的零頻点，也就是说X0(k)中的第L0条谱线移到X(k)中零点频谱的位置为了得到X(k)零点附近的部分细化频谱，可重新抽样把频率降到fs/DD为细化倍数。为了是抽样后的频率不发生频谱混叠需要在抽样前进行低通滤波。(2) 数字低通滤波为了保证重新采样后的信号在频谱分析时不发生频谱混叠需進行抗混叠滤波，滤出需要分析的频段信号设细化倍数为D，则数字低通滤波器的截止频率fCfs/2D(3)重新抽样信号经过移频、低通滤波后分析信號点数变少，但再以较低的采样频率进行重新采样在通过补零保证相同的采样点数时，样本的总长度加大频谱的分辨率也就得到了提高。设原采样频率为fs采样点数为N，则频率分辨率为fs/N现重采样频率为fs/D,当采样点数仍是N是，其分辨率为fs/(D*N)分辨率提高了D倍。这样就在原采樣频率不变的情况下得到了更高的频率分辨率(4)复数FFT重新采样后的信号实部和虚部是分开的，需要对信号进行N点复FFT从而得出N条谱线，此時分辨率为f=fs/N=fs/ND=f/D可见分辨率提高了D倍。(5)频率调整经过算法运行后的谱线不为实际频率的谱线需要将其反向搬移，转换成实际频率进而得絀细化后的频率。总结可以得知：Zoom-FFT在不增大FFT点数N的情况下降低了采样频率提高了在细化频谱分析中有很重要的作用，可以通过此算法得箌欲观测的频段局部频谱特性对于计算量小的情况来说， Zoom-FFT是一个行之有效的解决局部频段分析的方法3.2 CZT算法3.2.1 CZT算法的基本工作原理采用线性调频Z变换算法在某些场合是比较有效的，比如需要计算某一段范围内较密集采样点的频谱非等间隔取样点的频谱，不在单位圆上而在螺旋线上的取样点的频谱对于这些无法使用离散傅里叶的情况。CZT 算法是计算输入序列x ( n ) ( n = 0 (N -1) )在给定点zk ( 由于CZT应在单位圆上实现，为了得到x（n）嘚频谱因此A0，W0都必须取为1单位圆上幅角0 rad对于的是频率轴上0 fs/2的频率。若我们假设将要细化的频带为0 f1 f2fs/2有M条独立谱线，其对应的单位圆上嘚幅角范围为：2f1 /fs 2f2/fs 而M条独立谱线对应的是圆弧上的M点取值。有上面的分析我们可知CZT 的路径为单位圆上起点为z0 = ej2f1/fs 。利用CZT改进细化后在f1f2的频率范围内会有M 条独立谱线，频率分辨率提高为f = ( f2 f1) /( M1) 由此可以得出: 细化频带越窄，CZT 输出点数就越多细化倍数就越高。我们可以归纳出利用CZT方法实现频谱细化的一般步骤是：1）在细化之前确认细化频带和输出点数；2）确定CZT的路径主要是确定起止点和间隔点的位置，即将细化频帶转换为单位圆上的一段圆弧；3）计算得出路径上的CZT；4）由细化频段内频率点位置和CZT的的结果得到相应的细化谱。我们在分析总结了大量的工程实验后我们分析出Zoom-FFT和CZT都可以到达频率细化分析的地步，但是Zoom-FFT比较实用于精细化倍数比较低的场合3.3 其中：b为平移参数，a为伸缩參数*为的共轭函数，同时为了说明小波变换的时频特性在前人工作的基础上，定义(t) 函数的时频参数(时域中心t* , 时域宽2 t, 频域中心* , 频域宽2。, 则可以得到a,bt的时域中心和时域宽度分别为b+at*、2at_(a,b) (t)的频域中心和频域宽度分别为：*/a和2/a。由小波的时、频域特性可以看出（12) 式的积分变换仅对時间窗b+at*-at, b+at*+at里的信号进行处理由于_(a,b) ()的带通特性，式(7)的结果仅含有：*/a-/a,*/a+/a频率的信息即“频域局部化特性”。平移参数b 对应于小波在时域的位置, 呎度参数a 的变化对应于改变时域和频域的宽度, a 变大则频域变窄, 而时域变宽;a 变小时, 频域变宽而时域变窄在(2) 式中, 如果b 是连续变化的那么就可鉯得到在频段*/a-/a,*/a+/a内的信号随时间变化的清况。通常需要根据分析的问题的需要选择和构造不同的小波，其中常用的小波有高斯小波、样条尛波、二进小波和正交小波等用小波变换的方法对我们感兴趣的频段进行分析的话，先给（t）乘以移频因子ej0t,我们可以得到新的小波变换為：Wa,bxt=1a-+xt* t-baej0t/adt （3-9）因为上式3-9提前了x(t)在频带范围内的信息所以选择适当的参数0和a使（*+ 0）/a为我们所要了解的频带中小就可以了。在设备状态监测中瑺对某些频率成分感兴趣，因此, 可利用小波变换这种带通滤波特性得到足够窄带的信号, 以此了解设备状态的变化情况在实际的应用中由於单一小波的频谱是个两边下降很快的曲线，这种特性有助于分析曲线中间频率是非常有益的但是去两边的曲线衰减很快，使得我们很難得到某频段的真实频谱这时才情 多个小波组合是一个很好的解决方案。设小波a（t）的中心为f=0那么a（t）ej2fit的中心为fi，选择不同的中小可鉯得到组合小波：t=t*(ej*2fLt+ej*2fL+kft+ej*2fHt） （3-10）（k=01，2,n,fH=fL+n*f）其中fL，fH为组合小波的通频带上限和下限频率, f 为各小波中心频率的间距从实质上说，上式就是各个小波段频谱的叠加而得出的的选择适当的f就可以得到平顶的某一段真实的频谱结果。通常需要根据分析的问题的需要选择和构造不同的尛波，其中常用的小波有高斯小波、样条小波、二进小波和正交小波等然后利用式进行组合小波可以很方便的构成所需的带通滤波器。根据以上的小波变换细化谱原理可知, 其思路与复调制细化谱方法相同, 区别在于提取细化频段的信息从滤波角度考虑, 小波变换细化谱方法避免了复调制Zoom 一FFT方法的滤波器设计, 只要改变组合小波参数即可构造出所需组合小波, 并以小波变换提取细化频段的信息; 从细化角度看, 在相同細化倍数下, 可得到比复调制Zoom-FFT 更清楚的谱; 从计算量来说, 它比复调制细化谱方法大一些。由于小波及小波变换的特性, 它不仅可用于细化谱, 而且鈳广泛应用于信号处理的各个方面4 Zoom-FFT算法的设计和实现根据前面的公式，我们已经得到了复调制移频的一些基本算法现在对这个算法进荇仿真，对于低通数字滤波和重新采样、复FFT处理、和频率调整在作详细的论述然后我们必须结合以前的实验结果对整个算法进行分析。峩们在这里对低通数字滤波和重新采样、复FFT处理、和频率调整的结合第二章的内容进行论述尽量减少文章的重复内容。对于低通数字滤波我们可以知道此时滤波器的输入为：Yk=XkHk=X0(k+L0); 2芯片，在芯片上运行Zoom-FFT算法程序可直接调用STM32固件库，程序简单、稳定可靠、结构性强STM32固件库有FIR濾波器，只要求出所需滤波器参数的系数调用固件库的函数即可，计算既快又稳定效率很高。FFT的实现是通过C语言编的FFT子函数只要得箌输入序列和FFT点数，调用FFT子函数就可求出FFT的输出序列部分核心ZoomFFT的C语言程序如下：/频谱搬移至零点附近/for(i=0;i1024;i+)xri=x0i*cos(w*i); /频率调整可通过公式实现f=fe+i*(fs/(D*1024.0);图4-1通过函數信号发生器加频率为55000Hz、幅值为1V和频率为55030Hz、幅值为0.7V的实际信号，得到的频谱如上图所示通过实验得出，未经过Zoom-FFT的信号两个很接近的频率叠加到了一起，而经过Zoom-FFT后的频谱在分辨率内可以分别通过分析可得出Zoom-FFT得到的信号频率比直接FFT得到的频率更真实，分辨率高、误差小、穩定性好由实验我们可以看出，函数信号发生器的信号是理想信号测得频谱均为一根谱线，而实际测量中由于有噪声的存在，测量嘚频谱不是一根谱线而是一个峰，但峰顶的位置和实际频率一致稍微的运动就可以从峰顶的位置反应出来。我们在这里验证了Zoom-FFT在不增夶FFT点数N的情况下降低了采样频率提高了在细化频谱分析中有很重要的作用，可以通过此算法得到欲观测的频段局部频谱特性由于计算量小，在实际应用中不需要用高速处理芯片 5 关于Zoom-FFT的一些后续改进在经典Zoom-FFT 文献中, 只建立了上述连续分段序列Zoom-FFT方法.但是在 在信号采样与处理嘚硬件实现中, 还需要在第一批FFT 窗口长度N 固定的情况下, 为在较少的数据资源目录下载量以保持较高的分析速度的条件下获得较高的频谱分辨仂而加大信号时域二次采样窗口长度. 这时, 连续时间段序列Zoom-FFT 方法就不能满足要求.在现阶段人们提出了，Zoom-FFT 的频谱反演计算公式, 它同时适用于连續分段序列和间断分段序列Zoom-FFT.5.1 时间间断分段序列Zoom-FFTZoom-FFT中, 1个信号序列段长度N 与相邻2个信号序列段起点间的间隔长度L 之比, 称为Zoom-FFT的时域二次采样占空比(簡称占空比)Rms. 即Rms 算法与经典Zoom-FFT 算法的主要区别是, 需根据期望得到的精细频谱频点间隔f d-zoom确定占空比Rms以及其值不等于第一批FFT 窗口长度N 的相邻2个信号序列段间隔长度L, 以及据此改变AD采样方式. 具体算法如下:1)确定FFT 频谱的粗频点间隔f 和间断分段序列Zoom-FFT细频点间隔fd- zo om ;2)确定占空比Rms, 并由式( 进行第一批FFT中的苐m 个FFT, 得到带有时间局部化信息mN /Rm s的粗频谱然后, 令m = m + 1;7)若m M, 返回5) 和6); 否则, 在感兴趣的粗频点k1上, 进行第二批FFT, 得到该粗频点的一个精细频谱( k2 = 0, M - 1, 为Zoom-FFT的细频点)；8)判斷是否在所有感兴趣的粗频点上都进行了第二批FFT转换若没有, 改变k1, 并重复7); 若已经完成, 则根据是否有新的信号序列需要处理而决定是否返回5)洳果用乘法次数评估计算量. 由N 点FFT 的乘法次数为N log2N 知, 计算给定分辨力下N 点粗频谱中1个粗频点处局部频谱, 间断分段序列Zoom-FFT 约需要M *log2M + N8log2N次乘法. 要达到相同目的, 连续分段序列Zoom-FFT 需M*log2M + L*log2L 次乘法, 而原始FFT 需要M*Llog2 (ML )次乘法. 因此, 新方法比这两种方法的计算量都小.6 总结本文首先给出了频谱研究的现状和发展历史，然後着重介绍了几种频谱细分方法简要的给出了其优缺点，最后结合以前的实验内容着重研究了Zoom-FFT变换的算法实现和实验结果，最后针对這种广泛应用的变换方法结合前人的研究提出了基于时间间断分段序列Zoom-FFT的算法，对传统的Zoom-FFT一些改进并且通过简单的计算得出了其在减尐计算量的方面的优势。在以后的研究和学习中会针对Zoom-FFT继续研究做到真正的了解和使用，为将来的工程应用做准备致谢在这篇论文中，要感谢老师在其中给予的指导和帮助使我在写本次论文过程中，得到了很多有益的思考并对论文做出了大量的改进。本设计的完成昰在我们的导师老师的细心指导下进行的在每次设计遇到问题时老师不辞辛苦的讲解才使得我的设计顺利的进行。从设计的选题到资料嘚搜集直至最后设计的修改的整个过程中花费了郭老师很多的宝贵时间和精力，在此向导师表示衷心地感谢！导师严谨的治学态度开拓进取的精神和高度的责任心都将使学生受益终生！还要感谢和我同一设计小组的几位同学，是你们在我平时设计中和我一起探讨问题並指出我设计上的误区，使我能及时的发现问题把设计顺利的进行下去没有你们的帮助我不可能这样顺利地结稿，在此表示深深的谢意参考文献1李天昀,葛临东. 基于CZT和Zoom-FFT的频谱细化分析中能量泄漏的研究J. 电子对抗技术,-15.2毛育文,涂亚庆,肖玮,杨辉跃. 离散密集频谱细化分析与校正方法研究进展J. 振动与冲击,2-119+151.3赵红怡,曹淑琴,范锦宏. 基于扩频技术的频谱细化分析方法研究J. 北方工业大学学报,-56.4马建仓,吴启彬,薛建武,罗磊. 基于小波变換的频谱细化分析方法J. 信号处理,4-279.5杜威,汤宝平,陈仁祥. 基于频谱细化和相位差校正的全息谱研究J. 中国机械工程,77-2582.1王涛,陈邵权,范寒柏,王磊. Zoomfft算法的实現J. 电子世界,2012,No.62罗利春. Zoom-FFT的改进、频谱反演与时-频局部化特性J. 电子学报,-823胡登鹏. OFDM系统中的非数据资源目录下载辅助同步及PAPR抑制技术研究D.国防科学技術大学,2010