[B, A] = butter( 2, 3000/16000*2, 'high' );

fvtool

plot(w)

[H,W] = freqz(B,A,N) returns the N-point complex frequency response
vector H and the N-point frequency vector W in radians/sample of
the filter.

plt(w, angle(h));