资料参考文稿03讲义dsp c4.pptVIP

重叠保留法例2 若选L＝2，M＝3，N＝M+L-1＝4，则有 * 大连海事大学信息学院电子信息基础教研室 * ④ 丢掉2点 四 线性相关的FFT算法 自相关与互相关在统计通信与数字信号处理中应用较广，同利用圆周卷积计算线性卷积一样，相关的计算也可以利用圆周相关来代替线性相关，称快速相关。 * 大连海事大学信息学院电子信息基础教研室 * 线性相关： 圆周相关： 快速相关的算法步骤： 设x(n)长为L，y(n)长为M： * 大连海事大学信息学院电子信息基础教研室 * 选 求 求 求 §4.6 线性调频z变换（Chirp-z变换）算法 * 大连海事大学信息学院电子信息基础教研室 * 一 DFT算法及存在问题 DFT算法：DFT是序列x(n)的z变换X(z)在z平面单位圆上N个等间隔点的抽样值，即整个单位圆上(0～2π)的频谱，可用FFT计算； DFT算法存在的问题： FFT要求N＝2m或高度复合数(N=r1 ? r2 ? r3 …)； 无法求解到单位圆外的变换值； 若信号为窄带信号，则不得不算许多冗余值。 线性调频z变换是解决这一问题的适用方法，它是在Z平面上采用螺旋抽样。 Chirp－z Transform----CZT 二 Chirp－z变换的原理 设x(n)为有限长(N点)序列，其z变换为： * 大连海事大学信息学院电子信息基础教研室 * 序列x(n)长度为N，要分析z平面上M点频谱采样值，分析点为zk，k＝0,1,2，…，M－1。 设　　zk＝AW－k，　0≤k≤M－1 式中A和W为复数，用极坐标形式表示为　　　A＝A0ej?0　　　W＝W0e－j?0　　∴zk＝A0W0－kej( ?0+k?0 ) 式中，A0和W0为实数。当k＝0，有 　　　z0＝A0ej?0 则zk表示z平面上起始半径为A0，角度为?0的螺旋线上角度间隔为?0的抽样点。如图所示： Chirp－Z变换分析频率点分布图 * 大连海事大学信息学院电子信息基础教研室 * 1 z0 A0 ?0 A0W0－1 z1 ?0 zM－1 当w01：螺线内缩当w01：螺线外展当φ00：逆时针转当φ00：顺时针转当A01：单位圆外当A01：单位圆内 Chirp－z变换的运算 * 大连海事大学信息学院电子信息基础教研室 * 利用Bluestein等式： Chirp－Z变换计算方框图 相位随时间n的变化率为角频率，由相位可知h(n)是一个角频率随时间线性增加的复振荡，故称“线性调频”， 国外称Chirp，相应的算法称为Chirp－z Transform(CZT)。 * 大连海事大学信息学院电子信息基础教研室 * h(n) v(n)=v(k) x(n) g(n) N点 M点 三 CZT的实现（用基2－FFT计算卷积） 用圆周卷积计算线性卷积 * 大连海事大学信息学院电子信息基础教研室 * 是N点序列 是偶对称无穷序列 所以求M点卷积和时，h(n)只取h(n)[u(n-N+1)-u(n-M)] 用圆周卷积计算线性卷积时，为保证不发生混叠失真，最少要求L＝N+M-1，考虑到要用基2－FFT，取 的最小数。 g(n) n N-1 h(n) n M-1 -(N-1) CZT算法计算步骤： 选择L，满足以上条件； g(n)补零至L点，求L点DFT，G(R)= FFT[g(n)]; 求L点DFT，H(R)= FFT[h(n)]; * 大连海事大学信息学院电子信息基础教研室 * ? 四 CZT运算量的估算（复乘次数） 3次FFT： * 大连海事大学信息学院电子信息基础教研室 * 有N次； 有2N次： ：当NM时有2N次，NM时有2M次； ：有L次； ：有M次； 共计： 次； 例L＝128，N＝50， M＝50： 则CZT：1772次 M·N：2500次 五 CZT的特点 Chirp－Z变换与DFT算法相比较，有以下特点: 输入序列长度N和输出序列长度不需要相等，