数值作业：Crout克洛特分解法之C语言代码

最新推荐文章于 2025-08-07 09:14:29 发布

原创

最新推荐文章于 2025-08-07 09:14:29 发布 · 5.5k 阅读

11 ·

CC 4.0 BY-SA版权

文章标签：

#c语言 #Crout分解法 #克洛特分解法 #算法

本文介绍了Crout分解法的矩阵分解过程，与Doolittle分解法的区别在于次序。在Crout分解中，A被分解为L（下三角矩阵）和U（单位上三角矩阵）。提供了C语言代码实现，并给出了一个示例矩阵的分解。运行代码得到的L和U矩阵可供验证正确性。

矩阵crout分解的次序与Doolittle分解的次序不同，其他差不多．在矩阵A=LU的crout分解形式中，L是下三角矩阵，U是单位上三角矩阵
这里写图片描述
　　　当A的各阶主子式不为0,依次计算L的第一列，U的第一行，L的第K列，Ｕ的第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 ]＝　Ｌ(单位下三角矩阵）　＊　　　Ｕ（上三角矩阵）
[ 4 14 9 -13]
下面上代码．
　　　

/********************************************
    > File Name: Crout.c
    > Author:chendiyang
    > School:WUST_CST_1501班
    > Myblog:www.chendsir.com
    > Mail:1441353519@qq.com 
    > Created Time: 2017年