The integration of different GNSS constellations offers considerable opportunities to improve Precise Point Positioning (PPP) performance. Being aware of the limited number of the alternatives that utilize the potential advantages of the multi-constellation and multi-frequency GNSS, we developed a MATLAB-based GNSS analysis software, named PPPH. PPPH is capable of processing GPS, GLONASS, Galileo and BeiDou data, and forming their different combinations depending on user's preference. Thanks to its user-friendly graphical interface, PPPH allows users to determine a variety of processing options and parameters. In addition to an output file including the estimated parameters for every single epoch, PPPH also presents several analyzing and plotting tools for evaluating the results, such as positioning error, tropospheric zenith total delay, receiver clock estimation, satellite number, dilution of precisions. On the other hand, we conducted experimental tests to both validate the performance of PPPH and assess the potential benefits of multi-GNSS on PPP. The results indicate that PPPH provides comparable PPP solution with the general standards and also contributes to the improvement of PPP performance with the integration of multi-GNSS. Consequently, we introduce a GNSS analysis software that is easy to use, has a robust performance and is open to progress with its modular structure.