スライド タイトルなし

5.連立一次方程式
1
連立方程式とその解
ここでは、連立方程式を解くということを再考する。
そもそも方程式を“解く”とは、
与えられた式を満たす全ての実数の集合を
求めることである。すなわち、方程式 f ( x)  0 の解とは、
f ( x)  0 を満たす実数 x の
x  R | f ( x)  0
集合(要素が一つの場合もあ
である。
る。)
同様に、連立方程式
 f1 ( x )  0


 f ( x)  0
 m
ここで、
x
t
 x1
x は
x2
xn 
なるベクトル。
を“解く”とは、与えられた複数の式の全てを満たすベクトルの
集合を求めることである。すなわち、以下が解である。
 x  Rn | f1 ( x)  0 AND f 2 ( x )  0 AND AND f m ( x )  0 2
連立方程式とその解1
a1 x  b1 y  c1

a2 x  b2 y  c2
y
a1x  b1 y  c1
この連立方程式を“解く”とは、
2つの方程式を同時に満たす
( x, y ) の組を見つけることである。
y
解く
a2 x  b2 y  c2
O
a1x  b1 y  c1 を
満たす点の集合
x
a2 x  b2 y  c2 を
満たす点の集合
a1x  b1 y  c1
a2 x  b2 y  c2
x
O
a1x  b1 y  c1 と
a2 x  b2 y  c2 を
同時に満たす
点の集合(1点)
3
連立方程式とその解2
 c1
a1 x  b1 y

2a1 x  2b1 y  2c1
この連立方程式の答えは、
a1x  b1 y  c1 上の点全てである。
解く
a1x  b1 y  c1
21 x  2b1 y  2c1
O
a1x  b1 y  c1
y
y
c
(0, 1 )
b1
x
O
解
b1  0 とし、 k を任意のスカラーとして、
 1  0
 x
 a    c  とも表せる。

k
 y
1
1



 
 b1   b1 
(1, 
a1
)
b1
x
(x, y) | a1x  b1 y  c1
このように、無限個の
ベクトルが式を満たす
ことを不定という。 4
連立方程式とその解3
 c1
a1 x  b1 y

2a1 x  2b1 y  c2 ( 2c1 )
この連立方程式を満たすものはない。
y
y
解く
a1x  b1 y  c1
解   {}
21 x  2b1 y  c2
O
O
x
x
このように、式を満たす
ベクトルがひとつもない
ことを不能という。
5
連立方程式とその解4
a1 x  b1 y  c1 z  d1

