From 654d9ee4aad5389798a37f521a7a4238be3fdc34 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sun, 21 Mar 2021 17:21:45 +0100 Subject: JK bilder --- vorlesungen/slides/5/beispiele/kernbild.m | 78 +++++++++++++++++++++++++++++++ 1 file changed, 78 insertions(+) create mode 100644 vorlesungen/slides/5/beispiele/kernbild.m (limited to 'vorlesungen/slides/5/beispiele/kernbild.m') diff --git a/vorlesungen/slides/5/beispiele/kernbild.m b/vorlesungen/slides/5/beispiele/kernbild.m new file mode 100644 index 0000000..ca6bd6a --- /dev/null +++ b/vorlesungen/slides/5/beispiele/kernbild.m @@ -0,0 +1,78 @@ +# +# kernbild.m +# +# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# + +rand("seed", 1291) +rand("seed", 4711) + +lambda1 = 1; +lambda2 = 1.8; + +A = [ + lambda1, 0, 0; + 0, lambda2, 1; + 0, 0, lambda2 +]; + +B = eye(3) + rand(3,3); +det(B) + + +C = B*A*inverse(B) +rank(C) + +# Eigenwert lambda1 +E2 = C - lambda1 * eye(3) +rref(E2) + +# Eigenwert lambda2, k = 1 +E1 = C - lambda2 * eye(3) +D = rref(E1); +K1 = [ + -D(1,3); + -D(2,3); + 1 +]; +K1(:,1) = K1(:,1) / norm(K1(:,1)); +K1 + +f = fopen("JK.inc", "w"); +fprintf(f, "//\n// JK.inc\n//\n// (c) 2021 Prof Dr Andreas Müller\n//\n\n"); +fprintf(f, "// Kern und Bild von C - %.3f I\n", lambda2); +fprintf(f, "#declare k11 = < %.5f, %.5f, %.5f>;\n", K1(1,1), K1(2,1), K1(3,1)); +fprintf(f, "#declare j11 = < %.5f, %.5f, %.5f>;\n", E1(1,1), E1(2,1), E1(3,1)); +fprintf(f, "#declare j12 = < %.5f, %.5f, %.5f>;\n", E1(1,2), E1(2,2), E1(3,2)); +fprintf(f, "\n"); + +# k = 2 +D = rref(E1 * E1); +K2 = [ + -D(1,2), -D(1,3); + 1, 0; + 0, 1 +] +K2(:,1) = K2(:,1) / norm(K2(:,1)); +K2(:,2) = K2(:,2) / norm(K2(:,2)); +K2 + +fprintf(f, "// Kern und Bild von (C - %.3f I)^2\n", lambda2); +fprintf(f, "#declare k21 = < %.5f, %.5f, %.5f>;\n", K2(1,1), K2(2,1), K2(3,1)); +fprintf(f, "#declare k22 = < %.5f, %.5f, %.5f>;\n", K2(1,2), K2(2,2), K2(3,2)); +fprintf(f, "#declare j21 = < %.5f, %.5f, %.5f>;\n", E2(1,1), E2(2,1), E2(3,1)); +fprintf(f, "\n"); + +fclose(f); + +# Verifikation +x = K2 \ K1 +K2 * x + +eig(C) + +[U, S, V] = svd(C) + + +s = rand("seed") + -- cgit v1.2.1 From e2737868b67b183152c4436ececf12e68b5a145b Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Andreas=20M=C3=BCller?= Date: Sun, 21 Mar 2021 20:37:06 +0100 Subject: new slides --- vorlesungen/slides/5/beispiele/kernbild.m | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) (limited to 'vorlesungen/slides/5/beispiele/kernbild.m') diff --git a/vorlesungen/slides/5/beispiele/kernbild.m b/vorlesungen/slides/5/beispiele/kernbild.m index ca6bd6a..28cd552 100644 --- a/vorlesungen/slides/5/beispiele/kernbild.m +++ b/vorlesungen/slides/5/beispiele/kernbild.m @@ -47,7 +47,8 @@ fprintf(f, "#declare j12 = < %.5f, %.5f, %.5f>;\n", E1(1,2), E1(2,2), E1(3,2)); fprintf(f, "\n"); # k = 2 -D = rref(E1 * E1); +E12 = E1 * E1 +D = rref(E12); K2 = [ -D(1,2), -D(1,3); 1, 0; @@ -60,7 +61,7 @@ K2 fprintf(f, "// Kern und Bild von (C - %.3f I)^2\n", lambda2); fprintf(f, "#declare k21 = < %.5f, %.5f, %.5f>;\n", K2(1,1), K2(2,1), K2(3,1)); fprintf(f, "#declare k22 = < %.5f, %.5f, %.5f>;\n", K2(1,2), K2(2,2), K2(3,2)); -fprintf(f, "#declare j21 = < %.5f, %.5f, %.5f>;\n", E2(1,1), E2(2,1), E2(3,1)); +fprintf(f, "#declare j21 = < %.5f, %.5f, %.5f>;\n", E12(1,1), E12(2,1), E12(3,1)); fprintf(f, "\n"); fclose(f); -- cgit v1.2.1