不选主元的三角分解法.docVIP

2012-2013（1）专业课程实践论文 不选主元的三角分解法 葛禹泽，0818180109，R数学08-1班 徐金龙，0818180215，R数学08-2班 一、算法理论 设为非奇异矩阵，且有分解式，其中为单位下三角阵，为上三角阵， 即 =…………………………(1) 下面说明的元素可以由步直接计算定出，其中第步定出的第行和的第行元素。由式（1），有 于是得的第一行元素； 于是得得第一行元素。 总结上述讨论，得到第行，用直接三角分解法解（要求所有顺序主子式都不为零）的计算公式: 二、算法框图 三、算法程序 #includestdio.h #includestdlib.h #includemath.h #define N 20 double A[N][N], x[N], y[N], b[N], L[N][N], U[N][N], f1=0, f2=0; int n; int i,j,k; void main() { printf(请输入方程的维数n!\n n = ); scanf(%d, n); getchar(); printf(\n输入%d行%d列矩阵\n, n, n); for (int i=0; in; i++) { for (int j=0; jn; j++) scanf(%lf, A[i][j]); getchar(); } printf(\n输入线性方程组右端项b