线性方程组的简单迭代法.pptVIP

3.1 简单迭代法 引言 两个简单的例子 例1 已知 迭代法的基本思想是构造一个向量序列{X(k)}，使其收敛到某个极限向量 X*，而X*就是 AX = b 的准确解。 问题：如何构造迭代序列？ 迭代序列在什么情况下收敛？ 简单迭代法的迭代格式 迭代矩阵 易知，雅各布（Jacobi）迭代有 简单迭代法的算法如下： 迭代法的收敛性 一般迭代法收敛性的基本定理 计算机实现程序 用雅各比迭代法下面线性方程组 #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