原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
一编写Gauss列主元消元法求解线性方程组的程序
数值计算报告书
日期：2014 年11月8日
问题
名称 编写 Gauss 列主元消元法求解线性方程组的程序，要求附有算例。 问题描述：
编写一MATLAB程序，算法为列主元消元法，功能为求解线性方程组。 MATLAB程序：
function[x,flag]=Gauss(A,b)
[n,m]=size(A);nb=length(b);
if n~=m
error(The rows and columns of matrix A must be equal!);
% 不是方阵，则显示错误。
end
if m~=nb
error(The columns of A must be equal the length of b!)
% 如果列数不等于b的个数，则显示错误。
end
flag=OK;x=zeros(n,1);
for k=1:n-1
max1=0;
for i=k:n
if abs(A(i,k))max1
max1=abs(A(i,k));r=i;
end
end
% 记max1为最大值，r为最大值所在的行数。
if max11e-10
flag=failure;return;
end
% 最大值太小，近似为0，则返回。
if rk
for j=k:n
z=A(k,j);A(k,j)=A(r,j);A(r,j)=z;
end
z=b(k);b(k)=b(r);b(r)=z;
end
% 两条方程交换位置。
for i=k+1:n
m=A(i,k)/A(k,k);
for j=k+1:n
A(i,j)=A(i,j)-m*A(k,j);
end
b(i)=b(i)-m*b(k);
end
% 当前方程已下全部方程消元。
end
if abs(A(n,n))1e-10
flag=failure;return;
end
x(n)=b(n)/A(n,n);
for k=n-1:-1:1
for j=n:-1:k+1
b(k)=b(k)-A(k,j)*x(j);
end
x(k)=b(k)/A(k,k);
end
% 求出方程组的解x(k)。
算例：
求解下列方程组：
在MATLAB的Command Window中输入：
A=[0.729 0.81 0.9;1 1 1;1.331 1.21 1.1];
b=[0.6867;0.8338;1.000];
Gauss(A,b)
计算结果与分析：
计算结果如下:
ans =
0.2245
0.2814
0.3279
数据分析：
将=0.2245，=0.2814，=0.3279代入，得结果为0.6867045；代入，得结果为0.8338；代入，得结果为0.9999935。可见，
Gauss 列主元消元法的求解精度较高。
文档评论（0）