#
Jacobi点迭代法该线性方程组的增广矩阵为:=include<> include<> include<> include<> define n 3define precision 1e-16static double aa[n][n1]={{2255}{3476}{1335}} void main() {int ijdetdouble a[n1][n2]x[n1] int GaussElimin
解(1):采用Jacobi迭代法时Matlab计算程序为:clearclci=1a=[5 2 1-1 4 22 -3 10]d=diag(diag(a))l=d-tril(a)u=d-triu(a)d0=inv(d)b=[-12203]x0=[111]B=d0(lu)f=d0bx=Bx0fwhile norm(x-x0inf)>=1e-4 x0=x x=Bx0f
:
: :问题:用Jacobi迭代法解线性方程 解线性方程有两种方法:直接法和间接法迭代法属于后者考虑方程组Ax=b其中A为非奇异矩阵当A为低阶稠密矩阵时想去发是解方程组的有效方法但是对于工程技术中产生的大型稀疏矩阵方程组(A阶数n很大但零元素较多)利用迭代法是合适的在计算机内存和运算两方面迭代法都有很大优势具体算法对方程组Ax=b其中A为非奇异矩阵设并将A写为三部分=D-L
2 迭代法 迭代法的一般概念迭代法是数值计算中一类典型方法不仅用于方程求根而且用于方程组求解矩阵求特征值等方面迭代法的基本思想是一种逐次逼近的方法首先取一个精糙的近似值然后用同一个递推公式反复校正这个初值直到满足预先给定的精度要求为止 对于迭代法一般需要讨论的基本问题是:迭代法的构造迭代序列的收敛性天收敛速度以及误差估计这里主要看看解方程迭代式的构造 对方程()在区间内可改写
#
#
#
1. Jacobi迭代法例1 用jacobi迭代法求解代数线性代数方程组 保留四位有效数字(err1e-4)其中A=[8 -1 12 10 -11 1 -5]b=[1 4 3]解:编写jacobi迭代法的函数文件保存为 [xk]=jacobi(Abx0epsN) 求解Ax=bx0为初始列向量eps为误差容限N为最大迭代次数 输出x为近似解k为迭代次数n=length(A)x=zeros(n1
违法有害信息,请在下方选择原因提交举报