如何使用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

我的第一个赌注是，因为频域(为了方便)最常(从您的图中也为您)显示为对数缩放。因此，如果你拟合，函数不会像你“想象”的那样拟合，而是缩放后再拟合。在对数尺度上，更高的值被表示得更多，->你的拟合度在那里更好。

因此，您应该做的是:找出应用了什么缩放，并尝试线性频率缩放。请记住，这也不是一个“好”的想法。因此，试着找到一个频率矢量，对于它，你需要接近并与之匹配。