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=(a11a12a1na21a22a2nam1am2amn),a=(a1a2am),x=(x1x2xn),y=(y1y2ym).

(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(aik0) 称为主元; 它位于主列和主行的交点处. 右边新格式中元素 αμναμ 按下列选主元法则计算:

(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 次选主元步骤后 (例如 y4x4,y2x1,y1x3 ) 下列左表变成右表:

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 次步骤就可得到逆矩阵 A1 .A=(351245122)重新排列元素后就得到 A1=(28211513012) . (按 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