如何使用Julia中的Gadfly.jl绘制线性函数？

Question

我想知道如何使用Julia的Gadfly包绘制线性函数图形（比如y = 3x + 2）。我想到的一种方法是在该直线上绘制两个点，并添加Geom.line。

using Gadfly

function f(x)
    3 * x + 2
end

domains = [1, 100]
values = [f(i) for i in domains]
p = plot(x = domains, y = values, Geom.line)
img = SVG("test.svg", 6inch, 4inch)
draw(img, p)

有更好的绘制线条的方法吗？我在Gadfly文档的Functions and Expressions部分找到了相关内容。但是对于线性函数，我不确定该如何使用它？

- Jay Wong

Answer 1

大多数绘图工具只需要接受一个x值列表和一个y值列表即可。这样就可以创建一系列(x,y)点，然后根据您的要求绘制成个别点、通过直线相连或者拟合最佳模型等形式。

线性函数仅由两个点完全定义，因此使用一对点以及Geom.line选项应该就足够了。通常情况下，您需要使用linspace生成域中的点。然后将这些点传入您的函数以生成范围值：

julia> using Gadfly

julia> f(x) = 3x + 2
f (generic function with 1 method)

julia> domain = linspace(-2,2,100);

julia> range1 = f(domain);

julia> p = plot(x=domain, y=range1, Geom.line)

julia> img = SVG("test.svg", 6inch, 4inch)

julia> draw(img, p)

就像之前所说的那样，一条线由两个点定义，因此在这种情况下，我添加的所有额外点都是多余的。然而，更一般地说，您需要更多的点来准确地捕捉函数的形状：

julia> range2 = sin(domain);

julia> p2 = plot(x=domain, y=range2, Geom.line)

julia> img2 = SVG("test2.svg", 6inch, 4inch)

julia> draw(img2, p2)

- evan.oman

---

使用推导式似乎过于冗长：[f(x) for x in domain]; 更好的方法是使用 f.(domain)。 - DNF

好的，问题已解决。 - evan.oman

Answer 2

您可以直接绘制函数，将其作为参数传递给 plot 命令（正如您在 文档 中所指出的），然后在 X 轴（定义域）上跟随起始点和结束点：

julia> using Gadfly

julia> f(x) = 3x + 2
f (generic function with 1 method)

julia> plot(f, 1.0, 100.0)

甚至更复杂的表达式也可以用同样的方式绘制：

julia> f3(x) = (x >= 2.0) ? -2.0x+1.0 : 2.0x+1.0
f3 (generic function with 1 method)

julia> plot(f3, -5, 10)

如果您想在同一张图中绘制两个（或更多）函数，可以将函数数组作为第一个参数传递：

julia> plot([f, f3], -5, 10)

_{测试版本：Julia版本0.5.0}

- Gomiero