杜立特分解法.docVIP

2013-2014（1）专业课程实践论文 题目：杜立特分解法 一、算法理论 阶线性方程组的系数矩阵非奇异且有分解式，其中为单位下三角矩阵，为上三角矩阵，即，当时，；，当时，，矩阵的这种分解方法为Doolittle的分解。 比较等号两边的第行和第列的元素，得。 因为，所以，， 从而当，，时，，从而， 于是就得到了计算LR分解的一般计算公式。 二、算法框图 三、算法程序 #includeiostream #includestdlib.h #define N 3 using namespace std; int main() { double A[N+1][N+1]={{0,0,0,0},{0,2,1,1},{0,1,3,2},{0,1,2,2}}; double L[N+1][N+1]={0}; double U[N+1][N+1]={0}; double b[N+1]={0,4,6,5}; double y[N+1]; double x[N+1]; int i,j,k,p; for(j=1;j=N;j++) U[1][j]=A[1][j]; L[1][1]=1; for(i=2;i=N;i++) {L[i][1]=A[i][1]/U[1][1]; L