A Missile-Borne Angular Velocity Sensor Based on Triaxial Electromagnetic Induction Coils

Aiming to solve the problem of the limited measuring range for angular motion parameters of high-speed rotating projectiles in the field of guidance and control, a self-adaptive measurement method for angular motion parameters based on the electromagnetic induction principle is proposed. First, a framework with type bent “I-shape” is used to design triaxial coils in a mutually orthogonal way. Under the condition of high rotational speed of a projectile, the induction signal of the projectile moving across a geomagnetic field is acquired by using coils. Second, the frequency of the pulse signal is adjusted self-adaptively. Angular velocity and angular displacement are calculated in the form of periodic pulse counting and pulse accumulation, respectively. Finally, on the basis of that principle prototype of the sensor is researched and developed, performance of measuring angular motion parameters are tested on the sensor by semi-physical and physical simulation experiments, respectively. Experimental results demonstrate that the sensor has a wide measuring range of angular velocity from 1 rps to 100 rps with a measurement error of less than 0.3%, and the angular displacement measurement error is lower than 0.2°. The proposed method satisfies measurement requirements for high-speed rotating projectiles with an extremely high dynamic range of rotational speed and high precision, and has definite value to engineering applications in the fields of attitude determination and geomagnetic navigation.


Introduction
High-speed rotating projectiles, with rotational speeds generally in a range from 2 rps to 50 rps, have become one of the most important kinds of ammunition among conventional ammunition. On the one hand, a high-speed rotating projectile ensures flight stability by way of rotation; on the other hand, this kind of projectile changes flight attitude and corrects its flight trajectory through adjusting its rotational speed, so as to accurately attack targets. In the course of structure improvement and performance testing of high-speed rotating projectiles, the accurate measurement of angular motion parameters, such as angular velocity and angular displacement, is key for optimizing the structure and improve the attacking accuracy. The precision guidance of high-speed rotating projectiles is a research hotspot in the field of international arms and ammunition at present.
Scientific researchers mainly use gyroscopes to obtain the angular velocity information of projectiles nowadays [1][2][3][4][5][6][7][8][9][10][11][12]. For example, the angular velocity sensor based on a magnetically suspended control moment gyroscope (MSCMG) [13], the capacitive angular velocity sensor and the angular velocity sensor based on magnetohydrodynamics at low frequency [14,15], etc. Although these sensors have certain advantages of small volume and high damage impact resistance, they still have problems of small measuring ranges, zero drift and high cumulative errors when measuring rotational problems of small measuring ranges, zero drift and high cumulative errors when measuring rotational speed. Though some of the errors can be corrected by some algorithms like Kalman prediction, the requirements of large dynamic range and high precision for measuring the rotational speed of high-speed rotating projectiles still cannot be satisfied under actual combat conditions. Aiming at addressing the problems described above, this paper presents a method for measuring angular motion parameters based on the electromagnetic induction principle. The specific working of this method is as follows: under different rotation speeds, an induction coil of a bent "I-shape" type is used to acquire signals of a high-speed rotating projectile moving across a geomagnetic field. A self-adaptive frequency tracking measurement algorithm is proposed to obtain the angular velocity and angular displacement information of the projectile rotation. The cumulative error during the measurement is eliminated by using the periodically pulse clear method, thereby realizing the measurement of the angular motion parameters of the high-speed rotating projectile, and satisfying its measurement requirements of wide range of rotational speed and high precision. The proposed method has definite value to engineering applications in the fields of attitude determination and geomagnetic navigation [16].

