数字信号处理 4人文科技.ppt

令0≤k≤N-1，则上式为： 这里的P’(k) 和G’(k) 的实部都是周期性偶序列，而它们的虚部都是周期性奇序列，此情况与（4.36）式中的复序列P(k)和G(k)的情况相同。因此有： 上两式中0≤k≤N-1。 4.7.2 一个2N点的实序列 将一个2N点的实序列x(n)按偶数点和奇数点分组形成两个N点实序列： 则有： k=0,1,┅,N-1 (4.48) 其中： 实序列p(n) 和g(n) 的DFT P(k) 和G(k) 可以采用4.7.1节所说的方法作一次N点复序列的FFT而同时得到，然后再按（4.48）式进行组合便得到了2N点实序列x(n) 的DFT。 4.8 Matlab方法 4.8.1 利用Matlab计算FFT 在第3章中介绍用Matlab方法计算信号的DFT时，提到了函数fft(x,N) 和ifft(x.N)。对于这两个函数，如果N为2的正整数幂，则可以得到本章中介绍的基2 FFT快速算法； 如果N既不是2的正整数幂，也不是质数，则函数将N分解成质数，得到较慢的混合基 FFT算法；如果 N 为质数，则fft函数采用原来的 DFT 算法。这里我们不再详细介绍，具体用法参见3.6节。 4.8.2 用Matlab实现有限长序列的Chirp z变换 在Matlab中实现有限长序列的CZT(Chirp z-transform)算法的函数为： y=czt(x,M,W,A) 该函数的用法。 （1）y=czt(x,M,W,A) 该函数返回信号x的线性调频z变换（CZT）的值y，它是信号x沿着W 和A定义的螺旋线进行的z 变换。其中M是信号x的长度（点数），W是z平面上感兴趣的那部分螺旋线上抽样点之间的比值，A是螺旋线上的复数起始点。 （2）y=czt(x) 同样返回信号x的线性调频z变换（CZT）的值y，其中M、W和A都取缺省值，即有： W=exp(2*j*pi/M), A=1 可以看出，对于这些缺省值，y返回 x 信号在单位园上等间隔的M个点上的z变换，也就是x的离散傅立叶变换(DFT)。 例4.1 已知序列x(n)=(0.5)n , 0≤n≤12，试计算序列x(n)在单位圆上的CZT, 并与该序列的DFT 进行比较。 解： 运行结果如图4.14所示，从图中可以看出，此时序列在单位圆上的CZT就等于该序列的DFT。 例4.1 已知序列x(n)=(0.5)n , 0≤n≤12，试计算序列x(n)在单位圆上的CZT, 并与该序列的DFT进行比较。 解： 运行结果如图4.14所示，从图中可以看出，此时序列在单位圆上的CZT就等于该序列的DFT。 图4.14 序列 在单位圆上的CZT及其DFT 例4.2 假设序列x(n)由4个频率分别为6 Hz、6.3 Hz、9 Hz和8 Hz的正弦序列组合而成，抽样频率为40 Hz，时域抽样200点。 （1）?????? 用 CZT 计算DFT； （2）?????? 直接计算DFT； （3）?????? 在5～10Hz的频段范围求CZT。 (a) CZT (b) DFT 图4.15 CZT的应用 在图4.15中，（a）图是利用CZT求取DFT，（b）图是直接求DFT，两者结果相同；图中频率为6Hz和6.3Hz的两个正弦信号的频谱不易分辩出。 图（c）是在5～10Hz这个频率范围内求出的CZT，它的起始点不在处而且，由于它的分点比较细，所以图中四个正弦信号的频谱都可以分辨出来。 与DFT的正变换式比较，可知完全可以利用FFT的计算程序，只需要将X*(k) 作为输入序列，并将最后结果取共轭，再除以N就得到x(n)。 4.6 线性调频z变换（CZT）算法 如果所需要的频率抽样点并不均匀地分布在单位圆上，此时，常常采用Chirp z变换（CZT）算法 4.6.1 基本原理 用CZT算法可以计算下列给定点zk上的X(z)（即在这些点处的复频谱）： (4.26) 式中， ， ，W0 、A0 为正实数。 这些z平面上的抽样点所沿的周线是一条螺旋线，参数W