Abstract: The equations for the solution of global point-mass model based on the complex combination of global gravity gradients constituents and the full tensor forms are established with the observational information output from the gradiometer equipped on the satellite. With the latitudinal and longitudinal partition arrangement of satellite gravity gradients data and point-mass model, a method is put forward to ensure the stability of the solutions of the large linear equations in the establishment of the global point-mass model using blocked cyclic matrices. Then, the point-mass harmonic analysis solutions of both single determined and over determined boundary value problems based on satellite gravity gradients are obtained. Finally, simulations experiments with EGM2008 are made to validate the efficiency and availability of the point-mass harmonic analysis of single or full gravity gradient tensor constituents. The results show that with proper depth of the point mass, the method considered in this thesis could attain a little higher precision gravity field than the general torus harmonic analysis.