离散时间信号与系统的频域概要.ppt

思考：当输入为 输出波形如何变化？ * b=[1,0,0,0,0,0,0,0,-1]；a=[1,-1] 例 求下图所示一阶离散系统的频率响应。 差分方程 系统函数 频响特性 幅度响应 相位响应 结论： （1）原点处的极点和零点对频率响应的幅度无影响， 它们只是在相位中引入一个线性分量； （2）极点主要影响频响的峰值，极点越靠近单位圆， 峰值就越尖锐，当极点处于单位圆上，该点的频响就出现 ∞，这相当于该频率处出现无耗谐振； （3）零点主要影响频响的谷值，零点越靠近单位圆， 谷值越小，当处于单位圆上时，幅度为0。 例2-21 已知 利用几何法分析系统的幅频特性。 解: 极点：z=0（N阶极点） 零点：令 则 N个零点等间隔分布在单位圆上。 取N=8时，极零点分布和幅频特性如图 例2-22 利用几何法分析矩形序列的幅频特性 解： 零点 极点 （N-1阶）， 设N=8，z=1处的极点和零点相互抵消。 利用滤波器分析设计工具——FDATool 例2-23 设一个因果系统的差分方程为 为实数 求系统的频率响应。 解 将差分方程等式两端取Z变换，可求得 单位脉冲响应为 该系统的频率响应为 幅度响应为 相位响应为 h(n)无限长 例2-24 设系统的差分方程为 试求其频率响应。 解 这是M-1个单元延时及M个抽头相加所组成的电路， 称之为横向滤波器。 令 将所给差分方程等式两端取Z变换，可得 系统函数为 零点满足 ，即 极点 （M-1阶极点） 其中 第一个零点 和单极点 相抵消。 当输入为 时，系统只延时（M-1）位就不存在了 故 只有M个值，即 M=6及 条件下 h(n)有限长 2.3.5 IIR和FIR系统 1.无限长单位冲激响应(IIR)系统 如果一个离散时间系统的单位抽样响应h(n)延伸到无穷 长，即n→∞时， h(n)仍有值，这样的系统称作无限长单位 脉冲响应系统，简称IIR（Infinite Impulse Response）系统。 一个线性时不变系统的系统函数可以表示为 只要有一个 不为零，则序列就是无限长的。 该系统的差分方程为 在任何时刻系统的输出响应不仅与此时刻和此时刻以前时刻的输入有关，而且与此时刻以前的输出有关。在由差分方程确定输出时，需要进行迭代运算。因而通常将这种差分方程称为递归方程，这种方程所描述的系统也称为递归系统。 2．有限长单位冲激响应（FIR）系统 如果一个离散时间系统的单位抽样响应h(n)是有限长序 列，这样的系统称作为有限长单位脉冲响应系统，简称FIR （Finite Impulse Response）系统。 ak全为零，则序列就是有限长的。 描述该系统的系统函数和差分方程分别为 在任何时刻系统的输出只与此时刻和此时刻以前的输入 有关。在由差分方程确定输出时，不需要进行迭代运算。因 而通常将这种差分方程称为非递归方程，这种方程所描述的 系统也称为非递归系统。 2.3.6 MATLAB实现 1.零极点图 在MATLAB中，可以用DSP工具箱中的zplane(b,a)函数 或pzplotz(b,a)函数，由给定的分子行向量和分母行向量绘制 成系统的零极点图，符号“o”表示零点，符号“?”表示极 点，图中还给出了用作参考的单位圆。 例2-25 已知某系统的系统函数为 求其零、极点并绘出零、极点图。 解 MATLAB实现程序： b=[0.3 0.1 0.3 0.1 0.2]; a=[1 -1.2 1.5 -0.8 0.3]; r1=roots(a) % 求极点 r2=roots(b) % 求零点 zplane(b,a) 运行结果为： r1=[ 0.1976 + 0.8796i 0.1976 - 0.8796i 0.4024 + 0.4552i 0.4024 - 0.4552i] r2=[0.3236 + 0.8660i 0.3236 - 0.8660i -0.4903 + 0.7345i -0.4903 - 0.7345i] 2.系统的频率响应 可以用freqz函数来求系统的频率响应。用法为： [H, w]= freqz(b,a,N)在上半单位圆（0??）的等间隔的N 个点上计算频率响应。 [H, w]= freqz(b,a,N,’whole’)在整个单位圆（0?2?）等间 隔的N个点上计算。 [H]= freqz(b, a, w)计算在矢量w中指定的频率处的频率 响应。 例2-26 已知因果系统