数值分析-迭代法.docxVIP

华北科技学院上机报告 系（部） 专业、班级 姓名 学号 课程名称 数值分析 上机题目 实验六，实验七 任课教师 指导教师 成绩（优、良、中、及格、不及格） 华北科技学院基础部 实验六 解线性方程组的迭代法 目的与要求： 熟悉求解线性方程组的有关理论哈方法。 会编制雅可比迭代和高斯 —塞得尔迭代法。 通过实际计算，进一步了解各算法的优缺点，选择合适的数值方法。 2．雅可比迭代法 算法 设方程组 AX=b的系数矩阵的对角元素 aii 0(i 1,2,L ,n),M 为迭代次数容许的最大值 为容许误差 . ① 取初始向量 x (x1(0),x(20),L ,x(n0))T,令 k=0; ② 对 i 1,2,L ,n 计算 xi(k 1)1 (b xi (k 1) 1 (bi aii aij x(jk)); j1 n ③ 如果 xi(k 1) xi(k) ,则输出 x(k 1) ,结束;否则执行④ i1 ④ 如果 k M ,则不收敛 ,终止程序 ;否则 k k 1,转② . 1.分别用雅可比迭代法与高斯 1.分别用雅可比迭代法与高斯 -塞德尔迭代法解下列方程组 2)RI V,其中 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 R 0 0 0 30 57 7 0 5 0 ,V 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 1.用雅可比迭代法计算： #include #include #include #include #include #include #include #define EPS 1e-8 #define MAX 100 float *Jacobi(float a[9][10],int n) { float *x,*y,s; double epsilon; int i,j,k=1; x=(float *)malloc(n*sizeof(float)); y=(float *)malloc(n*sizeof(float)); for(i=0;in;i++) x[i]=0; while(1) { k++; epsilon=0; for(i=0;in;i++) { s=0; for(j=0;jn;j++) { if(j==i) continue; s=s+a[i][j]*x[j]; } y[i]=(a[i][n]-s)/a[i][i]; epsilon=epsilon+fabs(y[i]-x[i]); } 高斯 -塞德尔迭代法： #include #include #include #include #include iostream #include #define N 100 void main() { } } } } int i; float *x; float c[90]={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 *GauseSeide(float *,int); x=GauseSeide(c,9); for(i=0;i9;i++) printf(x[%d]=%f\n,i,x[i]); } float *GauseSeide(float *a,int n) { int i,j,nu=0; float *x,dx,d,wucha; x=(float *)malloc(n*sizeof(float)); for(i=0;i=n-1;i++) x[i]=; while(fabs(wucha)1e-8) { for(i=0;i=n-1;i++) { d=; for(j=0;j=n-1;j++) d+=*(a+i*(n+1)+j)*x[j]; d=d-*(a+i*(n+1)+i)*x[i]; dx=(*(a+i*(n+1)+n)-d)/(*(a+i*(n+1)+i)); wc=x[i]-dx; x[i]=dx; } if(nu=N) { printf( 迭代发散 \n); exit(1); nu++; } r