解线性方程组（二）

实验类型：●验证性实验  ○综合性实验  ○设计性实验

实验目的：进一步熟练掌握用Jacobi迭代法和Gauss-Seidel法解线性方程组的算法，提高编程能力和解算线性方程组问题的实践技能。

实验内容：

1)取初值性x(0)=(0,0,0,0)T, 精度要求ε=10-6，用Jacobi迭代法解线性方程组

2)取初值性x(0)=(0,0,0,0)T, 精度要求ε=10-6，用Gaoss-Seidel迭代法解线性方程组

实验原理：Jacobi迭代算法，Gaoss-Seidel迭代算法。

实验任务及说明：

1. 从一般迭代法收敛性理论上分析并讨论1）、2）两个问题所选用的迭代方法是否收敛，若分析过程涉及到四次代数方程求根可用Matlab的roots函数求多项式的根，请将分析、讨论过程及收敛性结论写在规格为A4的打印纸上。   
2. 分别用Jacobi迭代算法，Gaoss-Seidel迭代算法手工解算1）、2）两个问题的线性方程组，清楚的写出迭代公式，要求严格按照算法步骤手工计算，其中间计算结果保留五位小数，手工迭代至少迭代三步，计算出X(1),X(2),X(3)，请将手工计算过程写在规格为A4的打印纸上，对计算过程进行总结和梳理后在A4纸上分别手工画出Jacobi迭代算法，Gaoss-Seidel迭代算法的程序设计流程图（提示：算法开始或结束--用椭圆框；类似赋值语句（算法公式）的执行处理--用矩形框，赋值语句（算法公式）写在矩形框内；条件判断---用菱形框，条件判断表达式写在菱形框里，条件判断必须有两个出口即两个流出的流程线，在流出的流程线傍边要清楚的标出Yes或No）。
3. 用Matlab编写程序上机计算并输出中间结果（可参考马昌凤柯艺芬编著的数值计算方法(MATLAB版) ，科学出版社2020.11.  Page49，Page51 ）；请将算法流程图（使用Word画图）、程序代码、调试截图记录、输入数据截图记录、输出的中间结果截图记录及最终计算的结果截图记录等作为实验报告内容编排成word文档，然在规格为A4的打印纸上双面打印，实验报告中的word文档的文字采用宋体5号，单倍行距，报告中截图上的文字不得小于宋体六号，截图文字的字迹必须清晰可见！。
4. 请在实验报告中验证计算机程序迭代计算的前三步近似解向量X(1),X(2),X(3)与手工迭代计算的前三步近似解向量X(1),X(2),X(3)是否一致，要求将手工计算的前三步结果和计算机程序计算的前三步结果列表对比，经过对比分析在列表后的文档中清楚给出是否一致的验证结论。
5. 纸质版实验报告装订顺序：①实验任务（本文档版式不得做任何改动）；②实验报告正文（Word文档）；③写在A4纸上的手工分析与计算（也可以将A4纸上详细手工分析与计算过程拍照成JPG图片插入在实验报告的word 电子版中随Word文档双面打印）。

实验步骤

1. 要求上机实验前先进行算法收敛性分析并手工计算、对手工计算进行总结后画出流程图，再编写出程序代码；
2. 编辑录入程序；
3. 记录运行时的输入（截图）和输出（截图）。调试程序并记录调试过程中出现的问题（截图）及修改程序的过程（截图）；
4. 经反复调试后，运行程序并验证程序运行是否正确，必须明确写出分析对比过程及验证结论。

  

实验报告：根据实验情况和结果撰写并递交实验报告（含手工分析与计算）。