Basic Principle of Measuring Angular Velocity
The measurement method of angular velocity proposed in this paper draws from the theory of digital frequency meters. The basic principle is shown in Figure 1. A measured signal x(t) is converted into a gate signal through a shaping circuit. Under the circumstances that the standard frequency signal xs(t) with higher frequency is taken as a periodic pulse signal, the measured signal is calculated by pulse counting within a gate time T. Assuming that the count value is N, the cycle of measured signal and standard signal are Tx and T0, respectively. Then, the cycle of the measured signal is: and the frequency of the measured signal is: It is known by observation of Equations (1) and (2) that when counting a standard frequency signal with periodic pulse counting, the frequency of the measured signal can be calculated accordingly. However, rotational angular motion parameters of the projectile measured in this paper include not only the angular velocity information of the projectile rotation, but also real-time angular displacement information of the rotation. For the angular displacement information, the rotation angle of projectile that one cycle Tx of the measured signal corresponds to 360°. The angular A measured signal x(t) is converted into a gate signal through a shaping circuit. Under the circumstances that the standard frequency signal x s (t) with higher frequency is taken as a periodic pulse signal, the measured signal is calculated by pulse counting within a gate time T. Assuming that the count value is N, the cycle of measured signal and standard signal are T x and T 0 , respectively. Then, the cycle of the measured signal is: and the frequency of the measured signal is: It is known by observation of Equations (1) and (2) that when counting a standard frequency signal with periodic pulse counting, the frequency of the measured signal can be calculated accordingly. However, rotational angular motion parameters of the projectile measured in this paper include not only the angular velocity information of the projectile rotation, but also real-time angular displacement information of the rotation. For the angular displacement information, the rotation angle of projectile that one cycle T x of the measured signal corresponds to 360 • . The angular displacement is output in a cumulative way by using standard frequency signal x s (t). The period of each pulse represents a fixed step-progressive accumulation angle θ, where θ = 360 • /N. In this way, θ accumulates to 360 • with N cycles. Because the angular displacement that each pulse represents is fixed, therefore, the count value N is a stable constant. Variable N stands for the resolution of the angular displacement. However, in the process of frequency measurement, the frequency of a standard signal is changeless, while the count value N varies with the difference of the measured signal. Consequently, angular displacement information cannot be output effectively. Aiming at the problems mentioned above, this paper presents a self-adaptive frequency tracking measurement algorithm. On the premise that the value of pulse counting is N, the algorithm realizes measurement of angular displacement and angular velocity by changing the frequency of the standard signal self-adaptively.

Principle of Self-Adaptive Frequency Tracking Measurement
The principle of the self-adaptive frequency tracking measurement method is shown in Figure 2. Supposing that measured signal is x(t), the standard signal is x s (t). When the pulse counting value is fixed at N, the cumulatively output angular displacement is 180 • . The frequencies of the first three cycles of the measured signal are f 1 , f 2 and f 3 , respectively. The corresponding frequencies of the first three cycles of the standard signal are fs 0 , fs 1 and fs 2 , respectively. The half period of the measured signal x(t) is taken as gate time after processing the signal with shaping circuit modulation.
Sensors 2016, 16,1625 3 of 12 displacement is output in a cumulative way by using standard frequency signal xs(t). The period of each pulse represents a fixed step-progressive accumulation angle , where  = 360/N. In this way,  accumulates to 360° with N cycles. Because the angular displacement that each pulse represents is fixed, therefore, the count value N is a stable constant. Variable N stands for the resolution of the angular displacement. However, in the process of frequency measurement, the frequency of a standard signal is changeless, while the count value N varies with the difference of the measured signal. Consequently, angular displacement information cannot be output effectively. Aiming at the problems mentioned above, this paper presents a self-adaptive frequency tracking measurement algorithm. On the premise that the value of pulse counting is N, the algorithm realizes measurement of angular displacement and angular velocity by changing the frequency of the standard signal selfadaptively.

