/*************************************************************************************************************************************/
/* */
/* C A L C U L D E L ' E N T R O P I E D ' U N F I C H I E R N U M E R I Q U E : */
/* */
/* */
/* Author of '$xrv/entropie.01$K' : */
/* */
/* Jean-Francois COLONNA (LACTAMME, 20041016213745). */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* I N T E R F A C E ' listG ' : */
/* */
/* */
/* :Debut_listG: */
/* :Fin_listG: */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D I R E C T I V E S S P E C I F I Q U E S D E C O M P I L A T I O N : */
/* */
/*************************************************************************************************************************************/
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* F I C H I E R S D ' I N C L U D E S : */
/* */
/*************************************************************************************************************************************/
#include INCLUDES_BASE
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* V A L E U R S I M P L I C I T E S D E S P A R A M E T R E S : */
/* */
/*************************************************************************************************************************************/
#define CALCULER_L_ENTROPIE \
VRAI \
/* Indique s'il faut calculer l'entropie ('VRAI') ou l'information ('FAUX'). */
#define UTILISER_SLOGX_PLUTOT_QUE_LOGX \
FAUX \
/* Indique s'il faut utiliser 'SLOGX(...)' ('VRAI') ou 'LOGX(...)' ('FAUX'). Ceci a ete */ \
/* introduit le 20160704154227... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D U F I C H I E R : */
/* */
/*************************************************************************************************************************************/
#include xrv/ARITHMET.1d.I"
/* Passage a l'allocation dynamique le 20060214190858... */
#define NOMBRE_DE_DECIMALES \
SEIZE
#include xrv/ARITHMET.21.I"
#include xrv/champs_5.41.I"
#define VALEUR_IMPLICITE \
FZERO
gGENERATION_D_UN_FICHIER(fichier_des_valeurs,liste_des_valeurs);
/* Definition en memoire du fichier. */
#define ELEMENT_DU_FICHIER(index) \
gELEMENT_DU_FICHIER(liste_des_valeurs,index) \
/* Acces a un element courant du fichier. */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* D E F I N I T I O N D E S E D I T I O N S : */
/* */
/*************************************************************************************************************************************/
#define FORMAT_EFFECTIF_D_EDITION(message) \
Cara(chain_Aconcaten5(message,INTRODUCTION_FORMAT,valeurs_signees,".*",format_d_edition)) \
/* Format effectif de certaines editions... */
/*===================================================================================================================================*/
/*************************************************************************************************************************************/
/* */
/* C A L C U L D E L ' E N T R O P I E D ' U N F I C H I E R N U M E R I Q U E : */
/* */
/*************************************************************************************************************************************/
BCommande(nombre_d_arguments,arguments)
/*-----------------------------------------------------------------------------------------------------------------------------------*/
Bblock
DEFV(Logical,INIT(calculer_l_entropie,CALCULER_L_ENTROPIE));
/* Indique s'il faut calculer l'entropie ('VRAI') ou l'information ('FAUX'). */
DEFV(Logical,INIT(utiliser_SLOGX_plutot_que_LOGX,UTILISER_SLOGX_PLUTOT_QUE_LOGX));
/* Indique s'il faut utiliser 'SLOGX(...)' ('VRAI') ou 'LOGX(...)' ('FAUX'). Ceci a ete */
/* introduit le 20160704154227... */
DEFV(Float,INIT(entropie,FZERO));
/* Entropie du fichier... */
#include xrv/ARITHMET.22.I"
#include xci/valeurs.03.I"
/*..............................................................................................................................*/
#include xrv/champs_5.1A.I"
/* Ceci fut introduit le 20070103173646... */
GET_ARGUMENTS_(nombre_d_arguments
,BLOC(PROCESS_ARGUMENT_I("nombre_elements=""ne=",nombre_d_elements
,BLOC(VIDE;)
,BLOC(Bblock
PRINT_AVERTISSEMENT("'ne=' doit etre defini avant toute entree de fichiers");
Eblock
)
);
PROCESS_ARGUMENTS_DE_DEFINITION_DES_FICHIERS_01;
PROKESF_ARGUMENT_FICHIER("fichier="
,fichier_des_valeurs
,liste_des_valeurs
,VALEUR_IMPLICITE
,lTRANSFORMAT_0d
,iGENERATION_D_UN_FICHIER
);
GET_ARGUMENT_L("entropie=""H=",calculer_l_entropie);
GET_ARGUMENT_N("information=""I=",calculer_l_entropie);
GET_ARGUMENT_L("SLOGX=",utiliser_SLOGX_plutot_que_LOGX);
GET_ARGUMENT_N("LOGX=",utiliser_SLOGX_plutot_que_LOGX);
/* Parametres introduits le 20160704154227... */
PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_3;
/* Cette procedure fut introduite le 20070103173646... */
PROCESS_ARGUMENTS_DE_PARAMETRAGE_DE_LA_GENERATION_DE_SUITE_DE_VALEURS_1;
/* Cette procedure fut introduite le 20061226193427... */
)
);
DoIn(index
,PREMIER_ELEMENT_D_UN_FICHIER
,DERNIER_ELEMENT_D_UN_FICHIER
,I
)
Bblock
DEFV(Float,INIT(valeur_courante,ELEMENT_DU_FICHIER(index)));
/* Recuperation de la valeur courante dans le fichier. */
Test(IZEQ(valeur_courante))
Bblock
/* Dans le cas de la nullite, on considere que : */
/* */
/* 0.log(0) = 0 */
/* */
/* d'apres le cours sur l'entropie... */
Eblock
ATes
Bblock
Test(IL_FAUT(utiliser_SLOGX_plutot_que_LOGX))
Bblock
INCR(entropie,MUL2(valeur_courante,SLOGX(INVE(valeur_courante))));
/* Possibilite introduite le 20160704154227... */
Eblock
ATes
Bblock
Test(IZGT(valeur_courante))
Bblock
INCR(entropie,MUL2(valeur_courante,LOGX(INVE(valeur_courante))));
/* Jusqu'au 20041021161307, l'entropie etait calculee par : */
/* */
/* NEGA(MUL2(valeur_courante,LOGX(valeur_courante))) */
/* */
/* mais il est plus logique d'utiliser l'inverse ainsi que cela est explique plus bas... */
Eblock
ATes
Bblock
PRINT_ERREUR("le logarithme neperien d'un nombre negatif est demandee");
Eblock
ETes
Eblock
ETes
Eblock
ETes
Eblock
EDoI
Test(IL_FAUT(calculer_l_entropie))
Bblock
EDITION_DANS_gOPERATION_SUR_LES_FICHIERS(entropie);
/* Edition de l'entropie 'H' : */
/* */
/* -- */
/* \ */
/* H = - / X .log(X ) */
/* -- i i */
/* */
/* ou mieux : */
/* */
/* -- */
/* \ 1 */
/* H = / X .log(----) */
/* -- i X */
/* i */
/* */
/* car, en effet, X(i) est en general une probabilite et alors 1/X(i) est un nombre */
/* d'evenements et ainsi, en utilisant log(1/X(i)) on est plus proche de la definition */
/* physique de l'entropie... */
Eblock
ATes
Bblock
EDITION_DANS_gOPERATION_SUR_LES_FICHIERS(SOUS(LOGX(FLOT(nombre_d_elements)),entropie));
/* Edition de l'information 'I' : */
/* */
/* I = log(Cardinal(fichier)) - H */
/* */
Eblock
ETes
CAL2(Prin0("\n"));
lGENERATION_D_UN_FICHIER(liste_des_valeurs,VALEUR_IMPLICITE);
RETU_Commande;
Eblock
ECommande