原创力文档- 1、本文档共11页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
数值分析上机实验3
昆 明 理 工 大 学 理 学 院
信息与计算科学专业 操作性实验报告
年级： 10级 姓名： 刘陈 学号： 201011101128 指导教师： 陈智斌
实验课程名称： 数值分析matlab程序设计 开课实验室： 理学院机房
实验成绩：
学风(5) 观察能力(15) 操作能力(30) 调试能力(50) 其它 总分
实验内容：
1．实验/作业题目:数值分析与matlab程序设计
2．实验/作业课时：2学时
3．实验过程(包括实验环境、实验内容的描述、完成实验要求的知识或技能)：
实验环境：matlab
实验内容：
1.列主元消去法求L矩阵，U矩阵，P矩阵，det(A)和解向量X的值。通过实验，掌握matlab列主元消去法操作。
2. Implement the Gaussian elimination with partial pivoting and solve the following systems of linear equations:
(i) =
(ii) =
(iii) =
Requirement:
(1) Output the LU decomposition of the coefficient matrix A, i.e., L and U such that PA = LU.
(2) Output the solution X = A?1b.
(3) Output det(A).
完成实验要求的知识或技能：通过这个实验，了解matlab开发环境的配置以及如何通过matlab来编写《数值分析》程序，掌握matlab程序的一些基本结构。
4．程序结构(程序中的函数调用关系图)
5．算法描述、流程图或操作步骤：
高斯列主元消去法M文件算法：
function [l,u,p]=Gauss_x(A,b)
[m,n]=size(A);
if m~=n
error(矩阵不是方阵)
return
end
if det(A)==0
error(矩阵不能被三角分解)
end
u=A;p=eye(m);l=eye(m);
for i=1:m
for j=i:m
t(j)=u(j,i);
for k=1:i-1
t(j)=t(j)-u(j,k)*u(k,i);
end
end
a=i;B=abs(t(i));
for j=i+1:m
if Babs(t(j))
B=abs(t(j));
a=j;
end
end
if a~=i
for j=1:m
c=u(i,j);
u(i,j)=u(a,j);
u(a,j)=c;
end
for j=1:m
c=p(i,j);
p(i,j)=p(a,j);
p(a,j)=c;
end
c=t(a);
t(a)=t(i);
t(i)=c;
end
u(i,i)=t(i);
for j=i+1:m
u(j,i)=t(j)
文档评论（0）