From 155ff4351fb8f64f5efc7f822f1babedac94062a Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Anton=20Luka=20=C5=A0ijanec?= Date: Mon, 9 May 2022 09:18:12 +0200 Subject: nemscina --- "fiz/naloga/numeri\304\215no.c" | 192 +++++++++++++++++++++++++++++++++------- fiz/naloga/prodani/filelist.txt | 2 + 2 files changed, 162 insertions(+), 32 deletions(-) create mode 100644 fiz/naloga/prodani/filelist.txt (limited to 'fiz') diff --git "a/fiz/naloga/numeri\304\215no.c" "b/fiz/naloga/numeri\304\215no.c" index e102899..c1ca08c 100644 --- "a/fiz/naloga/numeri\304\215no.c" +++ "b/fiz/naloga/numeri\304\215no.c" @@ -2,28 +2,45 @@ #include #include #include +#include +#include +#include #define UVOD "program za numerični izračun jakosti magnetnega polja okoli helmholtzove tuljave\n" \ "sem spisal anton luka šijanec za projektno nalogo pri fiziki v tretjem letniku gimb.\n" \ - "uporaba: %s in argumenti po vrsti:\n" \ + "uporaba: %s in nujni argumenti po vrsti:\n" \ " 1. radij enega navitja v metrih\n" \ " 2. tok, ki teče po vodniku v amperih\n" \ " 3. število navojev na enem navitju\n" \ " 4. razmak med merilnimi točkami v metrih\n" \ " 5. koliko meritev od središča v obe dimenziji naj napravimo\n" \ " 6. koliko kotov naj ima navitje - računamo, kot da je mnogokotnik\n" \ - " 7. tip izhodnih podatkov (pgm ali csv)\n" \ + "nenujni argumenti: prazen niz argument nastavi na privzeto vrednost.\n" \ + " 7. tip izhodnih podatkov (pgm, ppm ali tsv). privzeto je to pgm.\n" \ + " 8. razmak med navitjima. privzeto je to r/2.\n" \ + " 9. zamik enega navitja v metrih. privzeto sta osi navitij ista premica.\n" \ + " A. polmer druge tuljave. privzeto je enak polmeru prve tuljave.\n" \ + " B. tok druge tuljave v metrih. privzeto je enak toku prve tuljave.\n" \ "oblika izhodnih podatkov, če je 7. parameter pgm, so pgm slike z vrednostmi 0-255\n" \ " - slika je prerez tuljave. magnetno polje teče vodoravno.\n" \ " - vrednosti direktno korelirajo z izračunano jakostjo v decigaussih: 10e-5 tesla\n" \ " - slika je široka 1+2*koliko in visoka 1+2*koliko (5. argument) slikovnih točk\n" \ - "oblika izhodnih podatkov, če je 7. parameter csv, je csv, z naslednjimi stolpci:\n" \ - " 1. vodoravna komponenta oddaljenosti od središča tuljave v metrih\n" \ - " 2. navpična komponenta oddaljenosti od središča tuljave v metrih\n" \ - " 3. jakost magnetnega polja v teslah - tokrat ni v decigaussih!\n" \ - " 4. smer vektorja magnetnega polja v radianih. 0 radianov je smer v desno" + "oblika izhodnih podatkov, če je 7. parameter ppm, so barvne slike z 8 biti na barvo\n" \ + " - barvi rdeča in zelena predstavljata komponenti vektorjev polja i in j\n" \ + "oblika izhodnih podatkov, če je 7. parameter tsv, je tsv, z naslednjimi stolpci:\n" \ + " 1. komponenta i krajevnega vektorja točke meritve. 0,0 je v središču tuljave.\n" \ + " 2. komponenta j krajevnega vektorja točke meritve. 0,0 je v središču tuljave.\n" \ + " 3. komponenta i vektorja magnetnega polja v točki meritve.\n" \ + " 4. komponenta j vektorja magnetnega polja v točki meritve.\n" \ + " 5. jakost magnetnega polja v teslah - tokrat ni v decigaussih!\n" \ + "način izdelave animacije: podan naj bo samo en argument - animacija - TOLE NE DELA\n" \ + " - delajo se datoteke animacija0000.ppm, animacija0001.ppm, ...\n" \ + " - parametri se sinusoidno spreminjajo po vdelanih konstantah.\n" \ + " - slike lahko recimo s ffmpeg(1) nato pretvorite v videoposnetek." enum oblika { PGM, - CSV + PPM, + TSV, + ANIMACIJA }; struct vektor { long double i; // x - desno na sliki @@ -87,27 +104,29 @@ struct vektor tuljava (long double R, unsigned kotov, struct vektor m /* meritev B = množi(B, MU0/(4*M_PI)); return B; } // ena zanka ob toku 1 A. pomnoži s tokom in številom navitij. 0,0 je v sredini. B teče v desno. +struct vektor zavrti_okoli (struct vektor točka, struct vektor središče, long double kot) { + točka = seštej(točka, množi(središče, -1)); + struct vektor r = { + .i = cosl(kot)*točka.i + sin(kot)*točka.j, + .j = cosl(kot)*točka.j - sin(kot)*točka.i + }; + return seštej(r, središče); +} // TODO: implement void natisni (FILE * f, struct vektor v, const char * i) { fprintf(f, "vektor %s {\n\t.i = %Lf,\n\t.j = %Lf,\n\t.k = %Lf\n}\n", i, v.i, v.j, v.k); } -int main (int argc, char ** argv) { - if (argc != 1 + 7) - error(1, 0, UVOD, argv[0] ? argv[0] : "./numerično"); - long double R = strtold(argv[1], NULL); - long double I = strtold(argv[2], NULL); - unsigned n = strtol(argv[3], NULL, 10); - long double razmak = strtold(argv[4], NULL); - int koliko = strtold(argv[5], NULL); - unsigned kotov = strtol(argv[6], NULL, 10); - enum oblika oblika = argv[7][0] == 'p' || argv[7][0] == 'P' ? PGM : CSV; +int nariši (long double R, long double I, unsigned n, long double razmak, int koliko, + unsigned kotov, enum oblika oblika, long double med_tuljavama, + struct vektor zamik_izven_osi, long double R2, long double I2, FILE * out) { struct vektor merilno_mesto = { // krajevni vektor .k = 0 }; - if (oblika == CSV) - error(2, 0, "CSV oblika še ni implementirana."); - printf("P5 %u %u 255\n", koliko*2+1, koliko*2+1); + if (oblika == PGM) + fprintf(out, "P5 %u %u 255\n", koliko*2+1, koliko*2+1); + if (oblika == PPM) + fprintf(out, "P6 %u %u 255\n", koliko*2+1, koliko*2+1); struct vektor Rpolovic = { - .i = R/2, + .i = med_tuljavama, .j = 0, .k = 0 }; @@ -115,8 +134,12 @@ int main (int argc, char ** argv) { merilno_mesto.i = i*razmak; for (int j = -koliko; j <= koliko; j++) { merilno_mesto.j = j*razmak; - unsigned long long int out = - 1000*absolutno(seštej( + struct vektor merilno_mesto2 = seštej( + merilno_mesto, + zamik_izven_osi + ); + struct vektor B = + seštej( množi( množi( tuljava( @@ -134,10 +157,10 @@ int main (int argc, char ** argv) { množi( množi( tuljava( - R, + R2, kotov, seštej( - merilno_mesto, + merilno_mesto2, množi( Rpolovic, -1 @@ -146,14 +169,119 @@ int main (int argc, char ** argv) { ), n ), - I + I2 ) - )); - if (out > 255) - putchar(255); - else - putchar(out); + ); + switch (oblika) { + case PGM: + if (10000*absolutno(B) > 255) + putc(255, out); + else + putc(10000*absolutno(B), out); + break; + case PPM: +#define NATISNI_KOMPONENTO(x) if (10000*fabsl(B.x) > 255) \ + putc(255, out); \ + else \ + putc(10000*fabsl(B.x), out); + NATISNI_KOMPONENTO(i); + NATISNI_KOMPONENTO(j); + NATISNI_KOMPONENTO(k); + break; + case TSV: + fprintf(out, "%Lf\t%Lf\t%Lf\t%Lf\t%Lf\n", + merilno_mesto.i, merilno_mesto.j, + B.i, B.j, absolutno(B)); + break; + case ANIMACIJA: + abort(); // invalid + break; + } + } + } + return 0; +} +int shouldexit = 0; +void handler (int s __attribute__((unused))) { + shouldexit++; +} +int main (int argc, char ** argv) { + if (argv[1] && !strcmp(argv[1], "animacija")) + goto animacija; + if (argc < 8) + error(1, 0, UVOD, argv[0] ? argv[0] : "./numerično"); + long double R = strtold(argv[1], NULL); + long double I = strtold(argv[2], NULL); + unsigned n = strtol(argv[3], NULL, 10); + long double razmak = strtold(argv[4], NULL); + int koliko = strtold(argv[5], NULL); + unsigned kotov = strtol(argv[6], NULL, 10); + enum oblika oblika = PGM; + long double med_tuljavama = R/2; + struct vektor zamik_izven_osi = {0, 0, 0}; + long double R2 = R; + long double I2 = I; + if (argc > 7 && argv[7][0]) + switch (argv[7][1]) { + case 'G': + case 'g': + oblika = PGM; + break; + case 'P': + case 'p': + oblika = PPM; + break; + case 'S': + case 's': + oblika = TSV; + break; + case 'N': + case 'n': + oblika = ANIMACIJA; + break; } + if (argc > 8 && argv[8][0]) + med_tuljavama = strtold(argv[8], NULL); + if (argc > 9 && argv[9][0]) + zamik_izven_osi.j = strtold(argv[9], NULL); + if (argc > 0xA && argv[0xA][0]) + R2 = strtold(argv[0xA], NULL); + if (argc > 0xB && argv[0xB][0]) + I2 = strtold(argv[0xB], NULL); + if (oblika != ANIMACIJA) + return nariši(R, I, n, razmak, koliko, kotov, oblika, med_tuljavama, + zamik_izven_osi, R2, I2, stdout); +animacija: + signal(SIGTERM, handler); + signal(SIGINT, handler); + unsigned čas = 0; + while (!shouldexit && čas < 21) { + char fn[25] = "animacija"; + sprintf(fn+strlen(fn), "%04d.ppm", čas); + fprintf(stderr, "%s\n", fn); + FILE * f = fopen(fn, "w"); + if (!f) + error(2, errno, "fopen"); + long double R = 0.06; + long double I = 2; + unsigned n = 320; + long double razmak = 0.0007; + unsigned koliko = 350; + unsigned kotov = 30; + enum oblika oblika = PPM; + long double med_tuljavama = (R/2/10)*(čas); + struct vektor zamik = { + .i = 0, + .j = 0, // zamik osi + .k = 0 + }; + long double R2 = R; + long double I2 = I; + if (nariši(R, I, n, razmak, koliko, kotov, oblika, med_tuljavama, zamik, R2, I2, f)) + error(3, errno, "nariši"); + if (fclose(f)) + error(4, errno, "fclose"); + čas++; } return 0; } diff --git a/fiz/naloga/prodani/filelist.txt b/fiz/naloga/prodani/filelist.txt new file mode 100644 index 0000000..52fbd12 --- /dev/null +++ b/fiz/naloga/prodani/filelist.txt @@ -0,0 +1,2 @@ +file helmholtz.mkv +file helmholtz.počasnejši.mkv -- cgit v1.2.3