Accueil | Cast3M

   1 : $$$$ AGRE     NOTICE  FD218221  26/02/16    21:15:05     12436          
   2 :                                              DATE     26/02/16
   3 : 
   4 :       Objet :                                                            
   5 :       _______
   6 :                                                                          
   7 :       Cet operateur calcule une fonction d'agregation d'une liste de valeurs.
   8 : 
   9 : 
  10 :       Arguments d'entree :
  11 :       ____________________
  12 :                                                                          
  13 :       X      : Liste de nombres (type LISTREEL).
  14 : 
  15 :       MOT1   : Mot clef indiquant la fonction a appliquer, a choisir dans la liste suivante :
  16 : 
  17 :                -- Statistiques, normes, moyennes
  18 :                   'SOMM' : Somme
  19 :                            = Sum[i=1;n](x_i)
  20 :                   'PROD' : Produit
  21 :                            = Pro[i=1;n](x_i)
  22 :                   'MOYE' : Moyenne arithmetique (moment d'ordre 1)
  23 :                            = SOMM/n
  24 :                   'MOHA' : Moyenne harmonique
  25 :                            = n/(Sum[i=1;n](1/x_i))
  26 :                   'MOGE' : Moyenne geometrique
  27 :                            = PROD**(1/n)
  28 :                   'VARI' : Variance (moment centre d'ordre 2)
  29 :                            = (1/n) * Sum[i=1;n]((x_i-MOYE)**2)
  30 :                   'ECTY' : Ecart type
  31 :                            = VARI**0.5
  32 :                   'ASYM' : Coefficient d'asymetrie (moment centre reduit d'ordre 3)
  33 :                            = (1/n) * Sum[i=1;n](((x_i-MOYE)/ECTY)**3)
  34 :                   'KURT' : Kurtosis (moment centre reduit d'ordre 4)
  35 :                            = (1/n) * Sum[i=1;n](((x_i-MOYE)/ECTY)**4)
  36 :                   'MEDI' : Mediane
  37 :                   'PMOM' : Moment d'ordre P
  38 :                            = Sum[i=1;n](x_i**P)
  39 : 
  40 :                -- Normes/moyennes/maximums generalises (de parametre P)
  41 :                   (le mot clef 'ROBU' est disponible pour ces mots clefs)
  42 :                   'PMOY' : Moyenne generalisee d'ordre P
  43 :                            = (PMOM/n)**(1/P)
  44 :                   'PNOR' : Norme generalisee d'ordre P
  45 :                            = (Sum[i=1;n](|x_i|**P))**(1/P)
  46 :                   'LEHM' : Fonction de Lehmer d'ordre P
  47 :                            = Sum[i=1;n](x_i**P) / (Sum[i=1;n](x_i**(P-1)))
  48 :                   'KSL'  : Fonction de Kreisselmeir Steinhauser inferieure d'ordre P (ou MellowMax)
  49 :                            = (1/P) * ln( (1/n) * Sum[i=1;n](exp(P*x_i)) )
  50 :                   'KSU'  : Fonction de Kreisselmeir Steinhauser superieure d'ordre P (ou LogSumExp)
  51 :                            = (1/P) * ln(         Sum[i=1;n](exp(P*x_i)) )
  52 :                   'BOLT' : Fonction de Boltzmann d'ordre P
  53 :                            = (Sum[i=1;n](x_i*exp(P*x_i))) / (Sum[i=1;n](exp(P*x_i)))
  54 : 
  55 :       P      : Parametre d'exponentiation pour certaines fonctions (type FLOTTANT).
  56 : 
  57 :       'ROBU' : Mot clef, optionnel, pour utiliser une version "robuste" des fonctions parametrees.
  58 :                Ceci est utile lorsque les valeurs de X, ou de P, sont grandes et conduisent a des
  59 :                debordements (overflow) lors de l'application des fonctions puissance ou exponentielle.
  60 :                On procede en normalisant les valeurs x avec la norme infinie (VINF) ou bien avec le
  61 :                maximum VMAX :
  62 :                - P moyenne robuste :  PMOY(x,p) = VINF * PMOY(x/VINF,p)
  63 :                - P norme   robuste :  PNOR(x,p) = VINF * PNOR(x/VINF,p)
  64 :                - P LEHM    robuste :  LEHM(x,p) = VINF * LEHM(x/VINF,p)
  65 :                - KSL       robuste :  KSL(x,p)  = VMAX + KSL(x-VMAX,p)
  66 :                - KSU       robuste :  KSU(x,p)  = VMAX + KSU(x-VMAX,p)
  67 :                - BOLT      robuste :  BOLT(x,p) = VMAX + BOLT(x-VMAX,p)
  68 : 
  69 :       'DERI' : Mot clef, optionnel, pour calculer la liste des valeurs des derivees de la fonction
  70 :                par rapport a chaque variable de X.
  71 : 
  72 : 
  73 :       Arguments de sortie :
  74 :       _____________________
  75 : 
  76 :       F      : Valeur de la fonction d'agregation (type FLOTTANT).
  77 : 
  78 :       DFDX   : Liste des valeurs des derivees de F par rapport a chaque valeur de x (type LISTREEL).
  79 : 
  80 : 
  81 :       Remarques :
  82 :       __________
  83 : 
  84 :       1) La fonction PMOY (P moyenne) est prise egale a la moyenne geometrique MOGE si P=0.
  85 : 
  86 :       2) Certaines fonctions d'agregation necessittent de calculer de tres grands nombres,
  87 :          une prudence s'impose donc sur le choix de l'exposant P lorsque les valeurs de X
  88 :          sont elevees afin d'eviter les debordements de nombre flottant.
  89 :