jacobi迭代法线性方程组求解汇编.pdf

《MATLAB课程设计实践》课程考核

1、编程实现以下科学计算算法，并举一例应用之。（参考书籍《精通MATLAB科学计算》，

王正林等著，电子工业出版社，2009年）

“jacobi迭代法线性方程组求解”

算法说明：

如果系数矩阵的主对角元全不为0，在上节A的分解中取

Q=D

C=D-A

11

其中D是由A的主对角元素组成的对角阵，则有BIDA,rDb，迭代公式为：

11

x(IDA)xDb

k1k

这种迭代方法称为Jacobi迭代法。

在MATLAB中编程实现的Jacobi迭代法函数为：jacobi。

功能：用Jacobi迭代法求线性方程组Ax＝b的解

调用格式：〔x，n〕＝jacobi（A，b，x0，eps，varargin）

其中，A为线性方程组的系数矩阵；

b为线性方程组中常数向量；

x0为迭代初始向量；

eps为解的精度控制（此参数可选）；

varargin为迭代步数控制（此参数可选）；

x为线性方程组的解；

n为求出所需精度的解实际的迭代步数。

流程图：

Function[x,n]=jacobi(A,b,x0,eps,varargin)

N

nargin==3

Ｎ

nargin3

Y

Y

nargin==5

eps=1.0e-6;

error

M=200;

输入矩阵a,bD=diag(diag(A))L=-tril(A,-1)

L,U为上下三角矩阵

U=-triu(A,1)B=D\(L+U)f=D\

x=B*x0+fn=1

N

x-x0=

eps

YN

nM

x0=xx=B*x0+f

Y

输出x.

输出‘Warning:迭