线性调频Z变换及其应用(DOC)

1、分类号编号TP32015060101夭*呼朽字局毕业论文题目线性调频Z变换及其应用学院电子信息与电气工程学院姓名包亚飞专业班级11级电信一班学号 20111060101扌旨导教师 刘保童提交日期 2015.5.22天水师范学院2015届毕业生论文原创性声明本人郑重声明:本人所指导下独立进行研究所取得用他人已经发表或未经发表明确注明出处。除文中已经任何其/、他个人或集体已经发呈交的论文是在指导教师的 的成果。学位论文中凡是引 的成果、数据、观点等均已 注明引用的内容外，不包含 表或撰写过的科研成果。本声明的法律责任由本人承担论文作者签名：论文指导教师签名：目录1引言 12 傅立叶变换的应用 12

2、.1离散傅立叶变换（DFT 22.2快速傅里叶变换（FFT 33 CZT变换 33.1 CZT变换理论分析 33.2 CZT变换的实际应用 53.3 CZT变换的运算结果仿真 64结语 7参考文献 8致谢 9线性调频Z变换及其应用包亚飞（天水师范学院，电子信息与电气工程学院，甘肃天水741000）摘要：在频谱分析领域，有多种运算方法，主要有离散傅立叶变换（DFT）算法、 快速傅里叶变换（Fast Fourier Transform， FFT）算法、线性调频Z变换等。但 是，由于FFT算法反映不出精确的信号的频谱特性，对此，在这里我们主要讨 论一种建立在DSP上的，采用FFT算法的变换方法对实序

3、列进行离散傅里叶变 换（DFT）计算的方法，即线性调频Z变换（CZT）。对于一样的数据序列，使 用CZT运算的效率是FFT变换运算的23倍，其运算结果和FFT、DFT的一样。 线性调频Z变换（CZT）可以用任意长度的采样序列，并非一定要求基 -2FFT的 长度，从而，可以使得系统得到最有效的采样率和频谱分辨率。关键词：线性调频Z变换；傅立叶变换；频谱分辨率；数据处理Chirp Z tran sform and its applicati onBao Yafei(School of Electro nic in formatio n and electrical engin eeri ng, T

4、ian shui NormalUn iversity,Tia nshui 741000,Chi na)Abstract: There are many operation methods in the field of spectral analysis, mainly has discrete Fourier Tran sform (DFT) algorithm and Fast Fourier Tran sform, Fast Fourier Tran sform, FFT) algorithm, chirp Z Tran sform, etc. But, as a result of t

5、he FFT algorithm does not reflect the precise signal spectrum characteristics, and for this issue we mainly discuss a kind of based on DSP, the tran sform method by using FFT algorithm compaction sequenee of discrete Fourier transform (DFT) calculation method, namely the chirp Z transform (CZT). For

6、 the same data sequenee, using CZT operations the efficiency of the FFT transform operations of 2 3 times, its computational results is equal to FFT and DFT. Chirp Z transform (CZT) can use the sampli ng seque nee of arbitrary len gth, does not have to request the len gth of the base 2-FFT, thus, ca

7、n make the system to get the most effective sampli ng rate and the spectral resoluti on.Key Words: Chirp Z tran sform;Fourier tran sform; Spectral resoluti on; The data process ing天水师范学院2015届毕业生论文线性调频Z变换及其应用1引言随着世界电子领域的高速发展，集成电路的应用也越来越广，其主要在计算 机、自动控制到航空、通信等各种个人电子产品中。在这些应用中，为了确保产 品的质量，从而需要对其一些参数进行测试，

8、其中大多数参数都是混合信号，所 以要对这些混合信号的数据就行采样。 采样得到的信号为连续时域信号，在根据 傅里叶变换把采样得到的信号转换成频域表示。可以根据变换的结果得到原信号 的相位谱和幅度谱。傅里叶变换频谱分析，已经成为一个不可缺少的数据分析工 具，主要有快速傅里叶变换和离散傅里叶变换。 这两种变换得到的结果是一样的， 不同的是它们的变换速度。对于某一序列x（n），用离散傅立叶变换（DFT）运算，则需要大量的时间处理 数据，如果采用快速傅里叶变换（FFT）就可以避免处理时间过长的问题， 但是， 由于FFT对有限序列的长度有限制。对此，为了解决上面在时间和序列长度方 面的问题和限制，有人就提

9、出了一种新的运算法则，即线性调频Z变换（CZT）出现了，°FFT的频谱是等间隔抽样的离散序列。若采样频率为fs,抽样点数为N, 则频域抽样间隔F。二fs/N ,如果在两谱线间频谱有很大变化时，则无法将其检测 出来。当fs保持不变时，为了减小频率抽样间隔，提高频率分辨率，只能增加抽 样点数，这将使得运算量大为增加。基于上面讨论的问题，在本文中，我主要 学习了线性调频Z变换，并作出一定的理论和仿真。2傅立叶变换的应用离散傅立叶变换（DFT）是连续傅立叶变换在离散系统中的表示形式，因为离散傅立叶变换（DFT）的计算量大，所以其应用受到了很大的限制。在这些问题的 存在下，在1956年由库利（

