matlab - 如何使用 invfreqs 将低频拟合到 Octave 音阶？

Question

我正在尝试使用 Octave 音程中的 invfreqs 从频率响应数据中找到传递函数。原则上它是有效的，问题是生成的传递函数总是拟合最高频率，低频匹配得很差。

尝试权衡拟合误差与频率的关系是行不通的。我做错了什么吗？

Hg = 10.^(mg/20).*exp(i*pg*pi/180);
wt(fgrps>1500) = 0;
m = 44;
n = 52;
[Bg,Ag] = invfreqs(Hg,fgrps,m,n,wt);

这是我得到的结果:

分子和分母多项式的阶次不同，结果或多或少相同。高频匹配好，低频匹配差。

对此我能做什么？

提前非常感谢您!

亲切的问候斯特凡

Answer 1

我的第一个赌注是，因为频域(为了方便起见)通常(并且从您的图中也适合您)以对数缩放显示。因此，如果您适合，该函数并不像您“想象”的那样适合，而是先缩放然后适合。在对数标度上，更高的值表示得更频繁 -> 那里的拟合效果更好。

所以你应该做的是:找出应用了什么缩放，并尝试线性频率缩放。请记住，这也不是一个“好”主意。因此，尝试找到一个频率向量，您需要接近并适应该频率向量。