von Mises使用gnuplot进行分发

Question

我需要将一些数据放到von Mises发行版中：

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

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

在gnuplot中，可以用函数

绘制Jo（k）

> besj0(k)

但我不能：

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

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

Answer 1

您可以使用Gnuplot的开发版本，该版本目前supports修改了第一类贝塞尔函数。具体而言，I0可用besi0。