I want to get an average time series for a set of multivariate (2-3 coordinates) time series. My aim is finding the usual pattern of several processes.

I researched the literature a bit and I only reached this paper that showed a DTW based approach (DBA) and talked about PSA and NLAAF. I also looked at the links of this question but they are more related to clustering than averaging.

first question: What are the recommended books papers on the topic of time series averaging?

Regarding the implementation, I was able to find a R package that performs the DBA in the dtwclust package but I found nothing on PSA or NLAAF nor other techniques.

second question: What time series averaging software implementations are there available? (preferably open source, I do not care much about the programming language)

Finally, I tried R's DBA implementation in some of my training sets and found an unpleasant result when dealing with repetitive circular trajectories of different lengths and speed (next image).

Last questions: Is there some theoretical consideration I should bare in mind when averaging such time series? What are the recommended preprocessing techniques when dealing with time series?

Averaging of spins with DBA

I know the questions asked are a bit different from each other, but they cover the same topic so I thought they were better together.

Below, I post the code to reproduce the image. The example included in the code consists of 25 time series of different lengths, with a different starting point and a different advancing speed. Think of it as cars racing on a circuit track, the cars are tracked from a random starting point and they are monitored until they finish the race. I calculate the average time series using the DBA function. The problem I find is that the average (plotted in black) does not respect the circular nature of the track and the races performed by the cars (red lines). Note that this is just a particular case, the trajectories I want to track are slightly more complex and, hence, I am more interested in the theory of the field rather than solving this particular example (which I believe can be solved by switching to polar coordinates).

# Each element of the cars list will be a trajectory
# 25 time series are included in the list

# initialise the seed for reproducibility
# each i generates a random time series
for(i in 1:ncars){
  # a random starting point is set
  # the time series has the coordinates of a process
  # that moves in circles. This is just an example,
  # the actual processes I'd like to average are slightly
  # more complex and move in 3 axes
  # time counter, it indicates the time elapsed from
  # the beginning of the observations
  # the block simulates spins at a varying random speed
    # The time series advance randomly 
  print(paste0("Duration of observation ",i,": ",t," seconds"))
  # store the time series in cars list
# average time series with DBA
fit<-DBA(lapply(cars,as.matrix), centroid = NULL, window.size = NULL, norm = "L1",
         max.iter = 20L, delta = 0.001, error.check = TRUE, tlap = FALSE) 
#  plot the average trajectory of the observations
# original trajectories of the observations
for(i in cars){lines(i,col=alpha(2,0.05))}