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