This paper presents a method of ballistic missile launch point estimation using phase division least squares. The proposed algorithm employs smoothing to enhance estimation accuracy and generates functions of time for total velocity, flight path angle and heading angle, allowing extrapolation to estimate the launch point. Performance of the proposed algorithm is tested in conjunction with the extended Kalman filter and the Kalman filter.