运筹学线性规划问题——R语言代码
library(“lpSolve”)
f = c(3, 5)
a1 = c(1, 0)
a2 = c(0, 2)
a3 = c(3, 2)
a = rbind(a1,a2,a3)
a.dir = rep("<=", 3)
a.rhs = c(4, 12, 18)
solution = lp(“max”, f, a, a.dir, a.rhs, compute.sens = 1)
solution
s
o
l
u
t
i
o
n
s
o
l
u
t
i
o
n
solution solution
solutionsolutionduals #影子价格
solutionKaTeX parse error: Expected 'EOF', got '#' at position 19: …s.coef.from #̲f=c(3,5),在产品2的利…sens.coef.to #产品1的利润在(0,7.5)之间,最优解的答案不变
solutionKaTeX parse error: Expected 'EOF', got '#' at position 19: …ls.from #̲第二种产品的利润在(6,18)…duals.to
#all.bin = TRUE 0-1规划
#all.int = TRUE 整数规划
目标规划
p1 = 100000000000
p2 = 100000000
p3 = 10000
p4 = 1
library(“lpSolve”)
f = c(0, 0, p4, p1, p2, 0, 0, 5p3, 0, 3p3)
a1 = c(1, 2, -1, 1, 0, 0, 0, 0, 0, 0)
a2 = c(1, 2, 0, 0, -1, 1, 0, 0, 0, 0)
a3 = c(1, -2, 0, 0, 0, 0, -1, 1, 0, 0)
a4 = c(0, 1, 0, 0, 0, 0, 0, 0, -1, 1)
a = rbind(a1, a2, a3, a4)
a.dir = rep("=", 4)
a.rhs = c(6, 9, 4, 2)
solution_2 = lp(‘min’, f, a, a.dir, a.rhs)
solution_2$solution