24小时热门版块排行榜    

查看: 1682  |  回复: 10
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

aixilin

新虫 (小有名气)

[求助] 若已知上三角矩阵U和下三角矩阵L满足 L‘*L = U’*U,如何有效的通过L求U

如题,若已知上三角矩阵U和下三角矩阵L满足 L‘*L = U’*U,如何有效的通过L求出U?

P.S. 我已经知道可以通过MATLAB中的QR分解来计算:U=QR(L),但是QR分解针对的一般矩阵,所以计算复杂度是O(N^3),这里N 是矩阵L的维数。由于L是三角矩阵,我在想是否存在一种更快速的计算方法,比如复杂度为O(N^2)的解决方法。

数学知识有限,敬请大家帮忙给点建议。。谢谢!
回复此楼
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aixilin

新虫 (小有名气)

引用回帖:
Originally posted by sskkyy at 2011-07-12 13:42:28:
可不可以这样考虑:两边同时乘以L'的逆,得到L=(L')^-1 U' U. 注意L和U‘都是下三角矩阵,从U’的最后一行开始考虑。假设他们都是n乘以n的矩阵。如果你把这些三角矩阵写出来,你会发现L和U'的(n,1)-th 元素是一样 ...

很感谢你的回答,为了说明问题,我用2*2 矩阵简单测试下您的方法:

L = [a1, 0; c1,d1];
U = [a2, b2; 0, d2];

D = (L')^-1 = [a3, b3; 0 ,d3]
L = (L')^-1*U'*U
  =  [a3, b3; 0 ,d3]*[a2, 0; b2, d2]*[a2, b2; 0, d2]
  = [a3, b3; 0, d3] * [a2^2, a2*b2; b2*a2, b2^2+d2^2 ]
  = [a3*a2^2+b3*b2*a2, a3*a2*b2+b3*(b2^2+d2^2); d3*b2*a2, d3*(b2^2+d2^2)]

在这里等号左侧L的(n,1)-th元素是c1,而右侧为d3*b2*a2,很明显他们不一样!
L和U'的(n,1)-th 元素是一样的如何得到的?
9楼2011-07-12 17:38:12
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 11 个回答

sskkyy

银虫 (正式写手)

【答案】应助回帖

L'是什么?是L的转置还是另外一个下三角矩阵?还有u‘呢?
2楼2011-07-12 03:43:02
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aixilin

新虫 (小有名气)

引用回帖:
Originally posted by sskkyy at 2011-07-12 03:43:02:
L'是什么?是L的转置还是另外一个下三角矩阵?还有u‘呢?

L'㎝U'常琌锣竚痻皚
3楼2011-07-12 06:50:23
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aixilin

新虫 (小有名气)

引用回帖:
Originally posted by sskkyy at 2011-07-12 03:43:02:
L'是什么?是L的转置还是另外一个下三角矩阵?还有u‘呢?

L'琌L锣竚. U'琌U锣竚.
4楼2011-07-12 06:55:04
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] (出分啦!!!)26 考研调剂 20w+调剂信息火速更新中!【捡漏】qq 群:1087823172 +3 db58 2026-02-28 133/6650 2026-02-28 18:25 by nbf68
[教师之家] 版面费该交吗 +15 苹果在哪里 2026-02-22 18/900 2026-02-28 18:20 by mibaomingg
[考研] 285求调剂 +5 满头大汗的学生 2026-02-28 5/250 2026-02-28 18:10 by 材料专硕调剂;
[考研] 材料学调剂 +3 提神豆沙包 2026-02-28 3/150 2026-02-28 18:10 by djennjx
[考研] 材料类求调剂 +4 wana_kiko 2026-02-28 4/200 2026-02-28 18:08 by djennjx
[考研] 298求调剂 +4 人间唯你是清欢 2026-02-28 4/200 2026-02-28 18:08 by njzyff
[考研] 材料调剂 +3 爱擦汗的可乐冰 2026-02-28 3/150 2026-02-28 18:06 by houyaoxu
[考研] 085600材料工程一志愿中科大总分312求调剂 +5 吃宵夜1 2026-02-28 7/350 2026-02-28 17:57 by njzyff
[考研] 311求调剂 +6 南迦720 2026-02-28 6/300 2026-02-28 17:17 by 刘茉莉
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
[考博] 博士推荐 +4 花儿笑? 2026-02-21 5/250 2026-02-28 16:55 by mumin1990
[考研] 265分求调剂不调专业和学校有行学上就 +4 礼堂丁真258 2026-02-28 6/300 2026-02-28 16:18 by 求调剂zz
[考博] 26申博 +3 想申博! 2026-02-26 3/150 2026-02-28 16:07 by nxgogo
[考研] 290求调剂 +4 材料专硕调剂; 2026-02-28 5/250 2026-02-28 13:32 by houyaoxu
[考研] 0856调剂 +3 刘梦微 2026-02-28 3/150 2026-02-28 13:22 by houyaoxu
[考研] 寻找调剂 +3 LYidhsjabdj 2026-02-28 3/150 2026-02-28 12:59 by miniwendy
[考研] 304求调剂 +5 曼殊2266 2026-02-28 6/300 2026-02-28 12:44 by 迷糊CCPs
[硕博家园] 博士自荐 +6 科研狗111 2026-02-26 9/450 2026-02-28 12:32 by seaskyy
[考研] 272求调剂 +3 田智友 2026-02-28 3/150 2026-02-28 12:31 by 王加浩to
[考研] 298求调剂 +4 axyz3 2026-02-28 4/200 2026-02-28 11:21 by wang_dand
信息提示
请填处理意见