RU2386141C1 - Method of determining parametres of calibration characteristics of magnetometre - Google Patents

Method of determining parametres of calibration characteristics of magnetometre Download PDF

Info

Publication number
RU2386141C1
RU2386141C1 RU2008127070/28A RU2008127070A RU2386141C1 RU 2386141 C1 RU2386141 C1 RU 2386141C1 RU 2008127070/28 A RU2008127070/28 A RU 2008127070/28A RU 2008127070 A RU2008127070 A RU 2008127070A RU 2386141 C1 RU2386141 C1 RU 2386141C1
Authority
RU
Russia
Prior art keywords
magnetometer
induction
readings
vector
axis
Prior art date
Application number
RU2008127070/28A
Other languages
Russian (ru)
Other versions
RU2008127070A (en
Inventor
Григорий Иванович Соборов (RU)
Григорий Иванович Соборов
Александр Николаевич Схоменко (RU)
Александр Николаевич Схоменко
Юрий Ромуальдович Линко (RU)
Юрий Ромуальдович Линко
Original Assignee
Открытое акционерное общество "Раменское приборостроительное конструкторское бюро"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Открытое акционерное общество "Раменское приборостроительное конструкторское бюро" filed Critical Открытое акционерное общество "Раменское приборостроительное конструкторское бюро"
Priority to RU2008127070/28A priority Critical patent/RU2386141C1/en
Publication of RU2008127070A publication Critical patent/RU2008127070A/en
Application granted granted Critical
Publication of RU2386141C1 publication Critical patent/RU2386141C1/en

Links

Landscapes

  • Measuring Magnetic Variables (AREA)

Abstract

FIELD: physics.
SUBSTANCE: method involves superposition of the magnetosensitive axis (MSA) of the magnetometre with the axis of the direction of the Earth's magnetic field (EMF) induction, determination of the result of converting this value using the magnetometre and then determination of results of converting a series of additional sample values of magnetic induction at the input of the magnetometre, where the said sample values are generated due to the recorded turns of the magnetosensitive axis of the magnetometre about the axis of the direction of the Earth's magnetic field induction vector in the plane of rotation, passing through this axis. Before taking measurements, the magnetosensitive axis of the magnetometre is superposed with the direction of the Earth's magnetic field induction vector through rough initial setup of the magnetosensitive axis in the direction of the Earth's magnetic field induction vector and then by turning the magnetosensitive axis about its initial position by angles +90° and -90° in the same plane of rotation. Readings of the magnetometre at the respective turning angles are taken. If the readings are different, initial setup is corrected, as well as the procedures for turning and taking readings.
EFFECT: wider technical capabilities of the procedure of determining parametres of calibration characteristics of a magnetometre, which enables use of simple equipment.

Description

Изобретение относится к области магнитных измерений, в частности к измерениям компонент и полного вектора индукции магнитного поля Земли, а также к средствам калибровки магнитометров.The invention relates to the field of magnetic measurements, in particular to measurements of the components and the full induction vector of the Earth's magnetic field, as well as to means for calibrating magnetometers.

Известен способ повышения точности магнитометра путем обеспечения точности калибровки, использования меры магнитной индукции, то есть точного определения параметров его статической характеристики преобразования.A known method of increasing the accuracy of the magnetometer by ensuring the accuracy of calibration, using a measure of magnetic induction, that is, accurately determining the parameters of its static conversion characteristics.

Известен способ повышения точности магнитометра путем обеспечения точности его калибровки использованием меры магнитной индукции или колец Гельмгольца и магнитного экрана [1]. В процессе калибровки магнитометра определяются параметры его статической характеристики преобразования или, что то же самое, градуировочной характеристики (ГХ). По результатам измерения выходных сигналов магнитометра, магниточувствительный элемент которого помещен в магнитный экран, определяется смещение нуля характеристики преобразования. При воздействии магнитной индукции образцового магнитного поля на магниточувствительный элемент, чувствительная ось которого совмещена с направлением вектора индукции, по результатам измерения выходных сигналов магнитометра определяется коэффициент преобразования характеристики преобразования магнитометра. Целесообразность данного способа калибровки объяснима применением магнитометров с высокой линейностью характеристики преобразования. Упрощение и снижение конструктивных качеств феррозондовых датчиков приводит к ограничению точностных параметров, вызванных в значительной степени нелинейностью характеристики преобразования. Обеспечение линейности характеристики вызывает потребность повышения технологических и конструктивных затрат, что повышает стоимость магнитометров.There is a method of increasing the accuracy of a magnetometer by ensuring the accuracy of its calibration using measures of magnetic induction or Helmholtz rings and a magnetic screen [1]. In the process of calibrating the magnetometer, the parameters of its static conversion characteristic or, equivalently, the calibration characteristic (GC) are determined. According to the results of measuring the output signals of the magnetometer, the magnetically sensitive element of which is placed in the magnetic screen, the zero shift of the conversion characteristic is determined. Under the influence of magnetic induction of an exemplary magnetic field on a magnetically sensitive element, the sensitive axis of which is aligned with the direction of the induction vector, the conversion coefficient of the magnetometer conversion characteristic is determined by measuring the output signals of the magnetometer. The feasibility of this calibration method is explained by the use of magnetometers with high linearity of the conversion characteristics. Simplification and reduction of the design qualities of flux-gate sensors leads to a limitation of the accuracy parameters, caused to a large extent by the nonlinearity of the conversion characteristics. Ensuring the linearity of the characteristic causes the need to increase technological and structural costs, which increases the cost of magnetometers.

Широкое внедрение высокотехнологичных датчиков низкой стоимости (например, магниторезистивных, датчиков Холла и др.) расширяет возможности применения магнитометров в различные области техники измерений, в то же время ограничивается возможность их применения в технике прецизионных измерений. Снижение точности преобразования подобных датчиков в значительной мере обусловлено нелинейностью их характеристик преобразования, выявление которой невозможно известным способом калибровки. Таким образом, недостатком указанного способа определения параметров ГХ является низкая точность калибровки, снижающая в свою очередь точность измерения параметров магнитного поля.The widespread introduction of high-tech low-cost sensors (for example, magnetoresistive, Hall sensors, etc.) expands the possibilities of using magnetometers in various fields of measurement technology, while at the same time, the possibility of their use in precision measurement techniques is limited. The decrease in the conversion accuracy of such sensors is largely due to the nonlinearity of their conversion characteristics, the detection of which is impossible by a known calibration method. Thus, the disadvantage of this method of determining the parameters of the GC is the low calibration accuracy, which in turn reduces the accuracy of the measurement of the magnetic field parameters.

В работе [2] предложен и показан способ повышения точности магнитометра путем коррекции его ГХ, аппроксимированной полиномом n-й степени. Позже в работе [3] показан пример практического применения данного способа в феррозондовом магнитометре малогабаритной системы контроля положения буксируемых объектов морской техники. В данном способе определение параметров ГХ осуществляется при калибровке методом эмпирических исследований математической модели характеристики преобразования, что усложняет определение параметров и затрудняет автоматизацию процесса калибровки.In [2], a method for improving the accuracy of a magnetometer by correcting its GC approximated by an nth degree polynomial is proposed and shown. Later in [3], an example of the practical application of this method in a flux-gate magnetometer of a small-sized system for monitoring the position of towed objects of marine equipment is shown. In this method, the determination of GC parameters is carried out during calibration by the method of empirical studies of a mathematical model of the conversion characteristics, which complicates the determination of parameters and complicates the automation of the calibration process.

Известен способ определения параметров аппроксимированной градуировочной характеристики при компенсации нелинейности [4], осуществляемый графоаналитическим методом исследования характеристики преобразования. В данном случае необходимость графического построения характеристики преобразования прибора требует большого количества измерений и больших временных затрат, приводящих к заметному усложнению способа.A known method for determining the parameters of the approximated calibration characteristics when compensating for nonlinearity [4], carried out by the grapho-analytical method for studying the conversion characteristics. In this case, the need for graphical construction of the conversion characteristics of the device requires a large number of measurements and a large time cost, leading to a significant complication of the method.

Показанный в работах [5] и [2] способ определения параметров ГХ, основанный на применении метода образцовых сигналов, исключает указанные недостатки аналогов, наиболее близок по технической сущности к предлагаемому и принят за прототип.The method for determining GC parameters, shown in [5] and [2], based on the application of the model signal method, eliminates the indicated drawbacks of analogues, is closest in technical essence to the proposed one and adopted as a prototype.

Указанный способ включает совмещение магниточувствительной оси (МЧО) магнитометра, помещенного в меру магнитной индукции, с направлением вектора магнитной индукции, задаваемой образцовым источником магнитного поля, последовательное воздействие на магниточувствительный элемент магнитометра образцового магнитного поля с различными значениями индукции и осуществление соответствующих этим значениям измерений выходных сигналов магнитометра, причем количество задаваемых воздействующих величин равно степени аппроксимирующего градуировочную характеристику магнитометра полинома, затем определение параметров ГХ в виде результатов решения системы уравнений, заданной в виде системы функций (результатов измерений) соответствующих значений индукции и параметров ГХ.The specified method includes combining the magnetically sensitive axis (MCO) of the magnetometer, placed to the extent of magnetic induction, with the direction of the magnetic induction vector specified by the exemplary source of the magnetic field, sequential exposure of the magnetically sensitive element of the magnetometer of the exemplary magnetic field with various values of induction and the implementation of the corresponding output signal measurements a magnetometer, and the number of specified acting quantities is equal to the degree of the approximating degree alignment characteristic of the polynomial magnetometer, then determining the GC parameters in the form of the results of solving a system of equations given in the form of a system of functions (measurement results) of the corresponding values of induction and GC parameters.

Математическое выражение ГХ магнитометра представим в виде следующего аппроксимирующего полинома n-й степениWe represent the mathematical expression of the GC magnetometer in the form of the following approximating polynomial of the nth degree

Figure 00000001
Figure 00000001

где N - выходная величина магнитометра, например, формируемая в виде кодового эквивалента;where N is the output value of the magnetometer, for example, formed in the form of a code equivalent;

В - измеряемая индукция магнитного поля, воздействующая на вход магнитометра;B is the measured magnetic field induction acting on the magnetometer input;

a1, a2, …, an - определяемые параметры ГХ, являющиеся постоянными коэффициентами.a 1 , a 2 , ..., a n are the determined parameters of the GC, which are constant coefficients.

Предполагая достаточно точное совмещение МЧО магнитометра, помещенного в меру магнитной индукции, с направлением вектора индукции задаваемого образцового магнитного поля, процесс определения параметров ГХ представим в виде последовательности n тактов, выражения результатов измерения которых образуют следующую систему уравненийAssuming a fairly accurate combination of the magnetic frequency response of the magnetometer placed to the measure of magnetic induction with the direction of the induction vector of the specified model magnetic field, the process of determining the GC parameters will be represented as a sequence of n clock cycles, the expressions of the measurement results of which form the following system of equations

Figure 00000002
Figure 00000002

где BOi (i=1, 2, …, n) - последовательно задаваемые в каждом такте значения индукции образцового магнитного поля;where B Oi (i = 1, 2, ..., n) are the values of the induction of an exemplary magnetic field sequentially set in each clock cycle;

N1, N2, …, Nn - показания магнитометра в соответствующих тактах преобразования.N 1 , N 2 , ..., N n - magnetometer readings in the corresponding conversion clocks.

Параметры аi определяются решением полученной системы уравнений (2) с помощью вычислительного устройства.The parameters a i are determined by solving the resulting system of equations (2) using a computing device.

В процессе измерения измеряемая индукция В исследуемого магнитного поля определяется решением уравнения (1) относительно измеряемой величины В после подстановки в данное уравнение найденных значений ai в результате решения системы уравнений (2).During the measurement, the measured induction B of the investigated magnetic field is determined by solving equation (1) with respect to the measured value B after substituting the found values of a i into this equation as a result of solving the system of equations (2).

Недостатками указанного способа являются необходимость обязательного использования уникальной конструкции дорогостоящих экранированной меры магнитной индукции или колец Гельмгольца, предполагающих наличия источника образцовых значений магнитной индукции и сложно настраиваемых катушек или колечной системы. Необходимость осуществления сложных регулировок и периодической аттестации этих устройств требует специального их обслуживания, а большие габаритные размеры вызывают необходимость специального помещения с большой площадью их размещения. Кроме того, неравномерность пространственного распределения формируемого магнитного поля внутри катушек магнитной меры вызывает потребность локализации ограниченных размеров рабочего объема. Для размещения калибруемых датчиков в этом объеме необходимы специальные механически перемещаемые поворотные устройства. Таким образом, уникальность и ограниченные возможности приобретения и содержания меры магнитной индукции или колечных систем значительно ограничивают возможности применения известного метода определения параметров ГХ магнитометра.The disadvantages of this method are the need to use the unique design of an expensive shielded measure of magnetic induction or Helmholtz rings, assuming a source of exemplary values of magnetic induction and difficult to configure coils or ring systems. The need for complex adjustments and periodic certification of these devices requires special maintenance, and large overall dimensions necessitate a special room with a large area of their placement. In addition, the uneven spatial distribution of the generated magnetic field inside the coils of the magnetic measure necessitates the localization of limited dimensions of the working volume. To place calibrated sensors in this volume, special mechanically movable rotary devices are required. Thus, the uniqueness and limited possibilities of acquiring and maintaining a magnetic induction measure or ring systems significantly limit the possibilities of using the known method for determining the parameters of a GC magnetometer.

Целью предлагаемого изобретения является расширение технической возможности способа определения параметров градуировочной характеристики магнитометра, позволяющей применение простых технических средств его осуществления.The aim of the invention is to expand the technical capabilities of the method for determining the parameters of the calibration characteristics of the magnetometer, allowing the use of simple technical means for its implementation.

Предлагаемый способ определения параметров ГХ магнитометра включает совмещение МЧО магнитометра с осью направления опорного вектора магнитной индукции известного значения, определение результата преобразования этого значения магнитометром, а затем определение результатов преобразования последовательности дополнительных воздействующих на вход магнитометра образцовых значений магнитной индукции, причем общее количество входных воздействий и соответственно преобразований их магнитометром равно степени аппроксимирующего ГХ полинома, затем определение параметров ГХ в виде результатов решения системы уравнений, заданной в виде системы функций входных воздействий и параметров ГХ магнитометра, при этом в качестве опорного вектора магнитной индукции выбран вектор индукции магнитного поля Земли (МПЗ), а формирование последовательности входных воздействий образцового поля осуществляется последовательностью фиксированных поворотов МЧО магнитометра относительно оси направления вектора индукции МПЗ вокруг начала этой оси в плоскости поворота, проходящей через эту ось, при этом повороты осуществляются на углы с разными значениями их модулей в диапазоне изменения модулей, не превышающего 180°, перед началом измерений совмещение МЧО магнитометра с направлением вектора индукции МПЗ осуществляется сначала приближенной (грубой) начальной установкой МЧО по направлению вектора индукции МПЗ, затем поворотами МЧО относительно ее начального положения на углы +90° и -90° в одной и той же плоскости поворота, определяются показания магнитометра при соответствующих углах поворота, при разных показаниях корректируется начальная установка, процедуры поворотов, определения показаний и по результатам показаний корректировки начального положения повторяются до полного установления равенства показаний магнитометра при поворотах на +90° и -90°, затем проведение процедур корректировки начальной установки МЧО продолжаются аналогичным образом поворотами МЧО относительно направления вектора индукции МПЗ в плоскости ее поворота, перпендикулярной первоначальной плоскости поворота, при необходимости процедуры начальной установки в разных плоскостях поворота повторяются до окончательного равенства нулю разности показаний магнитометра при углах поворота +90° и -90°.The proposed method for determining the parameters of a GC magnetometer includes combining the magneto-magnetic element of the magnetometer with the direction axis of the reference vector of magnetic induction of a known value, determining the result of converting this value to a magnetometer, and then determining the results of the conversion of the sequence of additional values of the magnetic induction acting on the magnetometer input, the total number of input actions and, accordingly transformations by their magnetometer is equal to the degree of approximating GC floor Nom, then the determination of GC parameters in the form of the results of solving a system of equations given in the form of a system of functions of input actions and GC parameters of a magnetometer, while the Earth’s magnetic field induction vector (MES) was chosen as the reference vector of magnetic induction, and the formation of a sequence of input actions of the reference field is carried out by a sequence of fixed rotations of the magneto-optical magnetometer magnetometer relative to the direction axis of the MPZ induction vector around the beginning of this axis in the rotation plane passing through axis, while rotations are carried out at angles with different values of their modules in the range of modules not exceeding 180 °, before the start of measurements, the MCH of the magnetometer is combined with the direction of the MPZ induction vector first by setting the rough (rough) initial setting of the MCH in the direction of the MPZ induction vector, then by turning the MCO relative to its initial position by angles + 90 ° and -90 ° in the same plane of rotation, the magnetometer readings are determined at the corresponding rotation angles, with different readings initial setup, turning procedures, determination of readings and, based on the results of readings, the initial position adjustments are repeated until the magnetometer readings are completely equal when turning + 90 ° and -90 °, then the procedures for adjusting the initial set-up of the magnetic mass spectrometer are continued in the same way by turning the mechanical mass spectrometer relative to the direction of the induction vector MPZ in the plane of its rotation perpendicular to the initial plane of rotation, if necessary, the initial installation procedure in different planes Rotating repeated until the final difference equal to zero magnetometer readings at rotation angles of + 90 ° and -90 °.

В предлагаемом способе определение параметров ГХ, представленной выражением (1), осуществляется решением следующей системы уравнений, образованной последовательными тактами преобразования воздействующих на вход магнитометра различных значений индукции магнитного поляIn the proposed method, the determination of GC parameters represented by the expression (1) is carried out by solving the following system of equations formed by successive clock ticks of the conversion of various values of the magnetic field induction acting on the magnetometer input

Figure 00000003
Figure 00000003

где В0 - значение индукции вектора МПЗ, предварительно измеренное прецизионным магнитометром;where B 0 - the value of the induction of the vector MPZ, previously measured with a precision magnetometer;

N1, N2, …, Nn - показания магнитометра в соответствующих тактах преобразования;N 1 , N 2 , ..., N n - magnetometer readings in the corresponding conversion clocks;

a1, a2, …, аn - параметры градуировочной характеристики;a 1 , a 2 , ..., and n are the parameters of the calibration characteristic;

k1, k2, …, kn-1 - коэффициенты преобразования индукции (В0) МПЗ при осуществлении поворотов МЧО.k 1 , k 2 , ..., k n-1 are the induction conversion coefficients (B 0 ) of the MPZ when cornering the MCO.

В данном случае воздействующее на магнитометр значение индукции в каждом такте преобразования является проекцией вектора индукции МПЗ на магниточувствительную ось при соответствующем ее наклоне относительно этого вектора, то естьIn this case, the induction value acting on the magnetometer in each conversion step is the projection of the MPF induction vector onto the magnetically sensitive axis with its corresponding slope relative to this vector, i.e.

Figure 00000004
Figure 00000004

где αi - угол поворота МЧО относительно вектора индукции МПЗ.where α i is the angle of rotation of the MCO relative to the induction vector

Процесс измерения индукции В исследуемого магнитного поля сопровождается решением уравнения (1) относительно измеряемой величины В после подстановки в данное уравнение найденных значений параметров αi в результате решения системы уравнений (3).The process of measuring induction B of the investigated magnetic field is accompanied by the solution of equation (1) with respect to the measured value B after substituting the found values of the parameters α i into this equation as a result of solving the system of equations (3).

Следует отметить, что параметр a1 можно считать уже определенным при ki=cos±90°=0 после осуществления начальной выставки МЧО.It should be noted that the parameter a 1 can be considered already determined at k i = cos ± 90 ° = 0 after the initial exhibition of the ICE.

Для формирования системы независимых уравнений (3) в каждом i-ом такте преобразования необходима реализация разных значений коэффициентов ki. Учитывая четность косинусной функции в выражении (4), необходимым условием формирования разных коэффициентов является неравенство модулей углов поворота в диапазоне изменения модулей, не превышающем 180°.For the formation of a system of independent equations (3) in each i-th transformation step, it is necessary to implement different values of the coefficients k i . Given the parity of the cosine function in expression (4), a necessary condition for the formation of different coefficients is the inequality of the moduli of rotation angles in the range of the moduli not exceeding 180 °.

Рассмотрим примеры решений простых ГХ.Consider examples of simple GC solutions.

В случае линейной ГХ магнитометра имеем систему двух уравнений с двумя неизвестными a1, a2 In the case of a linear GC magnetometer, we have a system of two equations with two unknowns a 1 , a 2

Figure 00000005
Figure 00000005

решения которой определяются выражениямиwhose decisions are determined by the expressions

Figure 00000006
Figure 00000007
Figure 00000006
Figure 00000007

При α1=90° (k=0) имеем a1=N2,

Figure 00000008
.When α 1 = 90 ° (k = 0) we have a 1 = N 2 ,
Figure 00000008
.

При квадратичной ГХ имеем систему трех уравнений с тремя неизвестнымиWith a quadratic GC, we have a system of three equations with three unknowns

Figure 00000009
Figure 00000009

решения данной системы при α1=90° (k=0), α2=180° (k=-1) имеют следующий видthe solutions of this system at α 1 = 90 ° (k = 0), α 2 = 180 ° (k = -1) have the following form

a1=N2, a 1 = N 2 ,

Figure 00000010
Figure 00000011
Figure 00000010
Figure 00000011

Таким образом можно определить параметры ГХ при ее аппроксимации полиномом любой степени n. Точность определения параметров ГХ данным способом зависит от точности измерения значения индукции В0 вектора МПЗ, точности задания углов поворота и точности совмещения МЧО с направлением вектора индукции МПЗ (начальной выставки МЧО). Измерение В0 легко осуществимо, например, протонным или квантовым прецизионным измерителем модуля магнитной индукции. Реально достижимой точностью осуществления поворотов практически реализуемыми поворотными средствами является величина, не превышающая одной угловой минуты, что заметно выше необходимой для данного способа точности осуществления поворотов. Достаточно простой и удобной конструкцией поворотного средства является, например, карданное поворотное устройство с тремя степенями свободы.Thus, it is possible to determine the GC parameters when it is approximated by a polynomial of any degree n. The accuracy of determining the GC parameters in this way depends on the accuracy of measuring the induction value B 0 of the MPZ vector, the accuracy of setting the rotation angles, and the accuracy of combining the MCH with the direction of the induction vector of the MMP (initial exhibition of MCH). The measurement of B 0 is easily feasible, for example, by a proton or quantum precision meter of the magnetic induction module. Really achievable accuracy of the implementation of turns practically feasible by turning means is a value not exceeding one angular minute, which is noticeably higher than the accuracy of the implementation of turns necessary for this method. A rather simple and convenient design of the rotary means is, for example, a universal joint rotary device with three degrees of freedom.

Простота изготовления, сравнительно низкая стоимость, простота контроля и периодической аттестации, сравнительно малые габариты и вес, обеспечивающие возможность переноса и размещения в различных местах рабочего помещения, характеризуемых отсутствием магнитных помех, а также малая площадь размещения определяют преимущественные особенности и доступность изготовления или приобретения применяемого в предлагаемом способе поворотного средства по сравнению с мерой магнитной индукции, применяемой в способе-прототипе.The simplicity of manufacture, relatively low cost, ease of monitoring and periodic certification, relatively small dimensions and weight, providing the possibility of transfer and placement in various places of the working room, characterized by the absence of magnetic interference, as well as a small area of placement determine the advantageous features and availability of manufacturing or acquisition used in the proposed method of rotary means in comparison with the measure of magnetic induction used in the prototype method.

Немаловажным преимущественным достоинством предлагаемого способа является возможность использования равномерно распределенной в пространственном объеме индукции МПЗ, используемой при проведении измерительных процедур.An important advantage of the proposed method is the possibility of using the MPZ induction evenly distributed in the spatial volume used in the measurement procedures.

Следовательно, предлагаемое изобретение, обладая новизной, полезностью и реализуемостью, может найти широкое применение в технике магнитных измерений.Therefore, the present invention, having novelty, utility and feasibility, can find wide application in the technique of magnetic measurements.

ЛитератураLiterature

1. Афанасьев Ю.В., Студенцов Н.В., Хорев В.Н. и др. Средства измерений параметров магнитного поля. - Л.: Энергия. Ленингр. отд-ние, 1979.1. Afanasyev Yu.V., Studentsov N.V., Khorev V.N. et al. Measuring instruments for magnetic field parameters. - L .: Energy. Leningra. Department, 1979.

2. Соборов Г.И. Разработка и исследование бортовых трехкомпонентных феррозондовых магнитометров моноблочного типа с жестко связанными осями. Диссертация на соискание ученой степени к.т.н. Московский энергетический институт (технический университет). Москва 2003.2. Soborov G.I. Development and research of onboard three-component flux-gate magnetometers of monoblock type with rigidly connected axes. Thesis for the degree of Ph.D. Moscow Power Engineering Institute (Technical University). Moscow 2003.

3. Андреев М.Я., Гаспаров П.М., Геркус А.А. и др. Малогабаритная система контроля положения буксируемых объектов морской техники. Ж. Датчики и системы № 5, 2006.3. Andreev M.Ya., Gasparov P.M., Gerkus A.A. and others. Small-sized system for monitoring the position of towed objects of marine equipment. G. Sensors and systems No. 5, 2006.

4. Акчурин Д.В., Богонин Б.В., Богонин М.Б. Применение полиномов Ньютона, Гаусса и Лагранжа при компенсации погрешности нелинейности. Ж. Датчики и системы № 12, 2002. - С.16.4. Akchurin D.V., Bogonin B.V., Bogonin M.B. The use of Newton, Gauss and Lagrange polynomials in compensating for the nonlinearity error. G. Sensors and systems No. 12, 2002. - P. 16.

5. Бромберг Э.М., Куликовский К.Л. Тестовые методы повышения точности измерений. - М.: Энергия, 1978. - С.20-22.5. Bromberg E.M., Kulikovsky K.L. Test methods to improve measurement accuracy. - M .: Energy, 1978. - S.20-22.

Claims (1)

Способ определения параметров градуировочной характеристики магнитометра включает совмещение магниточувствительной оси (МЧО) магнитометра с осью направления опорного вектора магнитной индукции известного значения, определение результата преобразования этого значения магнитометром, а затем определение результатов преобразования последовательности дополнительных воздействующих на вход магнитометра образцовых значений магнитной индукции, причем общее количество входных воздействий и соответственно преобразований их магнитометром равно степени аппроксимирующего градуировочную характеристику полинома, затем определение параметров градуировочной характеристики в виде результатов решения системы уравнений, заданной в виде системы функций входных воздействий и параметров градуировочной характеристики магнитометра, отличающийся тем, что в качестве опорного вектора магнитной индукции выбран вектор индукции магнитного поля Земли (МПЗ), а формирование последовательности дополнительных входных воздействий образцового поля осуществляется последовательностью фиксированных поворотов МЧО магнитометра относительно оси направления вектора индукции МПЗ в плоскости поворота, проходящей через эту ось, при этом повороты осуществляются на углы с разными значениями их модулей в диапазоне изменения модулей, не превышающего 180°, причем перед началом измерений совмещение МЧО магнитометра с направлением вектора индукции МПЗ осуществляется сначала приближенной (грубой) начальной установкой МЧО по направлению вектора индукции МПЗ, затем поворотами МЧО относительно ее начального положения на углы +90° и -90° в одной и той же плоскости поворота, определяются показания магнитометра при соответствующих углах поворота, при разных показаниях корректируется начальная установка, процедуры поворотов, определения показаний, и по результатам показаний корректировки начального положения повторяются до полного установления равенства показаний магнитометра при поворотах на +90° и -90°, затем проведение процедур корректировки начальной установки МЧО продолжаются аналогичным образом поворотами МЧО относительно направления вектора индукции МПЗ в плоскости ее поворота, перпендикулярной первоначальной плоскости поворота, при необходимости процедуры начальной установки в разных плоскостях поворота повторяются до окончательного равенства нулю разности показаний магнитометра при углах поворота +90° и -90°. The method for determining the calibration characteristics of the magnetometer includes combining the magnetically sensitive axis (MCH) of the magnetometer with the direction axis of the reference vector of magnetic induction of a known value, determining the result of the conversion of this value by a magnetometer, and then determining the results of the conversion of the sequence of additional values of the magnetic induction acting on the magnetometer input, the total number input actions and, accordingly, transformations of their magnetometer rum is equal to the degree of approximation of the calibration characteristic of the polynomial, then the determination of the calibration characteristic parameters in the form of the results of solving a system of equations defined as a system of input action functions and the calibration characteristics of the magnetometer, characterized in that the Earth’s magnetic field induction vector is chosen as the reference magnetic induction vector ( MPZ), and the formation of a sequence of additional input effects of the model field is carried out by the sequence fixed rotations of the magneto-magnetic resonator magnetometer relative to the direction axis of the MPZ induction vector in the rotation plane passing through this axis, while rotations are carried out at angles with different values of their modules in the range of the modules not exceeding 180 °, and before the start of the measurements, the magneto-magnetic magnetometer magnetometer alignment with the direction of the vector the induction of the MPZ is carried out first by an approximate (rough) initial installation of the MPM in the direction of the vector of the induction of the MPZ, then by turning the MCO relative to its initial position by angles + 90 ° and -9 0 ° in the same rotation plane, the magnetometer readings are determined at the corresponding rotation angles, the initial setting, the rotation procedures, determination of readings are corrected for different readings, and the initial position adjustments are repeated according to the readings until the magnetometer readings are equal when turning +90 ° and -90 °, then the procedures for adjusting the initial setup of the MCH continue in the same way by turning the MCH relative to the direction of the induction vector the speed of its rotation perpendicular to the initial plane of rotation, if necessary, the initial installation procedures in different rotation planes are repeated until the magnetometer readings are completely equal to zero at rotation angles + 90 ° and -90 °.
RU2008127070/28A 2008-07-04 2008-07-04 Method of determining parametres of calibration characteristics of magnetometre RU2386141C1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2008127070/28A RU2386141C1 (en) 2008-07-04 2008-07-04 Method of determining parametres of calibration characteristics of magnetometre

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2008127070/28A RU2386141C1 (en) 2008-07-04 2008-07-04 Method of determining parametres of calibration characteristics of magnetometre

Publications (2)

Publication Number Publication Date
RU2008127070A RU2008127070A (en) 2010-01-10
RU2386141C1 true RU2386141C1 (en) 2010-04-10

Family

ID=41643800

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2008127070/28A RU2386141C1 (en) 2008-07-04 2008-07-04 Method of determining parametres of calibration characteristics of magnetometre

Country Status (1)

Country Link
RU (1) RU2386141C1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2481593C9 (en) * 2011-06-03 2013-09-10 Открытое акционерное общество "Раменское приборостроительное конструкторское бюро " (ОАО "РПКБ") Method of determining parameters of transformation characteristics of three-component magnetometer
RU2497139C1 (en) * 2012-05-18 2013-10-27 Закрытое акционерное общество "МЕРА" Method to calibrate three-component magnetometer
RU2687170C1 (en) * 2018-06-06 2019-05-07 Акционерное общество "Научно-исследовательский институт командных приборов" Method for determining a characteristic of a ferroprobe during temperature tests
RU2724314C1 (en) * 2019-09-16 2020-06-22 Акционерное общество "Научно-исследовательский институт командных приборов" Method to control the ferroprobe conversion characteristic
RU2749303C1 (en) * 2020-09-09 2021-06-08 Акционерное общество "Научно-исследовательский институт командных приборов" Method for determining admissibility of using fluxgate in magnetometer

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
СОБОРОВ Г.И. Разработка и исследование бортовых трехкомпонентных феррозондовых магнитометров моноблочного типа с жестко связанными осями. Диссертация на соискание ученой степени к.т.н. (05.11.01) МЭИ, 2003. *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2481593C9 (en) * 2011-06-03 2013-09-10 Открытое акционерное общество "Раменское приборостроительное конструкторское бюро " (ОАО "РПКБ") Method of determining parameters of transformation characteristics of three-component magnetometer
RU2497139C1 (en) * 2012-05-18 2013-10-27 Закрытое акционерное общество "МЕРА" Method to calibrate three-component magnetometer
RU2687170C1 (en) * 2018-06-06 2019-05-07 Акционерное общество "Научно-исследовательский институт командных приборов" Method for determining a characteristic of a ferroprobe during temperature tests
RU2724314C1 (en) * 2019-09-16 2020-06-22 Акционерное общество "Научно-исследовательский институт командных приборов" Method to control the ferroprobe conversion characteristic
RU2749303C1 (en) * 2020-09-09 2021-06-08 Акционерное общество "Научно-исследовательский институт командных приборов" Method for determining admissibility of using fluxgate in magnetometer

Also Published As

Publication number Publication date
RU2008127070A (en) 2010-01-10

Similar Documents

Publication Publication Date Title
CN108692743B (en) Redundant fault detection apparatus and method
Hu et al. A cubic 3-axis magnetic sensor array for wirelessly tracking magnet position and orientation
RU2386141C1 (en) Method of determining parametres of calibration characteristics of magnetometre
CN101887068B (en) Calibration compensation method for triaxial vector sensor and biaxial vector sensor
CN103941309B (en) Geomagnetic sensor calibrator (-ter) unit and its method
CN113325353B (en) Magnetometer spatial attitude calibration method and system
CN102853760A (en) Method for calibrating verticality of magnetic shaft of three-shaft magnetic sensor
CN107544042A (en) A kind of magnetometer array correcting method
US11175160B2 (en) Magnetic field sensor and method with reduced distortion measurement in sideshaft applications
RU2481593C2 (en) Method of determining parameters of transformation characteristics of three-component magnetometer
CN102589716B (en) Calibration method used for fully polarimetric microwave radiometer
CN109489686B (en) Method for testing bandwidth and angular resolution of four-frequency differential laser gyroscope
RU2572109C1 (en) Method to calibrate electronic magnetic compass
RU2433421C1 (en) Method of determining magnetometer calibration curve parameters
US6714008B1 (en) Gradiometric measurement methodology for determining magnetic fields of large objects
CN115097370A (en) Steering difference calibration method for self-balancing vector magnetometer in large-plane magnetic measurement system
RU138023U1 (en) DEVICE FOR CALIBRATING A THREE-COMPONENT MAGNETOMETER
RU2623192C1 (en) Method for calibrating electronic magnetic compass
KR101815584B1 (en) Compensation method for bias errors in 3D intensity probe used for the source localization in precision
EP0135143B1 (en) Imaging method and apparatus using nuclear magnetic resonance
Bogatyrev et al. Technology for calibration of measuring instruments of samsat nanosatellites' family
JP2003254784A (en) Method and device for calibrating displacement
US11402439B2 (en) Method for calibrating the sensitivity of monoaxial or multiaxial magnetic field sensors
RU2235200C2 (en) Magnetometric sensor
CN114234951B (en) Magnetic field fluctuation testing method of SERF inertial device based on nuclear spin polarization suppression