Jacobi迭代法C程序方法.docVIP

一、算法理论 迭代格式的引出是依据迭代法的基本思想：构造一个向量系列，使其收敛至某个极限，则就是要求的方程组的准确解。 迭代 将方程组： 在假设，改写成 如果引用系数矩阵 ， 及向量，，， 方程组（1）和（2）分别可写为：及，这样就得到了迭代格式用迭代解方程组时，就可任意取初值带入迭代可知式，然后求。但是，比较大的时候，写方程组和是很麻烦的，如果直接由，能直接得到，就是矩阵与向量的运算了，那么如何得到，呢 实际上，如果引进非奇异对角矩阵 将分解成：要求的解，实质上就有而是非奇异的，所以存在，从而有我们在这里不妨令就得到迭代格式： 现在考虑迭代法的计算程序 分别代表的系数和等号右边的常数项 即 先输入方程，运行函数，如果不为，则执行括号里的，否则执行里面的，最后会调用方法。在中;，之后运行循环，最后输出结果，算法结束。 二、算法框图 三、算法程序 #include math.h #include stdio.h #define max 10 float a[3][3]={{10,-2,-1},{-2,10,-1},{-1,-2,5}}; float b[3]={3,15,10}; float x[3]={0,0,0}; float y[3]