我有一个大矩阵(时间 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]
假设您有一个包含 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