- 1、本文档共8页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
线性方程组的J-迭代,GS-迭代,SOR-迭代,SSOR-迭代方法7
西京学院数学软件实验任务书
课程名称 数学软件实验 班级 数0901 学号 0912020119 姓名 王震 实验课题
雅克比迭代、高斯—赛德尔迭代、超松弛迭代
实验目的
熟悉雅克比迭代、高斯—赛德尔迭代、超松弛迭代
实验要求
运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成 实验内容
雅克比迭代法
高斯—赛德尔迭代法、
超松弛迭代法
成绩 教师
【实验课题】
雅克比迭代、高斯—赛德尔迭代、超松弛迭代
【实验目的】
学习和掌握线性代数方程组的雅克比迭代、高斯—赛德尔迭代、超松弛迭代法,并且能够熟练运用这些迭代法对线性方程组进行求解。
【实验内容】
1、问题重述:
对于线性方程组,即:
(1),
其中,
如何运用雅克比迭代、高斯—赛德尔迭代、超松弛迭代法对线性方程组进行求解。
2、方法原理:
2.1雅克比迭代
迭代思想:首先通过构造形如的等式,然后给定一个初值,再通过进行迭代。
step1 :对(1)相应第行中的用其它元素表示为:
即:
Step 2 :进行迭代
,,
取它的判断条件为小于一个确定的误差值,跳出循环。其中,为循环次数,则为所求的近视解。(参考程序)
2.2、高斯—赛德尔迭代
迭代思想:以雅克比迭代法为基础,假定每一次迭代得到的数据都优越于,即中的每一个分量都优越于中对应的分量。由于在求解第个分量的过程中,前个分量已经求出,所以可以直接调用所求出来的前个分量。
考虑雅克比迭代: ,
即:
所以,
当可逆时,由于在求解第个分量的过程中,前个分量已经求出,将上式改写为:
即:
所以,
设,,上式改写为:
,
迭代算法:
得出计算步骤:
,
由于再代入数据时尽可能采用最新数值,因此,高斯—赛德尔迭代法的收敛速度大于雅克比迭代法。
2.3、逐次超松弛迭代法(SOR)
逐次超松弛迭代法是对高斯—赛德尔迭代的一种改进,在高斯—赛德尔迭代的基础上加入松弛因子,使得与的误差由控制。
由高斯—赛德尔迭代公式:
得出SOR迭代算法:
,其中,为松弛因子。迭代的控制条件为:
,
得出SOR的计算公式为:
,
则存在,
令,则,
2.4、对称超松弛迭代法
对称超松弛迭代法是对逐次超松弛迭代法的改进,在逐次超松弛迭代法的基础上,首先对线性方程组按照顺序方式依次求解,然后在此基础上对线性方程组逆序求解,得出SSOR迭代法。
求解步骤:
Step 1 顺序求解得出:
Step 2 逆序求解由得出
得出SSOR公式:
【实验总结】
由以上实验得出SSOR在求解线性方程时最优。当时,逐次超松弛迭代法(SOR)等价于高斯赛德尔迭代法。(参考程序)
【程序】
%%%%%%%%%%%%%%%%%%%%%% 程序 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function maintest1(A,b)
N=input(请输入最大迭代次数N(足够大):);
esp=input(请输入近视解的误差限:);
X0=input(请输入初始值:);
if any(diag(A))==0
error(系数矩阵错误,迭代终止!)
end
D=diag(diag(A));
U=-triu(A);
L=-tril(A);
%%%%%%%%%%%%%%%jacobi%%%%%%
fJ=inv(D)*b;
BJ=inv(D)*(L+U);
PJ=max(abs(eig(BJ)));
fprintf(用雅克比求解的谱半径为:%2d\n,PJ)
%%%%%%%%%%%%%%%GS%%%%%%
fG=inv(D-L)*b;
BG=inv(D-L)*U;
PG=max(abs(eig(BG)));
fprintf(用GS求解的谱半径为:%2d\n ,PG)
%%%%%%%%%%%%% SOR %%%%%%%%%%
omega=1;
fS=omega*inv(D-omega*L)*b;
BS=inv(D-omega*L)*((1-omega)*D+omega*U);
PS=max(abs(eig(BS)));
fprintf(用SOR求解的谱半径为:%2d\n ,PS)
%%%%%%%%%%%% SSOR %%%%%%%%%%%%%%%
Lw=inv(D-omega*L)*((1-omega)*D+omega*U);
Uw=inv(D-omega*U)*((1-omega)*D+omega*L);
Sw=Uw*Lw;
fw=omega*(2-omega)*inv(D-omega*inv(D)*U)*inv(D-omega*inv(D)*L)*inv(D)*b;
PSS=max(
您可能关注的文档
- 第十届地球小博士地理科普知识大赛高中组试题9.doc
- 第十章 诉讼保障制度和程序.ppt
- 第十次平面广告文案新闻式故事式3.ppt
- 第十四章质谱分析法2.ppt
- 第十章 圆柱体作画步骤与解析3.ppt
- 第十章 国际货物买卖合同的订立6.ppt
- 第十章 发展中国家的贸易政策 25.ppt
- 第十章 建设中国特色社会主义文化PPT3.ppt
- 第十章 水污染及其控制6.ppt
- 第十章 抗震设计基本概念5.ppt
- 专题17 中国特色社会主义道路(第01期)-中考历史模拟题分项汇编(全国通用)(原卷版).docx
- 专题17 中国特色社会主义道路(第01期)-中考历史模拟题分项汇编(全国通用)(解析版).docx
- 专题20 宾语从句---三年(2024-2024)中考模拟题英语分项汇编(全国通用).docx
- 专题18 民族团结与祖国统一国防建设与外交成就(第01期)-2024年中考历史模拟题分项汇编(原卷版).doc
- 专题18 民族团结与祖国统一国防建设与外交成就(第02期)-2024年中考历史模拟题分项汇编(全国通用)(解析版).docx
- 专题15 中华人民共和国的成立和巩固(第02期)-中考历史模拟题分项汇编(全国通用)(解析版).docx
- 专题12 中华民族的抗日战争(第01期)-中考历史模拟题分项汇编(原卷版).doc
- 专题12 中华民族的抗日战争(第01期)-中考模拟题历史试题分项汇编(解析版).doc
- 专题12 中华民族的抗日战争(第02期)-中考历史模拟题分项汇编(全国通用)(原卷版).docx
- 专题12 中华民族的抗日战争(第02期)-中考历史模拟题分项汇编(全国通用)(解析版).docx
文档评论(0)