- 1、本文档共5页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析线性方程组迭代法实验.doc
实验报告
一、实验目的
体会求解线性方程组的Jacobi迭代法和Gauss-Seidel迭代法的收敛速度。
二、实验题目
用迭代法分别对解方程组,其中
(1)选取不同的初值和不同的右端向量,给定迭代误差,用Jacobi迭代法和Gauss-Seidel迭代法计算,观测得到的迭代向量并分析计算结果给出结论;
(2)取定初值和右端向量,给定迭代误差,将的主对角元成倍放大,其余元素不变,用Jacobi迭代法计算多次,比较收敛速度,分析计算结果并给出结论。
三、实验原理
求解线性方程组的Jacobi迭代法和Gauss-Seidel迭代法
四、实验内容及结果
1.Jacobi迭代法:
1)先输入Jacobi迭代的M程序,在MATLAB的M文件窗口中输入以下内容并保存:
function x=majacobi(A,b,x0,ep,N)
n=length(b);
if nargin5,N=500;end
if nargin4,ep=1e-6;end
if nargin3,x0=zeros(n,1);end
x=zeros(n,1);k=0;
while kN
for i=1:n
x(i)=(b(i)-A(i,[1:i-1,i+1:n])*x0([1:i-1,i+1:n]))/A(i,i);
end
if norm(x-x0,inf)ep,break;end
x0=x;k=k+1;
end
if k==N,Warning(已达到迭代次数上限);end
disp([k=,num2str(k)])
2)再输入A,b,x0,在M文件窗口输入以下内容并保存:
n=20;
A=diag(4*ones(1,n));
B=diag(-1/3*ones(1,n-1),1);
C=diag(-1/5*ones(1,n-2),2);
A=A+B+C+B+C;
b=ones(n,1);
x0=zeros(n,1);
x= majacobi(A,b,x0)
x’
3)在MATLAB的主程序窗口输出以下结果:
k=10
ans =
Columns 1 through 10
0.2935 0.3195 0.3366 0.3394 0.3406 0.3408 0.3409 0.3409 0.3409 0.3409
Columns 11 through 20
0.3409 0.3409 0.3409 0.3409 0.3408 0.3406 0.3394 0.3366 0.3195 0.2935
2.Gauss-Seidel迭代法:
1)先输Gauss-Seidel迭代的M程序,在MATLAB的M文件窗口中输入以下内容并保存:
function x=maseidel(A,b,x0,ep,N)
n=length(b);
if nargin5,N=500;end
if nargin4,ep=1e-6;end
if nargin3,x0=zeros(n,1);end
x=zeros(n,1);k=0;
while kN
for i=1:n
if i==1
x(1)=(b(1)-A(1,2:n)*x0(2:n))/A(1,1);
else if 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
end
if norm(x-x0,inf)ep,break;end
x0=x;k=k+1;
end
if k==N,Warning(已达到迭代次数上限);end
disp([k=,num2str(k)])
2)再输入A,b,x0,在M文件窗口输入以下内容并保存:
n=20;
A=diag(4*ones(1,n));
B=diag(-1/3*ones(1,n-1),1);
C=diag(-1/5*ones(1,n-2),2);
A=A+B+C+B+C;
b=ones(n,1);
x0=zeros(n,1);
x=maseidel(A,b,x0);
x
3)在MATLAB的主程序窗口输出以下结果:
k=7
ans =
Columns 1 through 10
0.2935 0.3195 0.3366
您可能关注的文档
- 《西游记》知识竞赛1-25回.doc
- 2019高考专题: 名句名篇默写.ppt
- 《观察水》教学设计反思.doc
- 《说勤》课文.doc
- 22.1.3沪科版九年级《平行线分线段成比例定理》.ppt
- 传感器与传感器技术(何道清)课后答案.doc
- 《走进空竹》说课案例.doc
- 实用综合教程(第二版)外教社-第1册教案.doc
- 《跨越式跳高——直腿摆动过杆》教学设计.doc
- 实验A观察种子结构实.ppt
- 黑龙江齐齐哈尔市(新版)2024小学语文人教版小升初摸底(预测卷)完整试卷(含答案).docx
- 黑龙江哈尔滨市(新版)2024小学语文部编版小升初质量检测(自测卷)完整试卷(含答案).docx
- 黑龙江齐齐哈尔市(新版)2024小学语文苏教版小升初质量检测(评估卷)完整试卷(含答案).docx
- 青海省玉树藏族自治州(新版)2024小学语文苏教版小升初测试(评估卷)完整试卷(含答案).docx
- 黑龙江佳木斯市(新版)2024小学语文统编版(五四制)小升初摸底(评估卷)完整试卷(含答案).docx
- 黑龙江齐齐哈尔市(新版)2024小学语文统编版小升初模拟(综合卷)完整试卷(含答案).docx
- 青海省玉树藏族自治州(新版)2024小学语文部编版小升初质量检测(评估卷)完整试卷(含答案).docx
- 黑龙江双鸭山市(新版)2024小学语文苏教版小升初摸底(评估卷)完整试卷(含答案).docx
- 黑龙江伊春市(新版)2024小学语文苏教版小升初摸底(自测卷)完整试卷(含答案).docx
- 黑龙江哈尔滨市(新版)2024小学语文苏教版小升初模拟(自测卷)完整试卷(含答案).docx
文档评论(0)