《信号与系统分析》课件第6章.ppt

xlabel(′Frequency＼omega(unit:＼pi)′);subplot(2,2,2);plot(w/pi,angle(H));title(′Phaseresponse′);xlabel(′Frequency＼omega(unit:＼pi)′);h=impz(b,a,40);subplot(2,1,2);stem(h);title(′Impulseresponse′);xlabel(′n′);运行结果如图6-5所示。图6-5【例6-25】运行结果即将y(0)=2代入差分方程得y(0)-2y(-1)=f(0)而所以将Yx(z),Yf(z)进行逆z变换就得出yx(n)和yf(n)为得出y(n)=yx(n)+yf(n)=[3(3)n-2n]U(n)方法二先求出零输入响应yx(n),再利用系统函数求出yf(n)。(1)求yx(n)。求yx(n)可以用z变换法，也可以用时域法，下面用这两种方法分别求出yx(n)。①z变换法。将差分方程写为齐次差分方程，即y(n)-2y(n-1)=0两边取z变换得Yx(z)-2z-1[Yx(z)+yx(-1)z］=0整理得而由前知所以则②时域法。由差分方程知，特征值λ=2,则yx(n)=C2n,n≥0将代入上式得C=1所以yx(n)=2n,n≥0(2)求yf(n)。由差分方程y(n)-2y(n-1)=f(n)可得h(n)-2h(n-1)=δ(n)两边取z变换得H(z)-2z-1H(z)=1由式(6-30)得得出yf(n)=[-2(2)n+3(3)n]U(n)综上可得y(n)=yx(n)+yf(n)=[3(3)n-2n]U(n)6.5.3离散系统因果性、稳定性与H(z)的关系在第5章已经知道，一个离散LTI系统是因果系统的充分必要条件是h(n)=0,n0或h(n)=h(n)U(n)即h(n)为因果序列。由于因果序列z变换的收敛域是|z|R,因此，如果系统函数的收敛域具有|z|R的形式，则该系统是因果的；否则，系统是非因果的。这样系统因果性的充分必要条件可以用H(z)表示，即系统函数H(z)的收敛域为|z|R,R为某非负实数类似地，可以用系统函数来研究稳定性问题。已经知道离散系统为稳定系统的充分必要条件是上式表明，在单位圆|z|=1上是收敛的，根据收敛域的定义，单位圆在的收敛域内。因此，系统为稳定的充要条件可以表示为：系统函数H(z)的收敛域包含单位圆。如果系统是因果的，那么稳定性的条件是H(z)的收敛域是包含单位圆在内的某个圆的外部，由于收敛域中不能含有极点，故H(z)的所有极点均应在单位圆内。因此，因果系统稳定的充要条件是：H(z)的所有极点均在单位圆内。6.5.4离散系统z域分析的MATLAB实现【例6-19】一离散系统的差分方程为y(n)-by(n-1)=f(n),求激励f(n)=anU(n),y(-1)=0,求响应y(n)。解程序代码如下：%programch6-19symsnabzf=a^n;F=ztrans(f);H=1/(1-b*z^(-1));Y=H*F;y1=iztrans(Y);y=simplify(y1);运行结果如下：Y=(a^(1+n)-b^(1+n))/(a-b)即。【例6-20】条件同【例6-19】，改变初始条件为y(-1)=2,求系统的响应y(n)。解对差分方程两边取单边z变换，得Y(z)-bz-1Y(z)-by(-1)=F(z)代入y(-1)=2得到程序代码如下：%programch6-20symsnabzF=z/(z-a);Y=(F+2*b)/(1-b*z^(-1));即。y1=iztrans(Y);y=simplify(y1);运行结果如下：y=(a^(1+n)-b^(1+n)+2*b^(1+n)*a-2*b^(2+n))/(a-b)6.5.5利用MATLAB分析H(z)的零极点与系统特性系统函数的零点和极点可以通过MATLAB函数roots得到，也可借助函数tf2zp得到，它们的调用形式为p=roots(s)［z,p,k］=tf2zp(b,a)其中,s为多项式系统向量，p