a2 x  b2 y  c2 z  d 2
解
a2 x  b2 y  c2 z  d2
a1x  b1 y  c1z  d1
解く
z
O
z
y
x
O
y
x
6
連立一次方程式
これまでも、何度も扱ってきたが、
連立方程式を行列とベクトルを用いて表現できる。
ìï a11x 1 + a12x 2 + L + a1n x n
ïï
ïï a x + a x + L + a x
22 2
2n n
ï 21 1
í
ïï
M
ïï
ïï am 1x 1 + am 2x 2 + L + amn x n
î
éa11 a12 L
ê
êa 21 O
ê
ê
êM
ê
êêëam 1 L
a1n ù
úéx 1 ù
Múê ú
úê Mú=
úê ú
úê ú
úêëx n ú
û
amn ú
ú
û
Ax = b
= b1
= b2
= bm
éb1 ù
ê ú
êb ú
ê2ú
ê ú
ê Mú
ê ú
êbm ú
êë ú
û
7
ここで、
A = [aij ] = [a 1 a 2 L
a m ],
éb ù
ê1ú
ê ú
b = êMú
ê ú
êbn ú
êë úû
éx 1 ù
ê ú
x = êê Mú
ú,
ê ú
êëx n ú
û
このとき、各行列およびベクトルは、以下のような名称
で呼ばれる。
Ax = b
係数行列
定数項ベクトル
未知数ベクトル
通常の一次方程式と対応させてみるとよい。
3x = 2
係数
未知数
定数
8
さらに、拡大係数行列1つで連立方程式を表す。
[A | b ]
拡大係数行列
9
正則な係数行列を持つ連立方程式
(正則行列と連立一次方程式)
連立1次方程式 A x = b が一意の解を持つための
必要十分条件は、係数行列 A が正則行列である
ことである。(したがって、 A は正方行列)
証明:
十分性:(正則→一意)
A は正則なので、
逆行列A - 1 が存在する。
Ax = b
- 1
A
の両辺の左から
を乗じる。
- 1
- 1
A Ax = A b
\ Ix = A - 1b
\ x = A - 1b
積は一意なので、
解 x は一意である。
10
必要性:(一意→正則)
背理法による。
係数行列 A - 1 が正則でなくても、解が一意と仮定する。
(背理法の仮定)
このとき、 A を階段行列に変形したとき、段数が減少する。
そのときの変形行列(基本変形行列の積)を T
T を左から乗じて得られる拡大係数行列を éC
êë
T
éA b ù¾ ¾
®
êë
ú
û
éC
êë
とし、
d ù とする。
ú
û
dù
ú
û
すなわち、
C º TA
d º Tb
11
ìï a11x 1 +
ïï
ïï a x +
ïï 21 1
ïa x +
í 31 1
ïï
ïï M
ïï
ïï an 1x 1 +
ïî
a11x 1 +
L
+ a1n x n
= b1
a22x 2 +
L
+ a 2n x n
= b2
O
= b3
M
L
+ ann x n
= bn
T
ìï c11x 1 +
ïï
ïï
ïï
ïï
ïï
ïï
í
ïï
ïï
ïï
ïï
ïï
ïï
ïîï
c12x 2 +
L
+ c1n x n
= d1
c2 p x 2 +
L
+ c 2n x n
= d2
crn x r
= dr
0gx n
= dn - 1
0gx n
= dn
M
crq x r +
L
M
0gx r + 1 +
12
階段行列の段数が減少しているので、最後の式の係数は全て、
0である。よって、dn = 0
である。
しかし、最後の式の形より、 x n の値を任意な実数 k に
選んだとしても、
Cx = d
を満たす。すなわち、形が
éx 1 ù
ê ú
êx 2 ú
x = êê ú
ú
M
ê ú
ê ú
êëk ú
û
となる複数(無数)の解が存在する。
同値変形であるので、 x は
Ax = b
も満たす。これは、解の一意性に矛盾する。
QED
13
連立方程式と階数
(連立方程式と階数)
連立1次方程式 A x = b が解を持つための
必要十分条件は、
rankA = rank éêA
ë
bù
ú
û
が成り立つことである。
ここで、左辺は係数行列の階数、
右辺は拡大係数行列の階数である。
この命題は、解が一意に求まるときと、
不定のときの両方を表している。
(逆の言い方をすると、不能でない条件
を示している。)
証明略
14
連立方程式の解法
éC
拡大係数行列を階段行列化したものを ê
ë
すなわち、
階段行列化
éA b ù¾ ¾
¾ ¾¾
®
êë
ú
û
éC
êë
dù
ú
û とする。
dù
ú
û
行列 C の形により、2つの場合に分けて考える。
場合1:(簡単な場合)
C =
場合2:(複雑な場合)
C =
15
場合1:(簡単な場合)
ìï c11x 1 +
ïï
ïï
ïï
ïï
ïï
ïï
í
ïï
ïï
ïï
ïï
ïï
ïï
îï
c12x 2 +
L
+ c1n x n
= d1
c22x 2 +
L
+ c 2n x n
= d2
crn x r
= dr
0gx n
= 0
0gx n
= 0
M
crr x r +
L
M
解を持つとき
には、必ず0
になる。
この場合には、行基本変形を用いてさらに変形できる。
T'
C ¾ ¾¾
®C '
éE r
C 'r ,n - r ù
ê
ú
C '= ê
ú
O
O
n - r ,n - r ú
êë n - r ,r
û
16
é1
ê
ê0
ê
ê
êM
ê
C ' = êê0
ê
ê0
ê
êM
ê
ê0
êë
0
L
1
0 c '1,r + 1 L
M M
O
O
L
1 c 'r ,r + 1
L
0 0
L
M M
L
0 0
L
c '1,n ù
ú
ú
ú
ú
ú
ú
c 'r ,n ú
ú
ú
0 ú
ú
M ú
ú
0 ú
ú
û
この変形行列を用いて、次のように連立方程式が変形された
とする。
éC
êë
T '
dù
¾
¾
¾
®
ú
û
éC ' d 'ù
êë
ú
û
17
ìï x 1
ïï
ïï
x2
ïï
ïï
O
ïï
ïï
í
ïï
ïï
ïï
ïï
ïï
ïï
îï
+ c '1,r + 1 x r + 1
+ c '1,r + 2 x r + 2 +
L
+ c '1,n x n
=
d '1
+ c '2,r + 1 x r + 1 + c '2,r + 2 x r + 2 +
L
+ c '2,n x n
=
d '2
+ c 'r ,r + 1 x r + 1 + c 'r ,r + 2 x r + 2 +
L
+ c 'r ,n x n
=
d 'r
0gx r + 1
L
+ 0gx n
=
0
M
=
M
0gx n
=
0
M
xr
+ 0gx r + 2 +
rank A = rank éêA
ë
= rank[C d ]
bù
ú
û
= rank[C ' d ']
であることに注意する。
18
é1
ê
ê0
ê
ê
êM
ê
T
T '
A ¾ ¾ ® C ¾ ¾¾
® C ' = êê0
ê
ê0
ê
êM
ê
ê0
êë
éA
êë
T
bù
¾¾
®
ú
û
éC
êë
0
L
1
c '1,n ù
ú
ú
ú
ú
ú
ú
c 'r ,n ú
ú
ú
0 ú
ú
M ú
ú
0 ú
ú
û
0 c '1,r + 1 L
M M
O
O
L
1 c 'r ,r + 1
L
0 0
L
M M
L
é1
ê
ê0
ê
ê
êM
ê
T '
ù
é
ù
d ú¾ ¾ ¾
® êC ' d 'ú= êê0
û
ë
û ê
ê0
ê
êM
ê
ê0
êë
0 0
0
L
1
O
L
L
0 c '1,r + 1
c '1,n
c '2,r + 1
c '2,n
M
1 c '1,r + 1 L
c '1,r + 1
0 0
0
L
M M
L
0 0
M
L
0
d '1 ù
ú
d '2 ú
ú
ú
ú
ú
ú
ú
0 ú
ú
ú
M ú
ú
0 ú
ú
û
19
ìï x 1
ïï
ïï
x2
ïï
ïï
O
ïï
ïï
í
ïï
ïï
ïï
ïï
ïï
ïï
îï
+ c '1,r + 1 x r + 1
+ c '1,r + 2 x r + 2 +
L
+ c '1,n x n
=
d '1
+ c '2,r + 1 x r + 1 + c '2,r + 2 x r + 2 +
L
+ c '2,n x n
=
d '2
+ c 'r ,r + 1 x r + 1 + c 'r ,r + 2 x r + 2 +
L
+ c 'r ,n x n
=
d 'r
0gx r + 1
L
+ 0gx n
=
0
=
M
=
0
M
xr
+ 0gx r + 2 +
M
0gx r + 1
+ 0gx r + 2 +
L
+ 0gx n
このときは、連立方程式を考えると、
x r + 1, L , x n
の n - r 個の変数は任意の実数でよいことがわかる。
ìï x = k
1
ïï r + 1
ïï
M
í
ïï
ïï x n = kn - r
ïî
というように、 n - r 個任意定数を用いて、
連立方程式の解を表現できる。
20
ìï x 1
ïï
ïï
x2
ï
í
ïï
O
ïï
ïï
î
= d '1
- c '1,r + 1 k1 L
= d '2 - c '2,r + 1 k1
- c '1,n kn - r
- c '2,n kn - r
M
xr
= d 'r
- c 'r ,r + 1 k1
- c 'r ,n kn - r
このように、 x 1, L , x r は n - r 個の任意定数 k1, L , kn - r
用いて自動的に決定される。
を
21
自由度
定義(自由度)
未知数が n 個の連立一次方程式
( n 元一次連立方程式)
Ax = b
において、
rankA = rank éêA b ù
= r
ú
ë
û
とする。このとき、自由に定めることのできる未知数の数
n- r
を方程式の自由度という。
連立方程式の解は、自由度の数だけの
任意定数を用いて表現される。
22
場合2:(複雑な場合)
C =
ìï c1l1 x l1 + L +
ïï
ïï
ïï
ïï
ïï
ïï
í
ïï
ïï
ïï
ïï
ïï
ïï
îï
c1l2 x l2 +
L
+ c1n x n
= d1
c2l2 x l2 + L +
L
+ c 2n x n
= d2
crn x r
= dr
0gx r + 1
= 0
0gx n
= 0
M
crlr x lr +
L
M
23
この場合は、各階段の先頭に注目して、
場合1と同様に考えることができる。
C =
C =
ìï x l1 « x 1
ïï
ïï x l2 « x 2
ïí
M
ïï
ïï
ïïî x lr « x r
とすれば、ほぼ同様に議論される。
結局、次のように解くことができる。
24
ìï
ïï x 1
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
íï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
ïï
îï
ï
n- r
= d '1 -
å
c pqki
i= 1
x2
= k1
O
M
x l2 - 1
= kl2 - 1
n- r
x l2
= d 'l2 -
å
c 'rs ki
i= 1
O
M
n- r
x lr
L
xn
= d 'lr -
å
c 'uv ki
i= 1
x lr + 1
= kn - lr
M
xn
M
= kn - r
25
任意定数
n -r個
x2 x3
C =
x l1
x l2
一定値
x lr
r個
éx l1 ù
ê
ú
êx 2 ú
ê
ú
êM ú
ê
ú
êx
ú
ê l2 - 1 ú
êx
ú
ê l2 ú
ú
x = êx
ê l2 + 1 ú
ê
ú
êM ú
ê
ú
êx l ú
ê r ú
ê
ú
êM ú
ê
ú
xn ú
ê
ê
ú
ë
û
として、やはり
n - r 個の任意定数を用いて解くことができる。
26
例題1
次の連立方程式を解け。
ìï x - y - 4z
ï
í
ïï 2x - 2y - 8z
î
é1 - 1 - 4ù
ú
A º êê
ú
êë2 - 2 - 8ú
û
= 0
= 0
éx ù
êú
x º êêy ú
ú
êú
êëz ú
û
é0ù
b º êê úú とすると、
êë0úû
Ax = b
拡大係数行列は、次のようになる。
é1 - 1 - 4 0ù
éA b ù= ê
ú
êë
ú
û ê2 - 2 - 8 0ú
êë
ú
û
係数行列と、拡大係数行列の階数(rank)を求める。
27
é1 - 1 - 4ù
(2)- 2´ (1)
ú¾ ¾
A = êê
¾ ¾®
ú
2
2
8
êë
ú
û
なので、
rank A = 1
é1 - 1 - 4ù
ê
ú
ê0 0
0 ú
êë
ú
û
である。
また、
éA
êë
é1 - 1 - 4 0ù
(2)- 2´ (1)
ù
ú¾ ¾
b ú= êê
¾ ¾®
û 2 - 2 - 8 0ú
êë
ú
û
é1 - 1 - 4 0ù
ê
ú
ê0 0
0
0ú
êë
ú
û
rank éêA b ù
= 1
ú
ë
û
\ rankA = rank éêA b ù
ú
ë
û
よって、解が存在する。(不能ではない。)
係数行列の階段行列より、元の連立方程式は、次の
連立方程式と同値。
28
ìï x - y - 4z
= 0
ï
í
ïï 0x - 0y - 0z = 0
î
このとき、3-rankA=2に注意して、
階段の先頭以外を任意定数とする。
ìï y = k1
ï
í
ïï z = k2
î
これより、
x = k1 + 4k2
と表せる。
よって、任意定数 k1, k 2 を用いて、
ìï x = k + 4k
1
2
ïï
ïï
í y = k1
ïï
ïï z = k 2
ïî
(例題1終)
29
例題2
ìï x 1
ïï
ïï 2x 1
ï
í
ïï - x 1
ïï
ïï - 2x 1
î
+ 2x 2
- 2x 3
+ 3x 4
= 2
+ 4x 2 - 3x 3
- 2x 4
= 3
- 2x 2
+ 2x 3
+ x4
= 2
- 4x 2
+ 4x 3 - 6x 4
= - 4
解)
拡大係数行列を階段行列化する。
é1
ê
ê2
ê
ê
ê- 1
ê
ê- 2
êë
2 ù
ú
4
- 3 - 2 3 ú
ú (2)- 2´ (1)
ú¾ ¾ ¾ ¾®
- 2 2
1
2 ú
ú
- 4 4
- 6 - 4ú
ú
û
2
- 2 3
é1
ê
ê0
ê
ê
ê- 1
ê
ê- 2
êë
2
- 2 3
0
1
- 2 2
- 4 4
2 ù
ú
- 8 - 1ú
ú
ú
1
2 ú
ú
- 6 - 4ú
ú
û
30
é1
ê
ê0
ê
(3)+ 1´ (1)
¾¾
¾ ¾® ê
ê0
ê
ê- 2
êë
2
- 2 3
0
1
0
0
- 4 4
é1
ê
ê0
1
´ (3)
ê
4
¾¾
¾® ê
ê0
ê
ê0
êë
é1
ê
ê0
ê
(1)- 3´ (3)
¾¾
¾ ¾® ê
ê0
ê
ê0
êë
2 ù
ú
- 8 - 1ú
ú (4)+ 2´ (1)
ú¾ ¾ ¾ ¾®
4
4 ú
ú
- 6 - 4ú
ú
û
2 - 2 3
0 1
0 0
0 0
2 ù
ú
- 8 - 1ú
ú (2)+ 8´ (3)
ú¾ ¾ ¾ ¾®
1
1 ú
ú
0
0 ú
ú
û
2 - 2 0 - 1ù
ú
0 1
0 7 ú
ú (1)+ 2´ (2)
ú¾ ¾ ¾ ¾®
0 0
1 1 ú
ú
0 0
0 0 ú
ú
û
é1
ê
ê0
ê
ê
ê0
ê
ê0
êë
2 - 2 3
é1
ê
ê0
ê
ê
ê0
ê
ê0
êë
2 - 2 3 2ù
ú
0 1
0 7ú
ú
ú
0 0
1 1ú
ú
0 0
0 0ú
ú
û
é1
ê
ê0
ê
ê
ê0
ê
ê0
êë
0 1
0 0
0 0
2 ù
ú
- 8 - 1ú
ú
ú
4
4 ú
ú
0
0 ú
ú
û
2 0 0 13ù
ú
0 1 0 7 ú
ú
ú
0 0 1 1 ú
ú
0 0 0 0 ú
ú
û
31
この階段行列化より、係数行列と拡大係数行列の階数が等しい
ことがわかり解が存在することがわかる。
また、次のような同値な連立方程式が得られる。
ìï x + 2x
2
ïï 1
ï
í
ïï
ïï
ïî
= 13
x3
= 7
x4
= 1
よって、自由度が1であるので、任意定数を1つ用いて、
x2 = k
と表すことができる。
よって、
ìï x 1 = 13 - 2k
ï
ïï
ïï x 2 = k
í
ïï x 3 = 7
ïï
ïï x 4 = 1
î
(例題2終)
32
練習
次の連立方程式を解け。
(1)
ìï 2x
ï
í
ïï - x
î
- 6y
= 1
+ 3y
= 2
(2)
ìï x
ïï 1 - 3x 2
ï
í 3x 1 - 9x 2
ïï
ïï 2x 1 - 6x 2
ïî
+ x3
+ x4
+ 2x 5
= 3
+ 2x 3
+ 4x 4
+ 3x 5
= 9
+ x3
+ 2x 4
+ 4x 5 = 8
33
同次連立一次方程式
定数項ベクトルが零ベクトルであるような
連立一次方程式、すなわち
Ax = 0
を同次連立一次方程式という。
実は,一般の連立一次方程式:
Ax = b
は、対応する同次連立一次方程式:
Ax = 0
を利用して解くことができる。
34
同次連立一次方程式の自明な解
同次連立一次方程式
Ax = 0
では、零ベクトル 0
を必ず解に持つ。
この解 x = 0 、即ち
éx 1 ù
ê ú
ê Mú=
ê ú
ê ú
êëx n ú
û
é0ù
êú
êú
êMú
êú
ê0ú
êë ú
û
を自明な解という。
35
自由度のある同次方程式の解1
0
ax  by

