diff options
author | Nao Pross <np@0hm.ch> | 2021-04-13 19:48:07 +0200 |
---|---|---|
committer | Nao Pross <np@0hm.ch> | 2021-04-13 19:48:07 +0200 |
commit | d1b602b59a428bea7a59655cd5af34a919e7acf5 (patch) | |
tree | c9ad2469eb5c287d60179e4b57f78373e977a4dc /vorlesungen/slides/7/images/commutator.m | |
parent | Add outline (diff) | |
parent | typos (diff) | |
download | SeminarMatrizen-d1b602b59a428bea7a59655cd5af34a919e7acf5.tar.gz SeminarMatrizen-d1b602b59a428bea7a59655cd5af34a919e7acf5.zip |
Merge branch 'master' of https://github.com/AndreasFMueller/SeminarMatrizen
Diffstat (limited to '')
-rw-r--r-- | vorlesungen/slides/7/images/commutator.m | 111 |
1 files changed, 111 insertions, 0 deletions
diff --git a/vorlesungen/slides/7/images/commutator.m b/vorlesungen/slides/7/images/commutator.m new file mode 100644 index 0000000..5a448db --- /dev/null +++ b/vorlesungen/slides/7/images/commutator.m @@ -0,0 +1,111 @@ +# +# commutator.m +# +# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# + +X = [ + 0, 0, 0; + 0, 0, -1; + 0, 1, 0 +]; + +Y = [ + 0, 0, 1; + 0, 0, 0; + -1, 0, 0 +]; + +Z = [ + 0, -1, 0; + 1, 0, 0; + 0, 0, 0 +]; + +function retval = Dx(alpha) + retval = [ + 1, 0, 0 ; + 0, cos(alpha), -sin(alpha); + 0, sin(alpha), cos(alpha) + ]; +end + +function retval = Dy(beta) + retval = [ + cos(beta), 0, sin(beta); + 0, 1, 0 ; + -sin(beta), 0, cos(beta) + ]; +end + +t = 0.9; +P = Dx(t) * Dy(t) +Q = Dy(t) * Dx(t) +P - Q +(P - Q) * [0;0;1] + +function retval = kurven(filename, t) + retval = -1; + N = 20; + fn = fopen(filename, "w"); + fprintf(fn, "//\n"); + fprintf(fn, "// %s\n", filename); + fprintf(fn, "//\n"); + fprintf(fn, "#macro XYkurve()\n"); + for i = (0:N-1) + v1 = Dx(t * i / N) * [0;0;1]; + v2 = Dx(t * (i+1) / N) * [0;0;1]; + fprintf(fn, "sphere { <%.4f,%.4f,%.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1)); + fprintf(fn, "cylinder { <%.4f,%.4f,%.4f>, <%.4f, %.4f, %.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1), v2(1,1), v2(3,1), v2(2,1)); + end + for i = (0:N-1) + v1 = Dx(t) * Dy(t * i / N) * [0;0;1]; + v2 = Dx(t) * Dy(t * (i+1) / N) * [0;0;1]; + fprintf(fn, "sphere { <%.4f,%.4f,%.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1)); + fprintf(fn, "cylinder { <%.4f,%.4f,%.4f>, <%.4f, %.4f, %.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1), v2(1,1), v2(3,1), v2(2,1)); + end + fprintf(fn, "sphere { <%.4f,%.4f,%.4f>, at }\n", + v2(1,1), v2(3,1), v2(2,1)); + fprintf(fn, "#end\n"); + fprintf(fn, "#declare finalXY = <%.4f, %.4f, %.4f>;\n", + v2(1,1), v2(3,1), v2(2,1)); + fprintf(fn, "#macro YXkurve()\n"); + for i = (0:N-1) + v1 = Dy(t * i / N) * [0;0;1]; + v2 = Dy(t * (i+1) / N) * [0;0;1]; + fprintf(fn, "sphere { <%.4f,%.4f,%.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1)); + fprintf(fn, "cylinder { <%.4f,%.4f,%.4f>, <%.4f, %.4f, %.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1), v2(1,1), v2(3,1), v2(2,1)); + end + for i = (0:N-1) + v1 = Dy(t) * Dx(t * i / N) * [0;0;1]; + v2 = Dy(t) * Dx(t * (i+1) / N) * [0;0;1]; + fprintf(fn, "sphere { <%.4f,%.4f,%.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1)); + fprintf(fn, "cylinder { <%.4f,%.4f,%.4f>, <%.4f, %.4f, %.4f>, at }\n", + v1(1,1), v1(3,1), v1(2,1), v2(1,1), v2(3,1), v2(2,1)); + end + fprintf(fn, "sphere { <%.4f,%.4f,%.4f>, at }\n", + v2(1,1), v2(3,1), v2(2,1)); + fprintf(fn, "#end\n"); + fprintf(fn, "#declare finalYX = <%.4f, %.4f, %.4f>;\n", + v2(1,1), v2(3,1), v2(2,1)); + + fclose(fn); + retval = 0; +end + +function retval = kurve(i) + n = pi / 180; + filename = sprintf("f/%04d.inc", i); + kurven(filename, n * i); +end + +for i = (1:60) + kurve(i); +end |