Решить и построить неявное уравнение в R - PullRequest
0 голосов
/ 27 февраля 2020

Я пытаюсь решить это уравнение ниже (см. Изображение) для «t» в R, но застрял в сценариях. Я пробовал некоторые функции, такие как uni root, оптимизировать и оптимизировать, но все еще мог получить что-то. Любая помощь будет приветствоваться. введите описание изображения здесь

Here are the equation and my lines of codes.


#soil parameters
Ks= 17000;
Beta= 1.3;
Eta= 4.5;
Theta_s= 0.33;
Theta_r= 0.03;
Psi_Water= -80;
Delta=Eta/Beta

#crop parameters
Tp= 5.5;
EC50= 7.5505;
p= 3;
Psi_Root=-6000;
b= 10;
Yr0= 0;

#Water Salinity (EC)
EC <- seq(0.5,6, 1)

#Irrigigation water
I <-seq(0.5, 12, by =.1)
#the function created for the equation
fct <- function (t){-t + (min(Tp,(Psi_Root-(Psi_Water/((I-t)/Ks)^(1/Eta)))(I-t)b)/
       (1+((ECI(Theta_r+(Theta_s-Theta_r)((I-t)/Ks)^(1/Delta)))/(EC50(I-t)*Theta_s))^p))}

for (i in 1:length(EC)) {
for (j in 1:length(I)){
ans <- optim (fct( Tp, EC50, p, Psi_Root, b, Ks, Eta, Theta_s, Theta_r, Psi_Water, Delta, 
  I[j], EC[i],t), c(0.1, 13))
  }
}
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...