1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
|
#
# 4003.m
#
# (c) 2020 Prof Dr Andreas Müller, OST Ostschweizer Fachhochschule
#
A = [
-13, 5, -29, 29;
-27, 11, -51, 51;
-3, 1, -2, 5;
-6, 2, -10, 13
];
eig(A)
lambda = 2
B = A - lambda*eye(4)
rref(B)
D = B*B*B*B
lambda = 3
B = A - lambda*eye(4)
rref(B)
D = B*B*B*B
b1 = [0;0;1;1]
b2 = [1;0;0;0]
b3 = [0;1;0;0]
b4 = [0;0;1;2]
T = zeros(4,4);
T(:,1) = b1;
T(:,2) = b2;
T(:,3) = b3;
T(:,4) = b4;
AA = inverse(T)*A*T
A1 = AA(2:4,2:4)
B1 = A1 - 2*eye(3)
B1 * B1
B1 * B1 * B1
c30 = [ 0; 1; 3; 1 ]
c3 = T*c30
lambda=2
B=A-lambda*eye(4)
c2=B*c3
c1=B*c2
T = zeros(4,4);
T(:,1) = [0;0;1;1]
T(:,2) = c1;
T(:,3) = c2;
T(:,4) = c3
det(T)
inverse(T)
det(T)*inverse(T)
inverse(T)*A*T
|