10、Cooley）和图基（Tukey）年发现了快速傅立叶变换（FFT）算法。FFT证明是非常适合于高效的数字实现，并且它将计算变换所需的时间减少了几个数量级。尽管快速傅里叶变换（FFT）能减少数据处理时间， 但是速傅里叶变换（FFT）的结果只会取得取样点的频谱值，却取不到取样点之 间的频谱信息。当实际频谱的峰值落在频谱取样点之间时，从FFT计算结果中得不到该峰值的真实频率、幅值和相位。如果把FFT谱的峰值作为真实频谱的峰值，必然带来频率、幅值和相位误差。2.1离散傅立叶变换（DFT）般在傅立叶变换中，由频域表示的各个分量，主要由复指数函数和正弦和余弦函数交叉来表示频域中的各个分量。 假定任何一个数

11、据波形都可以用正余弦 交叉来表示，根据DFT的定义可以得到得：N 1XkX(n )(cosn卫2 二 knN-j sin2 二 knN)(k =0,1，N -1)11在目前的一些测试系统中，通过离散傅立叶变换（DFT）可以从上面的这些 正余弦分量中直接计算出有限序列的频率。 而这样一种可以在每一个点计算的方 法，就好比是一个对振幅和相位都可调的滤波器。 如果在某一个频带内，对有限 序列进行测试，则这个可调滤波器，就可以运算出在每一频率Xk时正余弦的输出。假定将所有的余弦输出相加定义为实部，所有的正弦输出相加定义为虚部， 则实部和虚部就可用下面的两个式子来表示：2 二 kn" NJRe

12、XkX(nHcok =0,1, ；N -1)n=0N2二kn、八N 4Im Xk 】=' X(n)(sinn =0这时，在频率Xk的采样处，其功率频谱可定义为：poweXk二ReXk2 I mXk2（4）(5)在频率Xk的采样处，其相位频谱定义为：但对每一个Xk，要得到DFT结果都需要进行N次复数乘运算和N -1次复数加运算，所以对于N个X分量的总的DFT运算数据的结果共需要 N2次的复数乘运算和N(N 一1)次的复数加运算。可见这样的运算量是很大的，尽管我们才用 数列x(n)的各种对称性，这样可以加快 DFT的运算效率，但是这样会花费更大 的代价，在现实生产中这种方法往往是不能采用的

