hieunvce

ww=-pi:(pi/500):pi;
bb=1/4*ones(1,4);
HH=freqz(bb,1,ww);
zeroIndex=find(HH<10^-6);
HHzero=HH(zeroIndex);
wwzero=ww(zeroIndex);

subplot(2,1,1);
plot(wwzero,abs(HHzero));
subplot(2,1,2);