gnuplot - 使用 gnuplot 的 von Mises 分布

Question

我需要将一些数据拟合到 von Mises 分布中：

但在此之前，我想先在 gnuplot 中绘制这个函数。我想重现不同 k 值的曲线，您可以在维基百科中看到它们：https ://en.wikipedia.org/wiki/Von_Mises_distribution 。

我的问题依赖于术语 Io(k)，它是 0 阶的修正贝塞尔函数。Io(k)=Jo(ik)，其中 Jo 是 n=0 的贝塞尔函数，基本上来自这个等式：。

在 gnuplot 中，可以使用函数绘制 Jo(k)

> besj0(k)

但是我不能：

besj0(i*k) when the argument is imaginary.

有人知道如何在 gnuplot 中绘制 von Mises 分布吗？

Answer 1

可以使用 Gnuplot 的开发版，目前支持第一类修改后的 Bessel 函数。具体来说，I0可作为besi0.