我目前正在尝试在软件包pracma中使用功能fsolve。我想知道您是否可以帮助我解决一个小问题。
我有N个方程式的系统,带有一个未知参数(超定系统)。
让y是内生变量,x是解释变量,eta是未知参数。
等式为:Y = exp (X+eta)/(1+exp(X+eta))
Length_Vector <- 8
功能
F <- function(eta) {
y <- c(1,1,0,0,0,1,0,1)
x <- c(5,4,6,8,4,6,15,14)
x_modified <- c(x) + eta
return(c(exp(x_modified)/(1+exp(x_modified)) - y))
}
矩阵初始化
eta0 <- matrix(0, Length_Vector , 1)
计算解决方案
Solution <- matrix(fsolve(F, eta0, maxiter=1000, tol=0.01)x, nrow = Length_Vector , ncol = 1)
当我计算此程序时,该程序运行,但是我获得了eta值不同的向量。该函数似乎可以解决由N个方程和N个未知参数(而不是1个未知参数)组成的系统。
如果您能帮助我解决这个问题,我将不胜感激。