12/05/2015

[線性代數] 線性算子 與 特徵值/特徵向量(0)

======================
Definition: Linear Operator
VVn 維 向量空間 且 L:VV 為線性轉換 (Linear transformation):亦即給定任意兩向量 u,vVcR (or cC)滿足
L(u+v)=L(u)+L(v)L(cu)=cL(u)則我們稱該 L:VV 為定義在 V 上的線性算子 (Linear Operator)
======================

那麼我們現在想問一個基本問題:是否可以找到 一組非零向量 v0 與 純量 λR (or C) 使得
L(v)=λv
此問題在工程領域有諸多應用,一般而言上述問題又稱為特徵值問題。

Comments:
1. 上述討論中所提及的 Linear Operator 僅僅表示 domain 與 codomain 都為同一個向量空間 V,其餘皆與線性轉換定義相同。也就是說若我們將 domain V 與 codomain W 設為不同的向量空間,且若 L:VW 滿足
L(u+v)=L(u)+L(v)L(cu)=cL(u)則我們稱 L 為線性轉換 (Linear Transformation)。

2. 若 v=0L(v)=λv 自動滿足故我們只需關心 v0 的情況
3. 注意到若 V:=Rn 或者複數向量空間 V=Cn 則我們可從幾何觀點來看上述特徵值問題,則此問題變成決定是否 L(v)v 平行。
4. 給定任意線性算子 L:VV ,其(實數)特徵值與其對應的特徵向量不一定存在,比如說考慮 旋轉轉換 L:R2R2 滿足
L([x1x2])=[cosθsinθsinθcosθ][x1x2]其中 0<θ<π 則此算子為線性算子但不存在(實數)特徵值與特徵向量。但存在 (複數)特徵值與複數特徵向量

Example 1:
令 Linear operator L:R2R2 滿足
L(v)=L([v1v2]):=[v1v2]試求 v0λ 使得 L(v)=λv?

Solution
觀察
L(v)=L([v1v2]):=[v1v2]因為我們要求 L(v)=λv
[v1v2]=λ[v1v2]λ[v1v2][v1v2]=0[λ100λ+1][v1v2]=[00]由上式可知我們在決定 λ 使得並求解 Null Space of [λ100(λ+1)] 故若我們選 λ:=λ1=1 則對應的 v1=[s0]T 其中 sR1

另外若選 λ:=λ2=1 則對應的 v2=[0t] 其中 tR1


由上述討論所求出的 λv 即為所謂特徵值與特徵向量,現在我們引入 eigenvalue 與 eigenvector 定義


======================
Definition:  Eigenvalue and Eigenvector
Vn 維度向量空間且 L:VV 為定義在 V 上的 線性算子。我們稱 λL 的特徵值( eigenvalue of L) 若存在一組非零向量 xV 使得
L(x)=λx且 任意非零向量 x 滿足上式稱為 L 對應於特徵值 λ特徵向量 (eigenvector of L associated with the eigenvalue λ)
======================

Remark: 
1. 上述定義中的純量 λ 與向量 x皆可為 實數或者複數。有興趣的讀者請看 Example 2
2. 若我們允許 x=0 則 任意 λ 都可為 eigenvalue 因為
L(x)=λxL(0)=λ0則任何 λ 都滿足上式。
3. 事實上特徵值與特徵向量在無窮維向量空間 亦可被定義,有興趣讀者請看 Example 3


===================
Fact: L:VV 為線性算子,若 x 為 eigenvector of L associated with the eigenvalue λ 則 對任意實數 cR1 我們有
L(cx)=λ(cx)===================

Proof:
由於 x 為 eigenvector of L associated with the eigenvalue λ 我們有 對任意非零向量 x0
L(x)=λ(x) 現在觀察
L(cx)=cL(x)=c(λx)=λ(cx)讀者應注意到上述第一條等式應用 L 為線性算子的性質。


以下我們用一個例子來說明何時會發生複數的 eigenvalue λ 與 負數特徵向量 x

Example 2:
考慮線性算子 L:R2R2 滿足
L([x1x2])=[x2x1]試求 eigenvalue of L 與對應的 eigenvector

Solution
回憶前述定義,線性算子 L 的特徵值 λ 必須滿足 L(x)=λx 故現在觀察
L(x)=λxL([x1x2])=[x2x1]=λ[x1x2]λ[x1x2][x2x1]=[00][λ11λ][x1x2]=[00]觀察上述線性系統方程,求解 Null Space of [λ11λ] 可得到以下結果:

λ=λ1=i 則對應的特徵向量可透過將 λ=λ1 帶回上述線性系統方程並求解 x1 如下 x1=s[i1]T,sR1
λ=λi=i 則對應的特徵向量可透過將 λ=λ2 帶回上述線性系統方程並求解 x2 如下
x2=t[i1]T,tR1
注意到上述結果之中,L 的 特徵值皆為複數:亦即 λiC 對任意 i=1,2 且對應的特徵向量亦為 C2 複數向量。且這表示原本題目之中要求 L:R2R2 並無法找到對應的複數向量除非我們更改 domain 與 codomain 使其變成 C2


Example 3: Eigenvalue Problem in Function Space
考慮 V:=C(R) 亦即 V 為所有 單變數實數函數所成的向量空間 且我們假設其上的函數任意階導數存在。現在令 L:VV 為線性算子滿足
L(f):=f
則我們想問是否可找到 常數 λ 與 函數 f0fV 使得 L(f)=λf 成立?

 Solution
fV 為單變數函數,注意到 L(f)=f 為線性算子(why?),現在觀察 L(f)=f=λf 亦即 我們要找出 λ 與對應的 f0,fV 滿足
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=λfdf(t)dt=λf(t)df(t)f(t)=λdt1f(t)df(t)=λdtlogf(t)+C=λt+Df(t)=eλteDC:=Keλt

1 則留言:

  1. https://userdisk.webry.biglobe.ne.jp/020/691/47/N000/000/007/155428072189143651179.gif

    回覆刪除

[人工智慧] 本地端 DeepSeek R1 快速安裝:以 Macbook Pro M4 Chip為例

最近火熱的 DeepSeek R1 模型由於採用了 distill 技術,可以大幅降低計算成本,使得一般人有機會在自家筆電上跑性能逼近 Open AI ChatGPT o1的大語言模型。本文簡單介紹一步安裝在 Macbook Pro 的方法以及使用方法,以下測試採用 Macboo...