如何在Julia中使用Gadfly.jl绘制线性函数?

时间:2016-12-27 22:41:17

标签: plot graph julia gadfly

我想知道如何使用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)

enter image description here

有没有更好的方法来绘制线条?我在Gadfly文档中找到了Functions and Expressions部分。我不确定如何将它与线性函数一起使用?

2 个答案:

答案 0 :(得分:3)

大多数绘图工具只接受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)

Linear Plot

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

julia> range2 = sin(domain);

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

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

julia> draw(img2, p2)

Non-Linear Plot

答案 1 :(得分:3)

您可以直接绘制函数,将其作为参数传递给plot命令(如您指向documentation),然后是X轴的起点和终点(域):

julia> using Gadfly

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

julia> plot(f, 1.0, 100.0)

enter image description here

更复杂的表达式可以用同样的方式绘制:

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)

enter image description here

如果要在同一个图形上绘制两个(或更多)函数,可以将一组函数作为第一个参数传递:

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

enter image description here

测试:Julia Version 0.5.0