Commit 2f2c45c5 authored by pfreon's avatar pfreon
Browse files

Improvement of transitional state estimation of V

Improvement of transitional state estimation of V for models with effect on catchability or catchability and abundance.
parent dd24534b
......@@ -91,7 +91,7 @@ public class EquationModele {
case 33: // CPUE = a V exp(bV)
puec = par[0] * vba * Math.exp(par[1] * vba);
break;
case 14: // CPUE=a+b.V+c.(a+b.V)^2.E MODELE CAPTURABILITE
case 14: // CPUE=a+b.V+c.(a+b.V)^2.E MODELE CAPTURABILITE #22
puec = par[0] + par[1] * vv + par[2] * (par[0] + par[1] * vv) * (par[0] + par[1] * vba) * ff;
break;
case 15: // CPUE=a.V^b+c.V^(2.b).E MODELE CAPTURABILITE
......@@ -151,12 +151,12 @@ public class EquationModele {
break;
case 26: // CPUE=a.V^(b+c)+d.V^(2.b).E MODELE CAPTURABILITE + ABONDANCE
tmp1 = Math.sqrt(vba * vv);
puec = par[0] * Math.pow(tmp1, par[1] + par[2]) + par[3] * Math.pow(vba, 2 * par[1]) * ff;
puec = par[0] * Math.pow(tmp1, par[1] + par[2]) + par[3] * Math.pow(vv*vba, par[1]) * ff;
break;
case 27: // CPUE=a.V^(1+b)+c.V^(2+b)+d.V^(2.b).E MODELE LINEAIRE QUADRATIQUE CAPTURABILITE + ABONDANCE.
tmp1 = Math.pow(vv, par[1]);
tmp2 = Math.pow(vba, par[1]);
puec = par[0] * tmp1 * vba + par[2] * (vba * vba) * tmp1 + par[3] * tmp1 * tmp2 * ff; // Formulation vérifiée (X^(a+b) = X^a . X^b)
puec = par[0] * tmp1 * vba + par[2] * (vv * vba) * tmp1 + par[3] * tmp1 * tmp2 * ff; // Formulation vérifiée (X^(a+b) = X^a . X^b)
break;
case 28: // CPUE=a.V^b.exp(c.V^d.E) MODELE ABONDANCE
puec = par[0] * Math.pow(vba, par[1]) * Math.exp(par[2] * Math.pow(vba, par[3]) * ff);
......@@ -165,11 +165,11 @@ public class EquationModele {
tmp1 = Math.sqrt(vba * vv);
puec = par[0] * Math.pow(tmp1, par[1]) * Math.exp(par[2] * Math.pow(vba, par[3]) * ff);
break;
case 30: // CPUE=(a.V^(1+b)+c.V^(2+b)).exp(d.V^b.E) MODELE CAPTURABILITE + ABONDANCE EXP. QUADRATIQUE
case 30: // CPUE=(a.V^(1+b)+c.V^(2+b)).exp(d.V^b.E) MODELE EXP. QUADRATIQUE-PUISSANCE CAPTURABILITE + ABONDANCE
tmp1 = Math.pow(vv, par[1]);
tmp2 = Math.pow(vba, par[1]);
// AVANT: puec = (par[1] * tmp1 * vba + par[0] * tmp1 * (vba * vba)) * Math.exp(par[3] * tmp2 * ff);
puec = (par[0] * tmp1 * vba + par[2] * (vba * vba) * tmp1) * Math.exp(par[3] * tmp2 * ff); // Formulation vérifiée (X^(a+b) = X^a . X^b)
puec = (par[0] * tmp1 * vba + par[2] * (vv * vba) * tmp1) * Math.exp(par[3] * tmp2 * ff); // Formulation vérifiée (X^(a+b) = X^a . X^b)
break;
}
return (puec);
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment