- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
习题二
1证明:三对角方程组 错误!未找到引用源。 如果满足条件 错误!未找到引用源。, 则Gauss顺序消去法一定可进行到底(等价地系数矩阵存在 LU分解),进一步说明系数矩
阵是可逆矩阵。
证明已知
_biaiC1b2
_bi
ai
C1
b2
an J
cn _J
bn
卜 崗0 满足条件]匕亠冋_|
Jbn|
* Cj , a^j, ci
0
= 0(i
= 2,3,| 山 n-1)
只需证明按式(1.19)的计算步骤:
PiPkqk
Pi
Pk
qk
Pn
二 D 4=4 / P1
=bk - Ck」q2 八
/ Kk Fk/ Pk
二 bn - Cn Jqn A
= 2,|||, n-1)
分解得到的Pk =0(k =1,2JH, n),这样必有LU
分解
■
a1 b2 \
q 1
P2
? Cn」
I- 1
h r
■ Cn」
an .A bn
[
qn」1 _
- Pn 一
bi
Pi
C1
G
从而也说明了 A是可逆矩阵。
由假设bj G
由假设bj G 0,得
g
Pi
P1 二 b = 0,再由
:::1 得
又知又知最后C2P2C3P3:::1
又知
又知
最后
C2
P2
C3
P3
:::1,从而
b3 - 2 a?
P2
占 - 2 |还 |鸟-比| ^|C3 0
P2
:::1,类似可证 p4 0,如此继续下去,可得
Pk 0,
Ck
Pk
:1,k=1,2,川,n_1
P2
b2 _Ciqi
=
b2 --CLa1
cl
ai a b2 - a
C2
Pi
Pn| = bn -CnAqn』纠加—学 a.」 bn — anj 0pn J
(i十,川,n)2 ?设A是严格对角占优矩阵,即 A二(aij )
(i十,川,n)
3n Z 3ij j 4 j:J:
又经过一步Gauss顺序消去之后, A化为
311_0a:A
311
_0
a:
A
证明A2也是严格对角占优矩阵,进一步说明
A是可逆矩阵。
证明记氏二a(2),
则有
3i(
3i(2)
ai1
=aj … a1 j
an
又A是对角占优矩阵,可知aiinE aj ,i =1,2」||,n。故nzj =2j止⑵aij3in=zj =2j Y313
又A是对角占优矩阵,可知
aii
n
E aj ,i =1,2」||,n。故
nzj =2j止
⑵
aij
3i
n
=z
j =2j Y
31
311
3i1aj … a1 j
an
an
n
_ 乞 a1jjf
ai1
(2)
ai 一 a1i
=
aii
a11
n
z
j=2j=i:
aij
J lai1
ai1
a11
311
aii
a1 j
n
=z
j=1
K-i
n
-z
j=2
aij
3j — 3i1
aii
(i =2,川,n)
a1an
aii
n
为a1j
ji2
ai1
a11
3ii
即A也是严格对角占优矩阵。
由A是对角占优矩阵知an = 0 ;又因为矩阵A2也是严格对角占优矩阵,以此类推可知, 经过n -1步消元后A所化成的上三角矩阵对角线上的元素都大于 0,故A是可逆矩阵。
3.证明如果 非奇异A有LDU分解错误!未找到引用源。,则分解式是唯一的。
证明 因为A非奇异,所以 L,D,U都非奇异,假设 A还存在另一个 LDU分解:
- L1 D1U1,这里L1, D1,U1也一定非奇异。
从
LDU = L1DU 1
可得
2-
2
-22
lFl 二 DU lU JD J
上式左端是单位下三角矩阵,右端是上三角矩阵,所以两端都应是单位矩阵,因此
TOC \o 1-5 \h \z 1 _1 1
即因此L.L = I, DU iU D
即
因此
Li = L, UiU 丄=Di」D = I
Ui 二 U,Di 二 D
0 5 8 i3
4?设 A= 6 2 2 ,b= i0 ,
]6 i2 4一 .22 _
5.设
5.设A是实对称正定矩阵,试推导按
列计算L的Cholesky分解A - LLT算法。
lii
l
ii
(1) 试对A进行PLU分解:PA = LU ;
(2) 根据PLU分解求解Ax二b。
q 、
■z6 2 2
‘0 i 0、
解(i) L =
i i
,U =
i0 2
,P =
0 0 i
e 0.5 i 』
7」
J 0 0」
(2) x = (i,i,i)T
hi2 ik
hi
j 4
|ij - ( aij ■- |ik 1 jk ) / 1 jj ( j - i 们山 n)
k=i
end
6.证明定理2.i。
证明
— 2—)^4:
— 2—)^4:
—
■T:
又因为|a _P『2 =(□ _P )=(口 —P)T (口 —P )=(ocT _^打(口 _ p)
又因为
=
您可能关注的文档
最近下载
- UbuntuLinux操作系统第2版(微课版)钟小平课后习题答案.doc VIP
- 输电线路三跨的经验分享.ppt VIP
- 电动汽车分布式驱动系统动力学分析与优化设计.docx VIP
- 专题1.4数轴与动点经典题型(四大题型)(原卷版+解析).docx VIP
- 幼儿小蝌蚪找妈妈绘本ppt课件.pptx VIP
- 磷酸铁行业现状.docx VIP
- 马克思主义与社会科学方法论课后习题答案(2018版).pdf VIP
- 6.1 树立法治观念 统编版道德与法治八年级上.pptx VIP
- GB_T 43933-2024 金属矿土地复垦与生态修复技术规范.pdf VIP
- 车牌识别系统软件使用说明8.1软件基本设置.pdf VIP
原创力文档


文档评论(0)