矩阵crout分解的次序与Doolittle分解的次序不同,其他差不多.在矩阵A=LU的crout分解形式中,L是下三角矩阵,U是单位上三角矩阵
当A的各阶主子式不为0,依次计算L的第一列,U的第一行,L的第K列,U的第K行,Doolitlle三角分解法跟这个恰好相反.依据Doolittle分解的手法,crout的算法可以比较容易实现,Crout分解题目: 不知道是分解形式Ax=b的形式还是,分解A=L*U的形式,我这里是根据实验报告书,分解A=L*U的形式,弄成Ax=b的形式修改下代码就行了.
[ 2 10 0 -3 ]
[-3 -4 -12 13]
[ 1 2 3 -4 ]= L(单位下三角矩阵) * U(上三角矩阵)
[ 4 14 9 -13]
下面上代码.
/********************************************
> File Name: Crout.c
> Author:chendiyang
> School:WUST_CST_1501班
> Myblog:www.chendsir.com
> Mail:1441353519