저역통과 베셀 필터의 주파수 응답

최대 $$10^4$$ rad/second까지 군지연이 거의 일정한 5차 아날로그 저역통과 베셀 필터를 설계합니다. freqs를 사용하여 필터의 크기 응답과 위상 응답을 플로팅합니다.

wc = 10000;
[b,a] = besself(5,wc);
freqs(b,a)

필터의 군지연 응답을 펼쳐진 위상 응답의 도함수의 음수로 계산합니다. 군지연을 플로팅하여 차단 주파수까지 군지연이 거의 일정한지 확인합니다.

[h,w] = freqs(b,a);
grpdel = -diff(unwrap(angle(h)))./diff(w);

clf
loglog(w(2:end),grpdel)
xlabel('Frequency (rad/s)')
ylabel('Group delay (s)')
xline(wc)
grid

참고 항목

besself | freqs