matlab中求解线性方程组的rref函数

时间:2024-02-22 22:10:40
摘自:http://www.maybe520.net/blog/987/


matlab中怎么求解线性方程组呢?
matlab中求解线性方程组可应用克拉默法则(Cramer’s Rule)即通过det()函数计算各个矩阵的行列式来求,也可以用高斯消元法来求解。
matlab中的rref()函数可以将矩阵化成行最简形式,用法如下:

假如有一线性方程组为:
16 x1 + 2 x2 + 3 x3 = 13
5 x1 + 11 x2 + 10 x3 = 8
9 x1 + 7 x2 + 6 x3 = 12
4 x1 + 14 x2 + 15 x3 = 1

令A =
16 2 3 13
5 11 10 8
9 7 6 12
4 14 15 1

然后使用命令 B = rref(A) 则可得到矩阵A的航最简形式
B=
1 0 0 1
0 1 0 3
0 0 1 -3
0 0 0 0

从而可以很方便得到原线性方程组的解了。

如果原线性方程组无解,那么返回的是一个单位矩阵。