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

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

使用第一种方法时，我得到的ii2在值和向量大小方面都非常不同（由于向量大小太大，所以无法发布结果）。

有人可以在这里帮助绝望的新手吗？我只需要将循环部分重新构造为更好，更简洁的版本。

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

Answer 1

您可以使用hist3：

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

使用矢量化或其他方法重新格式化循环-八度

1 个答案: