实验要求:
分别用雅可比算法和高斯-塞德尔算法求解给定的线性方程组。
函数接口定义:
int Jacobi( int n, double a[][MAX_SIZE], double b[], double x[], double TOL, int MAXN );
int Gauss_Seidel( int n, double a[][MAX_SIZE], double b[], double x[], double TOL, int MAXN );
两函数的接口定义是一致的:n为矩阵a的维数,MAX_SIZE是由裁判程序定义的矩阵最大维数,b是方程组中的常向量,x传入叠代初始向量,求得的解将存储在x中返回;TOL是需要达到的精度上限,MAXN是叠代最大次数。函数的返回值为整数,有以下4种可能:
这里必须注意到,如我们在第五章中讨论过的,在两种算法的公式中,A的对角元aii都出现在分母上,所以为了算法的稳定,我们必须采取措施使得对角元的绝对值尽可能大。因为叠代中A不改变,所以我们可以事先调整好A的对角元。
不同的调整策略可能导致不同的收敛性,本书裁判程序要求的调整策略是:对每
个aii,首先向下扫描该列元素(包括aii)中绝对值最大的那个元,若该元不为0,则通过行交换将该元换到aii的位置。若该元为0,则向上扫描该列元素中绝对值最大的那个元,若该元不为0,则将该元所在行的全部元加到第i行。
实验算法:
雅可比迭代:
因篇幅问题不能全部显示,请点此查看更多更全内容
Copyright © 2019- ovod.cn 版权所有 湘ICP备2023023988号-4
违法及侵权请联系:TEL:199 1889 7713 E-MAIL:2724546146@qq.com
本站由北京市万商天勤律师事务所王兴未律师提供法律服务