学习视频来源:麻省理工公开课_线性代数导论 讲师:Gilbert Strang
http://open.163.com/special/opencourse/daishu.html
Lecture 8 求解Ax=b:可解性和解的结构
已知矩阵 A=⎛⎝⎜⎜1232462682810⎞⎠⎟⎟ A = ( 1 2 2 2 2 4 6 8 3 6 8 10 ) ,解方程组 ⎧⎩⎨⎪⎪x1+2x2+2x3+2x4=b12x1+4x4+6x3+8x4=b23x1+6x2+8x3+10x4=b3 { x 1 + 2 x 2 + 2 x 3 + 2 x 4 = b 1 2 x 1 + 4 x 4 + 6 x 3 + 8 x 4 = b 2 3 x 1 + 6 x 2 + 8 x 3 + 10 x 4 = b 3 。关键在于第三行是前两行之和。根据消元,右侧常数的关联立马可见:如果方程组有解,右侧常数需要满足 b1+b2=b3 b 1 + b 2 = b 3 。
让我们写成增广矩阵形式 (Ab) ( A b ) 并进行消元:
⎛⎝⎜⎜⎜1232462682810b1b2b3⎞⎠⎟⎟⎟ ( 1 2 2 2 2 4 6 8 3 6 8 10 b 1 b 2 b 3 ) -> ⎛⎝⎜⎜⎜100200222244b1b2−2b1b3−3b1⎞⎠⎟⎟⎟ ( 1 2 2 2 0 0 2 4 0 0 2 4 b 1 b 2 − 2 b 1 b 3 − 3 b 1 ) -> ⎛⎝⎜⎜⎜100200220240b1b2−2b1b3−b2−b1⎞⎠⎟⎟⎟ ( 1 2 2 2 0 0 2 4 0 0 0 0 b 1 b 2 − 2 b 1 b 3 − b 2 − b 1 )
得出方程有解的条件为 b3−b1−b2=0 b 3 − b 1 − b 2 = 0 ,和我们之前的主观判断一致。这就是所谓的Solvability 可解性: Ax=b A x = b 有解当且仅当 b b 在中。<=> 如果 A A 各行的线性组合会得到零行,那么右侧向量中分量的同样组合也必须为零。
假设 b=⎛⎝⎜⎜156⎞⎠⎟⎟ b = ( 1 5 6 ) ,代入得: (1020222413) ( 1 2 2 2 0 0 2 4 1 3 ) 。现在我们的目标是求这个方程组的所有解。
step 1:求特定解(教授将 Ax=b A x = b 特解表示为particular solution, Ax=0 A x = 0 基为special solution)。将所有自由变量 x2 x 2 和 x4 x 4 设为 0(自由变量可以任取),然后解出 Ax=b A x = b 的主变量 x1 x 1 和 x3 x 3 。
方程组变为 {x1+2x2=12x3=3 { x 1 + 2 x 2 = 1 2 x 3 = 3 ,解得 xparticular=⎛⎝⎜⎜⎜⎜−203/20⎞⎠⎟⎟⎟⎟ x p a r t i c u l a r = ( − 2 0 3 / 2 0 )
step 2:求零空间内的任何解 xnullspace=c1⎛⎝⎜⎜⎜⎜−2100⎞⎠⎟⎟⎟⎟+c2⎛⎝⎜⎜⎜⎜20−21⎞⎠⎟⎟⎟⎟ x n u l l s p a c e = c 1 ( − 2 1 0 0 ) + c 2 ( 2 0 − 2 1 ) ,也被称为基础解系,解法详见前文。
step 3: xcomplete=xparticular+xnullspace x c o m p l e t e = x p a r t i c u l a r + x n u l l s p a c e 即为 Ax=b A x = b 的所有解。( Axparticular=b A x p a r t i c u l a r = b , Axnullspace=0 A x n u l l s p a c e = 0 , Axparticular+Axnullspace=b A x p a r t i c u l a r + A x n u l l s p a c e = b 。对于方程组某解,其与零空间内任意向量之和仍为解。)图像上应该是一个由子空间从原点平移上来得到的平面。
推广到一般性的 m×n m × n ,秩为 r r 的矩阵,始终有 ,r≤m,r≤n , r ≤ m , r ≤ n 。矩阵的秩决定了方程组解的数目。
- full column rank 列满秩的情况( r=n<m r = n < m ):主变量有 n n 个,自由变量有 0 个。零空间内只有一个零向量 0⃗ 0 → 。而 Ax=b A x = b 的全部解(如果解存在),只有特定解 xparticular x p a r t i c u l a r 。这种情况下 Ax=b A x = b 只有 0 或 1 个解。
- row column rank 行满秩的情况( r=m<n r = m < n ):主变量有 m m 个,自由变量有个。对任意 b b ,都有解。
- row column rank 方阵满秩的情况( r=m=n r = m = n ):可逆,自由变量有 0 个。零空间 N(A) N ( A ) 内只有一个零向量 0⃗ 0 → 。而 Ax=b A x = b 的全部解,只有特定解 xparticular x p a r t i c u l a r 。这种情况下 Ax=b A x = b 只有 1 个解。
- 都不满秩的情况( r<m,r<n r < m , r < n ):这种情况下 Ax=b A x = b 无解或者无穷解。