如何使用invfreqs在八度音阶中拟合低频？

Question

我试图使用八度音阶中的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

我的第一个赌注是，因为频域（为方便起见）最常见（并且也是你的情节），以对数缩放显示。因此，如果你适合，那么这个功能并不像你想象的那样适合，而是缩放然后适合。在对数标度上，更频繁地表示更高的值 - ＆gt;你的适合度更高。

所以你应该做的是：找出应用的缩放比例，并尝试线性频率缩放。请记住，这也不是一个好的＆＃34;理念。因此，请尝试找到一个频率向量，您需要将其关闭并适合该频率向量。