Definition: Linear Operator
令 VV 為 n 維 向量空間 且 L:V→V 為線性轉換 (Linear transformation):亦即給定任意兩向量 u,v∈V 與 c∈R (or c∈C)滿足
L(u+v)=L(u)+L(v)L(cu)=cL(u)則我們稱該 L:V→V 為定義在 V 上的線性算子 (Linear Operator)
======================
那麼我們現在想問一個基本問題:是否可以找到 一組非零向量 v≠0 與 純量 λ∈R (or ∈C) 使得
L(v)=λv
此問題在工程領域有諸多應用,一般而言上述問題又稱為特徵值問題。
Comments:
1. 上述討論中所提及的 Linear Operator 僅僅表示 domain 與 codomain 都為同一個向量空間 V,其餘皆與線性轉換定義相同。也就是說若我們將 domain V 與 codomain W 設為不同的向量空間,且若 L:V→W 滿足
L(u+v)=L(u)+L(v)L(cu)=cL(u)則我們稱 L 為線性轉換 (Linear Transformation)。
2. 若 v=0 則 L(v)=λv 自動滿足故我們只需關心 v≠0 的情況
3. 注意到若 V:=Rn 或者複數向量空間 V=Cn 則我們可從幾何觀點來看上述特徵值問題,則此問題變成決定是否 L(v) 與 v 平行。
4. 給定任意線性算子 L:V→V ,其(實數)特徵值與其對應的特徵向量不一定存在,比如說考慮 旋轉轉換 L:R2→R2 滿足
L([x1x2])=[cosθ−sinθsinθcosθ][x1x2]其中 0<θ<π 則此算子為線性算子但不存在(實數)特徵值與特徵向量。但存在 (複數)特徵值與複數特徵向量
Example 1:
令 Linear operator L:R2→R2 滿足
L(v)=L([v1v2]):=[v1−v2]試求 v≠0 與 λ 使得 L(v)=λv?
Solution
觀察
L(v)=L([v1v2]):=[v1−v2]因為我們要求 L(v)=λv 故
[v1−v2]=λ[v1v2]⇒λ[v1v2]−[v1−v2]=0⇒[λ−100λ+1][v1v2]=[00]由上式可知我們在決定 λ 使得並求解 Null Space of [λ−100(λ+1)] 故若我們選 λ:=λ1=1 則對應的 v1=[s0]T 其中 s∈R1。
另外若選 λ:=λ2=−1 則對應的 v2=[0t] 其中 t∈R1。 ◻
由上述討論所求出的 λ 與 v 即為所謂特徵值與特徵向量,現在我們引入 eigenvalue 與 eigenvector 定義
======================
Definition: Eigenvalue and Eigenvector
令 V 為 n 維度向量空間且 L:V→V 為定義在 V 上的 線性算子。我們稱 λ 為 L 的特徵值( eigenvalue of L) 若存在一組非零向量 x∈V 使得
L(x)=λx且 任意非零向量 x 滿足上式稱為 L 對應於特徵值 λ 的特徵向量 (eigenvector of L associated with the eigenvalue λ)
======================
Remark:
1. 上述定義中的純量 λ 與向量 x皆可為 實數或者複數。有興趣的讀者請看 Example 2
2. 若我們允許 x=0 則 任意 λ 都可為 eigenvalue 因為
L(x)=λx⇒L(0)=λ0則任何 λ 都滿足上式。
3. 事實上特徵值與特徵向量在無窮維向量空間 亦可被定義,有興趣讀者請看 Example 3
===================
Fact: 令 L:V→V 為線性算子,若 x 為 eigenvector of L associated with the eigenvalue λ 則 對任意實數 c∈R1 我們有
L(cx)=λ(cx)===================
由於 x 為 eigenvector of L associated with the eigenvalue λ 我們有 對任意非零向量 x≠0
L(x)=λ(x) 現在觀察
L(cx)=cL(x)=c(λx)=λ(cx)讀者應注意到上述第一條等式應用 L 為線性算子的性質。◻
以下我們用一個例子來說明何時會發生複數的 eigenvalue λ 與 負數特徵向量 x
Example 2:
考慮線性算子 L:R2→R2 滿足
L([x1x2])=[−x2x1]試求 eigenvalue of L 與對應的 eigenvector
Solution
回憶前述定義,線性算子 L 的特徵值 λ 必須滿足 L(x)=λx 故現在觀察
L(x)=λx⇒L([x1x2])=[−x2x1]=λ[x1x2]⇒λ[x1x2]−[−x2x1]=[00]⇒[λ1−1λ][x1x2]=[00]觀察上述線性系統方程,求解 Null Space of [λ1−1λ] 可得到以下結果:
若 λ=λ1=i 則對應的特徵向量可透過將 λ=λ1 帶回上述線性系統方程並求解 x1 如下 x1=s[i1]T,∀s∈R1
若 λ=λi=−i 則對應的特徵向量可透過將 λ=λ2 帶回上述線性系統方程並求解 x2 如下
x2=t[−i1]T,∀t∈R1
注意到上述結果之中,L 的 特徵值皆為複數:亦即 λi∈C 對任意 i=1,2 且對應的特徵向量亦為 C2 複數向量。且這表示原本題目之中要求 L:R2→R2 並無法找到對應的複數向量除非我們更改 domain 與 codomain 使其變成 C2
Example 3: Eigenvalue Problem in Function Space
考慮 V:=C∞(R) 亦即 V 為所有 單變數實數函數所成的向量空間 且我們假設其上的函數任意階導數存在。現在令 L:V→V 為線性算子滿足
L(f):=f′
則我們想問是否可找到 常數 λ 與 函數 f≠0 且 f∈V 使得 L(f)=λf 成立?
Solution
令 f∈V 為單變數函數,注意到 L(f)=f′ 為線性算子(why?),現在觀察 L(f)=f′=λf 亦即 我們要找出 λ 與對應的 f≠0,f∈V 滿足
f′=λf 由於我們要找的 f 必須無窮維導數存在 且一次導數必須滿足 f′=λf,故我們猜 f(t)=eλt 則
ddt(eλt)=λeλt=λf(t)故現在我們找到一組 f 滿足該方程,但是否有其他人選?答案是肯定的,比如說我們改令
f(t)=Keλt仍為該方程 f′=λf 的解 (讀者可自行驗證),因此我們有以下結果,對任意 特徵值 λ∈R1,其特徵向量 f(t)=Keλt 其中 K 為任意非零常數。◻
Comment:
1. 上述例子中顯示微分方程 f′=λf 的解 f(t)=exp(λt) 剛好為該 L(f)=f′ 的對應於特徵值 λ 的特徵向量。
2. 若我們只關心 f>0 上述微分方程可直接求解不必猜測,解法如下:
f′=λf⇒df(t)dt=λf(t)⇒df(t)f(t)=λdt⇒∫1f(t)df(t)=∫λdt⇒logf(t)+C=λt+D⇒f(t)=eλteD−C:=Keλt
https://userdisk.webry.biglobe.ne.jp/020/691/47/N000/000/007/155428072189143651179.gif
回覆刪除