24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1764  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研二轮调剂 +8 故人?? 2026-04-09 8/400 2026-04-10 09:44 by 青梅duoduo
[考研] 一志愿郑州大学 22408 305分求调剂 +3 安小满zzz 2026-04-08 3/150 2026-04-09 16:16 by wp06
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 265求调剂 +4 风说她早忘了 2026-04-07 4/200 2026-04-09 13:59 by only周
[考研] 085801 总分275 本科新能源 求调剂 +8 bradoner 2026-04-08 9/450 2026-04-09 13:43 by only周
[考研] 0703总分331求调剂 +18 ZY-05 2026-04-04 22/1100 2026-04-09 09:09 by 哦哦123
[考研] 一志愿985初试354分生物调剂 +3 031001 2026-04-06 3/150 2026-04-09 00:30 by Evan_Liu
[考研] 材料调剂 +14 一样YWY 2026-04-06 14/700 2026-04-08 23:00 by 猪会飞
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 调剂 +18 不逢春 2026-04-05 19/950 2026-04-07 22:04 by lijunpoly
[考研] 302分求调剂 一志愿安徽大学085601 +12 zyx上岸! 2026-04-04 12/600 2026-04-07 02:09 by BruceLiu320
[考研] 化学357分,考研调剂 +11 .Starry. 2026-04-04 12/600 2026-04-06 06:28 by houyaoxu
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 328分调剂 +6 门men 2026-04-04 6/300 2026-04-05 13:40 by imissbao
[考研] 083200 333求调剂 +3 十二!! 2026-04-04 3/150 2026-04-05 08:28 by barlinike
[考研] 一志愿沪985,326分求调剂 +3 刘墨墨 2026-04-03 3/150 2026-04-04 11:16 by 悲伤的芋头
[考研] 考研调剂 +3 15615482637 2026-04-03 3/150 2026-04-03 22:50 by ms629
[考研] 求调剂 +8 akdhjs 2026-04-03 8/400 2026-04-03 18:17 by 戴维ING
[考研] 372分材料与化工(085600)一志愿湖南大学求调剂 +3 蓝笺片 2026-04-03 4/200 2026-04-03 17:58 by Jimmyandyou
信息提示
请填处理意见