Julia中实现MATLAB freqz功能？解决freqresp调用报错问题

你的错误在于直接将系数数组组合后传给freqresp，该函数需要接收构造好的滤波器对象，而非原始系数数组。以下是两种可行的实现方式：

方法1：使用freqz（与MATLAB用法最接近） #

DSP包提供了和MATLAB同名的freqz函数，可直接传入分子b和分母a系数，用法完全对应：

using DSP

# 定义系统系数
a = [1, -1.1, 0.46]
b = [1, 0.04, 0.76]

# 获取频率响应h和对应归一化角频率w（范围0到π）
h, w = freqz(b, a)

方法2：构造滤波器对象后调用freqresp #

若需使用freqresp，需先通过系数创建IIR滤波器实例，再传入函数：

using DSP.Filters

a = [1, -1.1, 0.46]
b = [1, 0.04, 0.76]

# 构造IIR滤波器对象
filter_obj = IIRFilter(b, a)

# 获取频率响应，默认生成512个频率点
h, w = freqresp(filter_obj)

# 也可手动指定频率点，例如生成1024个从0到π的点
custom_w = range(0, π, length=1024)
h_custom = freqresp(filter_obj, custom_w)

错误原因说明 #

你传入的[b,a]被解析为一个嵌套数组，但freqresp没有针对这种输入的方法定义，因此抛出MethodError。必须先将系数转换为DSP包认可的滤波器对象，或直接使用更易用的freqz函数。

内容的提问来源于stack exchange，提问作者t.o.