如果我有矩阵
A=rand(50,4);
和范围值矩阵,例如
range=[1 10 15 30 45 50];
,我想计算范围矩阵中所有A列中的最小值及其索引。例如
for i=1 :numel(range)-1
[value index]= find(min(A(range(i):range(i+1),[],1)
if i>1
index=index+ range(i)+1; % indx# was calculated onlyin a
% range of array and not for
%whole array.Updated here
end
B_ind(i,:)=index;
B_val(i,:)=value;
end
如何在不使用循环方法的情况下为min(A)
获得range(i:i+1)
的索引?
简单地说,矩阵“ B”应类似于
B(1,column 1:4)=index of min(A (1:10)) in every column
B(2,column 1:4)=index of min(A (11:15)) in every column
B(3,column 1:4)=index of min(A (16:30)) in every column
以此类推...
我的问题是如何在没有“ for循环”的情况下做到这一点?
答案 0 :(得分:2)
一种可能性是:
% generate example data set
A=rand(50,4);
range=[1 10 15 30 45 50];
% generate indexes of interest
tmp_idx= arrayfun(@colon,range(1:end-1),range(2:end),'un',0);
% calculate the min function over the indexes
[min_value,min_idx]=cellfun(@(x) min(A(x,:)), tmp_idx, 'UniformOutput', false);
% get matrix B (with the offset of ranges)
B=cell2mat(min_idx')+(range(1:end-1)+1)';
B(1,:)=B(1,:)-2;