% SM18.M
% Example 1 of the article
clear y den num
i=sqrt(-1);
den(1,:)=conv([-1 .55-.75*i],[-1 .55+.75*i]);
den(2,:)=conv([-1 .55-.65*i],[-1 .55+.65*i]);
den(3,:)=conv([-1 .55-.55*i],[-1 .55+.55*i]);
num(1,:)=conv([.55-.75*i -1],[.55+.75*i -1]);
num(2,:)=conv([.55-.65*i -1],[.55+.65*i -1]);
num(3,:)=conv([.55-.55*i -1],[.55+.55*i -1]);
for i = 1:3,
 y(:,i)=dimpulse(num(i,:),den(i,:),100);
end
%y(:,1)=y(:,1)/norm(y(:,1))^2;
%y(:,2)=y(:,2)/norm(y(:,2))^2;
%y(:,3)=y(:,3)/norm(y(:,3))^2;
%zpplot(0,[roots(den1) roots(den2) roots(den3)])
%pause
%disp('.')
plot(y)
