24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1761  |  回复: 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的回帖

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的回帖

aixilin

新虫 (小有名气)

小雨萌萌: 在火星?表弄火星文啊~ 2011-07-12 08:53:09
L'琌L锣竚. U'琌U锣竚.
5楼2011-07-12 07:01:21
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

aixilin

新虫 (小有名气)

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

L'是L的转置, U'是U的转置
6楼2011-07-12 07:05:28
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sskkyy

银虫 (正式写手)

【答案】应助回帖


soliton923(金币+1): 谢谢参与讨论~~ 2011-07-12 21:43:22
可不可以这样考虑:两边同时乘以L'的逆,得到L=(L')^-1 U' U. 注意L和U‘都是下三角矩阵,从U’的最后一行开始考虑。假设他们都是n乘以n的矩阵。如果你把这些三角矩阵写出来,你会发现L和U'的(n,1)-th 元素是一样的,接着考虑(n,2)-th 元素,有个简单的关系,然后递归就可以考虑任意(n,i)-th 元素,再考虑n-1行......
7楼2011-07-12 13:42:28
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

sskkyy

银虫 (正式写手)

【答案】应助回帖

注意到L'和U都是上三角,你就可以考虑L=(L')^-1 U' U的两边,开始比较对应元素,从第n行开始。
8楼2011-07-12 13:45:07
已阅   关注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的回帖

sskkyy

银虫 (正式写手)

【答案】应助回帖

aixilin(金币+3): 谢谢你的回复和帮助 2011-07-14 16:57:47
不好意思,我误解了你的“上三角矩阵”的意思。我还以为你说的是严格上三角呢,也就对角线上全是1. 你只需要乘以一个对角矩阵把一般的上下角变成严格上三角就行。也就是一个上三角总是一个对角矩阵和严格上三角的乘积。下三角是类似的。对角矩阵的乘法是很简单的,我相信你应该能够自己搞定。
10楼2011-07-13 01:43:17
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 aixilin 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程302分求调剂 +17 zyx上岸! 2026-04-04 17/850 2026-04-10 09:50 by may_新宇
[考研] 344求调剂 +6 丶风雪夜归人丶 2026-04-09 6/300 2026-04-10 09:23 by sailor_y
[考研] 一志愿矿大,材料工程专硕314分,0856可调都可以 +14 无懈可击的巨人 2026-04-09 14/700 2026-04-10 08:10 by Sammy2
[考研] 288求调剂,一志愿华南理工大学071005 +16 ioodiiij 2026-04-08 16/800 2026-04-09 23:08 by wolf97
[考研] 求调剂材料科学与工程一志愿985初试365分 +5 材化李可 2026-04-08 5/250 2026-04-09 17:00 by Lilly_Li
[考研] 材料295 +15 小英11 2026-04-03 16/800 2026-04-09 15:44 by only周
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 296求调剂 +3 汪!?! 2026-04-08 3/150 2026-04-08 22:00 by zhouyuwinner
[考研] 327求调剂 +12 Xxjc1107. 2026-04-06 12/600 2026-04-08 16:46 by luoyongfeng
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 材料专硕(0856) 339分求调剂 +16 哈哈哈鹅哈哈哈 2026-04-05 16/800 2026-04-08 16:02 by luoyongfeng
[考研] 化学0703-一志愿211-338分求调剂 +10 vants 2026-04-05 11/550 2026-04-08 16:02 by screening
[考研] 304求调剂 +10 素年祭语 2026-04-06 17/850 2026-04-08 09:05 by 蓝云思雨
[考研] 283分求调剂 +14 试试看呗 2026-04-04 14/700 2026-04-08 07:03 by lijunpoly
[考研] 307求调剂 +3 Youth@@ 2026-04-07 3/150 2026-04-07 22:00 by hemengdong
[考研] 生物调剂 +5 橙子橙子橙子啊 2026-04-05 9/450 2026-04-07 15:31 by 上岸快快
[考研] 生物学学硕求调剂:351分一志愿南京师范大学生物学专业 +6 …~、王…~ 2026-04-06 7/350 2026-04-06 18:54 by macy2011
[考研] 22408 331分求调剂 +4 y__1 2026-04-06 4/200 2026-04-06 17:26 by 土木硕士招生
[考研] 359求调剂 +7 hhhhaaaa$ 2026-04-04 7/350 2026-04-04 18:49 by imissbao
[考研] 266求调剂 +8 学员97LZgn 2026-04-03 8/400 2026-04-04 09:02 by 20021109
信息提示
请填处理意见