- 1、本文档共5页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
消元法实验报告4
西京学院数学软件实验任务书
课程名称数学软件实验班级***学号***姓名***实验课题线性方程组高斯消去法,高斯列主元消去法,高斯全主元消去法实验目的熟悉线性代数方程组高斯消去法,高斯列主元消去法,高斯全主元消去法实验要求运用Matlab/C/C++/Java/Maple/Mathematica等其中一种语言完成实验内容线性方程组高斯消去法
线性方程组高斯列主元消去法
线性方程组高斯全主元消去法成绩教师
《数值分析》实验报告
实验一
一、实验目的与要求
1.掌握高斯列主元消去法解线性方程组的基本思路;
2.了解一些计算机的算法,会以某种汇编语言实现算法结果(本实验主要用matlab编程)
二、实验内容
1.编写用高斯列主元消元法解线性方程组的MATLAB程序,并求解下面的线性方程组,然后用逆矩阵解方程组的方法验证.
(1) (2)
2.列主元消元法及其matlab程序
function [Ra,Rb,n,X]=GaussXQLineMain(A,b)
%高斯列主元消元法,其中B为增广矩阵
B=[A b];
%读入b的长度
n=length(b);
%读出矩阵a,b秩
Ra=rank(A);
Rb=rank(B);
if (Rb-Ra)0
disp(因为Ra不等于Rb,所以此方程组无解.)
return
end
if Ra==Rb
if Ra==n
disp(因为Ra=Rb=n,所以此方程组有唯一解.)
X=zeros(n,1);
C=zeros(1,n+1);
for p= 1:n-1
%找出列中最大的元素并指出他的位置
[Y,j]=max(abs(B(p:n,p)));
C=B(p,:);
B(p,:)= B(j+p-1,:);
B(j+p-1,:)=C;
for k=p+1:n
m= B(k,p)/ B(p,p);
B(k,p:n+1)= B(k,p:n+1)-m* B(p,p:n+1);
end
end
b=B(1:n,n+1);A=B(1:n,1:n);
X(n)=b(n)/A(n,n);
for q=n-1:-1:1
X(q)=(b(q)-sum(A(q,q+1:n)*X(q+1:n)))/A(q,q);
end
else
disp(因为Ra=Rbn,所以此方程组有无穷多解.)
end
end
%调用格式
% clear;
% A=[1 2 3;4 5 6;7 8 9 ];
% b=[21;22;23];
% [Ra,Rb,n,X] =GaussXQLineMain(A,b)
三、实验过程
方程组(1)过程
clear;
A=[1 2 -2;1 1 1;2 2 1 ];
b=[1;1;1];
[Ra,Rb,n,X] =GaussXQLineMain(A,b)
因为Ra=Rb=n,所以此方程组有唯一解.
Ra =
3
Rb =
3
n =
3
X =
-3.0000
3.0000
1.0000
方程组(2)过程
clear;
A=[2 -1 1;1 1 1;1 1 -2 ];
b=[1;1;1];
[Ra,Rb,n,X] =GaussXQLineMain(A,b)
因为Ra=Rb=n,所以此方程组有唯一解.
Ra =
3
Rb =
3
n =
3
X =
0.6666
0.3333
0
在MATLAB中利用逆矩阵法检验结果:
(1) 在command windows中直接运行命令:
A=[1 2 -2;1 1 1;2 2 1 ];
b=[1;1;1];
X=A\b
结果
X =
-3.0000
3.0000
1.0000
(2) 在command windows中直接运行命令:
A=[2 -1 1;-1 1 1;1 1 -2 ];
b=[1;1;1];
X=A\b
结果
X =
0.6666
0.3333
0
四、实验总结
通过本次实验再次熟悉了高斯列主元消元法的思想,加深了对matlab语言的理解,简洁明了,在实验过程中函数编写实现不了,最后
文档评论(0)