Aiming to address poor self-consistency of the satellite clock and orbit interpolations at the day boundary, which is caused by the discontinuities of their International GNSS Service (IGS) products for two consecutive days, an effective method is proposed to improve the precise point positioning (PPP) performances at the day boundary. According to different orders of Lagrange interpolator and different IGS Analysis Centers products of Center for Orbit Determination in Europe (CODE), GeoForschungsZentrum (GFZ) and Wuhan University (WUH), biases at the day boundaries are estimated and analyzed using a 4 day (DOY 94-97, 2022) data set of GPS, BDS-3 and Galileo from 123 IGS stations. These estimated biases show the time-varying characteristics. The differences in biases across 9th-11th orders of Lagrange interpolator are minimal, and these variances have a negligible impact on positioning. The results show that this poor self-consistency at the day boundary has an obvious influence on the kinematic PPP positioning, especially there is a centimeter-level variation at time of 24:00:00/00:00:00. This influence on the Up direction of kinematic PPP positioning is more serious than other directions. When the bias is estimated and corrected, the kinematic PPP positioning accuracies at the day boundary have a mean improvement of 0.043, 0.064 and 0.027 m for WUH, GFZ and CODE, respectively. The mean improvements for GPS, BDS-3 and Galileo are 0.021, 0.062 and 0.051 m. Additionally, the static PPP performances at the day boundary show the convergence times are shortened by 3.2, 6.2, and 2.5 min for WUH, GFZ and CODE, respectively, when the poor self-consistency of the satellite clock and orbit interpolations is estimated and corrected. Meanwhile, its 0.5 and 1 h positioning accuracies are improved.