Principle of Self-Adaptive Frequency Tracking Measurement
The principle of the self-adaptive frequency tracking measurement method is shown in Figure 2. Supposing that measured signal is x(t), the standard signal is xs(t). When the pulse counting value is fixed at N, the cumulatively output angular displacement is 180°. The frequencies of the first three cycles of the measured signal are f1, f2 and f3, respectively. The corresponding frequencies of the first three cycles of the standard signal are fs0, fs1 and fs2, respectively. The half period of the measured signal x(t) is taken as gate time after processing the signal with shaping circuit modulation. Presupposing that the initial frequency of standard signal is fs0, the count value is M1 during the gate time of the first cycle of the measured signal. Then, the current angular velocity can be calculated by using Equation (3): Presupposing that the initial frequency of standard signal is fs 0 , the count value is M 1 during the gate time of the first cycle of the measured signal. Then, the current angular velocity can be calculated by using Equation (3): In order that the angular displacement can be output effectively after calculating the angular velocity, when requiring that the value of pulse counting be N within the same time period, the frequency of the standard signal x s (t) should be updated by using Equation (4): Since the count value of standard signal is M 1 under the current cycle, the angular displacement information cannot be effectively output yet. From the second cycle on, the detailed process is as follows: (1) Updating the frequency of standard signal as fs 1 , outputting angular displacement information of the first cycle in a cumulative way. Two sawtooth waves (STWs) are output in one cycle. The value of pulse counting that each STW corresponds to is N. The first STW corresponds to an angular displacement from 0 • to 180 • , and the second STW corresponds to an angular displacement from 180 • to 360 • . (2) Outputting angular velocity f 1 of the first cycle.
(3) Measuring angular velocity of the current cycle by utilizing standard signal Js 1 and Equation (3), and updating the frequency of standard signal by using Equation (4).
In this way of combining step (1), (2) with (3), the frequency of the standard signal can be updated self-adaptively, and the angular velocity and angular displacement information can be output with a mode which lags one cycle behind the measured signal. Simultaneously, the designed system outputs the refresh pulse, which represents the initial zero phase of every cycle of the measured signal. The refresh pulse can provide input attitude parameters to a subsequent inertial navigation system (INS) by cooperating with angular displacement information.

Design of Angular Velocity Sensor
System architecture and a specific design scheme of the prototype of the sensor are given in Figure 3. The sensor consists of five parts, including a data acquisition module, standard signal frequency update module, self-adaptive frequency tracking measurement module and output module. 1 2 In order that the angular displacement can be output effectively after calculating the angular velocity, when requiring that the value of pulse counting be N within the same time period, the frequency of the standard signal xs(t) should be updated by using Equation (4): Since the count value of standard signal is M1 under the current cycle, the angular displacement information cannot be effectively output yet. From the second cycle on, the detailed process is as follows: (1) Updating the frequency of standard signal as fs1, outputting angular displacement information of the first cycle in a cumulative way. Two sawtooth waves (STWs) are output in one cycle. The value of pulse counting that each STW corresponds to is N. The first STW corresponds to an angular displacement from 0° to 180°, and the second STW corresponds to an angular displacement from 180° to 360°. (2) Outputting angular velocity f1 of the first cycle.
(3) Measuring angular velocity of the current cycle by utilizing standard signal Js1 and Equation (3), and updating the frequency of standard signal by using Equation (4).
In this way of combining step (1), (2) with (3), the frequency of the standard signal can be updated self-adaptively, and the angular velocity and angular displacement information can be output with a mode which lags one cycle behind the measured signal. Simultaneously, the designed system outputs the refresh pulse, which represents the initial zero phase of every cycle of the measured signal. The refresh pulse can provide input attitude parameters to a subsequent inertial navigation system (INS) by cooperating with angular displacement information.

Design of Angular Velocity Sensor
System architecture and a specific design scheme of the prototype of the sensor are given in Figure 3. The sensor consists of five parts, including a data acquisition module, standard signal frequency update module, self-adaptive frequency tracking measurement module and output module. Figure 3. Schematic diagram of the general design scheme of the sensor.

Data Acquisition Module
Under the conditions of different rotational speeds of a high-speed rotating projectile, an electromagnetic induction coil is employed to obtain the real-time information about the projectile movement across a geomagnetic field. According to Faraday's Law of Electromagnetic Induction, a framework of a bent "I-shape" type is used to wind an uniaxial coil, whose wire diameter is 0.09 mm,

