% coefficienti dell'equazione alle differenze
a = [1 -0.9];
b = [1];
M = numel(b)-1;
N = numel(a)-1;

% passo di campionamento in omega
L = 1000;
dw = 2*pi/L;

% vettore dei valori di omega
w = -pi+(0:L)*dw;
num = b(1); % b_0
for m = 1:M
    num = num + b(1+m)*exp(-1i*w*m);
end
den = a(1); % a_0
for n = 1:N
    den = den + a(1+n)*exp(-1i*w*n);
end
H = num./den;

subplot(2,1,1)
hold off
plot(w,abs(H))
% hold on
% plot(w,abs(1./(-log(-a(2))+1i*w)),'r')
xlabel('\omega')
ylabel('|H(e^{j\omega})|')
subplot(2,1,2)
plot(w,unwrap(angle(H)))
xlabel('\omega')
ylabel('arg(H(e^{j\omega}))')
