离散傅里叶变换(DFT)幻灯片.pptVIP

图 3.4.6 用DFT计算理想低通滤波器频响曲线 现在用DFT来分析ha(t)的频率响应特性。 由于ha(t)的持续时间为无穷长， 所以要截取一段Tp， 假设Tp=8 s， 采样间隔T=0.25 s(即采样速度fs=4 Hz)， 采样点数N=Tp/T=32。 此时频域采样间隔F=1/NT=0.125 Hz。 则 H(k)=T·DFT［h(n)］, 0≤k≤31 其中 h(n)=ha(nT)R32(n) 在已知信号的最高频率fc(即谱分析范围时)， 为了避免在DFT运算中发生频率混叠现象， 要求采样速率fs满足下式 fs＞2fc (3.4.9) 按照(3.4.5)式， 谱分辨率F=fs/N， 如果保持采样点数N不变， 要提高谱的分辨率(F减小)， 必须降低采样速率， 采样速率的降低会引起谱分析范围减少。 如维持fs不变， 为提高分辨率可以增加采样点数N， 因为NT=Tp，　T=f-1s， 只有增加对信号的观察时间Tp， 才能增加N。 Tp和N可以按照下式进行选择: (3.4.10) (3.4.11) 例 3.4.1 对实信号进行谱分析， 要求谱分辨率F≤10 Hz，信号最高频率fc=2.5 kHz， 试确定最小记录时间TPmin， 最大的采样间隔Tmax， 最少的采样点数Nmin。 如果fc不变， 要求谱分辨率增加一倍， 最少的采样点数和最小的记录时间是多少? 解： 因此TP min=0.1 s， 因为要求fs≥2fc， 所以 N取512 2. 用DFT对序列进行谱分析 我们已知道单位圆上的Z变换就是序列傅里叶变换， 即 为使频率分辨率提高一倍， F=5 Hz， 要求 N取1024 对周期为N的周期序列 ， 由(2.3.10)式知道， 其频谱函数为 用DFT的隐含周期性知道， 截取 的主值序列x(n)= ， 并进行N点DFT得到 其中 如果截取长度M等于 (n)的整数个周期， 即M=mN， m为正整数， 则 令n=n′+rN, r=0, 1, …, m-1, n′=0, 1, …, N-1,则 因为 k/m=整数 k/m≠整数 k/m=整数 k/m≠整数 由此可见： 如果 的周期预先不知道， 可先截取M进行DFT， 即 再将截取长度扩大一倍， 截取 比较 , 如果二者的主谱差别满足分析误差要求，则以二者之一近似表示 的频谱。 图 3.4.7 单位圆与非单位圆采样 例如， 要求计算序列在半径为r的圆上的频谱， 那么N个等间隔采样点为 , k=0, 1, 2, …, N-1， zk点的频谱分量为 令 则 (3.4.12) 3. Chirp-Z变换 设序列x(n)长度为N， 要分析z平面上M点频谱采样值， 分析点为zk, k=0, 1, 2, …, M-1。 设 zk=AW-k, 0≤k≤M-1 式中A和W为复数， 用极坐标形式表示为 (3.4.13) 式中A0和W0为实数。 当k=0时有 图 3.4.8 Chrip-Z变换分析频率点分布图 将zk代入Z变换公式得到 利用下面的关系式： 得到: 令 (3.4.14) 图 3.4.9 Chirp z变换计算框图 图 3.4.10 Chirp-Z变换中hL(n)序列的形成 由(3.4.3)式知， y(n) h(n)是V(n)的周期延拓序列的主值序列， 延拓周期为L， 即 综上所述， 可归纳出具体计算步骤如下： (1) 形成hL(n)序列 (1) (2) (3) L L (4) (5) 计算 (6) (7) 与标准DFT(FFT)算法相比较， Chirp-Z变换有以下特点：