用矢量化或其他方法重新制定循环-倍频程

有没有任何方法可以向量化(或重新表述)这个代码中循环的每个主体：

col=load('col-deau'); %load data
h=col(:,8); % corresponding water column
dates=col(:,3); % and its dates
%removing out-of-bound data
days=days(h~=9999.000);
h=h(h~=9999.000);
dates=sort(dates(h~=9999.000));
[k,hcat]=hist(h,nbin); %making classes (k) and boundaries of classes (hcat) of water column automatically
dcat=1:15; % make boundaries for dates
for k=1:length(dcat)-1 % Loop for each date class
ii=find(dates>=dcat(k)&dates<dcat(k+1));% Counting dates corresponding to the boundaries of each date class
for j=1:length(hcat)-1                                % Loop over each class of water column
ij=find(h>=hcat(j)&h<hcat(j+1)); % Count water column corresponding to the boundaries of each water column class
obs(k,j)=length(intersect(ii,ij));               % Find the size of each intersecting matrix
end
end

例如，我尝试过使用矢量化来更改这一部分：

for k=1:length(dcat)-1        
ii=find(dates>=dcat(k)&dates<dcat(k+1))  
endfor

这个：

nk=1:length(dcat)-1;
ii2=find(dates>=dcat(nk)&dates<dcat(nk+1));

以及使用bsxfun:

ii2=find(bsxfun(@and,bsxfun(@ge,dates,nk),bsxfun(@lt,dates,nk+1)));

但无济于事。这两种方法都产生相同的输出，并且与使用for循环的方法(在元素和向量大小方面)不一致。

作为信息，h是一个包含以米为单位的水柱的矢量，date是一个矢量(带两位数字的整数)，它包含对相应水柱进行测量的日期。输入文件可在此处找到：https://drive.google.com/open?id=1EomLGYleaNtiGG2iV_9LRt425blxdIsm

至于输出，我想有这样的ii：

相反，使用第一种方法，我得到了ii2，它在值和向量大小方面非常不同(我无法发布结果，因为向量大小太大)。

这里有人能帮助一个绝望的新手吗？我只需要把循环部分重新表述成一个更好、更简洁的版本。

如果需要添加更多细节，请随时询问我。

您可以使用hist3:

pkg load statistics
[obs, ~] = hist3([dates(:) h(:)] ,'Edges', {dcat,hcat});

相关内容

最新更新

热门标签：