http://setosa.io/ev/eigenvectors-and-eigenvalues/

Explained Visually

By Victor Powell and Lewis Lehe

Eigenvalues/vectors are instrumental to understanding electrical circuits, mechanical systems, ecology and even Google's PageRank algorithm.
Let's see if visualization can make these ideas more intuitive.

To begin, let v be
a 2-dimensional vector (shown as a point) and A be
a matrix with columns a1 and a2 (shown
as arrows). If we multiply v by A,
then A sends v to
a new vector Av.

012345012345xyva₁a₂Αv

Α
=
a₁,x
1.00
a₂,x
0.50
a₁,y
0.50
a₂,y
1.00
=
1.00
0.50
0.50
1.00
v
=
2.00
v, x
3.00
v, y
Αv
=
3.50
v, x
4.00
v, y

If you can draw a line through (0,0), v and Av,
then Av is
just v multiplied
by a number λ;
that is, Av=λv.
In this case, we call λ an eigenvalue and v an eigenvector.
For example, here (1,2) is
an eigvector and 5 an
eigenvalue.

Av=(1821)⋅(12)=5(12)=λv.

Below, change the bases of A and
drag v to
be its eigenvector. Note two facts: First, every point on the same line as an eigenvector is another eigenvector. That line is an eigenspace. Second, when λ<1, Av is
closer to (0,0) than v;
and when λ>1,
it's farther away.

012345012345xyva₁a₂Αvλ₁ = 1.5λ₂ = 0.5s₁s₂

What are eigenvalues/vectors good for?

Eigenvalues/vectors explain the behavior of systems that evolve step-by-step, where each step occurs as multiplication by a matrix A.
If you keep multiplying v by A,
you get a sequence v,Av,A2v, etc.
As you can see below, eigenspaces attract this sequence and draw it toward (0,0) or
farther away, depending on their eigenvalues.

02004006008001,00002004006008001,000xyvΑvΑ²va₁a₂λ₁ = 1.1λ₂ = 0.5s₁s₂

Let's explore some applications and properties of these sequences.

Fibonacci Sequence

Suppose you have some amoebas in a petri dish. Every minute, all adult amoebas produce one child amoeba, and all child amoebas grow into adults (Note: this is not really how amoebas reproduce.). So if t is
a minute, the equation of this system is

adultst+1childrent+1==adultst+childrentadultst

which we can rewrite in matrix form like

vt+1(adultst+1childrent+1)==A(1110)⋅⋅vt(adultstchildrent)

Below, press "Forward" to step ahead a minute. The total population is the Fibonacci Sequence.

childrenadults

012012v₀childrenadults
reset 
forward
1 child + 0 adults = 11123581321345589144233

As you can see, the system goes toward the grey line, which is an eigenspace with λ=(1+5√)/2>1.

Steady States

Suppose that, every year, a fraction p of
New Yorkers move to California and a fraction q of
Californians move to New York. Drag the circles to decide these fractions and the number starting in each state.

New YorkCalifornia1 − p = 0.7p = 0.3q = 0.11 − q = 0.938.33m19.65m

To understand the system better, we can start by writing it in matrix terms like:

vt+1(New
Yorkt+1Californiat+1)==Avt(1−pqp1−q)⋅(New
YorktCaliforniat)

It turns out that a matrix like A,
whose rows add up to zero (try it!), is called a Markov matrix, and it always has λ=1 as
an eigenvalue. That means there's a value of vt for
which Avt=λvt=1vt=vt.
At this "steady state," the same number of people move in each direction, and the populations stay the same forever. Hover over the animation to see the system go to the steady state.

90% stay70% stayHover over to play/restart0m10m20m30m40m50m0m10m20m30m40m50mCaliforniaNew Yorkv₀

Complex eigenvalues

So far we've only looked at systems with real eigenvalues. But looking at the equation Av=λv,
who's to say λand v can't
have some imaginary part? That it can't be a complex number? For example,

(1−111)⋅(1i)=(1+i)⋅(1i).

Here, 1+i is
an eigenvalue and (1,i) is
an eigenvector.

If a matrix has complex eigenvalues, its sequence spirals around (0,0).
To see this, drag A's
columns (the arrows) around until you get a spiral. The eigenvalues are plotted in the real/imaginary plane to the right. You'll see that whenever the eigenvalues have an imaginary part, the system spirals, no matter where you start things off.

steps: -3-2-1123-3-2-1123-33-33-33-33realimrealimλ₀λ₁

Learning more

We've really only scratched the surface of what linear algebra is all about. To learn more, check out the legendary Gilbert Strang's Linear
Algebra
 course at MIT's Open Courseware site. To get more practice with applications of eigenvalues/vectors, also ceck out the excellent Differential