实验总结：(学会了......; 掌握了......; 训练了......; 发现了......; 今后学习中......有待提高。）

程序代码

电子报告word文件命名规则：专业班级-完整学号-实验X-姓名.doc, 如信息123班学号为201212030315的郭海涛同学实验2报告word文件命名则应是：信息123-201212030315-实验2-郭海涛.doc, 其中 .doc是Word文件扩展名。

特别提醒：学委负责检查电子文件命名规范，命名不规的将不接收。每人必需提交电子报告和纸质报告各一份(两者内容一致，纸质版实验报告要求用A4纸张双面打印，按规定的装订顺序装订)。本实验任务书版式不得做任何改动。电子报告以班为单位提交压缩包.

问题一：

手写算法迭代过程及判断是否收敛过程：

流程图：

算法流程图：

程序设计流程图：

程序代码：

function [y,n]=Jacobi(A,b,x0,ep)

A=[14 4 4 4;4 14 4 4;4 4 14 4;4 4 4 14];

b=[-4;16;36;56];

x0=[0;0;0;0];

% epslion为误差值

epslion=1.0000e-06;

% N为自己设置的迭代次数

N=100;

L=tril(A,-1); %L为A的单位下三角矩阵

D=diag(diag(A)); %D为A的对角矩阵

U=triu(A,1); %U为A的上三角矩阵

B=-D\(L+U); %B为迭代矩阵

g=D\b;error=1;n=0;

for i=1:N

while error>epslion && n<N

y=x0;x0=B*x0+g;error=norm(y-x0,inf);

if error>epslion

n=n+1;

fprintf('第%d次jacobi迭代结果为：',n)

disp(x0);

fprintf('第%d次迭代结果的误差为：\n',n);disp(error)

break

end

if error<=epslion

n=n+1;fprintf('第%d次jacobi迭代结果为：',n)

disp(x0);fprintf('第%d次迭代结果的误差为：\n',n);disp(error);

fprintf('该方程组的jacobi迭代法的最小迭代次数为：')

disp(n);break

end

end

end

end

运行时的输入和输出：

手工计算结果与程序运行结果对比：

变量	X1(k)		X2(k)		X3(k)		X4(k)
迭代次数	手算	Matlab	手算	Matlab	手算	Matlab	手算	Matlab
0	0	0	0	0	0	0	0	0
1	-0.28571	-0.2857	1.14286	1.1429	2.57143	2.5714	4.00000	4.0000
2	-2.48980	-2.4898	-0.65306	-0.6531	1.18367	1.1837	3.02041	3.0204
3	-1.30029	-1.3003	0.65306	0.6531	2.60641	2.6064	4.55977	4.5598

手工计算过程中保留五位小数，matlab程序中保留的是四位小数，将我手工计算的结果保留四位小数和matlab的结果比对后完全一致。

问题2：

手写算法迭代过程及判断是否收敛过程
省略

程序代码：

function x=Gauss_Seidel_fun(A,b,n,x0,tol,N)

x=zeros(n,1);  % 给x赋值

k=-1;

while k<N

for i=1:n

if i==1

x(1)=(b(1)-A(1,2:n)*x0(2:n))/A(1,1);

elseif i==n

x(n)=(b(n)-A(n,1:n-1)*x(1:n-1))/A(n,n);

else

x(i)=(b(i)-A(i,1:i-1)*x(1:i-1)-A(i,i+1:n)*x0(i+1:n))/A(i,i);

end

end

if norm(x-x0)<tol

break;

end

x0=x;k=k+1;

disp(['when k=',num2str(k)]);

disp('x=');disp(x);                       %输出中间结果

end

if k==N

disp('迭代次数已到达上限!');

end

disp(['迭代次数 k=',num2str(k)])

end

运行时的输入和输出：

手工计算结果与程序运行结果对比：

变量	X1(k)		X2(k)		X3(k)		X4(k)
迭代次数	手算	Matlab	手算	Matlab	手算	Matlab	手算	Matlab
0	0	0	0	0	0	0	0	0
1	-2.66667	-2.6667	0.05556	0.0556	0.89815	0.8981	4.35725	4.3573
2	-1.77225	-1.7722	0.64258	0.6426	2.59747	2.5975	4.62232	4.6223
3	-1.35627	-1.2492	0.87560	0.8935	2.79399	2.8725	4.69278	4.7248

手工计算过程中保留五位小数，matlab程序中保留的是四位小数，将我手工计算的结果保留四位小数和matlab的结果比对后在k=0和k=1时完全一致，在k=3时有出入。

实验总结：

学会了另外两种迭代方法：Jacobi迭代法和Gauss-Seidel法求解线性方程组，以及如何判断算法收敛; 掌握了这两种迭代方法的求解迭代过程以及word画流程图，还有这两种算法在面对同一个线性方程组的优势; 训练了function函数的使用以及在纸上解题的能力和独立思考以及自学能力，还有独立编程并实现这两种算法的能力; 今后学习中对于算法框图的构建以及程序设计框图的展现还有待提高。