Advertisement
Bisus

Untitled

Oct 22nd, 2019
121
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
C 1.72 KB | None | 0 0
  1. include <stdio.h>
  2. #include <math.h>
  3. #include <time.h>
  4. int f(const char *name, double *d);
  5.  
  6. /* Ф-ия, получающая имя файла с вещественной последовательностью неизвестной длины.
  7.  * Она записывает по адресу d среднее квадратическое отклонения чисел этого файла от их среднего арифметического.
  8.  * Возвращает:
  9.  * длину последовательности, в случае успешного завершения;
  10.  * -1, если не удалось открыть файл;
  11.  * -2, если не удалось прочитать элемент.
  12.  * */
  13. int f(const char *name, double *d)
  14. {
  15.     FILE *input;
  16.     double x, arith = 0, arith_of_squares = 0;
  17.     int n = 0;
  18.  
  19.     if( !(input = fopen(name, "r")) )
  20.         return -1;
  21.     // Файл открыт
  22.  
  23.     // Подсчет количества элементов последовательности
  24.     while( fscanf(input, "%lf", &x)==1 ) { n++; }
  25.     if( !feof(input) ) {
  26.         fclose(input);
  27.         return -2;
  28.     }
  29.     // Все элементы можно прочитать
  30.  
  31.     rewind(input);
  32.  
  33.     while( fscanf(input, "%lf", &x)==1 )
  34.     {
  35.         arith += x/n;
  36.         arith_of_squares += x*x/n;
  37.     }
  38.     *d = sqrt(arith_of_squares - arith*arith);
  39.     return n;
  40. }
  41.  
  42. int main(int argc, const char *argv[])
  43. {
  44.     double d;
  45.     clock_t begin_time, end_time;
  46.  
  47.     if(argc!=2)
  48.     {
  49.         fprintf(stderr, "Usage: %s name_of_file\n", argv[0]);
  50.         return -1;
  51.     }
  52.  
  53.     begin_time = clock();
  54.     printf("f(%s, &d) returned %d\n", argv[1], f(argv[1], &d));
  55.     end_time = clock();
  56.     printf("d==%.16lf\ntime==%lf\n", d, (double)(end_time - begin_time)/CLOCKS_PER_SEC);
  57.  
  58.     return 0;
  59. }
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement