The time series with periodic behavior, such as the periodic autoregressive (PAR) models belonging to the class of the periodically correlated processes, are present in various real applications. In the literature, such processes were considered in different directions, especially with the Gaussian-distributed noise. However, in most of the applications, the assumption of the finite-variance distribution seems to be too simplified. Thus, one can consider the extensions of the classical PAR model where the non-Gaussian distribution is applied. In particular, the Gaussian distribution can be replaced by the infinite-variance distribution, e.g. by the $α-$stable distribution. In this paper, we focus on the multidimensional $α-$stable PAR time series models. For such models, we propose a new estimation method based on the Yule-Walker equations. However, since for the infinite-variance case the covariance does not exist, thus it is replaced by another measure, namely the covariation. In this paper we propose to apply two estimators of the covariation measure. The first one is based on moment representation (moment-based) while the second one - on the spectral measure representation (spectral-based). The validity of the new approaches are verified using the Monte Carlo simulations in different contexts, including the sample size and the index of stability of the noise. Moreover, we compare the moment-based covariation-based method with spectral-based covariation-based technique. Finally, the real data analysis is presented.