24小时热门版块排行榜    

查看: 1714  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 323求调剂 +3 洼小桶 2026-03-18 3/150 2026-03-20 22:54 by JourneyLucky
[考研] 317求调剂 +5 申子申申 2026-03-19 9/450 2026-03-20 22:26 by JourneyLucky
[考研] 350求调剂 +5 weudhdk 2026-03-19 5/250 2026-03-20 22:04 by luoyongfeng
[考研] 279分求调剂 一志愿211 +10 chaojifeixia 2026-03-19 11/550 2026-03-20 20:32 by JourneyLucky
[考研] 一志愿西安交通大学 学硕 354求调剂 +3 我想要读研究生 2026-03-20 3/150 2026-03-20 20:13 by JourneyLucky
[考研] 353求调剂 +3 拉钩不许变 2026-03-20 3/150 2026-03-20 19:56 by JourneyLucky
[考研] 环境工程调剂 +9 大可digkids 2026-03-16 9/450 2026-03-20 17:38 by 醉在风里
[考研] 085410人工智能专硕317求调剂(0854都可以) +4 xbxudjdn 2026-03-18 4/200 2026-03-20 09:07 by 不168
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 085600材料与化工求调剂 +6 绪幸与子 2026-03-17 6/300 2026-03-19 13:27 by houyaoxu
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 070300化学319求调剂 +6 锦鲤0909 2026-03-17 6/300 2026-03-18 13:22 by Iveryant
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 机械专硕325,寻找调剂院校 +3 y9999 2026-03-15 5/250 2026-03-16 19:58 by y9999
[考研] 0703化学调剂 290分有科研经历,论文在投 +7 腻腻gk 2026-03-14 7/350 2026-03-16 10:12 by houyaoxu
[考研] 0856求调剂 +3 刘梦微 2026-03-15 3/150 2026-03-16 10:00 by houyaoxu
信息提示
请填处理意见