首页 > 解决方案 > Octave中的高斯平滑

问题描述

我有大量的谱线,我想模拟这些谱线的高斯展宽谱。这些线是一个 N×2 矩阵,第一列是 x 坐标,y 是强度。所以最后我想在每个位置 x 放置一个带有比例因子 y 的高斯钟形曲线,并将所有这些的总和作为一个 M×2 矩阵,其中 M 除以范围就是分辨率。

N 大约为 1e6,数据点在 x 中不等距并且在 0:4000 的范围内,y 始终为正。

由于输入矩阵是在 Octave 中创建的,我也想在 Octave 中处理它,但我还没有找到合适的函数(或者我不明白如何正确使用它们)。

但是,我可以通过 gnuplot 获得所需的结果

plot "matrix.dat" smooth kdensity bandwidth 0.8

在我的帮助下,set table我还可以将其导出到文件中。不幸的是,这已经花费了大量时间来绘制一个情节,而且我有很多这样的矩阵。

有谁知道如何在 Octave 中加快这个过程?

标签: gnuplotoctavegaussiansmoothing

解决方案


推荐阅读