According to the method, a polar
decomposition technology is adopted, a measurement matrix of each strapdown triaxial
magnetometer is decomposed into a unit
orthogonal matrix and a symmetric
positive definite matrix, and an
ellipsoid equation is adopted for iterative linear least square fitting to obtain equivalent zero offset and the symmetric
positive definite matrix of each strapdown triaxial
magnetometer; and introducing a virtual coordinate
system, calculating a measured data correction value of each strapdown three-axis
magnetometer by using orthogonal Prusk
decomposition, and further calculating a correction value of a
magnetic gradient tensor measured by the strapdown three-axis cross-shaped magnetic array, thereby completing calibration of the strapdown three-axis cross-shaped magnetic array. According to the method, the iterative
linear least squares method is adopted, the
ellipsoid equation can be solved with high precision, the calibration precision is remarkably improved, the direction of the rotating shaft does not need to be kept stable in the rotating process, and the iterative
linear least squares algorithm and the orthogonal Prushing analysis
algorithm are small in calculated amount and easy to achieve in
engineering. Under the condition that the geomagnetic field model and the magnetic array attitude are provided, calibration can be completed quickly, simply and conveniently, and the method is suitable for practical
engineering application.