第4章线性方程组的数值解法.ppt

2) 回代过程 方程组(4.6)实际上是一组递推关系式，求出 就能求出 ，由 可求出 ，直至求出 即 (4.7) 实现追赶法的N—S图 例 4.6 用追赶法求解线性方程组，式中 #define N 4 /* N为方程组系数矩阵的阶数 */ void zhuigan(a[], b[], c[], d[], r[], y[]) float a[], b[], c[], d[], r[],t y[]; { float q; int i; r[0]=c[0]/b[0]; y[0]=d[0]/b[0]; for(i=1;iN;i++) { q=b[i]-r[i-1]*a[i]; r[i]=c[i]/q; y[i]=(d[i]-y[i-1]*a[i])/q; } } void solve(float r[N],float y[N],float x[N]) { int i; x[N-1]=y[N-1]; for(i=N-2;i=0;i--) x[i]=y[i]-r[i]*x[i+1]; } main() { float a[N]={0.0,90.860, -67.590,46.260};/*主对角线左下方元素*/ float b[N]={136.01,98.81,132.01,177.17};/*主对角线元素*/ float c[N]={90.860,-67.590,46.260,0.0};/*主对角线右上方元素 */ float d[N]={-33.254,49.709,28.067,-7.3244};/*右端常数项元素*/ float r[N],y[N]; /* 存放中间变量数组 */ float x[N]; /* 存放方程组的解向量的数组 */ int i; clrscr(); zhuigan(a,b,c,d,r,y); solve(r,y,x); for(i=0;iN;i++) printf(x%d=%13.6f\n,i+1,x[i]); } 程序运行结果： x1= -2954.868164 x2= 4422.830078 x3= 2492.836426 x4= -650.933777 4.5.8 对称高斯法（分解法） 许多实际问题中，往往遇到求解系数矩阵 为对称的线性方程组 ，写成矩阵元素的形式如下： 为此，这里介绍系数矩阵为对称的方程组的解法—对称高斯法。对称高斯法的基本思想是将分解为的形式，即 L 是对角元素为1的单位下三角阵，然后再进行分解和回代过程求解。 4.5.9 实现对称高斯法的基本步骤 （1）输入方程组的阶数N，系数矩阵A和右端常数矩阵b。 （2）分解过程 对 对 （3）回代过程 对 （4）结束。 对称高斯法的N-S图描述 例 4.7 用对称高斯法求解线性方程组 ，式中 #include stdio.h #define N 5 /* N为方程组系数矩阵的阶数 */ void sgaus(float a[N][N],float b[N]) { float q; int i,j,k; for(k=0;kN-1;k++) { for(i=k+1;iN;i++) { q=a[i][k]/a[k][k]; for(j=k+1;j=i;j++) a[i][j]-=q*a[j][k]; b[i]-=q*b[k]; } } } void solve(float a[N][N],float b[N],float x[N]) { float s; int j,k; x[N-1]=b[N-1]/a[N-1][N-1]; for(k=N-2;k=0;k--) { s=b[k]; for(j=k+1;jN;j++) s-=a[j][k]*x[j]; x[k]=s/a[k][k]; } } void zg_matric(float a[N][N