Scribd
Upload
68 views6 pages

MATH 20D MATLAB HW4

The document provides a series of mathematical computations involving matrices, eigenvalues, and eigenvectors. It includes inputs and outputs for various matrices, detailing their eigenvalues and eigenvectors, as well as stability analysis based on the eigenvalues. The document also discusses the behavior of systems based on the nature of the eigenvalues, including stability and growth patterns.

Uploaded by

bjonathan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
68 views6 pages

MATH 20D MATLAB HW4

The document provides a series of mathematical computations involving matrices, eigenvalues, and eigenvectors. It includes inputs and outputs for various matrices, detailing their eigenvalues and eigenvectors, as well as stability analysis based on the eigenvalues. The document also discusses the behavior of systems based on the nature of the eigenvalues, including stability and growth patterns.

Uploaded by

bjonathan
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 6

4.

1
a) Input:
B = [1.2 2.5;4 0.7]

Output:
B=
1.2000 2.5000
4.0000 0.7000

b) Input
[P,D] = eig(B)

Output → P is eigenvector, D is eigenvalues


P=
0.6501 -0.5899
0.7599 0.8075
D=
4.1221 0
0 -2.2221

4.2
a) Input
A = [2 4;-1 -3]

Output
A=
2 4
-1 -3

b) Input
[P,D] = eig(A)

Output → P is eigenvector, D is eigenvalues


P=

0.9701 -0.7071
-0.2425 0.7071
D=
1 0
0 -2

c) v(t) = c1et [0.9701;-0.2425] + c2e-2t [-0.7071; 0.7071]


As t gets large, As x approaches positive infinity in, the y reaches negative infinity if c1 is
positive and vice versa if c1 is negative
d)
It is similar to my answer because there is a downward almost linear slope as x
approaches positive infinity.

4.3
a) Input
C = [4.1 -3.7;2.7 1]
[Pc,Dc] = eig(C)

Output
C=
4.1000 -3.7000
2.7000 1.0000
Pc =
0.7603 + 0.0000i 0.7603 + 0.0000i
0.3185 - 0.5661i 0.3185 + 0.5661i
Dc =
2.5500 + 2.7545i 0.0000 + 0.0000i
0.0000 + 0.0000i 2.5500 - 2.7545i

b) c1e(2.5500 + 2.7545i)t[0.7603; 0.3185 - 0.5661i] + c2e(2.5500 - 2.7545i)t[0.7603; 0.3185 + 0.5661i]


c)
If the eigenvalues are complex, it exhibits a spiral pattern
If the eigenvalues are real and positive, it grows and reaches infinity
If the eigenvalues are real and negative, it will tend to 0

4.4
a) Input
A = [1.25, -0.97, 4.6; -2.6, -5.2, -0.31; 1.18, 10.3, 1.12]
[Pa,Da] = eig(A)

Output
A=
1.2500 -0.9700 4.6000
-2.6000 -5.2000 -0.3100
1.1800 10.3000 1.1200
Pa =
-0.4591 + 0.0000i 0.7862 + 0.0000i 0.7862 + 0.0000i
-0.5199 + 0.0000i -0.2481 + 0.0826i -0.2481 - 0.0826i
0.7204 + 0.0000i 0.0961 + 0.5516i 0.0961 - 0.5516i
Da =
-7.0663 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 2.1181 + 3.1253i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 2.1181 - 3.1253i

b) It is an unstable system because at least one of the real parts of the eigenvalues is
positive.
4.5
a) Input
A = [-0.0558 -0.9968 0.0802 0.0415;
0.598 -0.115 -0.0318 0;
-3.05 0.388 -0.465 0;
0 0.0805 1 0]
B = [0.01; -0.175; 0.153; 0]
[Pa,Da] = eig(A)
Output
B=

0.0100
-0.1750
0.1530
0
Pa =

0.1994 - 0.1063i 0.1994 + 0.1063i -0.0172 + 0.0000i 0.0067 + 0.0000i


-0.0780 - 0.1333i -0.0780 + 0.1333i -0.0118 + 0.0000i 0.0404 + 0.0000i
-0.0165 + 0.6668i -0.0165 - 0.6668i -0.4895 + 0.0000i -0.0105 + 0.0000i
0.6930 + 0.0000i 0.6930 + 0.0000i 0.8717 + 0.0000i 0.9991 + 0.0000i
Da =
-0.0329 + 0.9467i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -0.0329 - 0.9467i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.5627 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0073 + 0.0000i

b) It is asymptotically stable
c) The 4th component is the largest, the type of rotation this eigenvector is most closely associated
to pitch

4.6
a) Input
B = [0.01;-0.175;0.153;0]
F = [0 7 0 -1]
BF = B*F
Output
BF =
0 0.0700 0 -0.0100
0 -1.2250 0 0.1750
0 1.0710 0 -0.1530
0 0 0 0
b) Input
B = [0.01;-0.175;0.153;0]
F = [0 5 0 -0.1]
BF = B*F
Output
BF =

0 0.0500 0 -0.0010
0 -0.8750 0 0.0175
0 0.7650 0 -0.0153
0 0 0 0

c) Input
A = [-0.0558 -0.9968 0.0802 0.0415;
0.598 -0.115 -0.0318 0;
-3.05 0.388 -0.465 0;
0 0.0805 1 0]
B = [0.01;-0.175;0.153;0]
F = [0 5 0 -0.1]
BF = B*F;
ABF = A + B*F
[Pa,Da] = eig(A)
[Pabf,Dabf] = eig(ABF)

Output
Da =

-0.0329 + 0.9467i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i


0.0000 + 0.0000i -0.0329 - 0.9467i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.5627 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.0073 + 0.0000i

Dabf =
-0.3400 + 0.8104i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -0.3400 - 0.8104i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.0883 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i -0.7425 + 0.0000i

The eigenvalues are not the same


d) F = 1.708
eig(A+ BF)
D=
-0.1569 + 0.9441i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i -0.1569 - 0.9441i 0.0000 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i -0.6300 + 0.0000i 0.0000 + 0.0000i
0.0000 + 0.0000i 0.0000 + 0.0000i 0.0000 + 0.0000i 0.0091 + 0.0000i

All real parts of complex eigenvalues are less than -0.2, there is at least one real eigenvalue that
is within 0.02 or 0.

You might also like