例として、以下に正方行列「a」があります。行列 a は、nxn 正方行列です。
a = matrix(
c(1, 5 , 3, 7 , 3,
5, 1, 2, 2, 4,
3, 2 , 1, 2,4,
7, 2, 2,1,3,
2, 4,4 ,3 , 1
),ncol = 5,nrow =5)
これを最適化ルーチンに提供するために、以下のように R で function(x) を記述しようとしています。関数 (x) を最小化しようとしていますが、x は不明です。x はベクトルです。
sumx <- function(x) {
sum(((a[i,j]*a[j,k])-(x[i]/x[j]))^2) for all i,j,k such that i not eq to j not eq to k
}
このロジックと関数を R でプログラミングするのを手伝ってくれませんか?
とても有難い