我想找到很多低通滤波器的截止频率。因此,我想要一个可以为我做这件事的函数。我可以做一个波德图,找到-3dB的频率,但这很无聊,也很耗时。有人知道我怎样才能自动完成这个程序吗?我在尝试
[mag,phase] = bode(sys)
但失败了。我该怎么做?
我有可用的滤波器系数a
、b
。当我离散化变换函数时,我试图使用一些不同的技术,因此我希望从EKV:中获得经验截止频率
y_k = b(1)*x_k + b(2)*x_{k-1} – a(2)*y_{k-1}
您可以在LTI对象(传输函数等)上使用带宽命令:
G = tf(1, [1 1])
bandwidth(G)
ans =
0.9976
如果你没有控制系统工具箱,你可以这样做:
% some filter
[b, a] = butter(5, 0.6);
% Determine frequency response
[h, w] = freqz(b, a, 2048);
% linear approximation of 3 dB cutoff frequency
ind = find(abs(h) < sqrt(1/2), 1, 'first');
slope = (abs(h(ind)) - abs(h(ind - 1))) / (w(ind) - w(ind - 1));
w_3dB = ( sqrt(1/2) - abs(h(ind - 1)) + slope * w(ind - 1) ) / slope;
% check result
figure; plot(w,abs(h))
hold on;
plot(w_3dB, sqrt(1/2), 'rx');
此外,您还需要freqz()
的DSP工具箱。
matlab中有一个函数可以为您提供截止频率及其所谓的"freqz"
如果你已经知道你的滤波器系数"a"one_answers"b"。然后freqz提供两个可变输出。一个变量存储复数幅度增益,而另一个变量则存储归一化频率。基本上,你所需要做的就是运行一个循环,找到最接近通带增益减去3dB的增益。然后读出你的频率变量的相应指数,你就有了答案。