信号与系统课件第三章33频率抽样理论.ppt

§3.3 频率抽样理论 二，线性卷积和循环卷积的关系及循环卷积与线性卷积相等的条件 重叠相加法步骤： 重叠保留法（Overlap Save Method） 2 算出： 3.4.2 用DFT对信号进行谱分析 利用DFT对序列进行谱分析 利用DFT进行谱分析的误差 栅栏效应 Chirp－Z变换 CZT的定义 CZT计算方法 计算步骤： 分别做 和 的L点FFT得 和 ，通常，W预先确定，则 可先算好，无需实时运算。 求 ， 为L点序列。 做 的IFFT得v(n)，取v(n)的n∈[0，M-1]M点。 变换变量，令v(n)=V(k)，计算 ，k=0,1,…M-1，求得 。 CZT的优点 N与M无关，可以任意选取频域序列长度； 可随意选定 ，调整 的分辨率； 起始频率(z0)可以任意选定，即可考察任意一段频谱。 选W0≠1，CZT可适用于计算非单位园上频域点，W01，随着k增大分析曲线以步长 向外盘旋； W01，随着k增大分析曲线以步长 向内盘旋。 当A＝1，M=N， 时，zk均匀分布在Z平面的单位圆上，此时的CZT变换就是序列的DFT。亦即，DFT是CZT的特例。 3 但是在N一定条件下，旁瓣越小，其窗函数主瓣就越宽。即只能通过降低谱分辨率为代价，换取减小谱间干扰。 1 增大N可使窗函数主瓣变窄，提高频谱分辨率； 减小截断效应影响的措施 2 增大N对旁瓣无影响，为了减小谱间干扰，应选择其它旁瓣小而少的非矩形窗函数来截断。 *频域分辨率取决于N，分辨率与计算量是一对矛盾； *频域N点是均匀的，具有相同的分辨率。在有些应用中，可能要求一段频谱分辨率高，另一些频率范围要求低，DFT很难适应； *输入与输出点数相同，灵活性低； *DFT在整个谱中计算，有时是无必要的； DFT局限： *DFT在单位圆上计算取样，而有的应用需要在Z平面的其它区域取样分析。 若希望z可从任意点开始以螺旋线规律变化，可设： 式中，复变量z由极坐标可写作 Chirp-Z变换由雷达信号处理引入，在雷达中称“线性调频信号”为Chirp信号，故称Chirp-Z变换，简写CZT。 N点序列x(n)的Z变换定义： 式中，A0，W0为任意实数， 为起始角， 为在Z平面上相邻的zk（即zk与zk＋1，k=0,1,2,....）之间的夹角。A决定CZT的起始位置，W决定z域样点曲线的曲率。 于是， N点序列x(n)的CZT定义为： 与DFT不同，N点x(n)的CZT是M个频域点。它是分布在一段螺旋线上的频域点。 CZT定义中各参数的意义可由图3.4.8说明。 CZT可与线性系统联系起来，可由线性滤波来计算CZT。 由 CZT的定义可写作： CZT的线性滤波计算方案 若令 并代入上式，则有： 式中 由此易见，CZT可由图3.4.9所示系统实现计算。 4 对h(n)循环移位并取 点，以便与 对齐，即： 3 确定输出点数M，选择L≥N+M-1，而且L应适合作FFT运算，有 2 x(n)与 相乘，得y(n),n=0,1,...,N-1 ； 1 根据CZT分析初始点和分析曲线曲率确定A和W； 其移位情况可如下图示意： * 时间(频率)函数抽样 频率(时间)函数周期化 导 致 对Z变换在单位圆上等间隔抽样—即对 抽样 ，频率抽样点数为N，则： 的离散傅立叶级数 的主值序列，即 由DFT和DFS的关系知， 是 以N为周期的周期延拓序列 一，频域采样定理 ∵ 时域无混叠由N个 可以恢复?得到X(z) 若 MN 时域混叠 故： N? M 为频率抽样（不失真）条件 若 M?N 所以 可见， 在单位圆上的N点等间隔采样 的IDFT是原序列 以N 为周期的周期延拓序列的主值序列。 频域采样定理： 如果序列长度为M， 表示在区间[0，2π]上对 的N点等间隔采样. 则只有当 时，才能由X(k)恢复出 和 ，否则产生时域混叠现象。 且 二 ， 内插函数 通过内插函数恢复出 或 频率抽样 内