Data Acquisition Module
Under the conditions of different rotational speeds of a high-speed rotating projectile, an electromagnetic induction coil is employed to obtain the real-time information about the projectile movement across a geomagnetic field. According to Faraday's Law of Electromagnetic Induction, Sensors 2016, 16, 1625 5 of 12 a framework of a bent "I-shape" type is used to wind an uniaxial coil, whose wire diameter is 0.09 mm, number of turns is 9500, and resistance is about 665 Ω. Simultaneously, triaxial coils are installed in a mutually orthogonal way to realize the measurement of triaxial geomagnetic signals, as shown in Figure 4. number of turns is 9500, and resistance is about 665 Ω. Simultaneously, triaxial coils are installed in a mutually orthogonal way to realize the measurement of triaxial geomagnetic signals, as shown in Figure 4. The relationships between the test frequency and output voltage of the wound field coil are listed in Table 1. The linearity of the induction coil is calibrated by linear fitting. The characteristic output curve and fitting curve are shown in Figure 5.   The relationships between the test frequency and output voltage of the wound field coil are listed in Table 1. The linearity of the induction coil is calibrated by linear fitting. The characteristic output curve and fitting curve are shown in Figure 5. number of turns is 9500, and resistance is about 665 Ω. Simultaneously, triaxial coils are installed in a mutually orthogonal way to realize the measurement of triaxial geomagnetic signals, as shown in Figure 4. The relationships between the test frequency and output voltage of the wound field coil are listed in Table 1. The linearity of the induction coil is calibrated by linear fitting. The characteristic output curve and fitting curve are shown in Figure 5.    Figure 5 reveals that when employing an induction coil to cut a fixed magnetic field, the frequency of the induced magnetic signal is in a range from 4 Hz to 100 Hz. The cutting frequency and induced voltage generated show a good linear relationship with linearity of 0.5% F.S., thereby satisfying the measurement requirements of a large range of rotational speeds for a high-speed rotating projectile. Weak signals induced by the coils are amplified self-adaptively within an appropriate range by using a signal conditioning circuit, so as to output gate signals after half-wave shaping.

Self-Adaptive Frequency Tracking Measurement Module and Frequency Update Module
A digital to analog converter (DAC) is utilized to control voltage-frequency converter (VFC) and make it output a standard signal x s (t). The frequency of the standard signal is used to update Equation (4). The standard signal is updated under the control of the VFC at the initial time of every cycle of the measured signal.
The relationship between DAC and VFC is as follows. Assuming that output voltage of DAC is V i in the i-th cycle, the corresponding output frequency of VFC is f si . According to the transmission characteristics of the VFC, the connection between these two parameters above is: where a and b are linear constants of the VFC obtained through a multi-point fitting method. Combined with Equation (4), the corresponding output voltage of DAC of the next cycle is: Based on the approach that an angular velocity measurement model based on Equation (6) is built by taking advantage of FPGA, the value of DAC can be updated. Simultaneously, the count value M i is latched up, and then the angular velocity can be calculated by using Equation (3).

Output Module
The angular velocity and angular displacement information are both output in the form of voltage by utilizing the DAC. The angular displacement is output in the form of STW. Two STW are output in one cycle. When the pulse counting value reaches N, the DAC outputs a STW signal with presupposed amplitude. The full voltage of each STW corresponds to angular displacements of 180 • and 360 • , respectively. The corresponding angular displacement can be calculated by using the output voltage.

Linearity Calibration of the VFC
The sensitivity of the VFC determines the corresponding relation between the regulating voltage of the DAC and VFC. In order to realize a self-adaptive frequency tracking measurement algorithm, the sensitivity of the VFC is calibrated before experiments. An AD654 of Analog Device Instrument (ADI Company, Norwood, MA, USA) is chosen as VFC, whose relationships between input voltage and output frequency are listed in Table 2. The calibration curve of the VFC's linearity is obtained according to Table 2, as shown in Figure 6. Based on the linear fitting approach, the corresponding output characteristic equation of the VFC, which also implies the relationship between regulating voltage of the DAC and VFC is expressed as: It can be concluded that the linearity of the selected VFC is 0.052% F.S. This good linearity can satisfy the high precision demands for measuring the angular velocity of high-speed rotating projectiles.

