雅可比迭代实验报告c语言,雅可比迭代实验报告.doc

雅可比迭代法求解线性方程组的实验报告

一、实验题目

分别利用雅可比迭代法和高斯-塞德尔迭代法求解以下线性方程组：

使得误差不超过 0.00001。

二、实验引言

1.实验目的

= 1 \* GB3 ①掌握用迭代法求解线性方程组的基本思想和步骤，熟悉计算机fortran语言；

= 2 \* GB3 ②了解雅可比迭代法在求解方程组过程中的优缺点。

2.实验意义

雅克比迭代法就是众多迭代法中比较早且较简单的一种，求解方便实用。

三、算法设计

1.雅可比迭代法原理：设有线性方程组Ax=b 满足, 将方程组变形为: x=Bx+f, 则雅可比(Jacobi)迭代法是指，即 由初始解逐步迭代即可得到方程组的解。

算法步骤如下：

步骤1.给定初始值，精度e,最大容许迭代次数M，令k=1。

步骤2.对i=1,2,…,n依次计算

步骤3.求出，若，则输出结果，停止计算。否则执行步骤4.

步骤4.若转步骤2继续迭代。若表明迭代失败，停止计算。

2.算法流程图

四、程序设计

program jacobi

implicit none

integer::i,j

integer::k

save k

real,parameter::e=0.001

integer,parameter::n=3

real::x(n),y(n),b(n)

data b/7.2,8.3,4.2/

real::