options = bodeoptions;
options.FreqUnits = 'Hz';
options.FreqScale = 'linear'

fp=1000;
z=tf('z',1/fp);
n=1;
k=1;
r=0.98;

for fw=200:20:300
 b(n)=r*exp(j*2*pi*(fw/fp));
 n=n+1;
 b(n)=r*exp(-j*2*pi*(fw/fp));
 n=n+1;
end

for fz1=0:50:150
 z(n)=exp(j*2*pi*(fz1/fp));
 n=n+1;
 z(n)=exp(-j*2*pi*(fz1/fp));
 n=n+1;
end
for fz2=350:50:500
 z(n)=exp(j*2*pi*(fz2/fp));
 n=n+1;
 z(n)=exp(-j*2*pi*(fz2/fp));
 n=n+1;
end
Hz=zpk(z,bieg,k,1/fp)
bode(Hz)