复杂的最小二乘问题

计算科学 线性代数 最小二乘
2021-11-27 18:06:32

有一个矩阵ACm×n我解决了以下最小二乘问题

Re(AHA)x=Re(AHb).
如果矩阵A是一个实数矩阵,上式的解可以写成
x=i=1rank(A)uiTbsivi,
在哪里uivi是对应的左右奇异向量和si是一个i奇异值。

我的问题是第一个方程中所述的最小二乘问题的解是否可以用类似的方式写出来A=USVH?

我知道可以拆分矩阵A~=[Re(A); Im(A)]并且等效地解决了实际问题,但是保持在原始矩阵的复数 SVD 内的条件是这里的主要关注点。

谢谢你。

2个回答

我看不出它是如何做到的,但其他人可能会进一步帮助你。我明白了A=USVH,AH=VSHUH,

Re(AHA)x=Re(AHb)=>Re(VSHUHUSVH)x=Re(VSHUHb)=>Re(VSHSVH)x=Re(VSHUHb)

通常,你现在会减少 VSH离开,但因为它们在实部内部,并且因为Re(VSHSVH)Re(VSH)Re(SVH),我不知道如何进行。其他人可以得到一个想法吗?

查看 LAPACK 例程 ZGELLS 或 ZGELLD。他们使用 SVD 解决了 LS 问题。有关例程,请参阅官方 LAPACK 文档