對於可逆矩陣來說,Ax=0的解只有x=0,Ax=b對於任意b都有解。當nullspace裡出現了非零解或者A的列向量空間小於是我們這節要處理的情況。
梯隊形式(echelon form):
因為我們使用消元法獲得U,將之按照相反的過程還原,自然可以得到原來的矩陣A,因此:
L是方陣,行數等於A和U的行數
PA=LU是一般形式:
擷取R: 1.將U的pivots化成0 2.從下往上將pivots正上方的元素化成0,如下:
* matlab裡擷取R R=rref(A),對於可逆矩陣,R=I
上章中是PA=LDU,Gauss-Jordan方法中[A I]->[U L逆]->[I A逆],涉及A->U,pivots上方元素化0等,求R的方法和這兩個方法神似。
R still contains an identity matrix.in the pivot rows and pivot columns.
Rx=0 Ux=0 Ax=0 具有相同的解。
pivot variables: variables(in vector x) correspond to columns with pivots
Free variables: variables(in vector x) correspond to columns without pivots
從下面可以看出,在對x的表示中pivot varialbes全部由free variables表示了。而free variables之所以叫free variables,是因為這些變數是可以任意取值的。
運用R求Ax=0:
Special solutions for Ax=0:就是③中右側的每一列是一個special solutions.
A的nullspace可以看成special solutions的combination的集合
在上例中,A的nullspace是 2維 的向量空間。
注意:仔細觀察,special solutions的個數等於free variables的個數
對於A(m by n,n>m)來說,至少有n-m個free variables.
一個非常重要的結論:
從解Rx=0的過程可知,每多出一個free variables,就會多出一個special solution,對於A(m by n,n>m)來說,至少有n-m個free variables,這也是得出2C的原因。又因為free variables可以任意取值,所以解的個數是無限個。
Solving Ax = b, U x = c, and Rx = d
和Ax=0中不同的一點:即使n>m(more unkonwns than equations),也不一定有解。主要是方程右邊的b!=0造成了方程的平行(b=0的時候方程重合)。
講解這個問題的順序:先找到一個使方程有解的b,主要使用了b應該在列向量空間的方法。使用與上節相同的方法找到方程的解(利用free variables和pivot variables)得出結論。
求Ax=b:
結論:
Every solution to Ax = b is the sum of one particular solution and a solution to Ax = 0:
關於,可以將b看成0+b,按照上節方法b=0的部分的解為,b=b的部分的解為。。(由現象到結論,這是一種推廣)
當b=0時也有特解,只不過特解是0.
將Ax=b化成Rx=d之後就可以直接得出結論,注意這個方法,連同上一節的直接擷取結果的方法:
D中值對應到pivot variable對應位置的特解,特解中free variable對應位置為0.free variable對應到special solution裡。
Rank of the matrix: 矩陣的秩
You see how the rank r is crucial. It counts the pivot rows in the "row space" and the
pivot columns in the column space. There are n − r special solutions in the nullspace.
There are m − r solvability conditions on b or c or d(之所以這麼說,是因為剩下的m-r行決定方程是否有解).
問題集合:
1.
方程無解的時候仍可以有nullspace,此時有Xn但是無Xp.
2.
ans:
3.
也就是說pivot columns所在的列是線性無關的列
4.[??]
小於A好證(AB的結果是A的線性組合,逃不出A所在的‘plane‘),小於B呢?
5.注意若AB=I,在AB不為方陣的時候也可以成立,定義左逆右逆的公式中都是指的方陣,不是方陣的情況下左逆不等於右逆(維數都不相同),逆的定義中是同時滿足左逆右逆才叫逆。
6.
7.[!!]
8. A system Ax=b has at most one particular solution:錯。它的任意解都可以認做是一個特解。只要Xp+Xn=X成立就可以。
9.解的情況大概也就這麼多了
(a)rows>rank(A),columns=rank(A) (b).rows=rank(A),columns>rows (c)rows>rank(A),columns>rank(A) (d) invertible
10.
用上面的結論很容易就可以做出來:(a)秩為2且有解即是 (b)秩小於列數,無解或無限解。
11.
方程組相交成點,方程組中沒有未知量,相交成線,有一個未知量,因此秩是3.
(b)
12. 處於後列的pivot variable是無法使用處於前列的free variable的。
13.why does no 3 by 3 matrix have a nullspace that equals its column space
可逆時:nullspace只有一個元素,column space為
秩為2時:nullspace是線,column space是面
秩為1時:nullspace是面,column space是線