LU分解法数值实验报告MATLAB源码及运行结果.docxVIP

数值实验LU分解法的优点

【实验目的】

编写函数来实现矩阵的直接LU分解法,并利用它来解方程;

总结matlab解线性方程组的常用命令用实例来展示命令的用法。

【实验内容】

给定矩阵A与向量b

求A的LU分解

利用A的LU分解解下列方程：

=1\*GB3①A*x=b=2\*GB3②A2?x=b=3\*GB3③A3?x=b

对第=3\*GB3③题分析一下，如果先求M=A3,再解M*相比有何缺点？

利用A的LU分解法求A?1

提示：的第k列为方程组之解。的准确解为：

【解】：

由Cholesky分解定理：

设A是实对称正定矩阵，则存在唯一的可逆下三角矩阵L，使得。

并且由矩阵乘法规则，可以推出实对称正定矩阵LU分解的简便计算方法(平方根法)：由可知：，进而由(i≥2)可知；

当已求出时，由可知；

当时，由可知，。

【计算机求解】:

流程图：

输入方程阶数n,系数矩阵A,右端向量b

K=1,…,n(分解A=L*U)

u

ukk

=0T

输出失败信息，停

k