movingave#

Purpose#

Computes moving average of a series.

Format#

y = movingave(x, d)#

Parameters:

x (NxK matrix) – data
d (scalar) – order of moving average.

Returns:

y (NxK matrix) – filtered series. The first \(d-1\) rows of x are set to missing values.

Examples#

Smoothing of random walk data#

// Create random data, setting initial seed
// value to be the hundredths of a second
// since midnight
numPoints = 1000;
seed = hsec;
{ delta, state } = rndn(numPoints, 1, seed);

// Set mean of data to be > 0
// to give data long term positive bias
delta = delta + 0.01;

// Increase magnitude of delta to
// create desired level of volatiliy
delta = 2*delta;

// Instantiate y : index price data
y = 1000*ones(numPoints,1);

// Loop through y and add the cumulative
// sums of delta to create random walk
for i(2, numPoints, 1);
      y[i]= y[i-1] + delta[i];
endfor;

// Find moving average
twentyMA = movingave(y, 20);

Remarks#

movingave is essentially a smoothing time series filter. The moving average is performed by column and thus it treats the NxK matrix as K time series of length N.