Generate some sample data:
C = (1:3630).';
idxout = randperm(3630,500); % random index for outliers
C(idxout) = C(idxout)+2000*(1-2*rand(numel(idxout),1));
Using rmoutliers some of the values are removed:
window = 50;
[A, B] = rmoutliers(C, 'movmedian', window);
Seeing as you've indicated that you only have a column vector of data I assume you want to interpolate using index position as your input for x, you could try the following:
idx = find(~B); % C(idx) is the same as A
D = interp1(idx,A,1:size(C,1)); % Interpolated using index position as x values.
Let me know.
