假设的测量数据可能包含很宽的频率范围,我们通过它拟合阶多项式曲线。这样做会过滤掉某些本来会出现的频率。正如我可以想象的那样,它就像一个低通滤波器。1secnth
是的。和不。
我想您指的是非常、非常、非常低阶(相对于数据集的长度)的多项式。由于多项式拟合的工作方式,您很快就会发现自己在处理非常小或非常大的数字,并且数值不稳定性逐渐蔓延。44.1KHz 的一秒数据是 44100 个数据点。拟合 1024 阶多项式仅涵盖大约 2% 的原始数据点,并且仍需要将范围内的数字提高到 1024 次方。0…1
此外,请注意,当您转移到下一个插值窗口时(新的一秒数据进入),您现在必须考虑约束。也就是说,插值的开始现在不能随意移动。曲线必须从最后一条开始,以确保连续性……否则,当波形在过渡处跳跃时,您会得到“咔哒”声。事实上,就多项式拟合而言,您肯定会得到点击,因为曲线将以最小二乘法规定的任何斜率偏离,以最大限度地减少误差,因为您无法知道未来,所以很难强制约束(即“无论如何,我希望曲线在具有特定斜率的直线段处结束”)。
通过它拟合抛物线可以近似为正弦曲线的一半
不。
其余的正弦曲线在哪里?你将如何处理“中间”的抛物线?也就是说,随着峰值滑出窗口,我们现在进入低谷。抛物线可以拟合正弦曲线的一部分,但不能替代它。
这种曲线拟合充当低通滤波器。这种低通滤波器的特点是什么?
问这个,就是问,我如何用另一个来表达一个?也就是说,我可以在多项式拟合和低通滤波之间找到一个等价的三角函数之和吗?答案是否定的,因为这两者的结构截然不同。
要考虑的另一件事是最小二乘的工作方式,因为即使数据不存在,最小二乘也会努力将抛物线拟合到数据中。那么,如果我的没有像抛物线一样弯曲呢?事实上,如果您的信号不包含随多项式函数的某种组合而变化的分量,则拟合将失败。一个典型的例子是脉冲函数。拍一个包含低音鼓(低频)和高帽(高频)的节拍,并尝试使用多项式拟合“剪切”高帽。不可能。多项式拟合将尝试理解所有内容,包括节拍之间的静默。多项式必须适合。最小二乘必须找到局部最小值。n,x(n)
这里的一个例外似乎是分段样条插值,它将波形分成多个部分并拟合更小的多项式,在它们之间有约束,但同样,它们的定义不允许在它们的样条表示和傅里叶变换之间轻松转换,您可以通过它然后可以“跳”到脉冲响应。也就是说,给定拟合分段样条的系数,您可以找到一种方法来导出有限脉冲响应滤波器的脉冲响应(更不用说无限脉冲响应滤波器了)。
您始终可以尝试获取大量数据样本,拟合平滑样条曲线,然后获得其傅里叶变换,以查看哪种低通滤波器可以近似样条曲线的结果,但这不是推导等效滤波器的方法。
希望这可以帮助。