从高斯拟合的快速傅里叶变换中查找频率误差

计算科学 算法 傅立叶分析
2021-12-22 14:40:31

我在一个由频率范围内的频率扫描生成的盒子中对声音进行了 FFT,并有一组频率及其相应的 FFT 幅度。根据声波模型,某些频率应该是共振频率,这在数据中肯定会显示为类似高斯的峰值。

我想知道是否有一些程序或算法可以快速确定光谱中的峰值及其标准偏差(如果假设为高斯)。

我已经有一种蛮力的做法,那就是在 Igor Pro 上检查每个点及其值,并找到f使得频谱幅度,A=A0e1/2在哪里A0是高峰。

这是因为高斯方程是

A0ef22σf
在哪里σf是标准差f,所以在f=σf,指数为e1/2.

我想知道是否有一种方法可以自动执行此操作。寻找峰值本身可能涉及某种离散导数来找到极值(由于噪声波动,这将不得不忽略转折点)。然后也许我可以比较峰值范围内的值,但当然比较不会完全相等,而是找到最接近的值,因为我正在处理离散数据。

0个回答
没有发现任何回复~