1

只是一个简单的嵌套问题：

我有一个<100x100 double>矩阵mat_B，我cumsum。从生成的矩阵mat_A中，我只需要最后一行vec_C，我需要cumsum再次这样做。我的代码如下所示：

mat_A = cumsum(mat_B);
vec_C = cumsum(mat_A(end,:));

我的问题是，是否可以将所有这些都放在一行代码中。我知道它cumsum(mat_B)返回一个矩阵，但如果我把(end, :)表达式放在后面，它就行不通了。

我知道这听起来很傻，但我想知道在这种情况下嵌套是如何工作的。

4

1 回答 1

6

您可以跳过第一个 cumsum 而只使用sum，因为最后一行cumsum等于 sum 的结果：

>> mat_B=rand(5); 
>> cumsum(mat_B)

ans =

    0.2517    0.4522    0.8838    0.3751    0.2527
    0.6847    0.7778    1.3412    0.7487    0.8376
    1.5270    1.1579    2.1404    1.2327    1.3613
    1.7115    2.0444    2.2745    2.2021    1.5247
    2.2197    2.8056    2.3398    2.5442    2.0111

>> sum(mat_B)

ans =

    2.2197    2.8056    2.3398    2.5442    2.0111

所以

vec_C = cumsum(sum(mat_B));

应该做你想做的。

于 2013-07-22T16:56:19.240 回答