yulewalk
设计递归数字滤波器。
库::`工程师`
争论
输入参数
# f — 频率值
+
向量资料
Details
频率值设置为一个向量,其元素范围从 0 以前 1,在哪里 1 对应于奈奎斯特频率,或采样频率的一半。 第一个元素 f 必须等于 0,而最后一个 — 1. 所有中间元素应按升序排列。 向量资料 f 其可具有对应于频率响应的步骤的重复频率值。
| 数据类型 |
|
例子:
使用Yule-Walker方法设计低通滤波器
Details
让我们设计一个低通滤波器 8-具有归一化截止频率的阶数 0.6. 让我们构造它的频率响应,并在其上叠加相应理想滤波器的特性。
import EngeeDSP.Functions: yulewalk, freqz
f1= [0, 0.6, 0.6, 1]
m1= [1, 1, 0, 0]
b1,a1=yulewalk(8,f1,m1)
h1,w1=freqz(b1,a1,128)
p=plot(w1/π,mag2db.(abs。(h1)),
xlabel="ω/π",
ylabel="幅度(dB)",
网格=真,
传说=错误,
线色=:蓝色)
vline!(f1[2:3],linestyle=:dash,linecolor=:black)

通过指定更宽的过渡带来增加抑制带中的衰减。
f2 = [0, 0.55, 0.6, 0.65, 1]
m2 = [1, 1, 0.5, 0, 0]
b2,a2 = yulewalk(8,f2,m2);
h2,w2 = freqz(b2,a2,128)
plot!(w2/π, mag2db.(abs.(h2)),
linecolor = :red)
