实验5 LU分解法求解线性方程组.docxVIP

实验5LU分解法求解线性方程组

实验目的

了解LU分解的实现，并从矩阵分解的角度理解求解方程组的高斯消元过程。

实验内容

系数矩阵A的LU分解可用来求矩阵A的逆A-1。如果矩阵A非奇异，则A-1存在，且AA-1I。设X1,X2,…,Xn是A-1的列，而E1,E2,…,En是单位矩阵I的列。

方程AA-1I可表示为

则上式等价于n个线性方程组

这样，求A-1等价于求解n一个具有相同系数矩阵的线性方程组。使用LU分解法求解，只需完成一次三角分解，然后求解。n个下三角方程组和n一个上三角方程组，这种方法比用高斯消去法求解。n个线性方程组节省很多计算量，效率更高。

编写矩阵的带列选主元过程的LU分解程序。

实验要求

使用MATLAB的向量运算来实现算法。

实验习题

1）编写一个使用带列主元LU分解法求解A-1的MATLAB程序；利用MATLAB内置的lu函数，再编写前代法和后代法的MATLAB程序求解A。-1。

2）使用第1题中的两种方法，并结合MATLAB的反斜杠运算符来求解矩阵的逆。

3）（1）分别用以下3种方法求10阶希尔伯特（）矩阵的逆：第1）题中的两种方法和MATLAB内置的专门求希尔伯特矩阵的逆的函数invhilb；

（2）用MATLAB内置的cond函数计算10阶希尔伯特矩阵的条件数。讨论当A病态时，第1）题求矩阵的逆的方法可靠吗？通过计算AA-1验证方法的可靠性。

5、实验思考

LU分解法求线性方程组相比高斯消去法有哪些优势？