求解50个非线性方程组时遇fn(par,...)参数缺失错误求助

首先，咱们来分析你遇到的错误：Error in fn(par, ...) : argument "n" is missing, with no default，这是因为你定义的fn需要两个参数x和n，但调用nleqslv的时候只传了初始猜测和函数，没有传递n参数。更关键的是，你的函数定义逻辑有问题——在循环里反复重写fn，最后得到的fn只会对应最后一次循环的i=49，根本没法生成49个对数方程。

咱们重新梳理问题：你需要求解50个变量（49个x[i]和1个n），对应50个方程：

• 49个方程：log(x[i], y[i]) = n（等价于log(x[i])/log(y[i]) = n，R的log函数底层是这么计算的）
• 1个约束方程：sum(x[1:49]) = 1

正确的做法是把所有变量放进同一个参数向量里，让fn接收这个长度为50的向量，然后计算每个方程的残差。下面是修正后的完整代码：

library(xlsx)
library(nleqslv)  # 别忘了加载这个包！你原来的代码没写

rm(list=ls())
setwd("C:/Users/.../folder")
my_data <- read.xlsx("samplefile.xlsx", 1)

# 处理y向量，注意要避免y[i]为0或1（否则对数无意义）
y <- numeric(49)
for(i in 1:49) {
  if(my_data[i,1] != 0 && my_data[i,1] != 1) {  # 新增判断避免无效底数
    y[i] = 1/my_data[i,1]
  } else {
    stop(paste("my_data[",i,",1]是0或1，无法作为对数底数"))
  }
}

# 定义残差函数：par是长度为50的向量，par[1:49]是x，par[50]是n
fn <- function(par) {
  x <- par[1:49]
  n <- par[50]
  
  # 计算49个对数方程的残差
  log_residuals <- log(x, base = y) - n
  # 计算求和约束的残差
  sum_residual <- sum(x) - 1
  
  # 返回所有残差，长度为50
  c(log_residuals, sum_residual)
}

# 设置初始猜测：前49个x初始值设为1/49（满足sum=1），n设为一个合理值比如1
guess <- c(rep(1/49, 49), 1)

# 求解方程组
result <- nleqslv(guess, fn)

# 查看结果
print(result)
# 提取解
x_solution <- result$x[1:49]
n_solution <- result$x[50]
cat("求解得到的n值：", n_solution, "\n")
cat("x向量的和：", sum(x_solution), "\n")  # 验证约束是否满足

关键修正点说明： #

1. 函数结构调整：把x和n合并成一个参数向量par，这样nleqslv可以一次性优化所有变量，不用单独传递n参数。
2. 避免无效底数：新增了对my_data[i,1]是否为0或1的判断，因为对数的底数不能是0、1或负数，否则会报错。
3. 初始猜测合理性：前49个x的初始值设为1/49，确保初始状态下满足求和约束，n设为1是一个合理的起始值（你可以根据实际数据调整）。
4. 加载必要包：别忘了加载nleqslv包，否则函数会找不到。

如果运行时还出现收敛问题，可以尝试调整nleqslv的参数（比如method选项），或者根据y的实际值调整初始猜测的n值，让它更接近真实解。

内容的提问来源于stack exchange，提问作者user11523191