- 1、本文档共8页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
§3 LU分解法
——Gauss消去法的变形
知识预备:
1矩阵的初等行变换、初等矩阵及其逆、乘积
2矩阵的乘法
3上三角矩阵的乘积、单位下三角矩阵的乘积
4单位下三角矩阵的逆、可逆的上三角矩阵的逆
一、Gauss消去法的矩阵解释
Gauss消去法实质上是将矩阵A分解为两个三角矩阵相乘。
我们知道,矩阵的初等行变换实质就是左乘初等矩阵。
第一轮消元:相当于对A(1)左乘矩阵L1,即
其中
第二轮消元:对应于
一般地
……………(1)
其中
整个消元过程为
………(2)
从而
其中L是单位下三角矩阵,即
…(3)
【注】消元过程等价于A分解成LU的过程
回代过程是解上三角方程组的过程。
二、矩阵的三角分解
1、若将A分解成L?U,即A=L?U,其中L为单位下三角矩阵,U为非奇异上三角矩阵,则称之为对A的Doolittle分解。
当A的顺序主子式都不为零时,消元运算可进行,从而A存在唯一的Doolittle分解。
证明:若有两种分解,A=L1U1,A=L2U2,则必有L1=L2,U1=U2。
因为L1U1=L2U2,而且L1,L2都是单位下三角矩阵,U1,U2都是可逆上三角矩阵,所以有
因此
即
L1=L2,U1=U2、
2、若L是非奇异下三角矩阵,U是单位上三角矩阵时,A存在唯一的三角分解,A=LU,称其为A的Crout分解(对应于用列变换实施消元)
三、直接分解(LU分解)算法
LU分解算法公式——按矩阵乘法
第一步:利用A中第一行、第一列元素确定U的第一行、L的第一列元素。由
得 u1j=a1j )
li1=ai1/u11)
第r步:利用A中第r行、第r列剩下的元素确定U的第r行、L的第r列元素(r=2,3,…,n).由
得U的第r行元素为
由得
…………(4)
直接分解的紧凑格式:
u11 u12 u13 … … u1n 1
l21 u22 u23 … … u2n 2
l31 l32
ln1 ln2 unn n
方程组的三角分解算法(LU分解)
对于方程组Ax=b,设A=LU (Doolittle分解)。
由于
1、求解Ly=b:
…………………(5)
2、求解Ux=y:
…(6)
LU分解算法
步1,输入A,b;
步2,对j=1,2,…,n 求
对i=2,3,…,n 求
步3,对r=2,3,…,n 做(3.1)-(3.2):
(3.1)
(3.2)
步4,
步5,
步6,输出结束。
例子与程序:
【】
解:对系数矩阵A进行LU分解
因此
先解。
再解
程序:LU_factorization
%Not Select Column LU_factorization
clear all
n=3;a=[2 2 3;4 7 7;-2 4 5];b=[3;1;-7];
%n=3;a=[1 4 7;2 5 8;3 6 11];b=[1;1;1];
%LU_factorazation
for i=2:n
a(i,1)=a(i,1)/a(1,1);
end
a
for r=2:n
for j=r:n
s=0.;
for k=1:r-1
s=s+a(r,k)*a(k,j);
end
a(r,j)=a(r,j)-s;
end
for i=r+1:n
s=0.;
for k=1:r-1
s=s+a(i,k)*a(k,r);
end
a(i,r)=(a(i,r)-s)/a(r,r);
end
a
end
%Extract Lower/Upper Triangular Part
l=tril(a);
for i=1:n
l(i,i)=1;
end
u=triu(a);
l
u
%Linear Low
您可能关注的文档
- 正定矩阵的性质和判定方法及应用.doc
- 第四章 矩阵的相似标准形.doc
- 线性代数方程组求解.doc
- 考研数学三公式大全.doc
- 11111文献实验十二烷基苯磺酸钠的制备.doc
- 2008专升本数学与应用数学《专业综合》试卷.doc
- 中学数学教育专题选讲心得体会.doc
- 关于“矩阵的行列式不等式”的几点注记.doc
- 安徽大学2012—2013学年第二学期 高等代数.doc
- 教育基础知识和基本原理第一章.doc
- CN105452440B 细胞的培养方法、细胞培养装置及试剂盒.docx
- 2024年湖南省怀化市通道侗族自治县木脚乡招聘社区工作者真题参考答案详解.docx
- 2024年湖南省怀化市通道侗族自治县江口乡招聘社区工作者真题含答案详解.docx
- 2024年湖南省怀化市通道侗族自治县杉木桥乡招聘社区工作者真题及完整答案详解1套.docx
- 2024年湖南省怀化市通道侗族自治县播阳镇招聘社区工作者真题及完整答案详解1套.docx
- 高技术产业集聚对绿色全要素生产率影响研究——基于集聚外部性视角.docx
- 2024年湖南省怀化市通道侗族自治县戈冲苗族乡招聘社区工作者真题及答案详解1套.docx
- 2024年湖南省怀化市通道侗族自治县戈冲苗族乡招聘社区工作者真题附答案详解.docx
- 基于电子听诊的新生儿先天性心脏病智能听诊算法研究.docx
- 2024年湖南省怀化市通道侗族自治县戈冲苗族乡招聘社区工作者真题及完整答案详解1套.docx
最近下载
- 高中生物选择性必修二知识点总结 .pdf VIP
- 中小学学校食堂食品安全方面存在的问题及对策建议思考.docx VIP
- 笔墨丹青润中考书法鉴赏.ppt VIP
- 初中语文_《笔墨丹青润中考——书法欣赏专题》教学设计学情分析教材分析课后反思.doc
- 《干冰清洗技术应用》课件介绍.ppt
- 国际数据管理协会_dama_数据管理知识体系指南.pdf VIP
- 非金属材料试验检测含答案.doc VIP
- VTE风险评估工具及应用.pptx VIP
- 2024高考复习必背高中英语词汇3500单词完整打印版 .pdf VIP
- Schneider Electric施耐德HVX 40.5kV(K)中压真空断路器中文操作手册安装和用户指南(中文).pdf
文档评论(0)