官术网_书友最值得收藏!

第四節 克萊姆(Cramer)法則

本節討論有n個未知量也有n個方程的線性方程組

   (1.16)   

稱上式為n元線性方程組。這是我們前面討論的二元、三元線性方程組的更一般情形。

若右端的常數b1b2,…,bn不全為零,則稱式(1.16)為非齊次線性方程組;而當b1b2,…,bn全為零,則稱式(1.16)為齊次線性方程組。記

D為線性方程組的系數行列式,又記

Dj是用方程右端的常數項b1b2,…,bn來替換系數行列式D中的第j列的元素而得到的行列式。

定理1.2 [克萊姆(Cramer)法則]如果線性方程組(1.16)的系數行列式D≠0,則方程組有唯一解

   (1.17)   

證 先證式(1.17)是方程組(1.16)的解。為此只要驗證式(1.17)滿足方程組(1.16)中的每一個方程,即要證對任意i(1≤in)有

或  ai1D1+ai2D2+…+ainDn=bi·D  (1.18)

現將每個Djj=1,2,…,n)按第j列展開為

注意到系數行列式D與行列式Dj僅第j列的元素不相同,所以在上述展式中,Dj的第j列元素的n個代數余子式A1jA2j,…,Anj也就是系數行列式D的第j列對應元素的代數余子式。

現將這些展開式代入式(1.18)的左邊,即有

把上面的和式按照常數項b1b2,…,bn重新整理,得

式中最后一個等式應用了本章第二節中式(1.15)(關于行列式展開)的結論。

再證方程組解的唯一性。

x1=λ1x2=λ2,…,xn=λn為方程組(1.16)的任一解,我們證明必有

   (1.19)   

因為λ1λ2,…,λn是式(1.16)的一個解,所以它滿足式(1.16),即

ai1λ1+ai2λ2+…+ainλn=bii=1,2,…,n

從而由行列式的性質1.4和性質1.6知

再把從第2列開始的每一個第j列的λj倍(j=2,…,n)都加到第一列上去,有

同理可證j=2,…,n時,λj·D=Dj,即。所以方程組(1.16)的解是唯一的。證畢。

【例1.12】 求解線性方程組

 方程組的系數行列式

由克萊姆法則知它有唯一解。又因為

所以方程組的解是

在上面的解題過程中,為了避免四階行列式的復雜計算,也可以利用國際知名的數學與工程軟件MATLAB來計算這些行列式的值(關于MATLAB及其基本使用方法可以參見第二章的第八節)。例如對行列式D1,在MATLAB的命令窗口輸入:

?D1=det([1 1 2 3;-4 -1 -1 -2;-6 3 -1 -1;-4 2 3 -1])%矩陣記號[…]可參見第二章回車,即可得到其值。

線性方程組(1.16)每個未知數的值都等于零的解稱為零解。反之,至少有一個未知數不等于零的解稱為非零解

任何齊次線性方程組總是有解的,因為它至少有零解。那么齊次線性方程組在什么情況下才有非零解呢?

由克萊姆法則可知,方程個數與未知量個數相等的齊次線性方程組,當它的系數行列式D≠0時,方程組有唯一零解,所以齊次線性方程組有非零解的必要條件是系數行列式D=0。以后還可以證明D=0也是齊次線性方程組有非零解的充分條件。

【例1.13】 討論λ為何值時,線性方程組

有唯一解,并求出其解。

 方程組的系數行列式

由此可知,當λ≠1且λ≠-2時,方程組有唯一解,又

此時方程組的(唯一)解是

兩點注意:

(1)上面解題過程中要計算4個帶有符號參數λ的3階行列式的值,工作量也不小,對此同樣可以利用MATLAB來計算這些行列式的值。例如對系數行列式D,在MATLAB的命令窗口輸入

?syms t                % 將符號參數λ換為英文字母符號t

?D=det([t 1 1;1 t 1;1 1 t])   % 矩陣記號[…]可參見第二章

?D=factor(D)             % 因式分解

即得D=(t+2)*(t-1)^2。

(2)對于例1.13中的方程組,當λ=1時,方程組的系數行列式D=0,方程組與下列方程同解

x1+x2+x3=1

從而方程組有無窮多組解;而當λ=-2時,原方程組即為

這是個不相容方程組,即矛盾方程組(把第一、第二個方程加到第三個方程上去,得到的是一個矛盾方程),所以方程組此時無解。

值得指出的是,在實際生產決策過程中由于數據統計誤差等原因,使得某種決策方案所對應的線性方程組為矛盾方程組,這是經常會發生的。因而矛盾方程組往往也是需要求解的,這就是所謂的“最小二乘解”。需要了解這方面知識的讀者,可參見本書附錄部分的案例二。在該案例中,還同時介紹了“多項式曲線擬合”問題。

克萊姆法則是線性方程組理論中的一個重要結果,它不僅給出了線性方程組(1.16)有唯一解的條件,方程組解的具體表達式,并且還揭示了方程組的解與系數和常數項之間的關系,這在方程組的更一般問題的研究中將起到很重要的作用。

當然,從前面的例子我們也看到,用克萊姆法則求解有n個未知量n個方程的線性方程組時,需要計算n+1個n階行列式,這個計算量是很大的。所以在實際求解高階線性方程組時,一般不用克萊姆法則。

此外,在實際問題中還會碰到系數行列式D=0,以及方程個數m與未知數個數n不相等的線性方程組。在第四章中,我們將對最一般形式的線性方程組作進一步的研究。

主站蜘蛛池模板: 德州市| 祁东县| 海口市| 康乐县| 罗源县| 含山县| 台山市| 淅川县| 孟州市| 丁青县| 淮阳县| 天门市| 资兴市| 龙口市| 清原| 阿拉善右旗| 三台县| 舟山市| 达尔| 无极县| 栖霞市| 上饶市| 云梦县| 海城市| 安达市| 乌拉特前旗| 洪湖市| 玛多县| 姚安县| 琼中| 普宁市| 民勤县| 邯郸县| 奉化市| 和硕县| 尼木县| 额敏县| 拉孜县| 阳信县| 永修县| 富宁县|