New estimation method for periodic autoregressive time series of order 1 with additive noise

The periodic behavior of real data can be manifested in the time series or in its characteristics. One of the characteristics that often manifests the periodic behavior is the sample autocovariance function. In this case, the periodically correlated (PC) behavior is considered. One of the main models that exhibits PC property is the periodic autoregressive (PARMA) model that is considered as the generalization of the classical autoregressive moving average (ARMA) process. However, when one considers the real data, practically the observed trajectory corresponds to the "pure" model with the additional noise which is a result of the noise of the measurement device or other external forces. Thus, in this paper we consider the model that is a sum of the periodic autoregressive (PAR) time series and the additive noise with finite-variance distribution. We present the main properties of the considered model indicating its PC property. One of the main goals of this paper is to introduce the new estimation method for the considered model's parameters. The novel algorithm takes under consideration the additive noise in the model and can be considered as the modification of the classical Yule-Walker algorithm that utilizes the autocovariance function. Here, we propose two versions of the new method, namely the classical and the robust ones. The effectiveness of the proposed methodology is verified by Monte Carlo simulations. The comparison with the classical Yule-Walker method is presented. The approach proposed in this paper is universal and can be applied to any finite-variance models with the additive noise.