Semi-Physical Simulation Experiment-Alternating Magnetic Field
In order to calibrate the angular velocity measurement accuracy of the sensor, simulation experiments are carried out on the basis of using a magnetic shielding device and simulative alternating magnetic field generator of Helmholtz coils. The novel sensor researched and developed based on the designed system is gently put in a simulative alternating magnetic field of Helmholtz coils as shown in Figure 7. The inclination angle of the fixed table-board is set as 30° to ensure that under the action of the alternating magnetic field, the triaxial electromagnetic induction coils can all output cutting magnetic induction signals. Furthermore, the frequency of the signal generator, which represents the rotational speed frequency of a projectile under actual working conditions, is set at a range from 5 Hz to 100 Hz. The intensity of the applied magnetic field is ±1 Gauss. The count pulse value is 1024 per half cycle, which means that N = 1024.  Based on the linear fitting approach, the corresponding output characteristic equation of the VFC, which also implies the relationship between regulating voltage of the DAC and VFC is expressed as: It can be concluded that the linearity of the selected VFC is 0.052% F.S. This good linearity can satisfy the high precision demands for measuring the angular velocity of high-speed rotating projectiles.

Semi-Physical Simulation Experiment-Alternating Magnetic Field
In order to calibrate the angular velocity measurement accuracy of the sensor, simulation experiments are carried out on the basis of using a magnetic shielding device and simulative alternating magnetic field generator of Helmholtz coils. The novel sensor researched and developed based on the designed system is gently put in a simulative alternating magnetic field of Helmholtz coils as shown in Figure 7. The inclination angle of the fixed table-board is set as 30 • to ensure that under the action of the alternating magnetic field, the triaxial electromagnetic induction coils can all output cutting magnetic induction signals. Furthermore, the frequency of the signal generator, which represents the rotational speed frequency of a projectile under actual working conditions, is set at a range from 5 Hz to 100 Hz. The intensity of the applied magnetic field is ±1 Gauss. The count pulse value is 1024 per half cycle, which means that N = 1024. Based on the linear fitting approach, the corresponding output characteristic equation of the VFC, which also implies the relationship between regulating voltage of the DAC and VFC is expressed as: It can be concluded that the linearity of the selected VFC is 0.052% F.S. This good linearity can satisfy the high precision demands for measuring the angular velocity of high-speed rotating projectiles.

Semi-Physical Simulation Experiment-Alternating Magnetic Field
In order to calibrate the angular velocity measurement accuracy of the sensor, simulation experiments are carried out on the basis of using a magnetic shielding device and simulative alternating magnetic field generator of Helmholtz coils. The novel sensor researched and developed based on the designed system is gently put in a simulative alternating magnetic field of Helmholtz coils as shown in Figure 7. The inclination angle of the fixed table-board is set as 30° to ensure that under the action of the alternating magnetic field, the triaxial electromagnetic induction coils can all output cutting magnetic induction signals. Furthermore, the frequency of the signal generator, which represents the rotational speed frequency of a projectile under actual working conditions, is set at a range from 5 Hz to 100 Hz. The intensity of the applied magnetic field is ±1 Gauss. The count pulse value is 1024 per half cycle, which means that N = 1024.  Under the condition that the frequency of the alternating magnetic field is 100 Hz, the oscillogram of the output signal of the sensor is shown in Figure 8. Under the condition that the frequency of the alternating magnetic field is 100 Hz, the oscillogram of the output signal of the sensor is shown in Figure 8. As can be seen from Figure 8, there are three types of signal output by the sensor: angular velocity signal, angular displacement signal and the refresh pulse signal. When the presupposed frequency of the alternating magnetic field is 100 Hz, the angular velocity signal is stable with small ripples of less than 40 mV. The corresponding frequency of the refresh pulse stabilizes at 100 Hz, which illustrates that the designed sensor can accurately detect the initial zero during every cycle and achieve the function of automatic zero clearing during the process of measurement. Simultaneously, two STW representing angular displacement information from 0° to 360° are output between two refresh pulses. Experiments are conducted with the preset values listed in Table 3 one after another. Under different frequencies, the corresponding angular velocity parameter information and measurement errors are acquired as follows. As Table 3 shows, along with the change of alternating frequency of the magnetic field, the sensor can realize angular velocity measurements within an extremely high dynamic range from 5 Hz to 100 Hz, with a maximum measurement error of less than 0.3%. As can be seen from Figure 8, there are three types of signal output by the sensor: angular velocity signal, angular displacement signal and the refresh pulse signal. When the presupposed frequency of the alternating magnetic field is 100 Hz, the angular velocity signal is stable with small ripples of less than 40 mV. The corresponding frequency of the refresh pulse stabilizes at 100 Hz, which illustrates that the designed sensor can accurately detect the initial zero during every cycle and achieve the function of automatic zero clearing during the process of measurement. Simultaneously, two STW representing angular displacement information from 0 • to 360 • are output between two refresh pulses. Experiments are conducted with the preset values listed in Table 3 one after another. Under different frequencies, the corresponding angular velocity parameter information and measurement errors are acquired as follows. As Table 3 shows, along with the change of alternating frequency of the magnetic field, the sensor can realize angular velocity measurements within an extremely high dynamic range from 5 Hz to 100 Hz, with a maximum measurement error of less than 0.3%.

