数字信号处理第三章补.ppt

解 (1) 由分辨率的要求确定最小长度tp 所以记录长度为 (2) 从信号的最高频率确定最大可能的采样间隔T（即最小采样频率fs=1/T）。 按采样定理 即 (3) 最小记录点数N应满足 取 对于事先不知道信号的最高频率，可以根据信号的时域波形图来估计它。 例如， 某信号的波形如图 4-23 所示。 先找出相邻的波峰与波谷之间的距离，如图中t1，t2，t3，t4。 然后，选出其中最小的一个如t4。这里, t4可能就是由信号的最高频率分量形成的。 峰与谷之间的距离就是周期的一半。 因此，最高频率为 知道fh后就能确定采样频率 图 4-23 估算信号最高频率fh 3.6.5 用DFT对模拟信号作谱分析时的几个问题 利用FFT对连续信号进行傅里叶分析时可能造成的误差如下。 1. 频谱混叠失真 在前面的基本步骤中，A/D变换前利用前置低通滤波器进行预滤波，使xc(t)频谱中最高频率分量不超过fh。假设A/D变换器的采样频率为fs，按照奈奎斯特采样定理，为了不产生混叠， 必须满足 一般应取 fs=(2.5～3.0)fh （4-58） 如果不满足fs≥2fh，就会产生频谱混叠失真。 也就是采样间隔T满足 对于FFT来说，频率函数也要采样，变成离散的序列，其采样间隔为F(即频率分辨率)。 由式（4-55）可得 （4-59） 从以上Ｔ和tp两个公式来看，信号的最高频率分量fh与频率分辨率F存在矛盾关系，要想fh增加，则时域采样间隔T就一定减小， 而fs就增加，由式（4-57）可知，此时若是固定N，必然要增加F， 即分辨率下降。 * * 第3章 快速傅里叶变换 3.6 DFT的应用 3.6.1 利用DFT计算线性卷积 3.6.1 利用DFT计算线性卷积——快速卷积 1. 快速卷积法 以FIR滤波器为例，因为它的输出等于有限长单位脉冲响应h(n)与有限长输入信号x(n)的离散线性卷积。 设x(n)为L点，h(n)为M点， 输出y(n)为 用DFT算法也就是用圆周卷积来代替这一线性卷积时，为了不产生混叠，其必要条件是使x(n)，h(n)都补零值点，补到至少N=M+L-1， 即: 0≤n≤L-1 L≤n≤N-1 0≤n≤M-1 M≤n≤N-1 然后计算圆周卷积 N 这时，y(n)就能代表线性卷积的结果。  用FFT计算y(n)的步骤如下：  ① 求N点X(k)=DFT［x(n)］, N点； ② 求H(k)=DFT［h(n)］， N点； ③ 计算Y(k)=X(k)H(k); ④ 求y(n)=IDFT［Y(k)］，N点。 3.6.2 利用DFT计算线性相关 –快速相关 线性相关的定义由第一章已表明可用卷积运算来表示相关运算，重写如下： 就像线性卷积与圆周卷积一样，讨论圆周相关的目的也是用DFT来计算线性相关和，DFT运算所对应的圆周相关就是一种快速相关运算。 圆周相关如下： 若相关函数的DFT为： 则有： 圆周相关定理 利用DFT来求线性相关步骤如下： ③求L点X(k)=DFT［x(n)］, L点； 求L点Y(k)=DFT［y(n)］， L点； ④计算Rxy(k)=X(k)Y*(k); ⑤求 =IDFT［Rxy(k)］，L点。 ① 给定x(n),N1点，y(n),N2点； ②将x(n),y(n)补零补到 点； ⑥ 确定rxy（0）的定位，由于圆周相关定理求出的 的m全部是正值，而线性相关在m为正数及m为负数时皆有值，因而有m=0的定位问题。 3.6.3 利用DFT对模拟信号的傅里叶变换（级数） 对的逼近 讨论内容： 用DFT逼近连续非周期信号的傅里叶变换。 用DFT逼近连续周期信号的傅里叶级数。 用DFT逼近有限长信号的傅里叶变换。 用DFT做傅里叶变换(级数)的逼近时所产生的 问题。 1、用DFT逼近连续非周期信号的傅里叶变换 在信号与系统中详细讨论的连续非周期信 号的傅里叶变换是连续非周期性的频谱函数，数字计算机难于处理的，因而我们采 用 DFT对其进行逼近。 连续时间非周期信号x(t)的傅里叶变换为 设：对连续非周期信号进行时域抽样，抽样间隔 为 T (时域)；对其连续非周期性的频谱函数进行频域抽样，频域抽样周期为 F (频域). 又因时域抽样，频域必然周期延拓；且延拓周期为时域抽样的频率值，即频域周期 fs = 1/ T ; 从频域抽样理论知识可知：频域抽样后