使用eig函数，复模态正交的问题

peter戴

动力方程使用状态空间法求解复模态；(Ar+B)Q=0;A=[ 7.6623e+000            0 -1.1931e+000  8.6491e-001            0  2.0000e+000            0 -2.0000e+000  1.0000e+000            0
            0  7.6623e+000  2.7733e+000            0  8.6491e-001            0  2.0000e+000  3.0000e+000            0  1.0000e+000
-1.1931e+000  2.7733e+000  2.2559e+001 -2.2526e+000  3.3790e+000 -2.0000e+000  3.0000e+000  1.4000e+001 -2.0000e+000  3.0000e+000
  8.6491e-001            0 -2.2526e+000  1.2649e+000            0  1.0000e+000            0 -2.0000e+000  1.0000e+000            0
            0  8.6491e-001  3.3790e+000            0  1.2649e+000            0  1.0000e+000  3.0000e+000            0  1.0000e+000
  2.0000e+000            0 -2.0000e+000  1.0000e+000            0            0            0            0            0            0
            0  2.0000e+000  3.0000e+000            0  1.0000e+000            0            0            0            0            0
-2.0000e+000  3.0000e+000  1.4000e+001 -2.0000e+000  3.0000e+000            0            0            0            0            0
  1.0000e+000            0 -2.0000e+000  1.0000e+000            0            0            0            0            0            0
            0  1.0000e+000  3.0000e+000            0  1.0000e+000            0            0            0            0            0];
B=[ 1.5675e+003            0 -9.0588e+002  6.8377e+002            0            0            0            0            0            0
            0  1.5675e+003  2.2423e+003            0  6.8377e+002            0            0            0            0            0
-9.0588e+002  2.2423e+003  1.3609e+004 -1.7809e+003  2.6713e+003            0            0            0            0            0
  6.8377e+002            0 -1.7809e+003  1.0000e+003            0            0            0            0            0            0
            0  6.8377e+002  2.6713e+003            0  1.0000e+003            0            0            0            0            0
            0            0            0            0            0 -2.0000e+000            0  2.0000e+000 -1.0000e+000            0
            0            0            0            0            0            0 -2.0000e+000 -3.0000e+000            0 -1.0000e+000
            0            0            0            0            0  2.0000e+000 -3.0000e+000 -1.4000e+001  2.0000e+000 -3.0000e+000
            0            0            0            0            0 -1.0000e+000            0  2.0000e+000 -1.0000e+000            0
            0            0            0            0            0            0 -1.0000e+000 -3.0000e+000            0 -1.0000e+000
];
[x,d]=eig(B,-A);
x'*x
模态结果不正交，求解答，恳请大家能给出解释！

peter戴

求解答！。。。。

Spy

x是关于A和B正交，而不是直接正交。

peter戴

谢谢你的回复，按照v'*A*v 或者v'*B*v 也不正交。这个问题已解决了。 涉及到复数，其实是转置的问题，应该用非共轭转置v.'*A*v，而不是共轭转置v'*A*v。

peter戴

Spy 发表于 2014-8-31 22:02
x是关于A和B正交，而不是直接正交。

谢谢你的回复，按照v'*A*v 或者v'*B*v 也不正交。这个问题已解决了。 涉及到复数，其实是转置的问题，应该用非共轭转置v.'*A*v，而不是共轭转置v'*A*v。