Skip to content

4.5.1 线性系, 选主元法 ​

4.5.1.1 线性系 ​

一个线性系含 m 个线性型

y1=a11x1+a12x2+⋯+a1nxn+a1,y2=a21x1+a22x2+⋯+a2nxn+a2,

......

ym=am1x1+am2x2+⋯+amnxn+am

或

(4.174a)y―=Ax―+a―,

其中

A=(a11a12⋯a1na21a22⋯a2n⋮⋮⋮am1am2⋯amn),a―=(a1a2⋮am),x―=(x1x2⋮xn),y―=(y1y2⋮ym).

(4.174b)

大小为(m, n)的矩阵 A 的元素 aμν 及列向量 a― 的分量 aμ(μ=1,2,⋯,n) 是常数. 列向量 x― 的分量 xν(ν=1,2,⋯,n) 是独立变量,列向量 y― 的分量 xμ(μ= 1,2,⋯,m) 是相关变量.

4.5.1.2 选主元法 ​

1. 选主元法格式 ​

如果 (4.174a) 中元素 aik 不等于零,那么在所谓选主元步骤中,变量 yi 可调换为独立变量,而变量 xk 可调换为相关变量. 选主元步骤是主元法的基本组成部分, 应用它可以解 (例如) 线性方程组和线性最优化问题. 通过对应于方程组 (4.174a) 的格式:

x1

x2

...

xk

...

xn

1

y1

a11

a12

...

a1k

...

a1n

a1

y2

a21

a22

...

a2k

...

a2n

a2

⋮

⋮

⋮

⋮

⋮

⋮

yi

ai1

ai2

...

aik

...

ain

ai

⋮

⋮

⋮

⋮

⋮

⋮

ym

am1

am2

...

amk

...

amn

am

xk

αi1

αi2

αik

...

αin

αi

选主元步骤得到格式:(4.175)

x1

x2

...

xk

...

xn

1

y1

α11

α12

...

α1k

...

α1n

α1

y2

α21

α22

...

α2k

...

α2n

α2

⋮

⋮

⋮

⋮

⋮

xi

αi1

αi2

...

αik

...

αin

αi

⋮

⋮

⋮

⋮

⋮

⋮

ym

αm1

αm2

...

αmk

...

αmn

αm

2. 选主元法则 ​

在格式中画了框的元素 aik(aik≠0) 称为主元; 它位于主列和主行的交点处. 右边新格式中元素 αμν 和 αμ 按下列选主元法则计算:

(1) αik=1aik ;(4.176a)

(2) αμk=aμkaik(μ=1,⋯,m;μ≠i) ;(4.176b)

(3) αiν=−aiνaik,αi=−aiaik(ν=1,2,⋯,n;ν≠k) ;(4.176c)

(4) αμν=aμν−aμkaiνaik=aμν+aμkαiν,αμ=aμ+aμkαi (对每个 μ≠i 及每个 ν≠k ).(4.176d)

为使计算容易些 (法则 4),我们将元素 αiν 写在选主元法格式的第 m+1 行 (最低行). 应用这个选主元法则可以调换其他变量.

4.5.1.3 线性相关性 ​

如果每个 yμ 都可以调换为某个独立变量 xν ,那么线性型 (4.174a) 是线性无关的 (参见第 732 页 9.1.2.3, 2.). 线性无关性将用来 (例如) 确定矩阵的秩. 不然, 可直接由选主元法格式找出相关性关系.

◼ 3 次选主元步骤后 (例如 y4→x4,y2→x1,y1→x3 ) 下列左表变成右表:

x1

x2

x3

x4

1

y1

2

1

1

0

-2

y2

1

-1

0

0

2

y3

1

5

2

0

0

y4

0

2

0

1

0

y2

x2

y1

y4

1

x3

-2

-3

1

0

6

x1

1

1

0

0

-2

y3

-3

0

2

0

10

x4

0

-2

0

1

0

因为 α32=0 ,所以不可能作进一步的调换,并且可以看到相关关系: y3=2y1− 3y2+10 . 对于另外一个主元法序列,仍然有一对不可调换的变量.

4.5.1.4 逆矩阵的计算 ​

如果 A 是一个 n×n 正则矩阵,那么对于方程组 y―=Ax― 应用选主元法实施 n 次步骤就可得到逆矩阵 A−1 .A=(351245122)⇒重新排列元素后就得到 A−1=(28−21−1−51301−2) . (按 yi 的下标重排矩阵的列, 按 xk 的下标重排矩阵的行.)

x1

x2

x3

y1

3

5

1

y2

2

4

7 5

y3

2

2

y3

x2

x3

y1

3

-1

-5

y2

2

0

7

x1

-2

-2

y3

x2

y2

y1

13

-1

-5

x3

-2

0

1

x1

5

-2

-2

y3

y1

y2

x2

13

-1

-5

x3

-2

0

1

x1

-21

2

8

version 1.24.0