/* [wxMaxima batch file version 1] [ DO NOT EDIT BY HAND! ]*/ /* [ Created with wxMaxima version 19.05.7 ] */ /* [wxMaxima: title start ] Aufgabe gps_17: Skyplot [wxMaxima: title end ] */ /* [wxMaxima: input start ] */ grad(bogen):=bogen/%pi*180; /* [wxMaxima: input end ] */ /* [wxMaxima: section start ] Aufgabe 1: Ausgangswerte [wxMaxima: section end ] */ /* [wxMaxima: subsect start ] Daten [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ r:matrix( [4151529], [744103], [4769225] ); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ s:matrix( [-7209930], [14399010], [21785556] ); /* [wxMaxima: input end ] */ /* [wxMaxima: subsect start ] Vektor vom Empfänger zum Satelliten [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ v:s-r; /* [wxMaxima: input end ] */ /* [wxMaxima: section start ] Aufgabe 2: Elevation [wxMaxima: section end ] */ /* [wxMaxima: input start ] */ elevation(n,u):=block( [elev], elev:asin(abs(u.n)/(mat_norm(u,frobenius)*mat_norm(n,frobenius))), if sign(n.u)=neg then -elev else elev ); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ abs(r); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ r.v; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ mat_norm(r,frobenius); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ elev:elevation(r,v),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(elev),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: section start ] Aufgabe 3: Richtung zum Satelliten in der Ebene [wxMaxima: section end ] */ /* [wxMaxima: input start ] */ projektion(x,n,r0):=x-(((x-r0).n)/(n.n)).n; /* [wxMaxima: input end ] */ /* [wxMaxima: subsect start ] Projektion des Satelliten [wxMaxima: subsect end ] */ /* [wxMaxima: subsubsect start ] Projektion [wxMaxima: subsubsect end ] */ /* [wxMaxima: input start ] */ projektion_satvect:projektion(s,r,r),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: subsubsect start ] Richtung vom Empfänger zum projizierten Satelliten [wxMaxima: subsubsect end ] */ /* [wxMaxima: input start ] */ satrichtung:projektion_satvect-r; /* [wxMaxima: input end ] */ /* [wxMaxima: section start ] Aufgabe 4: Nordrichtung [wxMaxima: section end ] */ /* [wxMaxima: subsect start ] Projektion der z-Achse [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ nord_v:matrix([0],[0],[1]); /* [wxMaxima: input end ] */ /* [wxMaxima: subsubsect start ] Projektion [wxMaxima: subsubsect end ] */ /* [wxMaxima: input start ] */ projektion_nordvect:projektion(nord_v,r,r),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: subsubsect start ] Nordrichtung [wxMaxima: subsubsect end ] */ /* [wxMaxima: input start ] */ nordrichtung:projektion_nordvect-r; /* [wxMaxima: input end ] */ /* [wxMaxima: section start ] Aufgabe 5: Azimut [wxMaxima: section end ] */ /* [wxMaxima: subsect start ] Bestimmung des Azimuts [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ azimut(u,v):=block( acos((u.v)/(mat_norm(u,frobenius)*mat_norm(v,frobenius))) ); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ az:azimut(nordrichtung,satrichtung); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(az),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: section start ] Aufgabe 6: Weitere Satelliten berechnen [wxMaxima: section end ] */ /* [wxMaxima: subsect start ] Satellit 1 [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ s1:matrix( [16146937], [14754344], [14992035] ); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ v1:s1-r; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ elev1:elevation(r,v1),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(elev1),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ projektion_satvect1:projektion(s1,r,r),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ satrichtung1:projektion_satvect1-r; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ az1:azimut(nordrichtung,satrichtung1); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(az1),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: subsect start ] Satellit 2 [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ s2:matrix( [26124344], [4803272], [-1301059] ); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ v2:s2-r; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ elev2:elevation(r,v2),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(elev2),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ projektion_satvect2:projektion(s2,r,r),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ satrichtung2:projektion_satvect2-r; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ az2:azimut(nordrichtung,satrichtung2); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(az2),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: subsect start ] Satellit 3 (Azimut bei 270°!) [wxMaxima: subsect end ] */ /* [wxMaxima: input start ] */ s3:matrix( [17719707], [-8916564], [17798680] ); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ v3:s3-r; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ elev3:elevation(r,v3),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(elev3),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ projektion_satvect3:projektion(s3,r,r),numer; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ satrichtung3:projektion_satvect3-r; /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ az3:azimut(nordrichtung,satrichtung3); /* [wxMaxima: input end ] */ /* [wxMaxima: input start ] */ grad(az3),numer; /* [wxMaxima: input end ] */ /* Old versions of Maxima abort on loading files that end in a comment. */ "Created with wxMaxima 19.05.7"$