system.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)')
%%%系统的收敛域包含单位圆，故系统为稳定系统


%%零极图也可以这样画
%figure(2)
%zs=roots(B)
%ps=roots(A)
%plot(real(zs),imag(zs),'o')
%hold on 
%plot(real(ps),imag(ps),'x')
%grid
%%
%%绘制S域的零极图的方法
%sys=tf(B,A)
%[p,z]=pzmap(sys)
%pzmap(sys)
%sgrid
%%