用矢量化或其他方法重新制定 for 循环 - octave

Question

是否有任何方法可以向量化（或重新制定）此代码中循环的每个主体：

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 是一个向量，其中包含以米为单位的水柱，而 dates 是一个向量（具有两位数的整数），其中包含对相应水柱进行测量的日期。输入文件可以在这里找到：https://drive.google.com/open?id=1EomLGYleaNtiGG2iV_9LRt425blxdIsm

至于输出，我想要这样的 ii:

与第一种方法相反，我得到了 ii2，它在值和向量大小方面有很大不同（我不能 post 结果，因为向量大小太大）。

有人可以帮助绝望的新手吗？我只需要将循环部分重新表述为更好、更简洁的版本。

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

Answer 1

您可以使用 hist3:

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

用矢量化或其他方法重新制定 for 循环 - octave

reformulating for loop with vectorization or other approach - octave

for-loop

vectorization

octave

bsxfun