MATLAB: How to loop the date vector, that the restrictions are the months,years and hours?

datevectorfind positionsloopMATLABspecific positions in a loop

 %for k=xa:1:y2 xa is the first year of the data that i use a function to discover it and y2 is the last year
     ka = find(IntervaloHorarioDateVec(:,1) == 2015);     
     kb = find(IntervaloHorarioDateVec(:,1) == 2016);
     kc = find(IntervaloHorarioDateVec(:,1) == 2017);
     kd = find(IntervaloHorarioDateVec(:,1) == 2018);
     %1-12
        %for m = 1:12
            ke = find(IntervaloHorarioDateVec(:,2) == 1);
            kf = find(IntervaloHorarioDateVec(:,2) == 2);
            kg = find(IntervaloHorarioDateVec(:,2) == 3);
            kh = find(IntervaloHorarioDateVec(:,2) == 4);
            ki = find(IntervaloHorarioDateVec(:,2) == 5);
            kj = find(IntervaloHorarioDateVec(:,2) == 6);
            kk = find(IntervaloHorarioDateVec(:,2) == 7);
            kl = find(IntervaloHorarioDateVec(:,2) == 8);
            km = find(IntervaloHorarioDateVec(:,2) == 9);
            kn = find(IntervaloHorarioDateVec(:,2) == 10);
            ko = find(IntervaloHorarioDateVec(:,2) == 11);
            kp = find(IntervaloHorarioDateVec(:,2) == 12);            
            kq = ismember(ka,ke,'rows');
            kr = find(kq);
            ks = ismember(ka,kf,'rows');
            kt = find(ks);
            ku = ismember(ka,kg,'rows');
            kv = find(kv);
            kw = ismember(ka,kh,'rows');
            kx = find(kw);
            ky = ismember(ka,ki,'rows');
            kz = find(ky);
            kaa = ismember(ka,kj,'rows');
            kab = find(kaa);
            kac = ismember(ka,kk,'rows');
            kad = find(kac);
            kae = ismember(ka,kl,'rows');
            kaf = find(kae);
            kag = ismember(ka,km,'rows');
            kah = find(kag);
            kai = ismember(ka,kn,'rows');
            kaj = find(kai);
            kak = ismember(ka,ko,'rows');
            kal = find(kal);
            kam = ismember(ka,kp,'rows');
            kan = find(kan);
            ... doing this until 2018
                % and for the hours of the day...
                k1 = find(IntervaloHorarioDateVec(:,2) == 0);            
                k2 = find(IntervaloHorarioDateVec(:,2) == 1);
                k3 = find(IntervaloHorarioDateVec(:,2) == 2);
                k4 = find(IntervaloHorarioDateVec(:,2) == 3);
                k5 = find(IntervaloHorarioDateVec(:,2) == 4);
                k6 = find(IntervaloHorarioDateVec(:,2) == 5);
                k7 = find(IntervaloHorarioDateVec(:,2) == 6);
                k8 = find(IntervaloHorarioDateVec(:,2) == 7);
                k9 = find(IntervaloHorarioDateVec(:,2) == 8);
                k10 = find(IntervaloHorarioDateVec(:,2) == 9);
                k11 = find(IntervaloHorarioDateVec(:,2) == 10);
                k12 = find(IntervaloHorarioDateVec(:,2) == 11);
                k13 = find(IntervaloHorarioDateVec(:,2) == 12);
                k14 = find(IntervaloHorarioDateVec(:,2) == 13);
                k15 = find(IntervaloHorarioDateVec(:,2) == 14);
                k16 = find(IntervaloHorarioDateVec(:,2) == 15);
                k17 = find(IntervaloHorarioDateVec(:,2) == 16);
                k18 = find(IntervaloHorarioDateVec(:,2) == 17);
                k19 = find(IntervaloHorarioDateVec(:,2) == 18);
                k20 = find(IntervaloHorarioDateVec(:,2) == 19);
                k21 = find(IntervaloHorarioDateVec(:,2) == 20);
                k22 = find(IntervaloHorarioDateVec(:,2) == 21);
                k23 = find(IntervaloHorarioDateVec(:,2) == 22);
                k24 = find(IntervaloHorarioDateVec(:,2) == 23);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                Lia = ismember(k1,kr);
                Idrow = find(Lia);
                ...to 00:00 until 24:00
                    ... going thru all the years 2015-2018,months 1-12, hours 0-23h
                    %So i want to create a loop, that I can find 00 hours in the year, let's say 2015 of my data, or 2017, or both,
                    %and in that year i want to find the month, let's say January(01<-number of the month), or March(03),
                    %and in that YEAR and MONTH, i want to find the position of the 00:00 Hour or 12:00, in a data that is hourly
                    %so i want to find all the data that are in 00:00hour of all the first month of the first year, the 01:00 hour of the first month
                    %and so on til 23:00hour, the same with february to december until the end of y2(the last year).

Best Answer

Logical indexing:

for k = xa:y2
     ka = (IntervaloHorarioDateVec(:,1) == k);     % No FIND
     for m = 1:12
       kb = (IntervaloHorarioDateVec(:,2) == m);
       for n = 0:23	
         ke = (IntervaloHorarioDateVec(:,4) == n);
         match = find(ka & kb & ke);
         ...
       end
    end
end

Or:

for k = xa:y2
   for m = 1:12
      for n = 0:23	
         ke = (IntervaloHorarioDateVec(:,4) == n);
         match = ismember(IntervaloHorarioDateVec(:, [1,2,4]), [ka, kb, ke], 'rows');
         ...
       end
    end
end

Related Solutions

MATLAB: Nested loop will not enter second iteration

I get 3 tierations for i. Here is my code:

% here is some sample data:
CO2_FFD = [0  2000  1  1  0  0  0  340  ;
0  2000  1  2  0  0  0  380  ;
0  2000  2  3  0  0  0  280  ;
0  2000  2  4  0  0  0  450  ;
0  2000  3  5  0  0  0  330  ;
0  2000  3  6  0  0  0  450  ;
0  2000  4  7  0  0  0  380  ;
0  2000  4  8  0  0  0  290  ;
0  2000  5  9  0  0  0  340  ;
0  2000  5  10  0  0  0  380  ;
0  2000  6  11  0  0  0  280  ;
0  2000  6  12  0  0  0  450  ;
0  2000  7  13  0  0  0  330  ;
0  2000  7  14  0  0  0  450  ;
0  2000  8  15  0  0  0  380  ;
0  2000  8  16  0  0  0  290  ;
0  2000  9  17  0  0  0  340  ;
0  2000  9  18  0  0  0  380  ;
0  2000  10  19  0  0  0  280  ;
0  2000  10  20  0  0  0  450  ;
0  2000  11  21  0  0  0  330  ;
0  2000  11  22  0  0  0  450  ;
0  2000  12  23  0  0  0  380  ;
0  2000  12  24  0  0  0  290  ;
0  2001  1  25  0  0  0  340  ;
0  2001  1  26  0  0  0  380  ;
0  2001  2  27  0  0  0  280  ;
0  2001  2  28  0  0  0  450  ;
0  2001  3  29  0  0  0  330  ;
0  2001  3  30  0  0  0  450  ;
0  2001  4  1  0  0  0  380  ;
0  2001  4  2  0  0  0  290  ;
0  2001  5  3  0  0  0  340  ;
0  2001  5  4  0  0  0  380  ;
0  2001  6  5  0  0  0  280  ;
0  2001  6  6  0  0  0  450  ;
0  2001  7  7  0  0  0  330  ;
0  2001  7  8  0  0  0  450  ;
0  2001  8  9  0  0  0  380  ;
0  2001  8  10  0  0  0  290  ;
0  2001  9  11  0  0  0  340  ;
0  2001  9  12  0  0  0  380  ;
0  2001  10  13  0  0  0  280  ;
0  2001  10  14  0  0  0  450  ;
0  2001  11  15  0  0  0  330  ;
0  2001  11  16  0  0  0  450  ;
0  2001  12  17  0  0  0  380  ;
0  2001  12  18  0  0  0  290  ;
0  2002  1  19  0  0  0  340  ;
0  2002  1  20  0  0  0  380  ;
0  2002  2  21  0  0  0  280  ;
0  2002  2  22  0  0  0  450  ;
0  2002  3  23  0  0  0  330  ;
0  2002  3  24  0  0  0  450  ;
0  2002  4  25  0  0  0  380  ;
0  2002  4  26  0  0  0  290  ;
0  2002  5  27  0  0  0  340  ;
0  2002  5  28  0  0  0  380  ;
0  2002  6  29  0  0  0  280  ;
0  2002  6  30  0  0  0  450  ;
0  2002  7  31  0  0  0  330  ;
0  2002  7  32  0  0  0  450  ;
0  2002  8  33  0  0  0  380  ;
0  2002  8  34  0  0  0  290  ;
0  2002  9  35  0  0  0  340  ;
0  2002  9  36  0  0  0  380  ;
0  2002  10  37  0  0  0  280  ;
0  2002  10  38  0  0  0  450  ;
0  2002  11  39  0  0  0  330  ;
0  2002  11  40  0  0  0  450  ;
0  2002  12  41  0  0  0  380  ;
0  2002  12  42  0  0  0  290  ]
[M,ia,ic] = unique(CO2_FFD(:,2)); %find unique years, M = unique year, ia = row with unique year
m=[1;2;3;4;5;6;7;8;9;10;11;12];%set an array for 12 months
CO2_YM_averages = zeros(length(M)*12,2); %array for averages & std's, size of year*12months
for i=1:length(M)  
  fprintf('i = %d\n', i);
  yrows = find(M(i)==CO2_FFD(:,2)); %find rows where year is equal to M(i)  
  for j=1:length(m)    
  fprintf('    j = %d\n', j);
    mrows=find(CO2_FFD(yrows,3)== m(j));%find rows where month = m(i)    
    CO2_YM_averages(j,1) = mean(CO2_FFD(mrows,8)); %average data for each month    
    CO2_YM_averages(j,2) = std(CO2_FFD(mrows,8)); %std data for each month    
  end  
end

and here is the results:

i = 1
    j = 1
    j = 2
    j = 3
    j = 4
    j = 5
    j = 6
    j = 7
    j = 8
    j = 9
    j = 10
    j = 11
    j = 12
i = 2
    j = 1
    j = 2
    j = 3
    j = 4
    j = 5
    j = 6
    j = 7
    j = 8
    j = 9
    j = 10
    j = 11
    j = 12
i = 3
    j = 1
    j = 2
    j = 3
    j = 4
    j = 5
    j = 6
    j = 7
    j = 8
    j = 9
    j = 10
    j = 11
    j = 12

Why do you say there is only one iteration for i???

MATLAB: How to know the indices of multiple max values

One way:

A = [4 5 9 9];
maxval = max(A);
lia = ismember(A,maxval);
idx = find(lia);

or shorter

maxval = max(A);
idx = find(A == maxval);

Best Answer

Related Solutions

MATLAB: Nested loop will not enter second iteration

MATLAB: How to know the indices of multiple max values

Related Question