%rutina para obtener graficas de tiempo discreto a partir de una funcion en
%Z
k=[0:40];
A=2.8;
alpha=0.9;
tetha0=pi/12;
N0=2*pi/tetha0;
hk=alpha.^k;
fk=A*hk.*cos(tetha0*k);
figure(22),stem(k,fk,'b')
hold
numz=A*[1 -alpha*cos(tetha0) 0];
denz=[1 -2*alpha*cos(tetha0) alpha*alpha];
FFz=tf(numz,denz,1)
ffk=dimpulse(numz,denz,41);
M=size(ffk);
kk=[0:M-1]
figure(5),stem(k,ffk,'r')
%gk=hk.*fk;
%figure(2),stem(k,gk,'k')