systemz_2.m

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

%分析系统函数H(Z)=(2*z^(-1)-z^(-2)+z^(-3))/(1+z^(-1)+z^(-2)/2)零极图，稳定性，及系统单位抽样响应
syms z n
B=[0 2 -1 1];
A=[1 1 1/2 0];
sys=filt(B,A)
%%%求零极点，在命令窗口显示结果
sprintf('系统函数的零极点为:')
[z,p,k]=tf2zp(B,A);
z
p  
%绘制零极点图
zplane(B,A)
title('零极点图')
%求单位抽样响应
[r,p,k]=residue(B,A);
p1=abs(p');
a1=angle(p')/pi;    %angle()的结果是弧度
r1=abs(r');
b1=angle(r')/pi;
sprintf('系统单位抽样响应:\nh(n)=3e^(j*pi/2)*[0.7071e^(j*3*pi/4)]^n*u(n)+3e^(-j*pi/2)*[0.7071e^(-j*3*pi/4)]^n*u(n)-2Q(n)')