dm13402.m

the code of the book come form the book of marlix laboratory with the book of singal

%dm13402
%绘制系统(1)零极点分布图及系统单位序列响应
z=[0;-2];           %定义系统零点位置
p=0.25;             %定义系统极点位置
k=1;                %定义系统增益
p=[0.8*exp(pi*i/6);0.8*exp(-pi*i/6)]
subplot(2,2,1)
zplane(z,p)
[num,den]=zp2tf(z,p,k)
subplot(2,2,2)
impz(num,den,20)
title('h(n)')

%绘制系统(2)零极点分布图及系统单位序列响应
z=[0;2];           %定义系统零点位置
p=0.25;             %定义系统极点位置
k=1;                %定义系统增益
p=[0.8*exp(pi*i/6);0.8*exp(-pi*i/6)]
set(gcf,'color','w')
subplot(2,2,3)
zplane(z,p)
[num,den]=zp2tf(z,p,k)
subplot(2,2,4)
impz(num,den,20)
title('h(n)')