Question

我正在尝试将代码从matlab移植到C.为了将此行转换为C：

A = E*[SOLS' ; ones(1,10 ) ];

其中，

>>size(SOLS)

ans =

10     3

和

>> size(E)

ans =

 9     4

SOLS是一个复杂的单一矩阵，E是一个真正的单矩阵，A应该是一个大小为9x10的复杂单矩阵。

我用

替换了A = E*[SOLS' ; ones(1,10 ) ];

for i=1:9
  for j=1:10
    A1(i,j)=E(i,1)*SOLS(j,1))+E(i,2)*SOLS(j,2))+E(i,3)*SOLS(j,3))+E(i,4);
   end
end

复杂的结果矩阵元素具有与A相同的实部，但是具有不同的虚部。

>> real(A)=real(A1)
imag(A) and `imag(A1)` are different.

造成这种差异的原因是什么？如何将matlab命令正确转换为C？

以下是矩阵示例：

    E =

    0.2248         0         0         0
   -0.4487   -0.1632   -0.1955    0.6355
    0.4379   -0.0651   -0.1032   -0.0754
   -0.4008    0.3513    0.2707   -0.5936
   -0.2294   -0.7853   -0.3290   -0.4648
    0.0385    0.2623   -0.6363   -0.0978
   -0.5716    0.0851    0.0943    0.0587
    0.1160   -0.3911    0.5964    0.0947
    0.0363   -0.0039   -0.0092   -0.0018

和

    SOLS =

   1.0e+02 *

  -0.2410 + 0.0000i   2.3741 + 0.0000i  -0.0646 + 0.0000i
   0.0000 + 0.0000i   0.0000 + 0.0000i  -0.0113 - 0.0046i
   0.0000 + 0.0000i   0.0000 + 0.0000i  -0.0113 + 0.0046i
  -0.0028 + 0.0000i  -0.0114 + 0.0000i  -0.0038 + 0.0000i
   0.0000 + 0.0000i   0.0000 + 0.0000i  -0.0024 - 0.0043i
   0.0000 + 0.0000i   0.0000 + 0.0000i  -0.0024 + 0.0043i
   0.0000 + 0.0000i   0.0000 + 0.0000i  -0.0007 - 0.0191i
   0.0000 + 0.0000i   0.0000 + 0.0000i  -0.0007 + 0.0191i
  -0.0080 + 0.0000i   0.0064 + 0.0000i   0.0108 + 0.0000i
  -0.7289 + 0.0000i   4.9347 + 0.0000i   0.3841 + 0.0000i

Answer 1

在MATLAB中，SOLS'执行复共轭转置操作，即元素{i，j}变为元素{j，i}，其值转换为a + 1i*b --> a -1i*b。要保留复杂值的阶段，请使用SOLS.'，如下所示：

A = E*[SOLS.' ; ones(1,10 ) ];

此外，这是你想要执行循环的方式（当然要翻译成适当的C）：

for i=1:size(E,1)
  for j=1:size(SOLS,1)
        A1(i,j)=0;
        for k = 1:size(SOLS,2)
            A1(i,j)= A1(i,j) + E(i,k)*SOLS(j,k);
        end
        A1(i,j)= A1(i,j) + E(i,k+1);
   end
end

然后

A1 - A

ans =

   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0
   0   0   0   0   0   0   0   0   0   0

在Matlab中乘以复杂矩阵

1 个答案: