矛盾方程組及其求解
由線性代數理論知,求解線性方程組時,若方程式的個數多於未知數的個數,則方程組往往無解,此類方程組稱為 矛盾方程組(或 超定方程組)。最小二乘法是用來解矛盾方程組的一個常用方法 。
設有矛盾方程組
![矛盾方程組](/img/c/806/wZwpmL3QDNxgzMzkDOxYjN1UTM1QDN5MjM5ADMwAjMwUzL5gzL0IzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
或寫為
![矛盾方程組](/img/1/5a5/wZwpmLzQzMxQjNwcTNxYjN1UTM1QDN5MjM5ADMwAjMwUzL3UzLxUzLt92YucmbvRWdo5Cd0FmLxE2LvoDc0RHa.jpg)
通常找不到能同時滿足方程組(1)的解,因此我們轉而去尋求在某種意義下的近似解,這種近似解不是指對精確解的近似(因為精確解並不存在),而是指尋求各未知數的一組值,使方程組(1)中各式能近似相等。這就是用 最小二乘法 解矛盾方程組的基本思想。把近似解代入方程組(1)後,只能使各方程式的兩端近似相等,不妨記各個方程式兩端之差為
![矛盾方程組](/img/2/784/wZwpmLyYTNyUTMxkzNxYjN1UTM1QDN5MjM5ADMwAjMwUzL5czLxQzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
並稱該差為 偏差。按最小二乘原理,採用使偏差的平方和
![矛盾方程組](/img/d/441/wZwpmLyQDO3IjM2MTOxYjN1UTM1QDN5MjM5ADMwAjMwUzLzkzL2QzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
達到最小值來作為衡量一個近似解的近似程度的標誌 。
相關概念及定理
![矛盾方程組](/img/7/937/wZwpmL2EzM3YjNwIjN0MTN1UTM1QDN5MjM5ADMwAjMwUzLyYzL0MzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
定義1 如果 的取值使偏差平方和即式(2)達到最小,則稱這組值是 矛盾方程組(1)的最優近似解 。
預備知識:
(1) 矩陣的秩。設矩陣 A有一個n階子式D,其不等於0,而n+1階子式皆為0,那么n稱為 A的秩。
![矛盾方程組](/img/5/e01/wZwpmL3YzN2ITOwATOxYjN1UTM1QDN5MjM5ADMwAjMwUzLwkzL1QzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
(2) 對稱矩陣。如果n階方陣 A ,滿足 A=A,即
![矛盾方程組](/img/b/285/wZwpmL2cTO0ITN1cDOxYjN1UTM1QDN5MjM5ADMwAjMwUzL3gzLwYzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
則稱 A為對稱矩陣。
![矛盾方程組](/img/8/400/wZwpmL0ATMyIzNxITOxYjN1UTM1QDN5MjM5ADMwAjMwUzLykzL2gzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
![矛盾方程組](/img/f/673/wZwpmL1QTN1kzN1EDOxYjN1UTM1QDN5MjM5ADMwAjMwUzLxgzL4czLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![矛盾方程組](/img/b/6e5/wZwpmL2YTN2YjM1ITOxYjN1UTM1QDN5MjM5ADMwAjMwUzLykzL3gzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![矛盾方程組](/img/6/b76/wZwpmLyQzN0UjMxcDMxMzM1UTM1QDN5MjM5ADMwAjMwUzL3AzLzUzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
(3) 正定矩陣。設有實二次型 ,如果對於任何 ,都有 ,則稱 為 正定二次型,對稱矩陣 A是正定的。
如果 A為正定矩陣,則 A的特徵值皆為正的。且 A的各階主子式皆為正。
![矛盾方程組](/img/3/2f8/wZwpmLwIDNxAjMxMTMzEDN0UTMyITNykTO0EDMwAjMwUzLzEzLyIzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![矛盾方程組](/img/3/2f8/wZwpmLwIDNxAjMxMTMzEDN0UTMyITNykTO0EDMwAjMwUzLzEzLyIzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![矛盾方程組](/img/3/2f8/wZwpmLwIDNxAjMxMTMzEDN0UTMyITNykTO0EDMwAjMwUzLzEzLyIzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
![矛盾方程組](/img/3/2f8/wZwpmLwIDNxAjMxMTMzEDN0UTMyITNykTO0EDMwAjMwUzLzEzLyIzLt92YucmbvRWdo5Cd0FmL0E2LvoDc0RHa.jpg)
(4) 矩陣的特徵值。設 A為n階方陣,如果存在常數及非零的n維列向量 X,使 AX= X成立,則稱是方陣 A的特徵值。非零向量 X稱為方陣 A的屬於特徵值的特徵向量。
![矛盾方程組](/img/9/fd4/wZwpmLyUDM0ADM0IDO0YzM1UTM1QDN5MjM5ADMwAjMwUzLygzL4gzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
![矛盾方程組](/img/f/b11/wZwpmL3czNzEDMyYzNxYjN1UTM1QDN5MjM5ADMwAjMwUzL2czL0AzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
定理1 設n元實函式 在點 的某個鄰域內連續,且有一階及二階連續的偏導數,如果
![矛盾方程組](/img/f/00f/wZwpmL4YzN1QzN5ITOxYjN1UTM1QDN5MjM5ADMwAjMwUzLykzL2MzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
(1)
(2)矩陣
![矛盾方程組](/img/1/e87/wZwpmLwQDO0gjM3kzNxYjN1UTM1QDN5MjM5ADMwAjMwUzL5czLyYzLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![矛盾方程組](/img/9/1a6/wZwpmL3gTM3gzM0AzNxYjN1UTM1QDN5MjM5ADMwAjMwUzLwczL2czLt92YucmbvRWdo5Cd0FmLyE2LvoDc0RHa.jpg)
![矛盾方程組](/img/9/fd4/wZwpmLyUDM0ADM0IDO0YzM1UTM1QDN5MjM5ADMwAjMwUzLygzL4gzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
是正定矩陣,則 是n元實函式 的極值。
![矛盾方程組](/img/e/8f1/wZwpmLxMDO4cjN0YjNxYjN1UTM1QDN5MjM5ADMwAjMwUzL2YzL3MzLt92YucmbvRWdo5Cd0FmLzE2LvoDc0RHa.jpg)
定理2 設非齊次線性方程組 Ax=b的係數矩陣 A= ,若rank A=n,則
(1) 矩陣 AA是對稱正定矩陣;
(2) n階線性方程組 A Ax= Ab有唯一的解。
下面我們討論此二次函式是否存在最小值,若存在,如何求最小值?由高等數學可知,有以下定理。
定理3 設矛盾方程組(1)的係數矩陣 A的秩為n,則二次函式
![矛盾方程組](/img/0/17c/wZwpmLyQjN1czMwIzNxYjN1UTM1QDN5MjM5ADMwAjMwUzLyczL3UzLt92YucmbvRWdo5Cd0FmLwE2LvoDc0RHa.jpg)
一定存在最小值。
通常稱線性方程組 A Ax= Ab為 正則方程組。只要矛盾方程組(1)的係數矩陣 A的秩rank A=n,則由定理3可以得出:
(1)矛盾方程組(1)的最小二乘解存在;
(2)正則方程組 A Ax= Ab有唯一解,此解就是矛盾方程組(1)的最小二乘解。
用最小二乘法求解矛盾方程組 Ax= b的步驟歸納如下:
(1)計算 AA和 Ab,得正則方程組 A Ax= Ab;
(2)求解正則方程組,得出矛盾方程組的最優近似解 。