In this thesis, we modified an existing program, CMBFAST, to have the capability of considering the finite universes. The original CMBFAST can compute the power spectrum of the Cosmic Microwave Background (CMB) only in infinite universes. With this modified program, we simulated a set of power spectrum database in a six-dimensional parameter space that accounts for the six cosmological parameters including the size of the hyper-torus topology. By using the maximum-likelihood method, we then compared these theoretical predictions with the observed WMAP data, so as to estimate the six cosmological parameters. In this investigation, we consider both the CMB temperature and the polarization.