0

我解决了一个线性双目标混合整数问题，我想绘制结果。结果包括线和点。例如

list=[([0.0; 583000.0], 0), ([190670.0; 149600.0], 0), ([69686.0, 385000.0], 1), ([33296.0, 484000.0], 1), ([136554.0, 2.38075e5], 1), ([24556.0, 503800.0], 0), ([47462.0, 437800.0], 1), ([129686.0, 253000.0], 1), ([164278.0, 178200.0], 1)]

在这个列表中，第三个点的([69686.0, 385000.0], 1)第二个元素1被确定为通过一条线连接到先前点的这个点通过一条线([190670.0; 149600.0], 0)连接到第二个点。

我编码如下：

using  JuMP,Plots

list=[([0.0, 583000.0], 0), ([24556.0, 503800.0], 0), ([33296.0, 484000.0],1), ([47462.0, 437800.0], 1), ([69686.0, 385000.0], 1), ([129686.0, 253000.0], 1), ([136554.0, 23805.0], 1), ([164278.0, 178200.0], 1), ([190670.0, 149600.0], 0)]

x=zeros(1,1)
for i=1:size(list,1)
    x=[x;list[i][1][1]]
end
row=1
x = x[setdiff(1:end, row), :]

y=zeros(1,1)
for i=1:size(list,1)
    y=[y;list[i][1][2]]
end
row=1
y = y[setdiff(1:end, row), :]

for i=2:size(list,1)
    if list[i][2]==0
        plot(Int(x[i]),Int(y[i]),seriestype=:scatter)
        plot(Int(x[i+1]),Int(y[i+1]),seriestype=:scatter)
    end
    if list[i][2]==1
        plot(Int(x[i]),Int(y[i]))
        plot(Int(x[i+1]),Int(y[i+1]))
    end
end

但它不起作用。你能帮帮我吗？谢谢

4

1 回答 1

1

您可以简单地将每个线段的 x 和 y 值推送到两个单独的数组中，x并y在下面的代码中。在每个线段的值（即 x1 和 x2 或 y1 和 y2）之后，将 aNaN放入数组中。如果不应该有连接，这将防止将线段连接到下一个线段。（例如，您看到 1 然后是 0 的情况）。最后plot(x, y)。

下面的代码片段可以做到这一点。请注意，allx和ally用于保存所有点，无论连接状态如何。您可能希望从中排除连接点。x并y保持连接的线段。

using Plots
x, y, allx, ally = Float64[], Float64[], Float64[], Float64[]

# iterate through list
for i = 1:length(list)-1
    if list[i+1][2] == 1
        # push x1 from the first point, x2 from the second and a `NaN`
        push!(x, list[i][1][1], list[i+1][1][1], NaN)
        # push y1, y2, `NaN`
        push!(y, list[i][1][2], list[i+1][1][2], NaN)
    end
    push!(allx, list[i][1][1])
    push!(ally, list[i][1][2])
end
push!(allx, list[end][1][1])
push!(ally, list[end][1][2])
# scatter all points
scatter(allx, ally)
# plot connections with markers
plot!(x, y, linewidth=2, color=:red, marker=:circle)

这应该有希望给你你想要的情节。

如果你碰巧使用Gadfly.jl而不是Plots.jl，你可以得到一个类似的情节

using Gadfly
connectedpoints = layer(x=x, y=y, Geom.path, Geom.point, Theme(default_color="red"))
allpoints = layer(x=allx, y=ally, Geom.point)

plot(connectedpoints, allpoints)

附带说明，如果您打算在已创建的绘图对象之上绘制另一个系列，则应plot!使用plot.

于 2019-05-01T09:37:23.223 回答