This paper presents a new computational algorithm of phase-shifting interferometry which can effectively eliminate the uncertainty errors of the reference phases encountered in obtaining multiple interferograms. The algorithm treats the reference phases as additional unknowns and determines their exact values by analyzing interferograms using numerical least square technique. A series of simulations prove that this algorithm can improve measuring accuracy being unaffected by the nonlinear and random errors of phase-shifters.