IQA/funcoes.c at master · RodrigoAkiraMitani/IQA · 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
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
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
#include <stdio.h>
#include <math.h>

double qCF (double CF){

	double q1, A = 98.03, B = -36.45, C = 3.138, D = 0.06776;

	if (CF > 100000){

		return 3;
	}
	else{

		q1 = A + (B * log10(CF)) + C * (pow(log10(CF), 2)) + (D * pow(log10(CF), 3));

		return q1;

	}

}

double qpH (double pH){

	double q2, A = 0.05421, B = 1.23, C = -0.09873;

	if (pH >= 12){

		return 3;

	}else if (pH < 2){

		return 2;

	}else {

		q2 = A * pow(pH, (B * pH) + (C * pow(pH, 2))) + 5.213;

		return q2;

	}

}

double qDBO (double DBO){

	double q3, A = 102.6, B = -0.1101;

	if (DBO > 30){
		return 2;
	}else{

		q3 = A * exp(B * DBO);

		return q3;
	}
}

double qNT (double NT){

	double q4, A = 98.96, B = -0.2232, C = -0.006457;

	if (NT > 100){
		return 1;
	}else{

		q4 = A * pow(NT, (B + (C * NT)));

		return q4;

	}
}

double qFT (double FT){

	double q5, A= 213.7, B = -1.680, C = 0.3325;

	if(FT > 10){
		return 1;
	}else{
		q5 = A * exp(B * pow(FT, C));
		return q5;
	}
}

double qDT (double DT){

	double q6, A = 0.0003869, B= 0.1815, C = 0.01081;

	if(DT <= 5){
		return 1;
	}else if (DT > 15)
	{
		return 9;
	}else{

		q6 = 1 / (A * pow(DT + B, 2) + C);

		return q6;
	}
}

double qTU (double TU){

	double q7, A = 97.34, B= -0.01139, C = -0.04917;

	if(TU > 100){
		return 5;
	}else{

		q7 = A * exp((B * TU) + (C * sqrt(TU)));

		return q7;
	}
}

double qST (double ST){

	double q8, A = 80.26, B = -0.00107, C = 0.03009, D = -0.1185;

	if(ST > 500){
		return 32;
	}else{

		q8 = (A * exp((B * ST) + (C * sqrt(ST)))) + D * ST;

		return q8;
	}
}

double qOxi (double Oxi){

	double q9, A = 100.8, B = -106, C = -3745;

	if (Oxi == 140){
		return 47;
	}else{

		double q9 = A * exp(pow(Oxi + B, 2) / C);

		return q9;
	}
}