第2章离散时间信号与系统分析.pptVIP

第2章 离散时间信号与系统的频域分析 2.1引言 与连续时间信号频率域分析相对应，对于离散时间信号与系统，除了前面介绍的时域分析方法外，也有变换域的分析手段。就像用拉普拉斯变换和傅里叶变换将连续时间域函数f(t)转换到频率域成为F(s)或F(j?)一样，用Z变换和序列傅里叶变换把离散时间域的f(n)转换成 F(z)或 ，从而由另一个角度来观察及获取信号或系统的内在特性。 【例2.2.1】证明表2.1.1中的第12项的能量公式，亦即Parseval定理。 证： 【例2.2.2】求因果序列 的DTFT，a是实数。 解：按DTFT定义有 显然，无穷级数在|a|≥1情况下将不收敛，上式没有意义。而当|a|＜1时是一个收敛的等比级数。其和为 写成幅度与相位两部分 用matlab程序绘制信号频谱图 a=0.5；%取一个指数信号的衰减速度。 n=0：49；%观察50个采样周期的信号序列长度。 x=a.^n；%生成指数序列x(n)的值。 figure(1)； stem（n，x）；%绘制序列杆图 title(‘指数序列0.5n’)；xlabel(n)；ylabel(幅度); %标注坐标 figure(2); [X,w]=freqz(x,1,whole); %来计算序列x(n)的频谱 XA=abs(X); %求出幅度 XB=angle(X); %求出相位 subplot(2,1,1); plot(w/pi,XA); %以数字频率ω除π后为横坐标绘图。横坐标值0～2。 xlabel(数字频率ω/π); ylabel(幅度谱); subplot(2,1,2); plot(w/pi,XB); %绘制相位图 xlabel(数字频率ω/π); ylabel(相位谱); a=0.5 a=-0.5 2.2.2周期序列及其傅立叶级数表示 设 是一个以N为重复周期的序列，即满足 ，0≤n≤N－1，k是任意整数 比如，一个N=5的周期序列，主周期值｛1，2，-1，3，4｝，如图所示。 现在来考察一下复指数序列： 数字频率为 ，称它为基频分量。 k次谐波分量 第k+N次谐波 表明：所有谐波分量中，只有N个是独立的，其他谐波都是它们的重复，因而可以被替代。 将周期序列展开成这些谐波序列的无穷级数时，就可以只取k=0，1，2，…，N-1个独立谐波，从而把级数化简成有限项之和。即 ： 称周期序列的离散傅里叶级数展开，也称DFS逆变换 。 是k次谐波的系数 说明 也是以N为周期的频率域序列。周期序列在时域是周期性和离散性，到频率域对称地表现为频谱的离散性（谐波）和周期性。 k=0对应直流成分：频率=0； k=1是基频分量：频率 K次谐波：其究竟是多少呢？ 2.2.3周期序列的周期卷积（定理） 称为序列的周期卷积。根据DFS与IDFS的对偶特点，容易证明 【例2.2.3】图示如下两个N=6点的周期序列的周期卷积过程。 最后得到周期卷积输出序列 【例2.3.1】举例说明两个不同序列对应相同Z变换式的情况。 解：设 【例2.4.1】分析一阶系统的零极点位置和频率响应特性。 当0p1时，系统稳定。系统的极点z=p在单位圆内部，零点z=0位于原点。假设p=0.8，用matlab分析，程序如下 ： b=[1];% 给出分子多项式系数，按减幂排列。 a=[1,-0.8];% 给出分母多项式系数，按减幂排列。 [z,p,A]=tf2zp（b，a）；%多项式之比转换成因子形式，得到零极点和增益A。 [H,w]=freqz(b,a);%求出系统频率特性。 Ha=abs（H）；%求幅频特性 Hb=angle（H）；%求相频特性 [h,n]=impz（b,a,30）; %求单位脉冲响应,30个点 figure（1）； subplot（1，2，1）；%绘图 stem(n,h);xlabel(‘n’);ylabel(‘h(n)’)； subplot（1，2，2）； zplane(b,a); %自动绘零极点图 figure（2）； subplot（2，1，1）；%绘幅频特性图 plot（w/pi,Ha）;xlabel(‘w’);ylabel(‘H(w)’); subplot（2，1，2）；%绘相频特性图 plot（w/pi,(Hb.*180)./pi）;xlabel(‘w’);ylabel(‘φ(w)Deg’); 【例2.4.2】分析的梳状滤波器的零极点位置和频率响应特性。