clear;clc;
a0=[ 7.135916968387e-010,1.568395747901e-008,5.471539757694e-008,9.376656588003e-008,...
  9.376656588003e-008,5.471539757694e-008,1.568395747901e-008,7.135916968387e-010]*1e10;
n=1:8;
f=linspace(0.61,1.388,8)*pi;
fn=f'*n;
fnexp=exp(j*fn);
a08=repmat(a0,1,8);
a08=reshape(a08,8,8)';
b=a08.*fnexp;
for n=1:8
    freqz(b(n,:),1,4096,'whole',1000)
    hold on;
end
axis([0,1000,-50,100])