Physical Simulation Experiment-Three-Axis Flight Turntable
In order to calibrate the angular displacement measurement accuracy of the sensor, test experiments were conducted using a three-axis position-rate-swing temperature control flight turntable. The highest test frequency of the turntable is 1000 • /s, with a rate precision error of ±0.01 • . The angular velocity sensor is fixed at the heart of the turntable, as shown in Figure 9. The turntable is controlled and adjusted to skew 30 • . The rotational test frequency of the inner frame of the turntable is set at different angular rotation frequencies.

Physical Simulation Experiment-Three-Axis Flight Turntable
In order to calibrate the angular displacement measurement accuracy of the sensor, test experiments were conducted using a three-axis position-rate-swing temperature control flight turntable. The highest test frequency of the turntable is 1000°/s, with a rate precision error of ±0.01°. The angular velocity sensor is fixed at the heart of the turntable, as shown in Figure 9. The turntable is controlled and adjusted to skew 30°. The rotational test frequency of the inner frame of the turntable is set at different angular rotation frequencies. Under the condition that the rotational test frequency of the inner frame of the turntable is 1000°/s, an oscillogram of the sensor output signal is shown in Figure 10. It can be seen from Figure 10 that when rotational test frequency of the turntable is 1000°/s, the corresponding frequency of the refresh pulse stabilizes at 2.770 Hz, which shows that the sensor can accurately detect the initial zero during every cycle and achieve the function of automatic zero clearing during the measurement. Simultaneously, the output STW of the sensor is stable, the initial zero of the STW and the refresh pulse have a one-to-one correspondence. The output DC voltage of Under the condition that the rotational test frequency of the inner frame of the turntable is 1000 • /s, an oscillogram of the sensor output signal is shown in Figure 10.

