Matlab - contar solo los 1 de un vector

%%%%%%%%%%%%%% CONTAR NÚMERO DÍAS CUMPLEN CONDICIÓN%%%%%%

kk_mes=find(Pr<1); %Siendo  Pr la precipitación 

proba=diff(kk_mes);% con la función diff encuentro los días seguidos

%%%Este bucle que viene a continuación sirve para hacer cero calquier valor que sea mayor que 1

%%%es decir, los que no son seguidos, así obtenemos un proba2 que solo ten 0 o 1

[nc nf]=size(proba);

for i=1:nc

    xxx=proba(i);

    if xxx>1

        proba2(i)=0;

    else

        proba2(i)=1;

    end

end

clear all

clc

Pr=1.5*rand(1,15)

kk_mes=find(Pr<1) %Siendo  Pr la precipitación 

proba=diff(kk_mes) % con la función diff encuentro los días seguidos 

%%%Este bucle que viene a continuación sirve para hacer cero calquier valor que sea mayor que 1

%%%es decir, los que no son seguidos, así obtenemos un proba2 que solo ten 0 o 1

[nc, nf]=size(proba);

for i=1:nf

    xxx=proba(i);

    if xxx>1

        proba2(i)=0;

    else

        proba2(i)=1;

    end

end

proba2

dias=sum(proba2)

>> precipitaciones2

Pr =

    0.4289    0.8155    1.4772    1.0735    1.2585    0.6499    0.7059    0.8411    0.4036    1.1235    0.7558    0.9702    0.4616    0.2081    0.7134

kk_mes =

     1     2     6     7     8     9    11    12    13    14    15

proba =

     1     4     1     1     1     2     1     1     1     1

proba2 =

     1     0     1     1     1     0     1     1     1

dias =

     7

clear all,close all,clc

prec = round(1.3*rand(100,1));

nd = length(prec);

n = 0;

for nc= 1:nd

if nc ==1

    if prec(nc)== 1

      n = n +1;

      cuenta(nc) = n;

    else

        n = 0;

        cuenta(nc) = n;

    end

else

    if and(prec(nc)== 1 ,prec(nc-1)==1)

      n = n +1;

      cuenta(nc) = n;

    else

        n = 0;

        cuenta(nc) = n;

    end

end

end

%%

subplot(211)

plot(prec,'o-')

ylabel('Precipitación ')

subplot(212)

stem(cuenta)

ylabel('Días de precipitación')