Statistical inference for functions of the covariance matrix in the stationary Gaussian time-orthogonal principal components model