Householder正交矩阵的新观点.docVIP

Householder正交矩陣的新觀點 陳正宗 郭世榮 李慶鋒 國立臺灣海洋大學河海工程學系，基隆，台灣 摘要 文獻中有許多的方法可建立正交矩陣，Householder利用鏡射導得正交對稱矩陣。本文提出奇數階與偶數階Householder矩陣，均可分別由eAt與eiBt導得，其中A為反對稱實矩陣，B為對稱實矩陣，t為某特定時間。本文並分別舉一階到五階的例子，說明Householder矩陣均可由本文新觀點導得。 一. 前 言 正交矩陣在工程應用及數學、力學理論推導中佔有非常重要的角色。在數學與力學上，正交矩陣是在不同觀察座標系統下線性代數中矩陣轉換的基礎，也是求解特徵值問題最有力的工具。在工程應用上，藉由尤拉角描述剛體運動的轉換矩陣也具正交性。因此如何有效率地建立正交矩陣族是我們所關心的[1-10]。 函數矩陣在生物、經濟及物理上常應用的到，陳[11]發展矩陣的餘式定理，有效的利用Cayley-Hamilton定理來計算矩陣函數。在這篇文章中，我們將利用反對稱矩陣A，透過矩陣函數eAt的方法[2,18-21]來建立奇數階的Householder正交矩陣，及透過對稱矩陣B，利用矩陣函數eiBt的方法來建立偶數階的Householder矩陣，並證明出Householder矩陣均可以此矩陣函數表示之。亦即Householder矩陣可以說是建立正交矩陣族方法中的一個特例而已。 二. Householder矩陣回顧 如果為一非零向量，存在一個n(n的矩陣，使得，此時H稱為Householder 矩陣。從此定義中，我們可以發現許多有趣的性質，例如： (a) HT=H(對稱) (b) HHT=HTH=H2=I ，其中I為單位矩陣 (c) 其中為任一向量，經由H映射後之向量。 由(c)我們可從2(2的矩陣看出其映射關係： =,吾人可得到,其中(R2. x2 x2 (′ ( Hv() (′ Hv() 鏡射平面法向量，Hv()為鏡射面法向量 (x1cos(+x2sin(=0) 圖一： 鏡射面與鏡射 矩陣相似定理與矩陣餘式定理 3.1 矩陣相似定理 存在一個n階方陣A，特徵值分別為(1,(2,…,(n，其相對應的特徵向量分別為為(1,(2…,(n，則A與D矩陣滿足相似關係，亦即AC=CD，若C-1存在，則可寫成 A=CDC-1 (1) 其中 (2) (3) 3.2 矩陣餘式定理 在實數餘式定理中，給定一函數f(x)，若除以x-a，其餘式為f(a)，吾人可寫為 f(x)=(x-a)Q(x)+f(a) (4) 其中Q(x)為商式，同理當除式為n次式時，可表示為 f(x)=(anxn+an-1xn-1+…+a1x+a0)Q(x) +rn-1xn-1+…+r1x+r0 (5) 其中rn-1xn-1+…+r1x+r0為餘式。 若將實數改為方陣A，則亦有相同的特性，即 f(A)=rn-1An-1+…+r1A+r0I (6) 其中anAn+an-1An-1+…+a1A+a0I =0，係由 Cayley- Hamilton定理導得。亦即anAn+an-1An-1+…+a1A+a0I =0為A之特徵方程式。而當A方陣具有重根特徵值時，則由於方程式數目小於未定係數的數目。因此需將(5)式做微分處理，而n次重根，就需做n次微分，以此類推。 矩陣的餘式定理非常適合用在特徵值重根的情況，因無需要導得特徵向量即可求得矩陣函數。 我們用以下的兩個例子來說明： 例1：E=，利用矩陣的餘式定理，可導得 eEt=(-6te2t-3e2t+4e3t)I+(5te2t+4e2t-4e3t)E +(-te2t-e2t+e3t)E2 (7) 其中I=，E=，E2= 因此我們可得到 (8) 例2：E=，利用矩陣的餘式定理，可導得 eEt= (9) 其中 I=，E=，E2= 因此，我們可得到 (10) 由3(3方陣導eAt與Householder矩陣關係 令=，可導得3(3的Householder矩陣= 其中=1。 若令反對稱