clear;clc;
a0=[ 5.93141234807e-008,4.347755861717e-007,1.303747145161e-006,2.772196904973e-006,...
  4.74857800029e-006,6.910166596115e-006,8.774077218793e-006,9.856915099804e-006,...
  9.856915099804e-006,8.774077218793e-006,6.910166596115e-006,4.74857800029e-006,...
  2.772196904973e-006,1.303747145161e-006,4.347755861717e-007,5.93141234807e-008]*1e8;
n=1:16;
f=linspace(0.5,1.6,16)*pi;
fn=f'*n;
fnexp=exp(j*fn);
a016=repmat(a0,1,16);
a016=reshape(a016,16,16)';
b=a016.*fnexp;
for n=1:16
    freqz(b(n,:),1,4096,'whole',1000)
    hold on;
end