我一直在尝试解决R中使用constrOptim()(我的第一次)的约束优化问题,但我正在努力为我的问题设置约束。约束优化R建立约束
的问题是非常简单的,我可以设置功能确定,但在大约经过约束损失感到有点。
例如我已经定义的问题是(我打算以N固定在1000开始说,所以我只想解决X最终我想选择N和X的最大利润):
所以我可以设置功能为:
fun <- function(x, N, a, c, s) { ## a profit function
x1 <- x[1]
x2 <- x[2]
x3 <- x[3]
a1 <- a[1]
a2 <- a[2]
a3 <- a[3]
c1 <- c[1]
c2 <- c[2]
c3 <- c[3]
s1 <- s[1]
s2 <- s[2]
s3 <- s[3]
((N*x1*a1*s1)-(N*x1*c1))+((N*x2*a2*s2)-(N*x2*c2))+((N*x3*a3*s3)-(N*x3*c3))
}
我需要实现的约束是:
x1>=0.03
x1<=0.7
x2>=0.03
x2<=0.7
x3>=0.03
x2<=0.7
x1+x2+x3=1
的X这里表示成我需要最优地分配N个区段,所以N个X1 = pecent放置在桶1等,每个桶具有至少3%但不超过70%。
非常感谢任何帮助...
例如,这里是我用来测试功能的一个例子,我想要什么:
fun <- function(x, N, a, c, s) { ## profit function
x1 <- x[1]
x2 <- x[2]
x3 <- x[3]
a1 <- a[1]
a2 <- a[2]
a3 <- a[3]
c1 <- c[1]
c2 <- c[2]
c3 <- c[3]
s1 <- s[1]
s2 <- s[2]
s3 <- s[3]
((N*x1*a1*s1)-(N*x1*c1))+((N*x2*a2*s2)-(N*x2*c2))+((N*x3*a3*s3)-(N*x3*c3))
};
x <-matrix(c(0.5,0.25,0.25));
a <-matrix(c(0.2,0.15,0.1));
s <-matrix(c(100,75,50));
c <-matrix(c(10,8,7));
N <- 1000;
fun(x,N,a,c,s);
所以...'X [1:3],'是变量,而一个'[1:3]','S [1:3]'和'C [1:3]'中给出值?你能否详细说明你需要的配方?我不清楚... – digEmAll
是的 - a [1:3],s [1:3],c [1:3]代表历史平均激活率,支出,购置成本..我试图选择x [1:3]来最大限度地发挥 – andrewm4894