Я пытаюсь решить проблему линейного программирования в R
с использованием пакета lpsolve
.Линейное программирование в R с использованием lpsolve
Вот проблема:
Вот пример в R для воспроизводимой например:
library("lpSolve")
a <- matrix(c(1,2,5,
1/2,1,3,
1/5,1/3,1),nrow=3,byrow=T)
#
f.obj <- c(1,0,0,0)
f.con <- matrix (c(
1,1,-a[1,2],0, #Contraint 1 for a12
1,-1,a[1,2],0, #Contraint 2 for a12
1,1,0,-a[1,3], #Contraint 1 for a13
1,-1,0,a[1,3], #Contraint 2 for a13
1,0,1,-a[2,3], #Contraint 1 for a23
1,0,-1,a[2,3], #Contraint 2 for a23
0,1,1,1, #Contraint 3
0,1,0,0, #Constraint 4
0,0,1,0, #Constraint 4
0,0,0,1 #Constraint 4
), nrow=10, byrow=TRUE)
f.dir <- c(rep("<=",6), "=",rep(">",3))
f.rhs <- c(rep(1,6),1,rep(0,3))
g <- lp ("max", f.obj, f.con, f.dir, f.rhs)
g$solution
я могу решить эту проблему вручную для небольшой проблемы, что, если У меня была 7 X 7
или n x n
матрица a
. Как бы я указал ограничение 1
и 2
, особенно я изо всех сил пытаюсь определить ограничение, связанное с [i, j]?
a = matrix(
c(1,4,9,6,6,5,5,
1/4,1,7,5,5,3,4,
1/9,1/7,1,1/5,1/5,1/7,1/5,
1/6,1/5,5,1,1,1/3,1/3,
1/6,1/5,5,1,1,1/3,1/3,
1/5,1/3,7,3,3,1,2,
1/5,1/4,5,3,3,1/2,1
),nrow = 7,byrow =T)
решение вышеуказанной матрицы 0.986 0.501 0.160 0.043 0.060 0.060 0.1 0.075
Любая помощь будет высоко оценили.
Привет, Я думаю, вы последняя строка в определении f.con не представляющих ограничение 4: 0,1,1,1 соответствует w_1 + w_2 + w_3. вместо этого он должен быть 0,1,0,0, [w_1> 0] \\ 0,0,1,0 [w_2> 0] \\ 0,0,0,1 [w_3> 0] – Cettt
Hi @Cettt да, я получил свое условие 4 неправильно, спасибо, указав это. – forecaster