fastbsmooth.m

来自「Interactive smoothing for time-series si」· M 代码 · 共 16 行

function SmoothY=fastbsmooth(Y,smoothwidth)
% fastbsmooth(Y,w) smooths vector Y by unweighted sliding-average
% (boxcar) smooth of width = smoothwidth. Works well with signals up
% to 100,000 points in length and smooth widths up to 1000 points. 
% Faster than bsmooth for smooth widths above 500 points.
%  T. C. O'Haver, 2006.
w=round(smoothwidth);
SumPoints=sum(Y(1:w));
s=zeros(size(Y));
halfw=round(w/2);
for k=1:length(Y)-w,
   s(k+halfw)=SumPoints;
   SumPoints=SumPoints-Y(k);
   SumPoints=SumPoints+Y(k+w);
end
SmoothY=s./w;