13、。2.2快速傅里叶变换(FFT)快速傅立叶变换(FFT)，一般采用基-2 FFT分蝶形运算。对N个X分量，需 要完成N /2个蝶形运算，因此需要做N /2(log 2 N )复数乘法和N /(log 2 N)次复数加法。利用FFT算法和利用DFT算法直接计算的结果相比，它们所需的复数乘法的次数之比为:号 logN2Nlo gN而它们所需的复数加法的次数之比分别为:»兽 JlogNN2N从上面两个式子可以得出，采用 FFT运算要比采用DFT的运算减少很大的运算量。但是，我们知道FFT对有限长数列的长度有严格的要求，这也是它最大的局限性所在。它的采样数据的个数必须为2n个，即(64,28

14、,2n)。同时我 们又知道采样序列的长度与FFT谱谱分辨率成反比，即FFT频谱分辨率=采样频 率/采样数。从这个式子出发，在信号的检测中，我们必须要将信号的采样频率 做出一定的调整，才可以使采样数据满足上面的条件，但是这样的采样频率却会 影响混合信号电路的检测。3 CZT变换Chirp-z变换是一种可以有效计算数据序列的功率谱和相位谱的方法，它采用螺线抽样，可适用于更一般情况下有 x(n)到X(ZQ快速算法，这种用卷 积来计算DFT变换的方法称为线性调频Z变换，简称CZT3。在对一数据进 行相同的采样的情况下，运用CZT变换的运算结果和DFT、FFT变换的运算 结果是一致的，但是FFT对信号序

15、列的长度要求有很大的限制，而CZT变换对信号序列长度可以使任意的。3.1 CZT变换理论分析在这里我们将一长度为(N 一 1)的有限序列x(n)，使用X(z)表示它的Z变换,利用C Z T运算法则，然后计算给定点Zk上的X(Zk)，令：Zk =AW(k=0,1, . M厂1)，其中W =W0j，经过一定的推论可得：K2 N jn2(k)2X(Zk)=W x(n)AWTW(0,1, . .N-1)(8)n=0研究某一组采样数列x(n)，主要按照下面的步骤进行分析：(1)选择一个最小的数L ,使其满足L =2m且L (N M -1)，以便用基-2FFT算法来求得g(n)与h(n)的卷积4(2)令g

16、(n) =x(n)Aw'补上零值点变为L点序列，并利用FFT法求g(n)序列的L点DFT：N _1£G(r) = » g(n)em(0r -1)n=0(9)(3)对相关序列h(n)进行补零加长后，其周期延拓成L点的序列，即:h(n)=0(L)2_2-0乞n乞M -1M岂n空L - NL N 1n 岂 L1(10)(4)利用FFT方法求h(n)序列的DFT:L 4、_j 2 rnH (r)八 h(n)e L 0_r _ L-1(11)n出(5)将列长为L的二序列G(r)和H(r)逐点相乘得到的是列长任为L的频域离散序列 Q(r)二G(r)H(r)。(6 )将Q(r)做

17、L点离散傅立叶反变换，这样就得到q(k)二g(k) “ h(k)二g(k) ： h(k),由于只有M点相关，故在这里只对前 M点序列进行采样。(7)最后求X(ZQ :k2X(Zk)=w2q(k)(k=0,1, ,N1)(12)从上面的公式我们可以得出，对于一组任意长度的采样数据序列，CZT变换 运算法则可以通过补零点的方法将数据序列转换成长为L的序列后直接进行运算处理，这样便解决了 FFT变换对数据序列长度的限制，且 CZT变换也可以表 示为离散卷积，其根据是利用快速傅立叶变换和快速傅立叶逆变换实现的，从而能保证对数据序列的运算速度。CZT运算相对于FFT变换需要进行更多次的变 换，运算次数和

18、N“log2N成正比例，变换时间是FFT变换时间的23倍。对一输入信号波形(fi)进行采样，设采样周期为n，采样频率为fs，频率分辨率为fres，样本大小为m，则可得到下面的两个公式：fres = fs / m(13)fi = n fres(14)对上面的两个公式(13)、( 14)合并后可得：m 二 n* fs/ fi(15)对于在一定频率下的采样数据，其输入信号的最高频率必须满足是其采样频 率的一半，这称为奈奎斯特速率。对于那些高于奈奎斯特速率的数据频率都是采 样不到的数据频率，或者称为失真的数据频率，如果不进行特定的调整，这些采 样不到的信号将会使得DFT变换在低频下的测试结果误差偏大，

19、从而失去了采 样的意义。为了更好地运用FFT变换运算，我们必须要使采样数据的频率宽度小于奈奎 斯特速率。如果要测量频率高的数据，就必须要采用更高的采样频率。根据公式 (14)式可以得出：是由数据的数量多少决定了数据频率的分辨率，即为取如果分辨率越大，则其取样就越多。而由式(15)知：采样周期为n，采样频率为fs， 频率分辨率为fres，样本大小为m等都处在一个动态平衡中，如果想要得到准确 的采样数据，就必须要设置合理的参数。3.2 CZT变换的实际应用CZT变换运算法则在一些混合信号的测试系统中得到了深入的应用，其主要运算是基于DSP完成的。还有一些对于数据频谱的分析中，对于频率的提取是 基于

20、matlab软件实现的，本文就基于matlab软件对对CZT算法的高速性和准确 性来研究，其主要为通过比较采用FFT算法和CZT算法的得到的仿真图来得出本文的结论。CZT算法，它是在更大范围的z平面上(既可以是单位圆，也可以 不是单位圆，也可以是单位圆部分频段；更可以是 z平面任意螺旋线上)来计算 X(Zk)的值5 oMATLAB信号处理工具箱中，计算 CZT的函数为 czt，调用格式为 y=czt(x,m,w,a)。它可以计算信号x的m点线性调频z变换，m、w、a、皆为标量，此函数沿 z平面指定的弧线上计算z变换的抽样值，这个弧线由 m、w和a来确定，m是 变换的长度，w是弧线上各相邻点之间

21、的复数比值，a是复数的起点。3.3 CZT变换的运算结果仿真本文中，我采用CZT算法来得到FIR滤波器过渡带中相比采用FFT算法更 为详细的频率特性。设fp = 100Hz， fst二200Hz，抽样频率fs = 1000Hz。对以上举例的matlab程序如下：clc;clear allh=fir1(40,150/500);fs=1000;f1=100;f2=200;m=1024;w=exp(-j*2*pi*(f2-f1)/(m*fs);a=exp(j*2*pi*f1/fs);H=fft(h,1000);H1=czt(h,m,w,a);fH=(0:le ngth(H)-1)'*1000

22、/le ngth(H);fH1=(0:le ngth(H1)-1)'*(f2-f1)/le ngth(H1)+f1;figure(1)subplot(211)plot(fH(1:500),abs(H(1:500);axis(1,500,0,1.2);title('FFT');xlabel('f);ylabel('H(jf)');grid subplot(212) plot(fH1,abs(H1);axis(f1,f2,0,1.2);title('czt');xlabel('f);ylabel('H(jf)');grid在matlab下运行这些程序可以得到下面的仿真图,FFT1H 0.5050100150200250300350400450500fczt1H 0.5n1-100110120130140150160170180190200f图1-对相同信号分别用FFT和CZT得到的频率特性图从上面得到的结果图可以看到，采用CZT运算远比采用FFT运算得到的频率特性图得到的效果好很多，其结果很好的反映出了信号数据的频谱特性。