我有2个矩阵,我是单独计算的。
其中一个矩阵是Y,它的大小是2x25881,第二个是U,它由2个独立的向量组成,如下所示:
U = [ArrayAT ; ArrayAD]
并且U的大小是2x25881。
现在,我想计算(我的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
有没有更好的方法呢?
谢谢。
答案 0 :(得分:0)
这只是Y和U的元素乘法,然后是所有元素的总和:
Y_Int_Sum = dt * sum(sum(Y_int .* [ArrayAT; ArrayAD]))