线性代数作业2 Matlab计算计算矩阵最简形和计算一元一次矩阵方程（x = ax + b）

线性代数作业2

1、分别用笔算和MATLAB函数计算：利用初等行变换化矩阵A为行最简形，并求矩阵A的秩，其中A=

（1）笔算：(请用公式编辑器编辑数学公式呈现过程)

这个需要用到WPSword文档的公式编辑器（微软的word用不了）

（2）MATLAB代码：

>> a = [1 -1 3 -4 3; 3 -3 5 -4 1; 2 -2 3 -2 0; 3 -3 4 -2 -1]
a =

1    -1     3    -4     3
3    -3     5    -4     1
2    -2     3    -2     0
3    -3     4    -2    -1

>>rref(a)

ans =

1    -1     0     2    -3
0     0     1    -2     2
0     0     0     0     0
0     0     0     0     0
>> rank(a)

ans = 

	2

Matlab输入矩阵a = [1 -1 3 -4 3; 3 -3 5 -4 1; 2 -2 3 -2 0; 3 -3 4 -2 -1]他就会生成矩阵，然后再套用公式rref(a)就可以计算矩阵a的最简形。

2．分别用笔算和MATLAB函数计算：
设且X = AX + B,求矩阵X

（1）笔算：(请用公式编辑器编辑数学公式呈现过程)

（2）MATLAB代码：

E = [1 0 0;0 1 0;0 0 1]

>>E =

1     0     0
0     1     0
0     0     1

>>B = [1 -1; 2 0; 5 3]

B =

 1    -1
 2     0
 5     3

>> A = [0 1 0; -1 1 1; -1 0 -1]

A =

0     1     0
-1     1     1
-1     0    -1
>>C = E - A

C =

1    -1     0
1     0    -1
1     0     2

>> X = C\B

X =

3     1
2     2
1     1

在Matlab中是无法直接计算X = AX + B的，但是可以计算ax =b，我们先把X = AX + B变成 （E-A）X = B在这里E-A用C来代替，然后通过B左除C(左除是)，X = C\B,就可以计算出X = AX + B。

左除与右除的概念
Matlab提供了两种除法运算：左除（\）和右除（/）。一般情况下，x=a\b是方程ax =b的解，而x=b/a是方程xa=b的解。