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

4.5.6追赶法许多实际问题中，如在船体数学放样中建立的三次样条函数等，往往遇到求解形式的线性方程组，其中系数矩阵称为三对角阵，它的特点是其非零元素集中分布在主对角线及其相邻两条次对角线上，除此之外的其余元素全部为零。写成矩阵元素的形式如下：为此，这里介绍系数矩阵为三对角矩阵的方程组的解法—追赶法。实际上，追赶法是高斯消去法的一种简化形式，同样分为消元和回代两个过程。第62页，共119页，星期日，2025年，2月5日4.5.7实现追赶法的基本步骤(1)消去过程…………其中：第63页，共119页，星期日，2025年，2月5日2)回代过程方程组(4.6)实际上是一组递推关系式，求出就能求出，由可求出，直至求出即(4.7)第64页，共119页，星期日，2025年，2月5日实现追赶法的N—S图第65页，共119页，星期日，2025年，2月5日例4.6用追赶法求解线性方程组，式中第66页，共119页，星期日，2025年，2月5日#defineN4/*N为方程组系数矩阵的阶数*/voidzhuigan(a[],b[],c[],d[],r[],y[])floata[],b[],c[],d[],r[],ty[];{floatq;inti;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;}}第67页，共119页，星期日，2025年，2月5日voidsolve(floatr[N],floaty[N],floatx[N]){inti;x[N-1]=y[N-1];for(i=N-2;i=0;i--)x[i]=y[i]-r[i]*x[i+1];}第68页，共119页，星期日，2025年，2月5日main(){floata[N]={0.0,90.860,-67.590,46.260};/*主对角线左下方元素*/floatb[N]={136.01,98.81,132.01,177.17};/*主对角线元素*/floatc[N]={90.860,-67.590,46.260,0.0};/*主对角线右上方元素*/floatd[N]={-33.254,49.709,28.067,-7.3244};/*右端常数项元素*/floatr[N],y[N];/*存放中间变量数组*/floatx[N];/*存放方程组的解向量的数组*/inti;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]);}第69页，共119页，星期日，2025年，2月5日程序运行结果：x1=-2954.868164x2=4422.830078x3=2492.836426x4=-650.933777第70页，共119页，星期日，2025年，2月5日4.5.8对称高斯法（分解法）许多实际问题中，往往遇到求解系数矩阵为对称的线性方程组，写成矩阵元素的形式如下：为此，这里介绍系数矩阵为对称的方程组的解法—对称高斯法。对称高斯法的基本思想是将分解为的形式，即L是对角元素为1的单位下三角阵，然后再进行分解和回代过程求解。第71页，共119页，星期日，2025年，2月5日4.5.9实现对称高斯法的基本步骤（1）输入方程组的阶数N，系数矩阵A和右端常数矩阵b。（2）分解过程对对第72页，共119页，星期日，2025年，2月5日（3）回代过程对（4）结束。第73页，共119页，星期日，2025年，2月5日对称高斯法的N-S图描述第74页，共119页，星期日，20