数值计算
实验报告
化学实验报告单总流体力学实验报告观察种子结构实验报告观察种子结构实验报告单观察种子的结构实验报告单
数值计算
方法
快递客服问题件处理详细方法山木方法pdf计算方法pdf华与华方法下载八字理论方法下载
实验
实验名称: 解线性方程组的迭代法 班 级:
宿 舍:
成 员:
实验三 解线性方程组的迭代法
一、 实验目的
用雅克比和高斯-赛尔迭代法解线性方程组A,,b,式中A为非奇异矩阵。再给定迭
代初值的情况下,进行迭代,知道满足精度
要求
对教师党员的评价套管和固井爆破片与爆破装置仓库管理基本要求三甲医院都需要复审吗
。
二、 实验方法
(1)雅克比迭代法
、实验程序 1
雅克比迭代法的MATLAB函数文件agui_jacobi.m如下。
function x=agui_jacobi(a,b)
n=length(b);
N=100;
e=1e-4;
x0=zeros(n,1);
x=x0;
x0=x+2*e;
k=0;
d=diag(diag(a));
l=-tril(a,-1);
u=-triu(a,1);
while norm(x0-x,inf)>e&k
> a=[4 -1 0 -1 0 0;-1 4 -1 0 -1 0;0 -1 4 -1 0 -1;-1 0 -1 4 -1 0;0 -1 0 -1
4 -1;0 0 -1 0 -1 4]
>> b=[0;5;-2;5;-2;6]
>> x=agui_jacobi(a,b)
4、实验结果截图
该程序运行后得到一个错误的结果,如下图:
5、实验程序修改
将实验程序中x0=zeros(n,1);修改为x0=b; ,去掉n=length(b);后得下图源程序:
得到结果为下图,仍然为错误结果
(2)高斯-赛尔迭代法
1、实验程序
高斯-赛尔迭代法的MATLAB函数文件agui_ gauss.m如下。 function x=agui_gauss(a,b) n=length(b);
N=100;
e=1e-4;
x0=zeros(n,1);
x=x0;
x0=x+2*e;
k=0;
a1=tril(a);
a2=inv(a1);
while norm(x0-x,inf)>e&k> a=[4 -1 0 -1 0 0;-1 4 -1 0 -1 0;0 -1 4 -1 0 -1;-1 0 -1 4 -1 0;0 -1 0 -1
4 -1;0 0 -1 0 -1 4] >> b=[0;5;-2;5;-2;6] >> x=agui_gauss(a,b)
4、实验结果截图
三、 结果
分析
定性数据统计分析pdf销售业绩分析模板建筑结构震害分析销售进度分析表京东商城竞争战略分析
通过计算机运用高斯-赛尔迭代法求解方程,要经过14次迭代,从而简化求解过程,得到结果。
经过宿舍多次讨论并修改尝试,用雅克比迭代法求解方程,计算机编程计算中出现错误,得不到相应结果,因而无法比较迭代次数。