peter戴 发表于 2014-8-23 21:30

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

动力方程使用状态空间法求解复模态;(Ar+B)Q=0;A=[ 7.6623e+000            0 -1.1931e+0008.6491e-001            02.0000e+000            0 -2.0000e+0001.0000e+000            0
            07.6623e+0002.7733e+000            08.6491e-001            02.0000e+0003.0000e+000            01.0000e+000
-1.1931e+0002.7733e+0002.2559e+001 -2.2526e+0003.3790e+000 -2.0000e+0003.0000e+0001.4000e+001 -2.0000e+0003.0000e+000
8.6491e-001            0 -2.2526e+0001.2649e+000            01.0000e+000            0 -2.0000e+0001.0000e+000            0
            08.6491e-0013.3790e+000            01.2649e+000            01.0000e+0003.0000e+000            01.0000e+000
2.0000e+000            0 -2.0000e+0001.0000e+000            0            0            0            0            0            0
            02.0000e+0003.0000e+000            01.0000e+000            0            0            0            0            0
-2.0000e+0003.0000e+0001.4000e+001 -2.0000e+0003.0000e+000            0            0            0            0            0
1.0000e+000            0 -2.0000e+0001.0000e+000            0            0            0            0            0            0
            01.0000e+0003.0000e+000            01.0000e+000            0            0            0            0            0];
B=[ 1.5675e+003            0 -9.0588e+0026.8377e+002            0            0            0            0            0            0
            01.5675e+0032.2423e+003            06.8377e+002            0            0            0            0            0
-9.0588e+0022.2423e+0031.3609e+004 -1.7809e+0032.6713e+003            0            0            0            0            0
6.8377e+002            0 -1.7809e+0031.0000e+003            0            0            0            0            0            0
            06.8377e+0022.6713e+003            01.0000e+003            0            0            0            0            0
            0            0            0            0            0 -2.0000e+000            02.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            02.0000e+000 -3.0000e+000 -1.4000e+0012.0000e+000 -3.0000e+000
            0            0            0            0            0 -1.0000e+000            02.0000e+000 -1.0000e+000            0
            0            0            0            0            0            0 -1.0000e+000 -3.0000e+000            0 -1.0000e+000
];
=eig(B,-A);
x'*x
模态结果不正交,求解答,恳请大家能给出解释!

peter戴 发表于 2014-8-25 18:00

求解答!。。。。

Spy 发表于 2014-8-31 22:02

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

peter戴 发表于 2014-9-1 16:48

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

peter戴 发表于 2014-9-1 16:49

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


谢谢你的回复,按照v'*A*v 或者v'*B*v 也不正交。这个问题已解决了。 涉及到复数,其实是转置的问题,应该用非共轭转置v.'*A*v,而不是共轭转置v'*A*v。
页: [1]
查看完整版本: 使用eig函数,复模态正交的问题