Tags | |
UUID | b7ab0f8e-227c-11e6-9770-bc764e2038f2 |
The computation of eigenvalues and eigenvectors can serve many purposes; however, when it comes to differential equations eigenvalues and eigenvectors are most often used to find straight-line solutions of linear systems.
To find eigenvalues, we use the formula:
A→v=λ→v
where A=(abdc) and →v=(xy)
(abdc)(xy)=λ(xy), which can be written in components as
ax+by=λx
cx+dy=λy
We want to solve for non-zero solution, such that the system becomes
(a-λ)x + by=0
cx +(d-λ)y=0
We can prove that given a matrix A whose determinant is not equal to zero, the only equilibrium point for the linear system is the origin, meaning that to solve the system above we take the determinant and set it equal to zero.
det(a-λbcd-λ)=0
Every time we compute eigenvalues and eigenvectors we use this format, which can also be written as det(A-λ→I)=0, where I is the Identity matrix →I=(1001). Computation of det(A-λ→I)=0 leads to the Characteristic Polynomial, where the roots of this polynomial are the eigenvalues of the matrix A.
det(A-λ→I)=det(a-λbcd-λ)=(a-λ)(d-λ)-bc=0, which expands to the quadratic polynomial
λ2-(a+d)λ+(ad-bc)=0.
This is referred to as the characteristic polynomial, where the characteristic polynomial always has two roots. These roots can be real or complex, and they do not have to be distinct. If the roots are complex we say that the matrix has complex eigenvalues. Otherwise, we say that the matrix has real eigenvalues.
Here are examples of how to solve for both kinds of eigenvalues:
Let's begin with an example where we compute real eigenvalues:
Suppose we have the matrix:
A=(5432)
det(A-λI)=det(5-λ432-λ)=(5-λ)(2-λ)-4⋅3=0
(5-λ)(2-λ)-12=λ2-7λ+(-2)=0
The roots are:
λ=7±√49-482
λ=4,3
Now we will compute complex eigenvalues:
Before we start we should review what it means to have a complex number.
"Complex numbers are numbers of the form x + iy, where x and y are real numbers and I is the 'imaginary number' √-1 " (Blanchard, Devaney, Hall, 291).
Consider the system where A = (-2-33-2)
det(A-λI)=det(-2-λ-33-2-λ)=(-2-λ)(-2-λ)-(-3⋅3)=λ2+4λ+13=0.
The roots are:
λ=-4±√-362
We see that the √-36 is equal to 6i, such that the eigenvalues become:
λ=-4±6i2=-2±3i
Given a matrix A=(abcd) and we know that λ is an eigenvalue, we use the same equation from above A→v=λ→v to solve for →v of the form →v=(xy). We notice that A→v=λ→v turns into a system of linear equations:
ax+by=λx
cx+dy=λy
Because we have already solved for lambda, "we know that there is at least an entire line of eigenvectors (x, y) that satisfy this system of equations. This infinite number of eigenvectors means that the equations are redundant. That is, either the two equations are equivalent, or one of the equations is always satisfied" (Blanchard, Devaney, Hall, 266).
We will give an example to demonstrate what is meant by the statement above:
Suppose the matrix A = (2213)
det(A-λI)=(2-λ)(3-λ)-(2⋅1)=0
λ2-5λ+4=0
λ=1,4 or λ1=4,λ2=1
Let's use λ2 in the equation:
A(xy)=(2213)(xy)=1(xy)
Rewritten in terms of components, the equation becomes
2x+2y=x
1x+3y=y
or
x+2y=0
x+2y=0
It is obvious that -12x=y satisfies both equations, such that the eigenvector for λ2=(1-12)
Now let's view an example where there are complex eigenvalues and a complex eigenvector:
Let's begin where we left off in the example from before where A = (-2-33-2)
We found that eigenvalues were λ1=-2+3i,λ2=-2-3i
Let's take λ1 and plug it into the equation,
A(xy)=(2213)(xy)=(-2+3i)(xy)
As a system of equations we have
-2x-3y =(-2+3i)x
3x-2y=(-2+3i)y
which can be rewritten as
(-3i)x+3y=0
3x+(-3i)y=0.
Just as in the example above, the equations are redundant. We see that (i)x=y and →v=(1i)
https://youtu.be/bOreOaAjDno
http://tutorial.math.lamar.edu/Classes/DE/LA_Eigen.aspx
https://www.khanacademy.org/math/linear-algebra/alternate-bases/eigen-everything/v/linear-algebra-introduction-to-eigenvalues-and-eigenvectors
No comments |