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 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123
| #include <stdio.h>
#include <stdlib.h>
#include <math.h>
int writefich ( double *t, int n, char *Name )
{
FILE *file=(FILE *)NULL ;
int i ;
if ( (file = fopen(Name,"w")) != NULL )
{
for ( i = 0 ; i < n ; i++ )
fprintf ( file, "\n%d\t%g",i,t[i]);
fclose ( file );
}
return 0 ;
}
int main(int argc, char *argv[])
{
int i, i1, n, low, high, filtre=2, ntemp=20 ;
double temp[20]={20.0,48.0,42.3,39.9,39.0,41.4,44.0,52.0,53.1,49.7,42.0,40.1,100.0,45.0,50.0,47.0,43.1,41.0,39.0,42.2};
double moy[20]={0.0},ecart[20]={0.0},moyecart[20]={0.0}, moylocale ;
double tempresult[20], vrailiss[20]={0.0};
/*
* Calcul moyenne et ecart sur 5
*/
for ( i = 0 ; i < ntemp ; i++ ) {
low = i - filtre ;
high = i + filtre ;
if ( low < 0 ) low = 0 ;
if ( high >= ntemp ) high = ntemp ;
n = 0 ;
for ( i1 = low ; i1 < high ; i1++ ) {
moy[i] = moy[i] + temp[i1] ;
n = n + 1 ;
}
moy[i] = moy[i] / (double)n ;
ecart[i] = moy[i] - temp[i] ;
}
/*
* Calcul moyenne des ecarts et remplacement des pics aberrants
*/
for ( i = 0 ; i < ntemp ; i++ ) {
low = i - filtre ;
high = i + filtre ;
if ( low < 0 ) low = 0 ;
if ( high >= ntemp ) high = ntemp ;
n = 0 ;
for ( i1 = low ; i1 < high ; i1++ ) {
moyecart[i] = moyecart[i] + ecart[i1] ;
n = n + 1 ;
}
moyecart[i] = moyecart[i] / (double)n ;
/* Test en fonction de la distance entre l'ecart et la moyenne des ecarts */
if ( fabs(ecart[i]) < fabs(moyecart[i]) )
tempresult[i] = temp[i] ;
else {
/* Nouveau calcul de la moyenne sans le point central */
/* UNIQUEMENT SI LES ECARTS D'UN COTE ET DE L'AUTRE SONT DANS LE MEME SENS */
if ( i != 0 && i != (ntemp-1) ) {
if ( (ecart[i]*ecart[i-1] < 0.0) && (ecart[i-1]*ecart[i+1] > 0.0)) {
moylocale = 0.0 ;
n = 0 ;
for ( i1 = low ; i1 < high ; i1++ ) {
if ( i1 != i ) {
moylocale = moylocale + temp[i1] ;
n = n + 1 ;
}
}
moylocale = moylocale / (double)n ;
tempresult[i] = moylocale ;
}
else
tempresult[i] = temp[i] ;
}
else
tempresult[i] = temp[i] ;
}
}
/*
* Finalement calcul vrai lissage sur 3 en ayant elimine pics aberrants
*/
for ( i = 0 ; i < ntemp ; i++ ) {
low = i - filtre ;
high = i + filtre ;
if ( low < 0 ) low = 0 ;
if ( high >= ntemp ) high = ntemp ;
n = 0 ;
for ( i1 = low ; i1 < high ; i1++ ) {
vrailiss[i] = vrailiss[i] + tempresult[i1] ;
n = n + 1 ;
}
vrailiss[i] = vrailiss[i] / (double)n ;
}
/* Impression des resultats en table */
fprintf ( stdout, "\n T\tmoyenne\tecart\tmoyenne ecart\tT resultat\tLissee\n");
for ( i = 0 ; i < ntemp ; i++ )
fprintf ( stdout, "\n %g\t%g\t%g\t%g\t%g\t%g",temp[i],moy[i],ecart[i],moyecart[i],tempresult[i],vrailiss[i]);
/* Impression dans fichier */
writefich ( temp, 20, "Input.dat" );
writefich ( tempresult, 20, "Output.dat" );
writefich ( vrailiss, 20, "Lissage.dat" );
return 1 ;
} |
Partager