计算方法线性方程组直接求解.pptVIP

* 3.6 消元形式的追赶法 二、消元形式的追赶法的算法 输入原方程组的阶数n。 输入增广矩阵对应的数组a[],b[],c[],d[]。 for(i=1;i=n-1;i++)，循环1次消去1列。 （消元过程） a[i]/=-b[i-1]; b[i]+=c[i-1]*a[i]; d[i]+=d[i-1]*a[i]; x[n-1]=d[n-1]/b[n-1]; （回代过程） for(i=n-2;i=0;i--) x[i]=(d[i]-c[i]*x[i+1])/b[i]; 输出原方程组的解x[n]。 * 3.6 消元形式的追赶法 消元形式的追赶法的程序（1/2） #include stdio.h #include math.h #define MAXSIZE 50 void input(double a[],double b[],double c[],double d[],long n); void output(double x[],long n); void main(void) {double a[MAXSIZE],b[MAXSIZE],c[MAXSIZE],d[MAXSIZE],x[MAXSIZE]; long n,i; printf(\n请输入原方程组的阶数：); scanf(%ld,n); input(a,b,c,d,n); for(i=1;i=n-1;i++) { a[i]/=-b[i-1]; b[i]+=c[i-1]*a[i]; d[i]+=d[i-1]*a[i];} x[n-1]=d[n-1]/b[n-1]; for(i=n-2;i=0;i--) x[i]=(d[i]-c[i]*x[i+1])/b[i]; output(x,n);} * 3.6 消元形式的追赶法 消元形式的追赶法的程序 （2/2） void input(double a[],double b[],double c[],double d[],long n) { long i; printf(\n请输入原方程组的增广矩阵：\n); printf(\nb1,c1,d1:); scanf(%lf,%lf,%lf,b[0],c[0],d[0]); for(i=2;i=n-1;i++) { printf(\na%ld,b%ld,c%ld,d%ld:,i,i,i,i); scanf(%lf,%lf,%lf,%lf,a[i-1],b[i-1],c[i-1],d[i-1]); } printf(\na%ld,b%ld,d%ld:,n,n,n); scanf(%lf,%lf,%lf,a[n-1],b[n-1],d[n-1]);} void output(double x[],long n) { long i; printf(\n方程组的解为：\n); for(i=1;i=n;i++) printf( %lf,x[i-1]);} * 3.7 LU分解法 三角分解法求解线性方程组的方法，是指把线性方程组的系数矩阵进行三角分解，然后按照线性方程组的矩阵形式进行回代求解的方法。 三角分解法主要有LU分解法（又称直接三角分解法、Doolittle分解法）、LDR分解法、Crout分解法等。 LU分解法把系数矩阵分解为单位下三角阵L和上三角阵U的积； LDR分解法把系数矩阵分解为单位下三角阵L、对角阵D和单位上三角阵R的积，其中DR=U； Crout分解法把系数矩阵分解为下三角阵 和单位上三角阵R的积，其中 =LD。 除此之外，还有求解对称正定线性方程组的LLT分解法（又称平方根法），求解带状对角线性方程组的矩阵形式的追赶法等方法。 * 3.7 LU分解法 对任意矩阵Am×n做初等变换的结果，等于Am×n与对应初等方阵的积。 Em(i（k）,j)= Em(i（k）)= Em(i,j)= 由Em得到的初等方阵 数k乘以Em的j列，加到i列上去，得到 Em(i（k）,j)。Em(i（k）,j)右乘Bn×m的积，等于数k乘以Bn×m的j列，加到i列上去的结果。 数k乘以Em的i行，加到j行上去，得到 Em(i（k）,j)。Em(i（k）,j)左乘Am×n的积，等于数k乘以Am×n的i行，加到j行上去的结果。 数k乘以Em的i列得到Em(i（k）)。Em(i（k）)右乘Bn×m的积，等于数k乘以Bn×m的i列的结果。 数k乘以Em的i行得到Em(i（k）)。Em(i（k）)左乘Am×n的积，等于数k乘以Am×n的i行的结果。 Em交换i列和j列得到