基本介紹
- 中文名:吉文斯法
- 外文名:Given's method
- 所屬學科:數學
- 所屬問題:線性方程組數值解法
- 別名:平面旋轉變換
- 定義:正交分解法的一種
基本介紹,吉文斯變換,
基本介紹
雅可比法是利用矩陣對角化的方法得到實對稱矩陣A的特徵值問題的解,它是基於坐標旋轉概念的一種疊代法,且當A的對角元素為主元時十分有效。吉文斯法是基於坐標旋轉概念的另一種方法,二者相比,吉文斯法實際上不能得出特徵值問題的解而只是把實對稱矩陣A簡化成三對角形式。關於三對角矩陣特徵值問題的解必須用另一種方法得到。吉文斯法的顯著優點是簡化成三對角形式時不需要進行疊代,而只需要有限的幾步運算。
雅可比方法(Jacobi's method),它系統地消除矩陣的非對角項,將對稱矩陣轉化為對角陣。遺憾的是,由於每消去一個非零元素的同時,常會將之前為零的元素變成非零的,因此該方法所需要的運算元是無窮多次。儘管經過有限次數的運算之後,矩陣的非對角元素仍不可能完全為零,但矩陣最終將收斂到對角型。因而,雅可比方法需要反覆疊代,直到非對角項“足夠”少為止。
吉文斯法(Given's method)也是一種將對稱矩陣轉化為簡單形式的方法。與雅可比方法不同的是,該方法所對應的簡單形式是三對角陣。另一個不同之處在於,吉文斯法只在殘餘的非對角位置產生零元素。因此,它所需要的運算元是有限的,從而比雅可比方法的效率更高。
豪斯霍爾德法(Householder's method)也是一種將對稱矩陣轉化為三對角形式的方法。這種方法的運算元是有限的,並且,由於它同時將整行和整列的非對角元素變成零元素,所以它比吉文斯法更為高效。
一旦用吉文斯法或豪斯霍爾德法將矩陣化為三對角系統,那么剩下的工作就是求出特徵值。比較直接的做法是,將行列式展開,從而得到一連串的多項式,通過疊代可以求出特徵值。
除了對稱矩陣之外,一般矩陣的所有特徵值也可以通過許多方法求解。它們包括:Rutishauser 的LR 方法(LR method)和Francis的QR 方法(QR method)。雖然QR方法的效率相對低一點,但由於它更加穩定,所以通常是優先使用的。如此,該方法被認為是最通用的求解方法。
最後要提到的一點是,上述方法通常是結合它們各自的特點來使用的。例如,吉文斯法和豪斯霍爾德法也可用於非對稱矩陣,不過所得結果不是三對角陣,而是一類被稱為海森伯格型(Hessenberg form)的特殊形式。另有一種方法利用了豪斯霍爾德法的速度優勢,即先用豪斯霍爾德法將矩陣轉化為三對角形式,再根據穩定的QR算法來求解特徵值。
吉文斯變換
吉文斯變換(Givens transformation)亦稱平面旋轉變換,數值代數的基本工具之一,它是一種正交變換,最常用的吉文斯變換是使變換後的向量的某個分量(例如第k個)變為零。
設
,則變換
![](/img/7/c66/wZ2NnLihjYmdzNjhDMxcTOilDO3MGN1AjM4MDMlNjZ1UmM0Y2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/e/fa2/wZ2NnLiRzYjVDO3EWMyAjY2ImZyUWZ2EGZhRTOhhjZzQWZ5Y2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/c/67f/wZ2NnLyMGZ0QzMwMTNmhDNkRGNkdjY3kjM0cDNjNTMzADN3EzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/b/ce3/wZ2NnLxITOmdDZwImNkdTMzEmMkRDMlhjY3EjZ5MmYwgTN3EzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/5/80a/wZ2NnL4QzMmBTZmFmY2MjN1I2Y1YTOwYGO2Y2Y5QDZmFTO0AzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/f/5d7/wZ2NnLlZTN0YjYhNWY1cDZ5ATNwIGO5IWZ5UjMiRzNygDNxE2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/b/5a8/wZ2NnLzEWN2ATOjVDOjNTMlZTNldzYjRjNmRzMmV2NmVjYxIzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/4/c73/wZ2NnL4AzNmVmNkNjNjhDZlRDZmFWY3MTMyAjMiJGNzIjMykzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/7/91b/wZ2NnL3IGOlBDOiljZzgDN3QTMkdjYmRTYkVmYjZzN5UmN3E2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/5/5db/wZ2NnL1QTO1MmN0MGM3gTY2UDO0gzYykTOkhzN3ATY5YGO4Y2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/9/707/wZ2NnLyMDNhZzN4ImNwQzY1EDZkBDNjVTN1IWN2UmMyczN3E2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
顯然,
具有性質:
![](/img/c/3e9/wZ2NnLxcjN5ETM4gzYxIWOmFWNiFTMxQjNzMmZlhjZ5ATNjdzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
(1) P與單位陣
只是在
位置元素不一樣,其他相同。
![](/img/4/ea2/wZ2NnLyMDZmdjZ1EjM2MTM5QTOhZWYlBjNmNGMyUGM4YWM2QzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/1/c29/wZ2NnL0IWMycTOmBzMiZWYyUjZxgDN3ATO1QzN4EjY5UDZzczLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
(2) P為正交矩陣
。
![](/img/4/5e3/wZ2NnL1gjZ5YTYxYDMmNjY3MWMzgTM0MzY2kDZhFGZkZjZwY2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
(3)
(左乘)只需計算第
行與第
行元素,即對
有
![](/img/6/a5f/wZ2NnLkRDZiRzN2IjZkFzYjZjZ4IDN1YmMxETN2MTMlFWMxE2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/6/8f5/wZ2NnLxcTYiJjNxMGZmZjNyUWOkJDOhhDM4UzN1EWO4ITYjBzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/1/843/wZ2NnL2AjMzU2MiZmZ4MzM5kzM0QDMxITZ5UzNykzYmdTO5M2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/d/4fe/wZ2NnLwQjY4gTO2ETO0gzM3QzM0UWN1ADZ1Q2MhBTN4ATNkF2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/7/995/wZ2NnLhZDM1UWYhhTOzYmYxkjNxkjYzMjNiFGZhVDNwIWM2EzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/d/bff/wZ2NnLzkTY0E2N3UTMhV2MxYWZiVDZ5YTO5QmYmhjZmhzY2U2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
(4)
(右乘)只需計算第i列與第j列元素
![](/img/c/eae/wZ2NnL1MTY5QmZlBTYjhjM1AjN4AzMmZmMlJWYxcjN3UWY4gzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/0/cb3/wZ2NnLiN2YwAzYiVDZzQDNiZmZiljZ1U2M0UzM2ITZldTO1E2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/8/d90/wZ2NnL0Q2YyUzM2UDNxImNlFWZjNWZ4IDM0I2NkVDOxUzN2QzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
定理 (約化定理)設
,其中
不全為零,則可選擇平面旋轉陣
,使
![](/img/0/e7f/wZ2NnLjFGZwETNhRTMjFGZ0cDMjlDOzUDZwMjZzEjYxEjM2AzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/2/a2b/wZ2NnLkdjNlhzM5YmNyUjM1U2N5YGOwkzM1EGM0U2MkNzNwIzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/c/3e9/wZ2NnLxcjN5ETM4gzYxIWOmFWNiFTMxQjNzMmZlhjZ5ATNjdzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/9/88a/wZ2NnLlFWZidTYhJWOjNDOmJWY2cjN3EDZklTZwQGZ3EWZllzLhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)
![](/img/8/c72/wZ2NnLjVjNiFzNyQDMjVjN1MWY5EzNyYWM5Y2MxQTYkVTZmF2LhxWdtJ3bm9SbvNmLz9mYlNmYu4GZj5yZtl2ai9yL6MHc0RHa.jpg)