diff options
| author | LordMcFungus <mceagle117@gmail.com> | 2021-03-22 18:05:11 +0100 |
|---|---|---|
| committer | GitHub <noreply@github.com> | 2021-03-22 18:05:11 +0100 |
| commit | 76d2d77ddb2bed6b7c6b8ec56648d85da4103ab7 (patch) | |
| tree | 11b2d41955ee4bfa0ae5873307c143f6b4d55d26 /buch/chapters/40-eigenwerte/beispiele/n.m | |
| parent | more chapter structure (diff) | |
| parent | add title image (diff) | |
| download | SeminarMatrizen-76d2d77ddb2bed6b7c6b8ec56648d85da4103ab7.tar.gz SeminarMatrizen-76d2d77ddb2bed6b7c6b8ec56648d85da4103ab7.zip | |
Merge pull request #1 from AndreasFMueller/master
update
Diffstat (limited to 'buch/chapters/40-eigenwerte/beispiele/n.m')
| -rw-r--r-- | buch/chapters/40-eigenwerte/beispiele/n.m | 55 |
1 files changed, 55 insertions, 0 deletions
diff --git a/buch/chapters/40-eigenwerte/beispiele/n.m b/buch/chapters/40-eigenwerte/beispiele/n.m new file mode 100644 index 0000000..af0219b --- /dev/null +++ b/buch/chapters/40-eigenwerte/beispiele/n.m @@ -0,0 +1,55 @@ +# +# n.m -- Polynome mit dem gleichen Wert von p(A) +# +# (c) 2021 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule +# + +A0 = [ + 2, 1, 0; + 0, 2, 0; + 0, 0, 3 +]; + +# find a 3x3 matrix in SL(3,Z) + +function retval = zufallswert() + x = round(rand() * 10) - 2; + if (x >= 0) + x = x + 1; + endif + retval = x; +end + +function retval = zufallsmatrix(n) + retval = zeros(n, n); + for i = (1:n) + for j = (1:n) + retval(i,j) = zufallswert(); + end + end +end + +function retval = regulaer(n) + d = 0; + do + retval = zufallsmatrix(2); + d = det(retval); + until (d == 1); +end + +function retval = eingebettet(n,k) + retval = eye(n); + retval(k:k+1,k:k+1) = regulaer(2); +end + +format long + +B = eye(3); +B = B * eingebettet(3,2) +B = B * eingebettet(3,1) + +B +inverse(B) + +A = round(B * A0 * inverse(B)) + |