關(guān)于Cramer 法則的再推廣與應(yīng)用

（喀什大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院，新疆 喀什 844000）

設(shè)A=（aji）∈Rn×m行滿秩，x=[x1，x2…xm]T∈Rm，b=[b1，b2…bm]T∈Rn，m≥n.給定線性方程組

若令lj=[aj1aj2…ajm]T∈Rm，則線性方程組（2）可改寫為

其中〈·，·〉表示Rm中向量的內(nèi)積.

通常，利用初等行變換將（2）的增廣矩陣

化為行最簡(jiǎn)階梯形后，通過(guò)觀察可求得（1）的某些解.特別地，當(dāng)m=n 且detA≠0 時(shí)，（1）的唯一解x=A-1b可由著名的Cramer 公式給出.但是，當(dāng)m≠n 時(shí)，（1）一般不存在用行列式表示的求解公式.基于此，本文將Cramer 法則推廣到求（1）或（2）在R 中的公式解上.

Cramer法則 設(shè)A是一階實(shí)矩陣.若detA≠0，則線性方程組（1）的解由公式

給出，其中，Ai是用[b1b2…bn]T替換A的第i列后所得到的矩陣[1].

文獻(xiàn)[2]給出了Cramer 法則的一個(gè)簡(jiǎn)單而有趣的推廣，即

Burgstahler定理 若線性方程組

有唯一解x1，x2，…，xn，則對(duì)?λi∈R，有

1 主要結(jié)果

定理1對(duì)?b∈Rn，線性方程組（1）有解的充分必要條件是

在定理1 有解的情況下，（1）的解由公式

給出，其中A*是A的轉(zhuǎn)置（在復(fù)數(shù)域C 中是的共軛轉(zhuǎn)置）.

推論 若n=m，則線性方程組（1）、（2）和（3）由公式（9）給出的解可寫成

形式，并且它還是（1）、（2）和（3）的最小范數(shù)解.

證明 注意到，當(dāng)n=m 時(shí)，（9）與Cramer 公式一致.事實(shí)上，（9）可以由公式

給出，其中（AA*）j是將AA*的第j 個(gè)列向量用b=[b1b2…bm]T替換后所得到的矩陣.同時(shí)，（11）還是最小范數(shù)解[3]，即

特別地，‖x‖=‖w‖且Aw=b?x=w.

定理2對(duì)?b∈Rn，線性方程組（1）有解的充分必要條件是A的行向量組{l1，l2，…，ln}在Rm中線性無(wú)關(guān).在有解的情況下，（1）的解由公式

給出，其中向量集（v1，v2，…，vn）（vj=[vj1，vj2，…vjm]T，j=1，2，…，n）可由Gram-Schmidt 正交化過(guò)程得到[1]，系數(shù)c1，c2，…，cn由公式

給出.

2 結(jié)果證明

利用Moore-Penrose 求逆公式和Cramer 法則不難證明公式（11）和（13）.但是，為便于理解，這里給出其在現(xiàn)行的《線性代數(shù)》教科書和許多文獻(xiàn)中并不常見(jiàn)的直接證明.

用〈x，y〉表示Rm中向量x和y的Euclid 內(nèi)積，表示向量x的范數(shù)；Im（G）和Ker（G）分別表示集合G 的值域和核.

引理1[4]設(shè)W 和Z 是Hilbert 空間.若算子G∈L（W，Z）和G*∈L（Z，W）相伴，則下列陳述成立：

（i）Im（G）=Z?存在γ＞0，使得

定理1的證明 因矩陣A可看作線性算子A∶Rm→Rn，即A∈L（Rm，Rn），故A的相伴算子A*∶Rn→Rm是A的轉(zhuǎn)置.對(duì)?b∈Rn，（1）有解的充分必要條件是算子A為滿射[5].由引理1 可知，存在γ＞0，使得‖A*z‖Rm≥γ‖z‖Rn（?z∈Rn）.因此，〈AA*z，z〉≥γ2‖z‖2Rn（?z∈Rn），這表明AA*是雙射.又因AA*是一n×n階實(shí)矩陣，故det（AA*）≠0.反過(guò)來(lái)，若det（AA*）≠0，即（AA*）-1存在，則對(duì)給定的b∈Rn，x=A*（AA*）-1b就是Az=b的一個(gè)解.

因z=（AA*）-1b是方程組（AA*）w=b的唯一解，故由定理1 可得

其中（AA*）j是將AA*的第j 個(gè)列向量用b=[b1b2…bn]T替換后所得到的矩陣.因此，（1）的解x=A*（AA*）-1b可寫成

的形式.此外，它還是最小范數(shù)解.事實(shí)上，若取w∈Rm使得Aw=b，則有

另一方面，

因此，‖w‖2-‖x‖2=‖w-x‖2≥0，即‖x‖≤‖w‖；特別地，當(dāng)x＝w時(shí)，有‖x‖＝‖w‖.

的解.為此，令l1=[-1 -1 1 1]T，l2=[-1 1 -1 1]T，l3=[1 -1 -1 1]T，因 {l1，l2，l3} 在R3中正交，故（24）的解為x1=-1/4，x2=-1/4，x3=-1/4，x4=3/4.

3 變分解法

定理1 和定理2 雖然給出了（1）的最小范數(shù)解的公式，但它們并非求解的唯一途徑.下面介紹求解（1）的變分法，其解為二次函數(shù)j∶Rn→R，即

的最小值.

引理2[6]對(duì)給定的b∈Rn，線性方程組（1）有解x∈Rm的充分必要條件為

引理3若二次函數(shù)j 有最小值ξb∈Rn，則

是（1）的一個(gè)解.

注：實(shí)際上，（27）是j 有臨界點(diǎn)的最優(yōu)條件.同時(shí)，可以驗(yàn)證由公式（27）和（9）給出的解完全一致.

定理3線性方程組（1）有解的充分必要條件是對(duì)?b∈Rn，二次函數(shù)j 有最小值.

證明 設(shè)（1）有解.因A是Rm到Rn的滿射，故由引理1 可知?γ＞0，使得‖A*ξ‖2≥γ2‖ξ‖2（ξ∈Rn）.于是因此，.這是j 存在最小值的保證.

下面給出一個(gè)不能應(yīng)用定理1 和定理2，卻能應(yīng)用引理2 的例子.