离散时间系统的数字频域分析.doc

离散时间系统的数字频域分析 z平面单位圆代表了信号序列的频谱或系统函数的频响，研究它具有特殊重要的意义。Z平面单位圆可写作 (2.52) 其模恒为1，幅角是区间的模拟量。因不重叠转一圈是，一般取位于或作为主值区。以式（2.52）代入Z变换定义式和反Z变换式,可得 DTFT[x(n)] = (2.53) 及 IDTFT[X()]x(n)= (2.54) 我们把式（2.53）、(2.54)定义为序列傅氏变换及序列傅氏反变换，有的书上也叫“离散时间傅氏变换（DTFT）”及“离散时间傅氏反变换”，分别用DTFT和IDTFT来表示，把单位圆定义为数字频域，用来表征信号和系统的频率特性。从上面的定义可知：序列傅氏变换实质上就是单位圆上的Z变换。 特别指出的一点是：序列傅氏变换X()是的周期函数，周期为2。式（2.53）从另一角度看恰恰是周期函数X()的傅氏级数展开式（参见式（2.21）），而x(n)正是傅氏级数的系数。这个概念将来会有用。 由于序列傅氏变换是Z变换的特例，因此它除了具有Z变换的一切特性外，还具有一些特殊的性质。比如，收敛条件由式（2.35）变为 (2.55) 少乘的一个因子有可能恰好是一个衰减因子，可见，序列傅氏变换不存在不等于Z变换不存在。 表2.4列出了序列傅氏变换的一些特性，其中有两条新的。一条是说：时域序列的乘积等效于在数字频域的卷积 DTFT[x(n)y(n)]= (2.56) 另一条是帕赛伐定律 当y(n)=x(n)时写成 (2.57) 这条定理的物理意义是：时域序列的总能量等于频谱的总能量。 上述两条定理在Z变换中皆有相应的公式，只因计算复杂，物理意义不直观而在前面Z变换性质中没有提及。 表2.4 序列傅氏变换特性 序列 序列傅氏变换 1 ax(n)+by(n) 2 x(n) 3 X(e) 4 X(e) 5 x(n)y(n) 6 x(n) X(e) 7 x(-n) X(e) 8 9 若x(n)是实序列 共轭偶对称：X() 即|X()|=,幅度谱是偶函数 arg[]=-arg[],相位谱是奇函数 　 2.4.2 系统频响 单位抽样响应h(n)的序列傅氏变换，或者说系统函数在单位圆上的值，叫做系统频响。序列傅氏变换不但用于分析系统，也适用于分析信号，只不过分析结果称为信号频谱而不叫系统频响罢了。 求解系统频响可有三条途径：一是直接用序列傅氏变换；二是先求系统函数，再以代入；三是先求序列的快速傅氏变换，，再用插值公式求频响，这第三种途径将在下面章节中讲到。 系统频响可分解为模和幅角两部分，分别称为幅频特性和相频特性。 (2.58) 我们通过例题来看频响的计算方法。 例2.12 一阶因果系统 差分方程：y(n)=bx(n-1)+x(n),0<b<1 系统函数：， |z|>b 求其频响。 解 幅频、相频特性曲线如图2.6。从图上看，极大值出现在处，与极点b的相角一致。这不是偶然的，有一定规律性。 例2.13 二阶因果系统 差分方程：y(n)＝ 系统函数：,且,, 求系统频响。 解 先求系统极点。根据有理分式性质，二个极点（分母上二阶有理式的根）要么是重叠的实数（时），要么是二个不同实数()，要么是二个共轭复数（）。根据题意，本题的H(z)有二个共轭极点，设为 其中 于是 其中 A 以代入H(z)。H(z)两个项的模值分别在极点频率（取[]时）或处（取时）取得极大值，并随着r的趋于1而极值增大，见图2.7。 例2.14 若横向滤波器的差分方程： 系统函数: |z|>0 求系统频响。 解 先求系统零极点 零点 极点 即 () (N-1阶) 我们注意到：所有零点分布在以a为半径的圆的等分点上，两两夹一幅角。第0个零点（k=0时），而a点正好又是极点，这对零极点对消。原点处的（N-1）阶极点对频谱的形状没影响。将代入H(z)，分别可得a<1时的频幅、相频曲线如图2.8(a)(b),以及a=1时的幅频曲线图2.8(c)。图中显示：幅频曲线的谷值位置由零点相角决定，分别出现在()，间隔。零频()处由于零极点对消而不再是零点，使该处零点间的距离最远，反倒成了幅频曲线的最大值点。比较图2.8(a)和(c)，可知a的大小不改变谷点位置，却能改变谷点值的大小。a=1时零点正好在单位圆上，该点位置的频响为零即谷点值为零。在0<a<1范围内，a越小谷点值越大。 作为一个特例，当a=1且各项乘以时，原题的差分方程变为 这是一个求N长序列平均值的平均器，可用来计算信号序列的直流分量。 若令原题差分方程只留下头尾两项，即 y(n)=x(n)-x(n-N) 对应的系统函数是，这是一个FIR系统，原点处是N阶极点，而零点由关系式可算得是 , 幅频曲线如图2.8(b)。与图2.8(c)不同之处在于处得零点没有被对消，