【在 r****y 的大作中提到】 : : Are A and B square matrices? If not, may not have unique solutions.
f*r
4 楼
问题解决了,多谢rossby的提醒 解法如下(matlab) function x=normalize(A,B) x0=eye(size(A,1)); function F=myfun(x) F=A*x*x'*A'+B*x*x'*B'-eye(size(A,1)); end x=fsolve(@myfun,x0); end 只是似乎效率不是特别高,算一个4阶的矩阵有时需要0.5秒 有没有直接求解的方法呢? 还有,似乎不是所有的矩阵A,B都有解,有时解又依赖于初始的x0, 是不是不是所有的矩阵A,B都存在解呢?
【在 r****y 的大作中提到】 : : Are A and B square matrices? If not, may not have unique solutions.