Martin Gorjan

The above error is a flaw in the code. The programmer forgot to include a check for what happens during shifting of signal when the peak index is too small or too large (i.e., at the edges of the matrix).

One could include the following code just before the shifting to bypass all peaks near the edges. Just define number of elements to ignore on each edge.

    % check if index too big
    if ipeak >= length(y)-numElements
        break
    end

    % check if index too small
    if ipeak <= numElements
        break
    end

yy = smooth(y) smooths the data in the column vector y using a moving average filter. Results are returned in the column vector yy. The default span for the moving average is 5.

The first few elements of yy are given by

yy(1) = y(1)
yy(2) = (y(1) + y(2) + y(3))/3
yy(3) = (y(1) + y(2) + y(3) + y(4) + y(5))/5
yy(4) = (y(2) + y(3) + y(4) + y(5) + y(6))/5

So if you do not have the MATLAB "Curve Fitting Toolbox", you can make your own smooth function from the above description.

smooth works fine, it's a standard function in my matlab. maybe you don't have the right addons