
%rutina para obtener graficas de tiempo discreto a partir de una funcion en
%Z
tet=pi/8;
numz=[1 -cos(tet) 0];
denz=[1 -2*cos(tet) 1];
figure(11),dimpulse(numz,denz,51);
y=dimpulse(numz,denz,51);
N=length(y);
k=[0:N-1];
figure(12),stem(k,y);