数字信号处理第四章.pptVIP

* * 用表格法求解8点圆周卷积 1 2 3 0 0 0 0 0 y(n) 0 1 0 0 0 0 4 3 2 1 1 2 1 0 0 0 0 4 3 4 2 3 2 1 0 0 0 0 0 10 3 4 3 2 1 0 0 0 0 16 4 0 4 3 2 1 0 0 0 17 5 0 0 4 3 2 1 0 0 12 6 0 0 0 4 3 2 1 0 0 7 0 0 0 0 4 3 2 1 0 h(m) x(-m) m 信息与通信工程学院 * * 重叠保留法： 重叠保留法的出发点是由于圆周卷积中发生混叠时只影响序列的一部分： 如x1(n)、 h(n)的长度分别为L和M，则线性卷积的长度为L+M-1。 若L+M-1 N，则N点圆周卷积的前L+M－1－N个值不等于线性卷积，而后边余下的值相等。 若通过（在序列前边）补数使序列长度增加，并把卷积后（在序列前边）多余的序列值舍去，则可得线性卷积结果。 其它作法与重叠相加法相同。 (证明见教材p181-182.) 信息与通信工程学院 * * 重叠保留法的步骤： 设h(n)长为M，x(n)分段后每段长为L： 选L，使 x(n) L … L L x0(n) 补0 M-1 x0(n)是x(n)前边补M－1个零值，凑N =L+M-1个点； x1(n)是由x0(n)后边取M－1个数值，再取其后x(n)值凑够N点；如此反复取xi(n)。 x1(n) M-1 … x2(n) 求 ? 将yi(n)的前M－1个值舍去，余下的值顺序排列为 y(n)。 y0(n) y1(n) y2(n) 丢掉 丢掉 丢掉 信息与通信工程学院 * * 重叠保留法例1 x(n)分为两段，每段长为L＝6，M＝3，设N＝8，则 ⑧ ⑧ 丢掉M-1=2点 补零结果 信息与通信工程学院 * * 重叠保留法例2 若选L＝2，M＝3，N＝M+L-1＝4，则有 ④ 丢掉2点 信息与通信工程学院 * * 四 线性相关的FFT算法 自相关与互相关在统计通信与数字信号处理中应用较广，同利用圆周卷积计算线性卷积一样，相关的计算也可以利用圆周相关来代替线性相关，称快速相关。 线性相关： 圆周相关： 信息与通信工程学院 * * 快速相关的算法步骤： 设x(n)长为L，y(n)长为M： 选 求 求 求 信息与通信工程学院 * * §4.9 线性调频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)] 用圆周卷积计