1

我有一个大矩阵(时间 x 频率),我想部分减少它。我想将每 1000 行(时间样本)相加,保持频率信息,这是一种分割。

有没有什么方法可以在 MATLAB 中没有任何循环?

一个较小的例子:

M=[1 2 3; 2 3 4; 5 8 7; 5 6 7; 1 2 3; 1 2 4]; 

我想将每 2 行相加,所以我得到:

[3 5 7; 10 14 14; 2 4 7]
4

1 回答 1

2

假设您有一个包含 N 行和 M 列的矩阵,并且您想将每 R 行相加(其中 N 可以被 R 整除),

>> mat = [1 2 3; 2 3 4; 5 8 7; 5 6 7; 1 2 3; 1 2 4]
mat =
     1     2     3
     2     3     4
     5     8     7
     5     6     7
     1     2     3
     1     2     4

>> [N, M] = size(mat); %=> [6, 3]
>> R = 2;

以下将允许您对 R 行的组求和:

>> res = reshape(mat, R, [])
res = 
     1     5     1     2     8     2     3     7     3
     2     5     1     3     6     2     4     7     4

>> res = sum(res) 
res =
     3    10     2     5    14     4     7    14     7

>> res = reshape(res, [], M)
res =
     3     5     7
    10    14    14
     2     4     7

您还可以在一行中完成所有操作:

>> reshape(sum(reshape(mat, R, [])), [], M)
ans =
     3     5     7
    10    14    14
     2     4     7
于 2012-08-03T14:55:30.403 回答