Question

我想将我的excel求解器模型转换为R中的模型。我需要找到3组坐标，以最小化与其他5个给定坐标的距离。我编写了一个程序，该程序计算距离矩阵，该矩阵输出从每个输入到给定坐标的最小距离。我想通过更改输入来最小化此功能。同理，我想找到使最小距离之和最小的坐标。我尝试了几种方法来执行此操作，请参见下面的代码（是的，我的距离矩阵函数可能有些笨拙，但这是因为为了运行某些算法，例如nloprt，我必须将输入减少为1变量（否则将得到警告））。我还看到了其他一些问题（例如GRG Non-Linear Least Squares (Optimization)），但它们没有更改/改进解决方案。

# First half of p describes x coordinates, second half the y coordinates # yes thats cluncky
p<-c(2,4,6,5,3,2) # initial points

x_given <- c(2,2.5,4,4,5)
y_given <- c(9,5,7,1,2)

f <- function(Coordinates){ 

# Predining
Term_1                       <-     NULL
Term_2                       <-     NULL
x                            <-     NULL
Distance                     <-     NULL
min_prob                     <-     NULL
l                            <-     length(Coordinates)
l2                           <-     length(x_given)
half_length                  <-     l/2
s                            <-     l2*half_length
Distance_Matrix              <-     matrix(c(rep(1,s)), nrow=half_length)

# Creating the distance matrix
for (k in 1:half_length){
  for (i in 1:l2){
    Term_1[i]                <-     (Coordinates[k]-x_given[i])^2
    Term_2[i]                <-     (Coordinates[k+half_length]-y_given[i])^2
    Distance[i]              <-     sqrt(Term_1[i]+Term_2[i])
    Distance_Matrix[k,i]     <-     Distance[i]
  }
}
d                            <-     Distance_Matrix

# Find the minimum in each row, thats what we want to obtain ánd minimize
for (l in 1:nrow(d)){
  min_prob[l] <- min(d[l,])
}
som<-sum(min_prob)
return(som)
}

# Minimise
sol<-optim(p,f)
x<-sol$par[1:3]
y<-sol$par[4:6]
plot(x_given,y_given)
points(x,y,pch=19)

然而，解决方案显然不是那么理想。我尝试使用nloptr函数，但不确定使用哪种算法。我可以使用哪种算法，或者可以使用/编程解决该问题的另一个函数？在此先谢谢您（对冗长的详细问题深表歉意）

Answer 1

查看optim的输出。它达到了迭代限制，尚未收敛。

> optim(p, f)
$`par`
[1] 2.501441 5.002441 5.003209 5.001237 1.995857 2.000265

$value
[1] 0.009927249

$counts
function gradient 
     501       NA 

$convergence
[1] 1

$message
NULL

尽管结果并没有太大不同，但您将需要增加迭代次数才能收敛。如果仍然不能接受，请尝试使用其他起始值。

> optim(p, f, control = list(maxit = 1000))
$`par`
[1] 2.502806 4.999866 5.000000 5.003009 1.999112 2.000000

$value
[1] 0.005012449

$counts
function gradient 
     755       NA 

$convergence
[1] 0

$message
NULL

GRG非线性R

1 个答案: