线方程组简单迭代法.ppt

第三章 求解线性方程组的迭代方法;引言 ;两个简单的例子; 例1;（1.3）; 将这些值代入(1.3) 式右边 (若(1.3)式为等式即求得方程组的解，但一般不满足).;（1.4）; 从此例看出，由迭代法产生的向量序列 逐步逼近;迭代法的基本思想是构造一个向量序列{X(k)}，使其收敛到某个极限向量 X*，而X*就是 AX = b 的准确解。 问题：如何构造迭代序列？ 迭代序列在什么情况下收敛？ ;简单迭代法的迭代格式; 设 ,并将 写为三部分 ;易知，雅各布（Jacobi）迭代有;; 的雅可比(Jacobi)迭代公式如下: ;于是，解 的雅可比迭代法的分量计算公式为 ;;; 由可知计算过程可知，雅可比迭代法计算 公式简单，每迭代一次只需计算一次矩阵和 向量的乘法且计算过程中原始矩阵A始终不变. ;例1 用J法求解线性方程组;取初始向量x(0)=(0,0,0)T,迭代可得;k;简单迭代法的算法如下：; 对于任何由 变形得到的等价方程组 ，;一般迭代法收敛性的基本定理; 设有解 的迭代法 ;因此，研究迭代法(2.2)收敛性问题就是要研究迭代矩阵 满足什么条件时，有 ; 定理1; 推论;计算机实现程序;#includestdio.h #includemath.h #define eps 1e-3 #define max 100 void Jacobi(float *a,int n,float x[]) { int i,j,k=0; double epsilon,s; double *y=new double[n]; for(i=0;in;i++)x[i]=0; while(1) { epsilon=0; k++; for(i=0;in;i++) { s=0; for(j=0;jn;j++) { if(j==i)continue; s+=*(a+i*(n+1)+j)*x[j]; } y[i]=(*(a+i*(n+1)+n)-s)/(*(a+i*(n+1)+i)); epsilon += fabs(y[i]-x[i]); } for(i=0;in;i++)x[i]=y[i]; if(epsiloneps) {printf(迭代次数为:%d\n,k);return;} if(k=max) {printf(迭代发散);return;} }delete y;};void main() { int i; float a[4][5]={10,-1,2,0,-11,0,8,-1,3,-11,2,-1,10,0,6,-1,3,-1,11,25}; /*float a[9][10]={31,-13,0,0,0,-10,0,0,0,-15, -13,35,-9,0,-11,0,0,0,0,27, 0,-9,31,-10,0,0,0,0,0,-23, 0,0,-10,79,-30,0,0,0,-9,0, 0,0,0,-30,57,-7,0,-5,0,-20, 0,0,0,0,7,47,-30,0,0,12, 0,0,0,0,0,-30,41,0,0,-7, 0,0,0,0,-5,0,0,27,-2,7, 0,0,0,0,0,0,0,-2,29,-10};*/ float x[4]; //float x[9]; Jacobi(a[0],4,x); //Jacobi(a[0],9,x); for(i=0;i4/*9*/;i++)； printf(x[%d]=%f\n,i,x[i]); };作业题:P104 1