The conjugate gradient normal equations residual minimizing (CGNR) algorithm is a popular tool for solving large nonsymmetric linear systems. In this study, we propose the matrix form of the CGNR (MCGNR) algorithm to find the least squares solution group of the discrete-time periodic coupled matrix equations

{A1,tXtB1,t+C1,tXt+1D1,t+E1,tYtF1,t=G1,t,A2,tXtB2,t+C2,tXt+1D2,t+E2,tYtF2,t=G2,t,fort=1,2,...

with periodic matrices. We prove that the MCGNR algorithm converges in a finite number of steps in the absence of round-off errors, that is, this algorithm has the finite termination property. Also we show that the norms of the residual matrices of the MCGNR algorithm decrease monotonically during its iteration, that is, this algorithm has the residual reducing property. To show the efficiency of the MCGNR algorithm, two numerical examples are presented.