Question

我正在尝试设置“求解器”功能来查找正常的通道深度（yn）。参数在下面的代码中给出，我可以在其中估算方程式的一侧。所有其他参数都是yn的功能。我需要找到解决函数A *（R ^（2/3）= nQSo.yn的yn。

So=0.001                                       
n=0.013                                        
Q=30  
B=10 

nQSo=(n*Q)/(So^(1/2))

A=B*yn
P=B+2*yn
R=A/P

A*(R^(2/3)=nQSo

Answer 1

您可以看看optimize

So=0.001                                       
n=0.013                                        
Q=30  
B=10 

nQSo=(n*Q)/(So^(1/2))

error = function(yn,nQSo){
  A=B*yn
  P=B+2*yn
  R=A/P

  return(abs(A*(R^(2/3))-nQSo))
}

optimize(error,interval = c(0,2),nQSo = nQSo)

您看到的结果是yn = 1.239066

用规划求解函数查找变量值

1 个答案: