利用矩陣除法求線性方程組的特解(或一個(gè)解) 方程:AX=b 解法:X=A\b 例1-76 求方程組 的解。 解: >>A=[5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5 6 0 0 0 1 5]; B=[1 0 0 0 1]'; R_A=rank(A) %求秩 X=A\B %求解 運(yùn)行后結(jié)果如下 R_A = 5 X = 2.2662 -1.7218 1.0571 -0.5940 0.3188 這就是方程組的解。 或用函數(shù)rref求解: >> C=[A,B] %由系數(shù)矩陣和常數(shù)列構(gòu)成增廣矩陣C >> R=rref(C) %將C化成行最簡行 R = 1.0000 0 0 0 0 2.2662 0 1.0000 0 0 0 -1.7218 0 0 1.0000 0 0 1.0571 0 0 0 1.0000 0 -0.5940 0 0 0 0 1.0000 0.3188 則R的最后一列元素就是所求之解。 |
|