Equations
 course.

For more explanations, visit the Explained Visually project homepage.

Or subscribe to our mailing list

Eigenvectors and eigenvalues的更多相关文章

  1. A Beginner’s Guide to Eigenvectors, PCA, Covariance and Entropy

    A Beginner’s Guide to Eigenvectors, PCA, Covariance and Entropy Content: Linear Transformations Prin ...

  2. <<Numerical Analysis>>笔记

    2ed,  by Timothy Sauer DEFINITION 1.3A solution is correct within p decimal places if the error is l ...

  3. opencv 61篇

    (一)--安装配置.第一个程序 标签: imagebuildincludeinputpathcmd 2011-10-21 16:16 41132人阅读 评论(50) 收藏 举报  分类: OpenCV ...

  4. A geometric interpretation of the covariance matrix

    A geometric interpretation of the covariance matrix Contents [hide] 1 Introduction 2 Eigendecomposit ...

  5. OpenCV LDA(Linnear Discriminant analysis)类的使用---OpenCV LDA演示样例

    1.OpenCV中LDA类的声明 //contrib.hpp class CV_EXPORTS LDA { public: // Initializes a LDA with num_componen ...

  6. 三种方法实现PCA算法(Python)

    主成分分析,即Principal Component Analysis(PCA),是多元统计中的重要内容,也广泛应用于机器学习和其它领域.它的主要作用是对高维数据进行降维.PCA把原先的n个特征用数目 ...

  7. Oja’s rule

    目录 Oja's rule 背景 Hebbian learning 主要的一些理论 论文里面一些主要的假设 引理1 引理2 引理3 定理1 LEMMA 3(ALL) 引理 4 定理 2 定理 3(关于 ...

  8. <Numerical Analysis>(by Timothy Sauer) Notes

    2ed,  by Timothy Sauer DEFINITION 1.3A solution is correct within p decimal places if the error is l ...

  9. [BOOK] Applied Math and Machine Learning Basics

    <Deep Learning> Ian Goodfellow Yoshua Bengio Aaron Courvill 关于此书Part One重难点的个人阅读笔记. 2.7 Eigend ...

随机推荐

  1. C语言 链表的使用(链表的增删查改,链表逆转,链表排序)

    //链表的使用 #define _CRT_SECURE_NO_WARNINGS #include<stdio.h> #include<stdlib.h> #include< ...

  2. Boost_udp错误

      注意一点:当我们不同PC机间进行通信的时候,IP和端口号是不一样的.之前遇到的问题是,boost_system_error,这是因为我们在写程序的时候,发送和接收绑定了同一个端口,导致程序出错. ...

  3. js运动框架 step by step

    开启setInterval定时器之前,请先清除之前的定时器 window.onload = function() { var btn = document.getElementById('btn'); ...

  4. eclipse使用

    Eclipse 是一个开放源代码的.基于 Java 的可扩展开发平台. Eclipse 是 Java 的集成开发环境(IDE),当然 Eclipse 也可以作为其他开发语言的集成开发环境,如C,C++ ...

  5. python大数据工作流程

    本文作者:hhh5460 大数据分析,内存不够用怎么办? 当然,你可以升级你的电脑为超级电脑. 另外,你也可以采用硬盘操作. 本文示范了硬盘操作的一种可能的方式. 本文基于:win10(64) + p ...

  6. Linux 基础入门 第一周9.14~9.20

    第一节 Linux系统简介 Linux——操作系统 1.使多个用户从不同的终端同时操作主机(分时操作系统): 2.MINIX是一个功能有限的类似于UNIX的操作系统(UNIX 实现了 TCP/IP 协 ...

  7. 如何下载Hibernate

    官网: http://hibernate.org/ 打开hibernate官网,选择Hibernate ORM,点击左侧的Downloads 点击Downloads后,可以看到如下页面,右侧是各个版本 ...

  8. 安装.NET CORE

    需要安装两个包 https://github.com/dotnet/cli 1. .NET Core Installer 2. .NET Core SDK Installer

  9. 通通的最后一篇博客(附自制html5平面射击小游戏一枚)

    这是我最后一篇博客了,由于本人的人生规划吧,以后应该也写不出什么好的技术文章了,到现在在博客园写了2年, 今天一看,我也有了120个粉丝,好几万的浏览量,感谢大家的支持啊~~ 半年没有写博客了,由于半 ...

  10. Thrift搭建分布式微服务(二)

    第二篇 连接池  连接池配置,请前往Thrift搭建分布式微服务(一)  下面要介绍的其实不是单一的连接池,应该说是连接池集合.因为它要管理多个Tcp Socket连接节点,每个服务节点都有设置了自己 ...