The use of three-axis magnetic sensor array to locate underwater ferromagnetic target is a typical nonlinear least squares optimization issue. Traditional Gaussian-Newton(GN)method and Levenberg-Marquardt (LM)algorithm have an initial sensitivity problem when solving this issue. In this paper,the LM algorithm is improved by introducing the trust domain search technology,and magnetic target localization is realized based on the improved LM algorithm. And then,the distance between the iteration point and the optimal solution is evaluated by setting the threshold value. An improved LM-GN algorithm combining the characteristics of the improved LM algorithm and Gaussian-Newton method is proposed. It not only reduces the dependence of the algorithm on the initial value,but also improves the operation efficiency. The simulation results show that the proposed method can overcome the problem of the existing methods which are greatly affected by the initial value,and the estimation of the target characteristic parameters is more accurate and the convergence rate is fast,so it has certain practical application value.