我有一个矩阵A
1 1 0 0
0 1 0 0
1 0 0 1
0 0 1 0
0 0 0 0
0 1 1 1
1 1 0 0
1 0 0 0
0 0 0 1
我希望根据用户的输入说明d = [1 2 3].
for i=2:length(d)
d(i) = d(i) + d(i-1); % d = [1 3 6]
end
这会给d = [1 (1+2) (1+2+3)] = d[1 3 6]
。此矩阵中有9行,计算ceil
,[(1/6)*9]
和[(3/6)*9]
的{{1}}。因此,这给出了[(6/6)*9]
第一次拆分是前两行,第二次拆分是下一行(5-2 = 3)3行,第三次拆分是(9-5 = 4)4行。
输出应该是:
分手是:第一次分手 - >
[2 5 9].
第二次分手 - >
1 1 0 0 % first 2 rows in matrix A
0 1 0 0
第3次分手 - >
1 0 0 1 % next 3 rows
0 0 1 0
0 0 0 0
答案 0 :(得分:5)
您可以将mat2cell
与输入d = [1 2 3]
一起使用,将最终拆分存储在单独的单元格数组中
B = mat2cell(A, d+1, size(A,2));
或者,使其适应您对分割行大小的计算:
d = [1 2 3];
c = cumsum(d); % [1, 3, 6]
s = ceil(size(A,1)*c/c(end)); % [2, 5, 9]
n = [s(1) diff(s)]; % [2, 3, 4]
B = mat2cell(A, n, size(A,2));
要显示拆分,您可以添加类似于以下命令:
cellfun(@disp, B)