此题目的直接三角分解法,与上面的Guass列选主元消去法有着一些联系,至少前面的输入代码是差不多的:题目:
2、Doolittle三角分解法
[ 2 10 0 -3 ]
[-3 -4 -12 13]
[ 1 2 3 -4 ]
[ 4 14 9 -13]
有了前面的基础,这个题目就比较简单了,根据实验书给的算法,一步一步解决就OK了:
1.输入矩阵的维数,接着输入矩阵A,再进行一行一列的计算
2.计算u的第一行,再计算L的第一列....后面的就不多bb了,根据指导书上用一个循环就可以实现u的所有行的值,以及L的所有列的值.talk is cheap,let’s show my code. 下面上自己的代码(唯一的问题就是自己不会GUI图形编程,只怪当初懒,半途而废了..)
Doolittle三角分解法代码:
/********************************************
> File Name: Doolittle.c
> Author:chendiyang
> School:WUST_CST_1501班
> Myblog:www.chendsir.com
> Mail:1441353519@qq.com
> Created Time: 2017年04月18日 星期二 9时15分56秒
************************************************************************/
#include <stdio.h>
#include <math.h>
#define MAX 20
int main()
{