1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63
|
public static double [] exponential(double cor1,double cor2,
double cor3, double cor4, double temps1, double temps2, double temps3,
double temps4, boolean cpm1, boolean cpm2, boolean cpm3,
boolean cpm4, int nombrepoints)
{
double pond = FactorielleN(nombrepoints)/(FactorielleN(2)*FactorielleN(nombrepoints -2));
double b2= 0;
double b1 = 0;
double b3= 0;
double b4= 0;
double b5= 0;
double b6= 0;
double b;
double a2= 0;
double a1 = 0;
double a3= 0;
double a4= 0;
double lna;
double a;
if (cpm1 && cpm2 ){
b1 = (Math.log(cor2)-Math.log(cor1))/(temps2 - temps1);
}
if (cpm1 && cpm3 ){
b2 = (Math.log(cor3)-Math.log(cor1))/(temps3 - temps1);
}
if (cpm1 && cpm4 ){
b3 = (Math.log(cor4)-Math.log(cor1))/(temps4 - temps1);
}
if (cpm2 && cpm3){
b4 = (Math.log(cor3)-Math.log(cor2))/(temps3 - temps2);
}
if (cpm4 && cpm2 ){
b5 = (Math.log(cor4)-Math.log(cor2))/(temps4 - temps2);
}
if (cpm3 && cpm4 ){
b6 = (Math.log(cor4)-Math.log(cor3))/(temps4 - temps3);
}
b= (b1+b2+b3+b4+b5+b6)/pond;
if (cpm1 ){
a1 = Math.log(cor1)-b*temps1;
}
if (cpm2){
a2 = Math.log(cor2)-b*temps2;
}
if (cpm3){
a3 = Math.log(cor3)-b*temps3;
}
if (cpm4){
a4 = Math.log(cor4)-b*temps4;
}
lna= (a1+a2+a3+a4)/nombrepoints;
a = Math.exp(lna);
double [] exp = {b, a};
return exp;
} |
Partager