MATLAB线性方程组求解

MATLAB有多种求解线性方程组的方法，包括高斯消元法、LU分解法、追赶法等。下面我将为你介绍其中两种方法：MATLAB自带的求解函数和LU分解法。

1. MATLAB自带的求解函数

MATLAB自带的求解函数是最为简单的方法，只需要输入系数矩阵A和常数矩阵b，即可求解线性方程组Ax=b。下面是一个简单的示例：

% 定义系数矩阵A和常数矩阵b
A = [1, 2, 3; 4, 5, 6; 7, 8, 9];
b = [1; 2; 3];

% 求解线性方程组Ax=b
x = A\b;

% 输出解向量x
disp(x);

在这个例子中，我们定义了系数矩阵A和常数矩阵b，然后使用MATLAB自带的求解函数A\b来求解线性方程组Ax=b。最后，我们输出解向量x。

2. LU分解法

LU分解法是比较常用的求解线性方程组的方法之一。它将系数矩阵A分解为一个下三角矩阵L和一个上三角矩阵U，即A=LU。然后，将Ax=b转化为Ly=b和Ux=y两个方程组，分别使用前向替换和后向替换求解出y和x。下面是一个简单的示例：

% 定义系数矩阵A和常数矩阵b
A = [1, 2, 3; 4, 5, 6; 7, 8, 9];
b = [1; 2; 3];

% LU分解
[L, U] = lu(A);

% 前向替换求解Ly=b
y = L\b;

% 后向替换求解Ux=y
x = U\y;

% 输出解向量x
disp(x);

在这个例子中，我们首先使用MATLAB自带的lu函数对系数矩阵A进行LU分解，然后分别使用前向替换和后向替换求解Ly=b和Ux=y两个方程组，最后得到解向量x。