RbControl/mathTools.c at main · ahrendsen/RbControl · GitHub

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
#include <math.h>
float stdDeviation(float* value, int numValues); // Calculate std. deviation of the mean
float stdDeviationMean(float* value, int numValues); // Calculate std. deviation of the mean
float pu2(float value, float uncertainty); // Calculate the percent uncertainty squared of a value.

float stdDeviation(float* value, int numValues){
	float stdDev, sum, avg;

	// First Calculate the Average value
	sum = 0.0;
	int i;
	for(i=0; i < numValues;i++){
		sum += value[i];
	}
	avg = sum / (float) numValues;

	// Then calculate the Standard Deviation of the mean
	sum = 0.0;
	for(i=0; i < numValues;i++){
		sum += pow(avg - value[i],2);
	}
	stdDev = sqrt(sum/(numValues));

	return stdDev;
}

float stdDeviationMean(float* value, int numValues){
	float stdDev, sum, avg;

	// First Calculate the Average value
	sum = 0.0;
	int i;
	for(i=0; i < numValues;i++){
		sum += value[i];
	}
	avg = sum / (float) numValues;

	// Then calculate the Standard Deviation of the mean
	sum = 0.0;
	for(i=0; i < numValues;i++){
		sum += pow(avg - value[i],2);
	}
	stdDev = sqrt(sum/(float)(numValues*(numValues-1)));

	return stdDev;
}

// Calculate percent uncertainty squared
float pu2(float value, float uncertainty){
	return pow(uncertainty/value,2);
}