systemz_1.m

这是一个Z变换和Z反变换的matlab程序。

%分析系统函数H(Z)=(z^2+z)/(z^3+2*z^2-4*z-1)零极图，稳定性，及系统单位抽样响应
syms z n
B=[0 1 1 0];
A=[1 2 -4 1];
sys=filt(B,A)
[r,p,k]=residue(B,A);
%%%求零极点，在命令窗口显示结果
sprintf('系统函数的零极点为:')
[z,p,k]=tf2zp(B,A);
z
p  
%绘制零极点图
zplane(B,A)
%or:zplane(z,p)
title('零极点图')
%求单位抽样响应
sprintf('系统单位抽样响应:\nh(n)=[0.4902*(-3.30828)^n-0.1569*(-0.3028)^n+0.6667]*u(n)')
%%%系统的收敛域包含单位圆，故系统为稳定系统
%原则：若系统为因果系统，则收敛域在过离原点最远的极点的单位圆之外；
%      若此时系统也是稳定系统，则所有极点都在单位圆内。