Physical Simulation Experiment-Three-Axis Flight Turntable
In order to calibrate the angular displacement measurement accuracy of the sensor, test experiments were conducted using a three-axis position-rate-swing temperature control flight turntable. The highest test frequency of the turntable is 1000°/s, with a rate precision error of ±0.01°. The angular velocity sensor is fixed at the heart of the turntable, as shown in Figure 9. The turntable is controlled and adjusted to skew 30°. The rotational test frequency of the inner frame of the turntable is set at different angular rotation frequencies. Under the condition that the rotational test frequency of the inner frame of the turntable is 1000°/s, an oscillogram of the sensor output signal is shown in Figure 10. It can be seen from Figure 10 that when rotational test frequency of the turntable is 1000°/s, the corresponding frequency of the refresh pulse stabilizes at 2.770 Hz, which shows that the sensor can accurately detect the initial zero during every cycle and achieve the function of automatic zero clearing during the measurement. Simultaneously, the output STW of the sensor is stable, the initial zero of the STW and the refresh pulse have a one-to-one correspondence. The output DC voltage of It can be seen from Figure 10 that when rotational test frequency of the turntable is 1000 • /s, the corresponding frequency of the refresh pulse stabilizes at 2.770 Hz, which shows that the sensor can accurately detect the initial zero during every cycle and achieve the function of automatic zero clearing during the measurement. Simultaneously, the output STW of the sensor is stable, the initial zero of the STW and the refresh pulse have a one-to-one correspondence. The output DC voltage of the regulation DAC is steady with small ripples, which illustrates that the sensor has stable performance, and the obtained angular displacement is accurate.
The voltage information of the angular displacement signal, obtained through Figure 10, is converted into angle information (as shown in Figure 11a) by integrating with a time horizon of 1 s (as shown in Figure 11c). The rotational frequency of the turntable is 1000 • /s, which means that the angular displacement accumulates 1000 • in 1 s. When integrating the actual angular displacement signal, the error is 0.1 • , thus satisfying the measurement requirements. the regulation DAC is steady with small ripples, which illustrates that the sensor has stable performance, and the obtained angular displacement is accurate. The voltage information of the angular displacement signal, obtained through Figure 10, is converted into angle information (as shown in Figure 11a) by integrating with a time horizon of 1 s (as shown in Figure 11c). The rotational frequency of the turntable is 1000°/s, which means that the angular displacement accumulates 1000° in 1 s. When integrating the actual angular displacement signal, the error is 0.1°, thus satisfying the measurement requirements.  Table 4. On the basis of the analysis of Table 4, it is concluded that when the inner frame of the turntable rotates at 600, 900 and 1000°/s, respectively, the measurement errors of the angular displacement acquired through the sensor are all less than 0.2°, thereby satisfying the actual angular displacement measurement requirements for a high-speed rotating projectile.

Conclusions
Aiming at meeting the demands for measuring angular motion parameters of high-speed rotating projectiles, this paper puts forth a new method for measuring the angular velocity and angular displacement parameters of the rotation of a projectile based on the electromagnetic induction principle, and also designs the corresponding sensor. To begin with, triaxial electromagnetic induction coils are used to acquire the cutting geomagnetic signals of the projectile. Next, a self-adaptive frequency tracking measurement algorithm, which combines pulse counting with a way of updating the frequency of a standard signal self-adaptively, is proposed by utilizing DAC and VFC to calculate the angular velocity and angular displacement information, and a principle The corresponding experimental test results are given in Table 4. On the basis of the analysis of Table 4, it is concluded that when the inner frame of the turntable rotates at 600, 900 and 1000 • /s, respectively, the measurement errors of the angular displacement acquired through the sensor are all less than 0.2 • , thereby satisfying the actual angular displacement measurement requirements for a high-speed rotating projectile.

Conclusions
Aiming at meeting the demands for measuring angular motion parameters of high-speed rotating projectiles, this paper puts forth a new method for measuring the angular velocity and angular displacement parameters of the rotation of a projectile based on the electromagnetic induction principle, and also designs the corresponding sensor. To begin with, triaxial electromagnetic induction coils are used to acquire the cutting geomagnetic signals of the projectile. Next, a self-adaptive frequency tracking measurement algorithm, which combines pulse counting with a way of updating the frequency of a standard signal self-adaptively, is proposed by utilizing DAC and VFC to calculate the angular velocity and angular displacement information, and a principle prototype of the sensor is designed based on FPGA. Third, angular velocity performance tests are carried out by using a simulative alternating magnetic field generator of Helmholtz coils, within the alternative frequency range from 5 rps to 100 rps. The measurement errors are all less than 0.3%. Finally, angular displacement performance tests are conducted by utilizing a three-axis position-rate-swing temperature control turntable, in a range from 1 rps to 3 rps. The angular displacement measurement errors are all lower than 0.2 • . The experimental results both show that the sensor can measure angular motion information within a range from 1 rps to 100 rps, and can realize real-time measurements of the angular motion parameters of a high-speed rotating projectile with high precision.
Two research topics are scheduled to be developed in the near future. On the one hand, the final purpose of the sensor is to provide motion parameters for attitude determination of high-speed rotating projectiles and geomagnetic navigation. How to construct a navigation model and match with the interface of the current navigation system would be a direction for further research. On the other hand, the stability and balance of the self-adaptive frequency tracking measurement model should be analyzed to apply the available sensor for whole trajectory measurement of high-speed rotating projectiles.