原创力文档- 1、本文档共6页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
实验报告三
实验目的：
编写程序，完成实验题目
实验题目：
用LU分解和列主元消去法解线性方程组
输出Ax=b中系数A=LU分解的矩阵L和U，解向量x及detA；列主元法的行交换次序，解向量x及detA；比较两种方法所得的结果。
实验步骤：
1.新建函数lux.m
function [L,U,x]=lux(A,b)
[n,n]=size(A);
p=eye(n);//eye函数，生成一个单位矩阵
for k=1:n-1
[r,m]=max(abs(A(k:n,k)));
m=m+k-1;
if(A(m,k)~=0)
if(m~=k)
A([k m],:)=A([m k],:);
p([k m])=p([m k]);
end
for i=k+1:n
A(i,k)=A(i,k)/A(k,k);
j=k+1:n;
A(i,j)=A(i,j)-A(i,k)*A(k,j);
end
end
end
L=tril(A,-1)+eye(n,n);提取下三角矩阵
U=triu(A);提取上三角矩阵
newb=p*b;
y=zeros(n,1);全0
for k=1:n
j=1:k-1;
y(k)=(newb(k)-L(k,j)*y(j))/L(k,k);
end
x=zeros(n,1);
for k=n:-1:1
j=k+1:n;
x(k)=(y(k)-U(k,j)*x(j))/U(k,k);
end
2.输入矩阵A和b如下：
A=[10,-7,0,1;-3,2.099999,6,2;5,-1,5,-1;2,1,0,2]
b=[8;5.900001;5;1]
输入显示如下图：
3.调用函数lux计算结果
[l,u,x]=lux(A,b)
4.调用函数det求行列式
det(A)
计算结果如下图：
其中：lux计算结果如下:
l =
1.0000 0 0 0
0.5000 1.0000 0 0
-0.3000 -0.0000 1.0000 0
0.2000 0.9600 -0.8000 1.0000
u =
10.0000 -7.0000 0 1.0000
0 2.5000 5.0000 -1.5000
0 0 6.0000 2.3000
0 0 0 5.0800
x =
0.1949
-0.7661
0.9695
0.6882
调用函数det求行列式如下：
ans =
-762.0001
小组成员：4100311 王恒
4100320 周冀豪
4100325 郭静文
4100326 林佳钰
4100328 宁薇
文档评论（0）