积分矩阵乘法

Question

我有 2 个单独计算的矩阵。

其中一个矩阵是 Y，它的大小是 2x25881，第二个是 U，它由 2 个单独的向量组成，如下所示：

U = [ArrayAT ; ArrayAD]

U的大小为2x25881。

现在，我想计算 (my dt = 0.001) 的积分：

我的做法是：

 Y_Int_Sum = 0;

 for i = 1 : length(T)

     Y_Int_Sum = Y_Int_Sum + Y_Int(:,i)'*[ArrayAT(i) ; ArrayAD(i)]*dt;

 end

有没有更好的方法呢？

谢谢。

Answer 1

这只是 Y 和 U 的逐元素乘法，然后是所有元素的总和：

Y_Int_Sum = dt * sum(sum(Y_int .* [ArrayAT; ArrayAD]))