原创力文档- 1、本文档共4页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验2解线性方程组实验
实验2 解线性方程解线性方程实验内容：用高斯消去法解线性方程Uptrbk.m文件程序代码：
function X=uptrbk(A,B)
%Input -A is an N x N nonsingular matrix
% -B is an N x 1 matrix
%Output -X is an N x 1 matrix containing the solution to AX=B
%Initialize X and the temporary storage matrix C
n=length(B);
X=zeros(n,1);
Aug=[A B]
for p=1:n-1
for k=p+1:n
m=Aug(k,p)/Aug(p,p)
Aug(k,:)= Aug(k,:)-m*Aug(p,:)
end
end
X=backsub(Aug(:,1:n),Aug(:,n+1))
Backsub.m文件程序代码：
function X=backsub(A,B)
%Input - A is an n x n upper-triangular nonsingular matrix
% - B is an n x 1 matrix
%Output -X is the solution to the linear system AX=B
%Find the dimension of B and initialize X
n=length(B);
X=zeros(n,1);
X(n)=B(n)/A(n,n);
for k=n-1:-1:1
X(k)=(B(k)-A(k,k+1:n)*X(k+1:n))/A(k,k);
End
截图如右：
计算结果如下：
X=uptrbk(A,B)
Aug =
1 2 1 4 13
2 0 4 3 28
4 2 2 1 20
-3 1 3 2 6
m =
2
Aug =
1 2 1 4 13
0 -4 2 -5 2
4 2 2 1 20
-3 1 3 2 6
m =
4
Aug =
1 2 1 4 13
0 -4 2 -5 2
0 -6 -2 -15 -32
-3 1 3 2 6
m =
-3
Aug =
1 2 1 4 13
0 -4 2 -5 2
0 -6 -2 -15 -32
0 7 6 14 45
m =
1.5000
Aug =
1.0000 2.0000 1.0000 4.0000 13.0000
0 -4.0000 2.0000 -5.0000 2.0000
0 0 -5.0000 -7.5000 -35.0000
0 7.0000 6.0000 14.0000 45.0000
m =
-1.7500
Aug =
1.0000 2.0000 1.0000 4.0000 13.0000
0 -4.0000 2.0000 -5.0000 2.0000
0 0 -5.0000 -7.5000 -35.0000
0 0 9.5000 5.2500 48.5000
m =
-1.9000
Aug =
1.0000 2.0000 1.0000 4.0000 13.0000
0 -4.0000 2.0000 -5.0000 2.0000
0 0 -5.0000 -7.5000 -35.0000
0 0
文档评论（0）