如何normalize矩阵# Computation - 科学计算
f*r
1 楼
假设我有两个同样大小的square矩阵，A,B.我希望可以找到一个变换矩阵X，这样如果我做一个
变换
A->A*X,B->B*X,之后，就有A*A'+B*B'=I .
如何找到这样的变化矩阵呢？
忘大家不吝赐教，我在这里先谢谢了。
变换
A->A*X,B->B*X,之后，就有A*A'+B*B'=I .
如何找到这样的变化矩阵呢？
忘大家不吝赐教，我在这里先谢谢了。
r*y
2 楼
Are A and B square matrices? If not, may not have unique solutions.
【在 f********r 的大作中提到】
: 假设我有两个同样大小的square矩阵，A,B.我希望可以找到一个变换矩阵X，这样如果我做一个
: 变换
: A->A*X,B->B*X,之后，就有A*A'+B*B'=I .
: 如何找到这样的变化矩阵呢？
: 忘大家不吝赐教，我在这里先谢谢了。
f*r
3 楼
A和B都是方阵,只要有解就行，我并不需要unique的解
【在 r****y 的大作中提到】
:
: Are A and B square matrices? If not, may not have unique solutions.
【在 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.
解法如下（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.
f*t
5 楼
such X may not exist at all.
果我做一个
【在 f********r 的大作中提到】
: 假设我有两个同样大小的square矩阵，A,B.我希望可以找到一个变换矩阵X，这样如果我做一个
: 变换
: A->A*X,B->B*X,之后，就有A*A'+B*B'=I .
: 如何找到这样的变化矩阵呢？
: 忘大家不吝赐教，我在这里先谢谢了。
果我做一个
【在 f********r 的大作中提到】
: 假设我有两个同样大小的square矩阵，A,B.我希望可以找到一个变换矩阵X，这样如果我做一个
: 变换
: A->A*X,B->B*X,之后，就有A*A'+B*B'=I .
: 如何找到这样的变化矩阵呢？
: 忘大家不吝赐教，我在这里先谢谢了。
相关阅读
help: 怎么用mathematica 输出坐标轴标记?matlab question:dual variable求助: 解有mixed derivative的PDE时遇到麻烦问个蠢蠢的问题thank u, achillis, thinkme, hanzo....Re: 有个问题drive me crazyHelp change the code of optimization[转载] FFT里面的"码位倒置"用英文怎么说啊?[转载] 请maple 或者 mathematica 的高手帮个忙！有个问题drive me crazy为什么两台机器的计算结果不一样？Who has ever used the Tsyganenko code?2篇基本的matlab有限元编程的paper3-D dipole-dipole interaction[转载] 问个误差估计的问题向超星北美用户求助[转载] S-C-C bend force constants你的程序写法有问题 Re: 请教大家c＋＋ array一个用mathematica 求微分的问题做分类，要找数据