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
| #include <stdio.h>
#include <stdlib.h>
#include <math.h>
#define ARRSIZE 10
float avgsum(); // среднее арифметическое
float variance(); // дисперсия
float stddev(); // стандартное среднеквадратичное отклонение
float cov(); // ковариация
float pirson(); // коэффициент корреляции Пирсона
int main()
{
float arr_X[ARRSIZE]={2,9,13,12,4,6,14,10,6,7};
float arr_Y[ARRSIZE]={8,1,10,5,6,34,23,64,2,4};
float avg_sum,var,std_dev;
avg_sum=avgsum(arr_X,ARRSIZE);
printf("avgsum X = %f\n",avg_sum);
printf("avgsum Y = %f\n",avgsum(arr_Y,ARRSIZE));
var=variance(arr_X,ARRSIZE);
printf("variance X = %f\n",var);
printf("variance Y = %f\n",variance(arr_Y,ARRSIZE));
std_dev=stddev(arr_X,ARRSIZE);
printf("std_dev X = %f\n",std_dev);
printf("std_dev Y = %f\n",stddev(arr_Y,ARRSIZE));
printf("covariation (X,Y) = %f\n",cov(arr_X,arr_Y,ARRSIZE));
printf("correlation (X,Y) = %f\n",pirson(arr_X,arr_Y,ARRSIZE));
return 0;
}
float avgsum(float *array_adr, int array_size)
{
int i;
float sum=0;
for (i=0;i<array_size;i++)
{
sum+=array_adr[i];
}
return(sum / array_size);
}
float variance(float *array_adr, int array_size)
{
int i;
float sumsqare=0;
float mean=avgsum(array_adr,ARRSIZE);
for (i=0;i<array_size;i++)
{
sumsqare+=(array_adr[i]-mean)*(array_adr[i]-mean);
}
return(sumsqare / array_size);
}
float stddev(float *array_adr, int array_size)
{
return ( sqrt(variance(array_adr,ARRSIZE)) );
}
float cov(float *array_adr1, float *array_adr2, int array_size)
{
int i=0;
float mean1=avgsum(array_adr1,array_size);
float mean2=avgsum(array_adr2,array_size);
float producedelta=0;
for (i=0;i<array_size;i++)
{
producedelta+=(array_adr1[i]-mean1)*(array_adr2[i]-mean2);
}
return (producedelta/array_size);
}
float pirson(float *array_adr1, float *array_adr2, int array_size)
{
int i;
float covariat, deviation1, deviation2;
covariat = cov(array_adr1,array_adr2,array_size);
deviation1 = stddev(array_adr1,array_size);
deviation2 = stddev(array_adr2,array_size);
return (covariat/(deviation1*deviation2));
} |
(
