dm13401.m

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

%dm13401
%用MATLAB观察和分析系统零极点分布与系统时域特性之间的关系和规律
%绘制情况(a)系统零极点分布图及系统单位序列响应
z=0;                             %定义系统零点位置
p=0.25;                           %定义系统极点位置
k=1;                              %定义系统增益
zplane(z,p)                       %绘制系统零极点分布图
pause
set(gcf,'color','w')              %设置绘图窗口背景颜色
[num,den]=zp2tf(z,p,k);           %零极点模型转换为传递函数模型
impz(num,den)                     %绘制系统单位序列响应时域波形
title('h(n)')                     %定义标题
pause
%绘制情况(b)系统零极点分布图及系统单位序列响应
p=1;
zplane(z,p)
pause
[num,den]=zp2tf(z,p,k);
impz(num,den)
title('h(n)')
pause
%绘制情况(c)系统零极点分布图及系统单位序列响应
p=-1.25;
zplane(z,p)
pause
[num,den]=zp2tf(z,p,k);
impz(num,den,20)
title('h(n)')
pause
%绘制情况(d)系统零极点分布图及系统单位序列响应
p=[0.8*exp(pi*i/6);0.8*exp(-pi*i/6)]
zplane(z,p)
pause
[num,den]=zp2tf(z,p,k)
impz(num,den,20)
title('h(n)')
pause
%绘制情况(e)系统零极点分布图及系统单位序列响应
p=[exp(pi*i/8);exp(-pi*i/8)]
zplane(z,p)
pause
[num,den]=zp2tf(z,p,k)
impz(num,den,40)
title('h(n)')
pause
%绘制情况(f)系统零极点分布图及系统单位序列响应
p=[1.1*exp(3*pi*i/16);1.1*exp(-3*pi*i/16)]
zplane(z,p)
pause
[num,den]=zp2tf(z,p,k)
impz(num,den,40)
title('h(n)')