Matlab - Interpola en Matlab con varios NaN consecutivos

function x=interpola(x)

I=find(isnan(x));

m=length(I);

for i=1:m

    x(I(i))=(x(I(i)-1)+x(I(i)+1))/2;

end

end

>> x=[1,nan,3]

x =

     1   NaN     3

>> interpola(x)

ans =

     1     2     3

I=find(isnan(x)); % Encuentra los NaN

igrupos = [1 find(diff(I)~=1)+1 length(I)+1]; % Tiene los índices en donde empieza cada grupo de NaNs consecutivos y el número de elemento de I +1

for g=1:length(igrupos)-1

    grupo = I(igrupos(g)):I(igrupos(g+1)-1);

    paso = (x(grupo(end)+1) - x(grupo(1)-1))/(length(grupo)+1);

    for gg=1:length(grupo)

        x(grupo(gg)) = x(grupo(gg)-1)+paso;

    end

end