前面, 我们学习了线性方程组的一些基本的概念. 我们知道, 可用阵等式简单地写一个线性方程组 (或者, 形如 AX=B 的阵等式就是一个线性方程组). 本节, 我们讨论由 n 个 n 元 ⩽1 次方程作成的方程组. 用阵等式表示这类方程组, 就是 AX=B, 其中 A 是 n 级阵, B 是 n×1 阵, 且 X 是未知的 n×1 阵. A 是方阵, 故它有行列式. AX=B 的解跟 A 的行列式是否有关系? 此事的回答是 “是”. 我们讨论一个特别的情形. 这会是有用的.
设 A 是 n 级阵 (n⩾1). 设 B 是 n×1 阵. 设 X 是未知的 n×1 阵. 若 det(A)=0, 则线性方程组 AX=B 有唯一的解X=(det(A))−1adj(A)B.
在论证此事前, 我想用一个简单的例助您理解它.
设 a, b 是常数. 一个 1 元 ⩽1 次方程 ax=b 当然也是一个线性方程组: 这是方程数为 1 时的特别情形. 我们在中学就知道, 当 a=0 时, ax=b 有唯一的解 x=a−1b. 一方面, x=a−1b 是一个解: a(a−1b)=(aa−1)b=1b=b.另一方面, 若数 y 也适合 ay=b, 则y=1y=(a−1a)y=a−1(ay)=a−1b.
我们说, 此事是 Cramer 公式的一个特例. 我们可写 ax=b 为阵等式 [a][x]=[b]. (注意, [a], [x], [b] 都是 1 级阵.) Cramer 公式说, 若 det[a]=0, 则 [a][x]=[b] 有唯一的解 (注意, 1 级阵的古伴是 [1]) [x]=(det[a])−1adj([a])[b]=a−1[1][b]=a−1[b]=[a−1b].这跟我们已知的结论是一样的.
证. 我们先验证 C=(det(A))−1adj(A)B 适合 AC=B: AC=====A((det(A))−1adj(A)B)(det(A))−1(Aadj(A)B)(det(A))−1(det(A)InB)((det(A))−1det(A))(InB)B.
我们再证
AX=B 至多有一个解. 设
n×1 阵
D 也适合
AD=B. 则
D======((det(A))−1det(A))(InD)(det(A))−1(det(A)In)D(det(A))−1(adj(A)A)D(det(A))−1adj(A)(AD)(det(A))−1adj(A)BC. 考虑线性方程组 AX=B, 其中A=[1214],B=[3594],X=[x1x2].不难算出det(A)=1⋅4−2⋅1=2=0,所以, 此方程组有唯一的解X===(det(A))−1adj(A)B21[4−2−11][3594]=21[4⋅35−1⋅94−2⋅35+1⋅94]21[4624]=[2312].
我们可进一步地改写 Cramer 公式. 我们先具体地算出 adj(A)B 的每一个元:[adj(A)B]i,1==ℓ=1∑n[adj(A)]i,ℓ[B]ℓ,1ℓ=1∑n(−1)ℓ+idet(A(ℓ∣i))[B]ℓ,1.设 A{i,B} 是以 B 代阵 A 的列 i 后得到的阵. 则[A{i,B}]ℓ,k={[A]ℓ,k,[B]ℓ,1,k=i;k=i.由此可见, A(ℓ∣i)=(A{i,B})(ℓ∣i) (ℓ=1, 2, …, n). 所以==ℓ=1∑n(−1)ℓ+idet(A(ℓ∣i))[B]ℓ,1ℓ=1∑n(−1)ℓ+idet((A{i,B})(ℓ∣i))[A{i,B}]ℓ,1det(A{i,B}).从而[X]i,1====[(det(A))−1adj(A)B]i,1(det(A))−1[adj(A)B]i,1(det(A))−1det(A{i,B})det(A)det(A{i,B}).由此, 我们得到 Cramer 公式的另一个形式:
设线性方程组⎩⎨⎧a1,1x1+a1,2x2+⋯+a1,nxn=b1,a2,1x1+a2,2x2+⋯+a2,nxn=b2,…,an,1x1+an,2x2+⋯+an,nxn=bn.(注意, 方程数等于未知数数.) 记A=⎣⎡a1,1a2,1⋮an,1a1,2a2,2⋮an,2⋯⋯⋯a1,na2,n⋮an,n⎦⎤,B=⎣⎡b1b2⋮bn⎦⎤.若 det(A)=0, 则此线性方程组有唯一的解xi=det(A)det(A{i,B}),i=1, 2, …, n,其中 A{i,B} 是以 B 代阵 A 的列 i 后得到的阵.
考虑由 2 个 2 元 ⩽1 次方程作成的方程组{a1,1x1+a1,2x2=b1,a2,1x1+a2,2x2=b2.根据 Cramer 公式, 若 det[a1,1a2,1a1,2a2,2]=a1,1a2,2−a2,1a1,2=0, 则此方程组有唯一的解x1=det[a1,1a2,1a1,2a2,2]det[b1b2a1,2a2,2]=a1,1a2,2−a2,1a1,2b1a2,2−b2a1,2,x2=det[a1,1a2,1a1,2a2,2]det[a1,1a2,1b1b2]=a1,1a2,2−a2,1a1,2a1,1b2−a2,1b1.
考虑线性方程组{x1+x2=35,2x1+4x2=94.因为 det[1214]=1⋅4−2⋅1=2=0, 故, 此方程组有唯一的解x1=det[1214]det[359414]=235⋅4−94⋅1=23,x2=det[1214]det[123594]=21⋅94−2⋅35=12.
一般地, Cramer 公式只是一个理论的公式. 这是因为, 一般地, 计算行列式是复杂的事 (或许, 您还记得, 3 级阵的行列式的具体的公式含 6 项, 且 4 级阵的行列式的具体的公式含 24 项). 所以, 我们一般用别的方法 (如代入消元法、加减消元法) 解线性方程组.
还是以{x1+x2=35,2x1+4x2=94为例. 由方程 1, 有 x1=35−x2. 代入其到方程 2, 有 2(35−x2)+4x2=94, 即 2x2=24. 由此可知 x2=12. 代入其到 x1=35−x2, 有 x1=23. 最后, 经验证, (23,12) 的确是此方程组的一个解.