DFT与CZT的不同应用场合程序.pptVIP

DFT与CZT的不同应用场合 目录 一、CZT的引入 二、DFT与CZT实例讲解 三、总结 一、CZT的引入 1、DFT与Z变换的关系 序列的DFT是 的Z变换在单位圆上的N点等间隔采样。 采样间隔 。 一、CZT的引入 2、DFT的局限性 DFT 是均匀分布在 Z 平面单位圆上 N 点处的频谱，如果取样点不均匀时，则很麻烦。 只研究信号的某一频段，要求对该频段取样密集，提高分辨率(如窄带信号的频谱分析）； 研究非单位圆上的取样值（如频谱峰值探测）； 如果需要准确计算N点DFT，且N 为大的素数； 当x(n)是短时间序列时,则得到的频率分辨率2 /N 是很低的。提高频谱密度的办法：用补零的方法增加点数，但DFT的点数又大大增加，使计算工作量增大。 一、CZT的引入 3、CZT的引入 采用FFT算法可以很快计算出全部N点DFT值，即Z 变换在z平面单位圆上的全部等间隔取样值。而实际 中也许不需要计算整个单位圆上Z变换的取样值，实 际中有时也对非单位圆上的取样值感兴趣。 一、CZT的引入 CZT算法：对 z 变换采用螺旋线取样，chirp-z 变换(线性调频 z 变换，Chirp z-transform) CZT的原理 设 x(n) 为已知时间序列，其 Z 变换的形式为： 一、CZT的引入 按照上式计算 Z[x(n)] 必然是从 z 的实轴开始，为得到任意起始点和以 螺旋线规律变化的 z 值，做如下假设： 其中，A、W 为任意复数，θ0 为 A 的起始角（第1个取样点（k=0） ），A0为 A 起始半径，φ0 为在 Z 平面中相邻的 zk （即 zk 和 zk+1 ）之间的夹角，W0 为任意正数值。 M为要分析的点数，不一定等于 N。 一、CZT的引入 CZT 表达式 一、CZT的引入 一、CZT的引入 A = 0.8*exp(j*pi/6); %起始半径0.8,起始角 pi/6 W = 0.985*exp(-j*pi*0.05); %W01，外旋；取样间隔 0.05pi M = 100; %计算点数 z1 = A*(W.^(-(0:M-1))); zplane([ ],z1.) 一、CZT的引入 A = 0.8*exp(j*pi/6); %起始半径 0.8，起始角 pi/6 W = 1.031*exp(-j*pi*0.05); %W01，内旋；取样间隔 0.05pi M = 100; %计算点数 z2 = A*(W.^(-(0:M-1))); zplane([ ],z2.) 二、DFT与CZT的实例对比 CZT (Chirp z-transform) 变换的 Matlab 句法为： y = czt(x,M,W,A) y = czt(x) CZT 是 x 信号沿着 W和 A 定义的螺旋线进行的 z 变换。M 是说明了变换的长度，W 是 z 平面上感兴趣的那部分螺旋线上取样点之间的比值，A 是螺旋线上的复数起始点。如果 x 是矩阵，czt(x,M,W,A) 是 x 的列变换。 y = czt(x) 使用下列缺省值： M = length(x) W= exp(j*2*pi/m) A= 1 对于这些缺省值，czt 返回 x 信号在单位圆上 m 份等间隔的 z 变换，也就是 x 的离散傅立叶变换，或者 fft(x)。空矩阵[ ] 说明了这些参数的缺省值。 三、总结 CZT变换和DFT变换相比 （1）输入与输出序列的长度不需要相等，即不需要 M=N。 （2）N与M不必是合成数，即两者均可为素数。 （3） 点的间隔不必均匀相等，这样可以得到任意的分辨率 （4）不必要在Z平面上的同心圆的轮廓上来求z变换。 （5）起始点的选择可以是任意的，因此便于从任意频率上对输入数据开始进行高分辨率的分析。 （6）当A=1，N=M 时，利用Chirp-Z变换可以求得x(n) 的DFT. （1）取样点沿螺旋线按角度间隔φo分布，φ0 0 时，取样轨迹逆时针旋转；φ0 0 时，取样轨迹顺时针旋转。 （2）zk 周线是一条螺旋线：W0 表示螺旋线的伸展率；W01，随着 k 的增加，向内盘旋，朝向原点；W01，随着k 的增加，向外盘旋； （3）若W0=1，一段圆弧，（若同时A0=1，则为单位圆一部分，此时 CZT[x(n)] = DFT[x(n)] （M=N） 例：假设x(n)