在R中的另一個循環內部循環

Question

我在循環函數中存在循環結果問題。它只在循環內計數一次，爲第一個原始數據選擇最佳解決方案，然後停止。

我想記住矩陣zmienne的每一行的最佳解決方案。我究竟做錯了什麼？

schaffer <- function(xx) 
{x1 <- xx[1] 
x2 <- xx[2] 

fact1 <- (sin(x1^2-x2^2))^2 - 0.5 
fact2 <- (1 + 0.001*(x1^2+x2^2))^2 

y <- 0.5 + fact1/fact2 
return(y) 
} 

gradient_descent <- function(func, step, niter) { 

    N <- 3 #N- number of random points 
    zmienne <- matrix(runif(N*2, min = -100, max = 100), N, 2) 
    print(zmienne) 
    h = 0.001; 
    iter_count = 0; 

    for (i in 1:N) { 
    x_0 <- zmienne[i,] 
    x_n = x_0; 

    for (j in 1:niter) { 
     func_grad = (func(x_n+h) - func(x_n))/h; 
     if (abs(func_grad) < 0.0001) { break; } 
     x_n = x_n - step * func_grad; 
     print(x_n) 
     iter_count = iter_count + 1 
    } 
    } 
return(list(iterations = niter, best_value = func_grad, best_state = x_n, x0=x_0)) 
    } 

solution_m1 <- gradient_descent(schaffer, 0.1, 20) 
solution_m1 

Answer 1

返回不應該在內部循環內部，而是在函數結束時。

Answer 2

我想這是你想要什麼：

gradient_descent <- function(func, step, niter) { 

    N <- 3 #N- number of random points 
    zmienne <- matrix(runif(N*2, min = -100, max = 100), N, 2) 
    print(zmienne) 
    h = 0.001; 
    iter_count = 0; 
    best.vals <- NULL 
    for (i in 1:N) { 
    x_0 <- zmienne[i,] 
    x_n = x_0; 

    for (j in 1:niter) { 
     func_grad = (func(x_n+h) - func(x_n))/h; 
     if (abs(func_grad) < 0.0001) { break; } 
     x_n = x_n - step * func_grad; 
     print(x_n) 
     iter_count = iter_count + 1 
    } 
    best.vals <- c(best.vals, func_grad) 
    } 

    return(list(iterations = iter_count, best_value = best.vals, best_state = x_n, x0=x_0)) 
} 

solution_m1 <- gradient_descent(schaffer, 0.1, 20) 
solution_m1