如何计算Matlab中所有迭代的for循环输出的总和？

Question

我有两个不同的行向量 P 和 L。该算法对于每次迭代都有不同的 P 和 L，我的代码应该计算两个行向量在同一位置有多少对应的行向量，例如如果 P 和 L 是：

P=[1,2,2,1,1]
L=[2,2,1,1,1]

答案应该是 2。例如，如果我有另一个 P 和 L，答案是 3，我应该能够将 2 和 3 相加得到 5。

我的代码计算每次迭代的正确答案，但我不知道如何对这些答案求和以获得向量 P 和 L 在所有迭代中对应的总数。

这是我的代码：

     i = 1:numel(P); 
     j = 1:numel(L);
     Valuecompared=(P(i) == 1) & (L(j)==1);

     a=0;

   for k = 1:numel(Valuecompared)

   if Valuecompared(k) == 1
    a=a+1;
   end

   end

我在循环内部和外部尝试了sum(a), cumsum(a)，索引k 和没有，但它只是给了我相同的答案“a”，而没有添加先前计算的“a”答案。可能是因为它被覆盖并且每次都没有保存？

此外，我似乎无法在算法运行时存储“a”的值，我可以显示它，但如果我有 100 次不方便的迭代，我如何计算所有迭代的总和答案并存储像Atotal这样的一个变量的总和？

Answer 1

您必须从每个循环中存储a。不幸的是，您没有显示完整的循环，所以我不知道您正在迭代哪个变量。但是让我们假设它是q。然后你必须在循环外定义a，例如与

a = zeros(1,niter); % niter = number of times you run the loop

然后你必须用类似的东西替换最后几行

if Valuecompared(k) == 1
    a(q)=a(q)+1;
end

循环结束后，您可以使用sum(a)

或者，您可以将a=0;移出循环，然后您甚至不必再取总和。

有点离题，但您可以将目前显示的所有代码替换为

a = sum(P==1 & L==1)

做同样的事情。