2ax  2by  0
方向ベクトル
ax  by  0
y
y
解く
(b, a)
O
x
O
解
k
x
( x, y) | ax  by  0
を任意のスカラーとして、
 x   kb 
b


k
 y  ka 
 a 
  

 
と表せる。
36
自由度のある同次方程式の解2
0
x  y  z

2 x  2 y  2 z  0
z
 1
0
 
 1 
y
x yz 0
x
O
 1
1
 
 0 
k1 , k2  R を任意のスカラーとして、
 x    k1  k2 
 1
 1
 y   k
 k  1 k  0 
1
  
 1  2 
 z   k2

 0 
 1 
と表せる。
37
自明な解しかもたない条件
(同次方程式と自明な解)
係数行列 A を m ´ n 行列とし、
未知数ベクトル x を n ´ 1 の列ベクトルとする。
同次連立方程式 A x = 0 自明な解しか持たな
いための必要十分条件は、係数行列の階数が n
であることである。すなわち、
Ax = 0
の解が、 {0 }
rank( A )= n
38
自明な解と正則行列
(同次連立方程式と正則行列)
A を n 次の正方行列とし、 x を n 次元ベクトルとする。
このとき、同次連立一次方程式
Ax = 0
が自明な解しか持たないための必要十分条件は、
A が正則行列であることである。
証明
n = rank( A )
が必要十分条件であるが、
これは A が正則行列であることも意味している。
QED
39
練習
次の同次連立一次方程式が、自明な解以外を持つかどうかを
判定せよ。
ìï x
+ 2y - 2z = 0
ïï
(1) ìï 6x + 4y = 0
(2) ï 2x - y
+ 3z = 0
ï
ï
ï
í
í
ïï 9x + 6y = 0
ïï - x - 3y - z
= 0
î
ïï
4y
+ 3z = 0
ïï x
î
é1
ê
ê5
(3) êê
ê9
ê
ê13
êë
4 ùéx 1 ù
úê ú
6 7 8 ú
úêêx 2 ú
=
úêx ú
10 11 12ú 3 ú
úêêx ú
ú
4ú
14 15 16ú
ê
ú
ûë û
2
3
é0ù
êú
ê0ú
êú
ê0ú
êú
êú
êê0ú
ëú
û
40
非同次連立一次方程式の解1
Ax = 0
Ax = b
y
y
p
O
自明な解
x
x
O
Ap = b
41
非同次連立一次方程式の解2
c
ax  by

