summaryrefslogtreecommitdiffstats
path: root/fiz/naloga/podatki-obdelani/obdelaj.c
diff options
context:
space:
mode:
Diffstat (limited to 'fiz/naloga/podatki-obdelani/obdelaj.c')
-rw-r--r--fiz/naloga/podatki-obdelani/obdelaj.c14
1 files changed, 10 insertions, 4 deletions
diff --git a/fiz/naloga/podatki-obdelani/obdelaj.c b/fiz/naloga/podatki-obdelani/obdelaj.c
index 05241d7..05736d0 100644
--- a/fiz/naloga/podatki-obdelani/obdelaj.c
+++ b/fiz/naloga/podatki-obdelani/obdelaj.c
@@ -1,12 +1,13 @@
#define _XOPEN_SOURCE
#define PRAG 0.6
-#define KOLNAZ 10
+#define KOLNAZ 30
#include <stdio.h>
#include <dirent.h>
#include <time.h>
#include <unistd.h>
#include <string.h>
#include <stdlib.h>
+#include <math.h>
struct meritev {
int t_alpha; /* milisekunde */
long double F_m;
@@ -14,6 +15,7 @@ struct meritev {
long double F_alpha;
long double V_alpha; /* litri */
long double F_g;
+ long double F_teor;
};
struct poskus {
int id;
@@ -79,6 +81,7 @@ int main (int argc, char ** argv) {
ldbuf = strtold(cp, &cp);
nmer.F_m = ldbuf*9.81;
nmer.F_alpha = -1;
+ nmer.F_teor = -1;
fgets(buf, 256, vagaf);
poskus[i].meritevn++;
}
@@ -96,6 +99,8 @@ int main (int argc, char ** argv) {
omer.V_alpha = h2v(omer.h_alpha);
omer.F_g = 0.0005*9.81 /* teža plastenke */ + omer.V_alpha*9.81 /* teža vode */;
omer.F_alpha = omer.F_m > PRAG ? omer.F_m + omer.F_g : 0;
+ omer.F_teor = poskus[i].P_0*20000*pow((poskus[i].P_0*20000*((poskus[i].V_0+(poskus[i].V_0+0.005)*1)/poskus[i].V_0)), -1.4);
+ omer.F_teor = 1*3.141592*0.000016*(omer.F_teor-10000);
fgets(buf, 256, videof);
}
int intbuf = -1;
@@ -133,17 +138,18 @@ int main (int argc, char ** argv) {
if (poskus[i].meritev[j].F_m >= PRAG || (!odmislimo && (t0 != -1 && poskus[i].meritev[j].V_alpha > 0))) {
if (t0 == -1)
t0 = poskus[i].meritev[j].t_alpha;
- fprintf(outf, "%d,%Lf,%Lf,%Lf,%Lf,%Lf\n",
+ fprintf(outf, "%d,%Lf,%Lf,%Lf,%Lf,%Lf,%Lf\n",
poskus[i].meritev[j].t_alpha - t0,
poskus[i].meritev[j].F_m,
poskus[i].meritev[j].h_alpha,
poskus[i].meritev[j].F_alpha,
poskus[i].meritev[j].V_alpha,
- poskus[i].meritev[j].F_g
+ poskus[i].meritev[j].F_g,
+ poskus[i].meritev[j].F_teor
);
/*
else
- fprintf(outf, "%d,,,,,\n", poskus[i].meritev[j].t_alpha);
+ fprintf(outf, "%d,,,,,,\n", poskus[i].meritev[j].t_alpha);
*/
}
}