2019年1月3日 星期四

GeoGebra -- 求矩陣的特徵值(eigen value)及特徵向量(eigen vector)




GeoGebra 只能用運算區求特徵值及特徵向量,代數區是不能使用 Eigenvalues() 及Eigenvectors() 這兩個函式的。


因此我們直接進入運算區,關閉指令列及其他畫面。

↓運算區畫面

利用GGB模擬手算過程,同時也瞭解運算區的操作。

↓輸入 A:={{-1,0},{1,-5}}

還是提醒一下,在運算區不是等號(=),而是冒號加等號(:=)。
在手算求特徵值的過程中,會先求 det(AλI)=0det(AλI)=0。因此我們也先輸入一個單位矩陣I。

↓再定義一個單位矩陣I

↓再輸入 AλIAλI

也許你會問: 如何輸入 λλ

↓先點選視窗右角αα,會出現一個符號表,再點選λλ

↓輸入 Determinant($3),求第3式的行列式

在運算區的視窗最左邊有一排數字,錢符號($)加上數字就代表該運算式。因此Determinant($3)就是計算第3式的行列式值。

接下來就是求取行列式為零時的λλ值,也就是特徵值。

↓輸入solutions($4=0),求λ2+6λ+5=0λ2+6λ+5=0的解,其解就是特徵值


↓當然也可以直接使用eigenvalues(A)取代上述步驟求出特徵值


接下來由特徵值來求特徵向量,在此之前,先瞭解如何取出向量(或list)中的個別元素。

↓Elements($6,1)取出式6中第1個元素,Elements($6,2)取出第2個元素

手算特徵向量就是要求Ax=λxAx=λx中的向量xx,因此我們需要一個向量(x,y)(x,y)來代表x。

↓第9式定義一個向量、第10式就是AxAx、第11式就是λxλx

其實在前面我們已經算出λλ值,因此修改第11式。

↓第12式: 將第1個特徵值取代第11式的λλ;第13式: 求解Ax=λxAx=λx

上述第13式的輸入變數有兩個,第一個是方程式,第二個指定哪些是變數。結果是(4y,y)(4y,y),取y=1y=1可得特徵向量為(4,1)(4,1);同理可輸入第2個特徵值,λ=5λ=5,再求其特徵向量。

↓第14式: 代入第2個特徵值求特徵向量,取y=1y=1可得特徵向量為(0,1)(0,1)


↓當然也可以一步到位,直接使用Eigenvectors(A)求矩陣A的特徵向量


-- END --

沒有留言:

張貼留言