GNUPlot 中的第二类修正贝塞尔函数答案

【问题标题】：Modified Bessel function of second kind in GNUPlotGNUPlot 中的第二类修正贝塞尔函数
【发布时间】：2013-09-02 15:17:39
【问题描述】：

我必须在 GNUPlot 中绘制K-Distribution，为此我需要 GNUPlot 中的第二类伽玛函数和修改后的贝塞尔函数。 GNUPlot 有一个内置的 gamma 函数，但似乎没有第二类修改后的 bessel 函数。如果可能的话，我希望它的功能尽可能类似于MATLAB 中的解决方案。我已经尝试寻找解决方案，但至今仍找不到它（或贝塞尔函数的 GNUplot 版本）。

如果这是错误的地方，请引导我到更合适的地方。

【问题讨论】：

不，无法直接在gnuplot 中计算它们，但您可以使用任何外部工具生成数据文件并绘制它。 Octave 有例如besselk function。然后做plot '< octave script.oct'。
非常感谢。我会这样做的！

标签： matlab gnuplot

【解决方案1】：

无法直接在gnuplot中计算，但可以使用外部工具生成数据文件并绘制。使用例如以下脚本callbesselk.oct:

# callbesselk.oct start inc stop alpha
#
x_start = str2num(argv(){1});
x_inc = str2num(argv(){2});
x_end = str2num(argv(){3});
alpha = str2num(argv(){4});
x = x_start:x_inc:x_end;
[k, kerr] = besselk(alpha, x);
out = [x; k];
printf("%f %f\n", out);

并用gnuplot绘制数据：

plot '< octave -q callbesselk.oct 0 0.01 10 2'

【讨论】：

【解决方案2】：

这些特定的贝塞尔函数在 gnuplot 中不可用； gnuplot 有两个patches，允许您使用来自外部库的函数。

【讨论】：