% ARMA_GENE_RI.M
%
% P. Flandrin, Sept. 25, 2003
%
% calcule et trace une rˇalisation la rˇponse
% impulsionnelle d'un filtre ARMA
% d'ordres 0, 1, 2, 3 ou 4
%
% calls  : - ARMA_coeff.m

N = 100; % nombre de points de la rˇponse impulsionnelle

p = 4; % ordre de la partie AR du mod¸le (<= 4)
q = 2; % ordre de la partie MA du mod¸le (<= 4)

n = 1; % numˇro de la figure de placement des p™les

[A,B] = ARMA_coeff(p,q,n);

% pause

figure(n+1)
clf

ind0 = 5;
imp = zeros(1,N);
imp(ind0) = 1;
x = filter(B,A,imp);

plot([1 N],[0 0],'k')
hold on
plot([ind0 ind0],[1.1*min(x) 1.1*max(x)],'k')
stem(x,'r')
title(['filtre ARMA(',int2str(p),',',int2str(q),')'])
axis([1 N 1.1*min(x) 1.1*max(x)])
xlabel('temps')
set(gca,'YTick',[])
set(gca,'XTick',[])
hold off