f=200:200:3400;

for i=1:16


h=sin(pi*f./8000)./(pi*f./8000);
g=-20.*log10(h);
end
plot(f,g)