复数矩阵乘法本征与 Matlab

Question

谁能给我解释一下为什么结果不同。

C++ 代码：

MatrixXcd testTest;
testTest.resize(3,3);
testTest.real()(0,0) = 1;
testTest.real()(0,1) = 2;
testTest.real()(0,2) = 3;
testTest.real()(1,0) = 1;
testTest.real()(1,1) = 2;
testTest.real()(1,2) = 3;
testTest.real()(2,0) = 1;
testTest.real()(2,1) = 2;
testTest.real()(2,2) = 3;

testTest.imag()(0,0) = 1;
testTest.imag()(0,1) = 2;
testTest.imag()(0,2) = 3;
testTest.imag()(1,0) = 1;
testTest.imag()(1,1) = 2;
testTest.imag()(1,2) = 3;
testTest.imag()(2,0) = 1;
testTest.imag()(2,1) = 2;
testTest.imag()(2,2) = 3;

cout<< endl << testTest << endl;
cout<< endl << testTest.transpose() << endl;
cout<< endl << testTest*testTest.transpose() << endl;
cout<< endl << testTest << endl;

来自 C++ 的结果：

(1,1) (2,2) (3,3)
(1,1) (2,2) (3,3)
(1,1) (2,2) (3,3)

(1,1) (1,1) (1,1)
(2,2) (2,2) (2,2)
(3,3) (3,3) (3,3)

(0,28) (0,28) (0,28)
(0,28) (0,28) (0,28)
(0,28) (0,28) (0,28)

(1,1) (2,2) (3,3)
(1,1) (2,2) (3,3)
(1,1) (2,2) (3,3)

用 Matlab 写的同样的东西：

testTest = [ complex(1,1) complex(2,2) complex(3,3); 
             complex(1,1) complex(2,2) complex(3,3); 
             complex(1,1) complex(2,2) complex(3,3)];

testTest
testTest'
testTest*testTest'
testTest

Matlab 结果：

testTest =

1.0000 + 1.0000i   2.0000 + 2.0000i   3.0000 + 3.0000i
1.0000 + 1.0000i   2.0000 + 2.0000i   3.0000 + 3.0000i
1.0000 + 1.0000i   2.0000 + 2.0000i   3.0000 + 3.0000i


ans =

1.0000 - 1.0000i   1.0000 - 1.0000i   1.0000 - 1.0000i
2.0000 - 2.0000i   2.0000 - 2.0000i   2.0000 - 2.0000i
3.0000 - 3.0000i   3.0000 - 3.0000i   3.0000 - 3.0000i


ans =

28    28    28
28    28    28
28    28    28


testTest =

1.0000 + 1.0000i   2.0000 + 2.0000i   3.0000 + 3.0000i
1.0000 + 1.0000i   2.0000 + 2.0000i   3.0000 + 3.0000i
1.0000 + 1.0000i   2.0000 + 2.0000i   3.0000 + 3.0000i

C 中 testTest * testTest' 的乘法 returns returns 实数部分 0 和图像部分 28 的复数。Matlab returns 只需将值 28 加倍。

Answer 1

Matlab 中的

' 进行转置并采用复共轭 (http://uk.mathworks.com/help/matlab/ref/ctranspose.html)。如果您只想进行转置，请使用 .'（前面有一个点）。

因此，如果您将 MATLAB 测试更改为

testTest*testTest.'

结果应该是一样的。

如果你想在 eigen 中进行复杂的转置，那么你可以去 matrix.adjoint()（或 matrix.conjugate().transpose()）

复数矩阵乘法本征与 Matlab

Complex Number Matrix multiplication Eigen vs Matlab

c++

matlab

eigen