低通贝塞尔滤波器的频率响应

设计一个五阶模拟低通贝塞尔滤波器，其群延迟在频率达到 $$10^4$$ 弧度/秒前大致恒定。使用 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