- 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 宁薇
您可能关注的文档
- 操作系统模拟unix文件系统实验.doc
- 操作系统进出口预算表的核算过程.doc
- 攀枝花市块菌产业建设项目建议书攀西特色生物资源开发项目.doc
- 支付结算征文.doc
- 支出费用成本之间的关系研究.doc
- 支持向量机(SVM)及支持向量回归(SVR)简介.doc
- 支票、本票、汇票的联系和区别(完整列表).doc
- 支行货币押运突发事.doc
- 收购股权中介服务合同(样本).doc
- 改进型控制棒驱动机构电源研究.docx
- 移动机器人技术-课件 项目4-任务1-STM32集成开发环境的安装与应用.pptx
- 移动机器人技术-课件 项目4-任务4-舵机的控制.pptx
- 《茶文化与茶艺》 知识点汇总 重大 不同茶类冲泡水温建议 ----中国五大名泉.pdf
- 《WPS办公应用项目化教程(中级)》教案 项目1 长文档管理_教案.docx
- 大数据挖掘导论与案例课件-第5章 分类概念与方法.pptx
- 大数据挖掘导论与案例课件-第2章 数据分析和可视化分析.pptx
- 《工程与社会》 课件 第一章绪论.pptx
- 项目1-任务1-移动机器人分类及行业应用.pptx
- 商务统计学(第7版)英文ppt课件11 Chi-Square Tests、13 Multiple Regression.ppt
- 移动机器人技术-课件 项目5-任务4-超声波传感器的应用.pptx
文档评论(0)