2ax  2by  2c
0
ax  by

2ax  2by  0
y
y
解く
é0 ù
ê ú
p º ê cú
ê- ú
êë b úû
(b, a)
p
O
 x
b

k
 y
 a 
 
 
x
O
x
 0 
x
b
b
 
 
   

k

p

k
 y
 a 
 a     c 
 
 
   
 b
42
自由度のある同次方程式の解3
0
x  y  z

2 x  2 y  2 z  0
z
y
 1
0

 1
 1
1
 
 0 



O
2
x  y  z

2 x  2 y  2 z  4
z
y
 1
x yz 0
x



x yz 2
と表せる。
x
O
 1
1
 
 0 
k1 , k2  R を任意のスカラーとして、
 x    k1  k2 
 1
 1
 y   k
 k  1 k  0 
1
  
 1  2 
 z   k2

 0 
 1 
0

 1
(2, 0, 0)
 x   k1  k2  2 
 1
 1  2 
 y   k
  k  1   k  0   0
1
  
 1  2   
 z   k2

 0 
 1   0 
43
練習
次の同じ係数を持つ同次方程式と非同次方程式を解け。
(1)
 y z  0
x

2 x 3 y  z  0
  x 4 y 2 z  0

(2)
 y z  3
x

2 x 3 y  z  3
  x 4 y 2 z  6

44