24小时热门版块排行榜    

查看: 2590  |  回复: 10

anionxt

铁杆木虫 (著名写手)

[求助] 请教 AXA^T = B 这类方程怎么数值求解呢? 已有4人参与

请教各位朋友:

AXA^T = B

这种方程怎么样数值方法求解呢?
A^T 表示转置, A 和 B 都是实数矩阵, X为未知数矩阵。
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

bluesine

至尊木虫 (职业作家)

科苑小木虫

【答案】应助回帖

★ ★ ★ ★ ★
感谢参与,应助指数 +1
anionxt: 金币+5, ★★★很有帮助 2016-05-18 08:36:00
AXA^T = B
得到A^T AXA^T A = A^TBA
A^T AX  =  A^TBA*inv(A^TA)

X = inv(A^T A) A^T BA*inv(A^TA)
inv是求逆
板凳要做十年冷文章不发一个字
4楼2016-05-17 08:46:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

lukebest

禁虫 (正式写手)

本帖内容被屏蔽

2楼2016-05-16 18:54:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lukebest

禁虫 (正式写手)

★ ★ ★ ★ ★
感谢参与,应助指数 +1
anionxt: 金币+5, ★★★很有帮助, 谢谢! 2016-05-22 10:45:42
本帖内容被屏蔽

3楼2016-05-16 18:56:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cfc823

铁杆木虫 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
引用回帖:
4楼: Originally posted by bluesine at 2016-05-17 08:46:51
AXA^T = B
得到A^T AXA^T A = A^TBA
A^T AX  =  A^TBA*inv(A^TA)

X = inv(A^T A) A^T BA*inv(A^TA)
inv是求逆

用求逆来解答,前提需要确保A,B是满秩矩阵吧

发自小木虫Android客户端
5楼2016-05-17 09:08:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

bluesine

至尊木虫 (职业作家)

科苑小木虫

引用回帖:
5楼: Originally posted by cfc823 at 2016-05-17 09:08:14
用求逆来解答,前提需要确保A,B是满秩矩阵吧
...

B不需要满秩,A行满秩就够了。。。
板凳要做十年冷文章不发一个字
6楼2016-05-17 09:13:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
祝福
7楼2016-05-17 21:30:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

anionxt

铁杆木虫 (著名写手)

引用回帖:
4楼: Originally posted by bluesine at 2016-05-17 08:46:51
AXA^T = B
得到A^T AXA^T A = A^TBA
A^T AX  =  A^TBA*inv(A^TA)

X = inv(A^T A) A^T BA*inv(A^TA)
inv是求逆

谢谢你,但是矩阵的维度都是200x200左右,求逆可能速度较慢,有什么迭代的方法来求解吗

发自小木虫Android客户端
8楼2016-05-18 08:38:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
祝福
9楼2016-05-18 10:26:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

rock_faith

金虫 (小有名气)

【答案】应助回帖

引用回帖:
8楼: Originally posted by anionxt at 2016-05-18 08:38:23
谢谢你,但是矩阵的维度都是200x200左右,求逆可能速度较慢,有什么迭代的方法来求解吗
...

象楼上说的那样,将矩阵方程两边,按行拉直后,变成Ax =b的形式,然后就可以用高斯消元法,即Ax=LUx=b (三角分解),不用求逆,效率很高。

发自小木虫Android客户端
太初有道,道与 神同在,道就是 神(约1:1)
10楼2016-05-19 22:26:18
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 anionxt 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-04 5/250 2026-08-04 13:31 by jzdxtiUmJ6PL
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +4 Ocx3OKOzGU4q 2026-08-03 7/350 2026-08-04 13:22 by jzdxtiUmJ6PL
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 Ocx3OKOzGU4q 2026-08-03 8/400 2026-08-04 13:15 by jzdxtiUmJ6PL
[教师之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 jRl3mE6ddZGq 2026-08-03 11/550 2026-08-04 12:58 by jzdxtiUmJ6PL
[基金申请] 娱乐 +4 Tide man 2026-08-03 4/200 2026-08-04 11:51 by wgch518
[基金申请] 纯娱乐,不喜欢勿喷 +5 Tide man 2026-08-04 8/400 2026-08-04 11:43 by everxc
[基金申请] 面上提前没消息,有中的吗 +14 archvillain 2026-08-02 17/850 2026-08-04 09:46 by archvillain
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 f1dmPx051StA 2026-08-03 7/350 2026-08-04 09:31 by v1TJC4qJauAB
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 09:17 by v1TJC4qJauAB
[基金申请] 系统今天又提示维护了,估计离放榜不远了 +12 winnerche 2026-07-29 16/800 2026-08-04 08:18 by Equinoxhua
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 9/450 2026-08-04 07:37 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +4 cu9Nq1xK233Z 2026-08-03 13/650 2026-08-04 07:30 by k7OM8YghWbkC
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 jRl3mE6ddZGq 2026-08-03 6/300 2026-08-04 04:47 by k7OM8YghWbkC
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 cu9Nq1xK233Z 2026-08-03 7/350 2026-08-04 04:46 by k7OM8YghWbkC
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
[高分子] UV压敏胶开发 +3 ichall 2026-07-30 5/250 2026-08-03 14:40 by Sunrisepay
[基金申请] 面上再次挂了,太难了,躺也躺不了,倦也卷不过,小学校之殇! +19 低垂的野花 2026-07-31 25/1250 2026-08-03 09:25 by gy116024
[基金申请] 2026年国自然面上资助率 +16 布布和一二 2026-07-30 22/1100 2026-08-03 09:20 by gy116024
[高分子] HXDI做水性聚氨酯乳液,是不是特别容易出渣 15+3 yuyusuv 2026-07-29 3/150 2026-07-31 09:21 by huizingga
信息提示
请填处理意见