CN115824198A - Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model - Google Patents

Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model Download PDF

Info

Publication number
CN115824198A
CN115824198A CN202211288342.9A CN202211288342A CN115824198A CN 115824198 A CN115824198 A CN 115824198A CN 202211288342 A CN202211288342 A CN 202211288342A CN 115824198 A CN115824198 A CN 115824198A
Authority
CN
China
Prior art keywords
temperature
cold
state
model
hot
Prior art date
Legal status (The legal status 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 status listed.)
Pending
Application number
CN202211288342.9A
Other languages
Chinese (zh)
Inventor
卓超
任磊
刘晴晴
郝芮
张潇
杜建邦
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Aerospace Automatic Control Research Institute
Original Assignee
Beijing Aerospace Automatic Control Research Institute
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 Beijing Aerospace Automatic Control Research Institute filed Critical Beijing Aerospace Automatic Control Research Institute
Priority to CN202211288342.9A priority Critical patent/CN115824198A/en
Publication of CN115824198A publication Critical patent/CN115824198A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Gyroscopes (AREA)

Abstract

The invention provides a method for selecting and switching an optical gyro cold and hot start temperature error compensation model, which obtains the temperature space gradient among different positions of a gyro at the moment of starting the inertial navigation through a plurality of temperature sensors distributed in a gyro instrument area in the inertial navigation. The inertial navigation 'cold start' and 'hot start' are considered as two populations, and the 'mahalanobis distance' from the sample set formed by the spatial gradient to the two populations is calculated. The method can judge the temperature state of inertial navigation in a very short time when the inertial navigation is started, select a correct error model to compensate the output of the gyroscope, can complete the transition switching of the compensation model when the gyroscope enters a hot-state long-time stable working stage from a cold-state starting stage, and is compatible with different working states of cold-state starting and hot-state starting. The method has the advantages of small calculated amount, simplicity, practicability and convenience for engineering realization.

Description

Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model
Technical Field
The application relates to the technical field of inertial navigation, in particular to a method for selecting and switching an optical gyroscope cold and hot start temperature error compensation model.
Background
The optical gyroscope mainly comprises an optical fiber gyroscope and a laser gyroscope, has the outstanding advantages of wide angular velocity measurement range, overload impact resistance, short starting time and the like compared with a classical rotor gyroscope, and is particularly suitable for a strapdown inertial navigation system. As a core instrument of the inertial navigation system, the performance of a gyroscope directly determines the navigation precision of the inertial navigation system. However, the output of the optical gyro is affected by temperature, and the generated instrument temperature error seriously restricts the improvement of system alignment, navigation precision and quick response performance. The solution commonly adopted at present is temperature error compensation. Temperature error compensation is essentially a mathematical modeling based method that is excited by a temperature test and separates the temperature error from the gyro output, and models the error and temperature using their correlation and deducts them from the system output. However, in the specific implementation process, the temperature error behavior mode of the gyro is greatly different under two conditions of starting from a cold environment (the gyro is powered off and cooled for a long time before the start and is approximately in thermal balance with the environment) and starting from a hot environment (the gyro is powered on for a long time before the start and the internal temperature is higher than the ambient temperature).
The gyro temperature error model obtained through a long-time full-temperature range temperature test is basically a hot-state model, and cannot accurately describe the error mode in the cold-state starting stage of the starting, so that the targeted modeling needs to be carried out on the short process of the cold-state starting of the gyro. The gyro temperature error models established under the two conditions cannot be used universally, so that the temperature state of inertial navigation needs to be accurately judged in a very short time when the inertial navigation is started, so that the correct error model can be selected to compensate gyro output, and the model transition switching from a cold-state starting stage to a hot-state long-time stable working stage of the gyro can be completed.
At present, most of public researches are focused on gyro temperature error compensation in a long-time working stage of inertial navigation, the attention on the gyro cold-state starting condition is less, and a selection and switching method of cold and hot starting models for different working stages is hardly involved.
Disclosure of Invention
In order to solve the technical problem, the invention provides a method for selecting and switching an optical gyroscope cold and hot start temperature error compensation model. The technical scheme adopted by the invention is as follows:
a method for selecting and switching an optical gyroscope cold and hot start temperature error compensation model comprises the following steps:
step 1, distributing a plurality of temperature sensors at different positions of different gyros of inertial navigation, and setting the total distribution number to be n;
step 2, calculating the temperature space gradient by using the absolute temperature difference values measured by the temperature sensors at different positions;
step 3, respectively carrying out n 1 Group inertial navigation cold environment starting test and n 2 Taking the group inertial navigation hot-state starting test as a training data set, taking two environments of a cold state and a hot state as two universes, and recording the universes as G 1 And G 2 (ii) a Calculating estimates of ensemble mean of temperature gradients of two ensembles separately
Figure BDA0003900282370000021
And covariance estimation
Figure BDA0003900282370000022
Step 4, judging cold state and hot state starting based on a judgment threshold TSH according to the Mahalanobis distance overall judgment criterion;
step 5, if the hot start is judged through the step 4, directly selecting a hot start temperature error compensation model Mod 2 Compensating gyro output, mod 2 The error compensation amount calculated in real time for the pre-established thermal state temperature error model can be expressed as
Figure BDA0003900282370000031
Wherein T is g (t)=[T g1 (t)T g2 (t)...T gn (t)] T Is the real-time temperature collected by the gyro temperature sensor,
Figure BDA0003900282370000032
a thermal state temperature error compensation quantity function established according to the correlation between the error and the temperature;
if the cold start is judged through the step 4, the inertial navigation working time t is less than t G Internal cold-state starting temperature error compensation model Mod 1 Compensating gyro output, mod 1 The error compensation amount calculated in real time for the pre-established cold state temperature error model can be expressed as
Figure BDA0003900282370000033
Figure BDA0003900282370000034
Representing a cold temperature error compensation quantity function, t, established on the basis of the error-temperature dependence G The applicable duration of the model is started in a cold state; at t G ≤t≤t T Internal arrangement of the transition process, t T For the arranged model transition duration, smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using a transition function; when t is>t T Then, a thermal state starting temperature error compensation model Mod is adopted 2 And compensating the output of the gyroscope.
Further, in step 1, the disposing a plurality of temperature sensors at different positions of different gyros in inertial navigation includes:
1 temperature sensor is respectively arranged on the gyroscope in the X, Y, Z three directions, namely n =3; or
For the laser gyroscope, a temperature sensor is respectively arranged on each gyroscope anode and the shell, namely n =6; or
For the fiber-optic gyroscope, a temperature sensor is respectively arranged at the position of each gyroscope close to the optical fiber and close to the surface of the shell, namely n =6.
Further, in step 2, the calculating the temperature spatial gradient includes:
let the real-time temperature measurement values of n temperature sensors be T gi (t), i =1, …, n, the temperature gradient between the two temperature sensors is expressed as:
Figure BDA0003900282370000041
through the above calculation formula, the maximum value is obtained by calculation
Figure BDA0003900282370000042
A spatial gradient of temperature.
Further, in step 2, after the calculating the temperature spatial gradient, the method further includes:
according to the actual situation, the calculation results
Figure BDA0003900282370000043
K of the temperature space gradients are selected
Figure BDA0003900282370000044
Typical gradient composition measurement sample Δ T g (T), the measurement sample Δ T g (t) is expressed as:
Figure BDA0003900282370000045
recording inertial navigation starting time T =0 and corresponding to temperature gradient delta T g (T) smoothing to obtain a temperature gradient smoothing value Delta T gs Let the smoothing time be Δ t, expressed as:
Figure BDA0003900282370000046
the value Δ t =1s since the absolute temperature noise is small.
Further, the estimation of the ensemble mean of the temperature gradients of the two ensembles is calculated separately
Figure BDA0003900282370000047
And
Figure BDA0003900282370000048
the method comprises the following steps:
smoothing the value Δ T according to the test results and the temperature gradient gs Under the two population conditions, the temperature gradient elements belonging to the two populations are respectively expressed as
Figure BDA0003900282370000049
Figure BDA00039002823700000410
Estimation of the ensemble mean of the temperature gradients of two ensembles
Figure BDA00039002823700000411
And
Figure BDA00039002823700000412
calculated by the following formula:
Figure BDA0003900282370000051
Figure BDA0003900282370000052
further, the covariance estimation
Figure BDA0003900282370000053
Calculated by the following formula:
Figure BDA0003900282370000054
Figure BDA0003900282370000055
Figure BDA0003900282370000056
further, the discrimination threshold TSH is expressed as
Figure BDA0003900282370000057
Wherein the content of the first and second substances,
Figure BDA0003900282370000058
further, in step 4, the determining cold and hot starts based on the determination threshold TSH includes:
real-time temperature space gradient delta T obtained by using the discrimination threshold value to inertial navigation starting time gs Performing on-line discrimination if the function omega (delta T) gs ):
ω(ΔT gs )=a T ·ΔT gs ≥TSH
Judging that the inertial navigation is in a cold starting state; otherwise, judging that the inertial navigation is in a hot starting state.
Further, the smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using the transition function includes:
at t G ≤t≤t T Simultaneous cold state model Mod calculation in time period 1 And thermal state model Mod 2 Real-time compensation amount Δ M 1 (t) and Δ M 2 (t); smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using a transition function F (t), wherein the compensation quantity is represented as follows:
Figure BDA0003900282370000061
ΔM(t)=[1-F(t)]ΔM 1 (t)+F(t)ΔM 2 (t)t G ≤t≤t T
f (t) is a transition function, and delta M (t) is a transition process temperature error compensation quantity fusing a cold state model and a hot state model.
Further, in step 5, t G And t T The selection needs to be carried out according to the duration of the modeling data set of the cold-state starting model of the inertial navigation system and the temperature characteristic of the inertial navigation.
Through the embodiment of the application, the following technical effects can be obtained:
(1) The method can accurately judge the cold state and the hot state starting of the inertial navigation in a very short time when the inertial navigation is started, select a correct error model to compensate the output of the gyroscope, can complete the natural smooth switching of the compensation model from the cold state starting to the hot state stable working stage of the gyroscope, and is compatible with different working states of the cold state starting and the hot state starting.
(2) The method adopts the spatial gradient distributed in the gyro temperature field as the distinguishing mark of the cold and hot states, does not need to calculate the temperature time change rate, and avoids the influence of noise introduced by obtaining the first-order derivative of the temperature on the judgment accuracy and the starting delay introduced by smooth differential noise.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings needed to be used in the embodiments or the prior art descriptions will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present application, and those skilled in the art can also obtain other drawings according to the drawings without inventive labor.
FIG. 1 is a schematic diagram of a temperature error compensation model selection and switching process according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are some embodiments of the present application, but not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
FIG. 1 is a schematic diagram of a temperature error compensation model selection and switching process according to the present invention. The technology of the invention solves the problem that the temperature space gradient between different positions of the gyroscope at the starting moment of inertial navigation is obtained by a plurality of temperature sensors distributed in a gyroscope instrument area in the inertial navigation. The inertial navigation "cold start" and "hot start" are considered as two populations, and the mahalanobis distance of the set of samples formed by the spatial gradient to the two populations is calculated. The method comprises the steps of training cold-state starting and hot-state starting temperature tests in advance for multiple times to generate a distance discrimination threshold value, taking the distance discrimination threshold value as a discrimination basis for cold-state and hot-state starting, carrying out online discrimination on a temperature space gradient obtained at the moment of inertial navigation starting by using the criterion, and directly selecting a hot-state starting temperature error compensation model to compensate gyro output if the judgment result is that the gyro is hot-state starting; if the starting is judged to be cold, the working time t of inertial navigation is less than t G Compensating the gyro output by using a cold start temperature error compensation model internally (t) G The applicable duration of the cold start model). Then, at t G ≤t≤t T Arranging the transition process inside (t) T For a scheduled model transition duration), smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using a transition function. When t is>t T And then, compensating the output of the gyroscope by adopting a thermal state starting temperature error compensation model. The method for selecting and switching the optical gyro cold and hot start temperature error compensation model based on spatial gradient discrimination utilizes the temperature spatial gradient measured by different gyro temperature sensors to calculate the Mahalanobis distance so as to realize accurate discrimination of inertial navigation cold state and hot state start, and arranges a transition process to realize natural switching of cold state and hot state models, as shown in FIG. 1.
The temperature error compensation model selecting and switching method comprises the following steps:
step 1, distributing a plurality of temperature sensors at different positions of different gyros of inertial navigation, and setting the total distribution number to be n.
In step 1, the arranging a plurality of temperature sensors at different positions of different gyros of inertial navigation includes:
1 temperature sensor is respectively arranged on the gyroscope in the X, Y, Z three directions, namely n =3;
for the laser gyroscope, a temperature sensor is respectively arranged on each gyroscope anode and the shell, namely n =6;
for the fiber-optic gyroscope, a temperature sensor is respectively arranged at the position of each gyroscope close to the optical fiber and close to the surface of the shell, namely n =6;
step 2, calculating the temperature space gradient by using the absolute temperature difference values measured by the temperature sensors at different positions;
in step 2, the calculating the temperature spatial gradient includes:
let the real-time temperature measurement values of n temperature sensors be T gi (t), i =1, …, n, the temperature gradient between the two temperature sensors is expressed as:
Figure BDA0003900282370000081
through the above calculation formula, the maximum value is obtained by calculation
Figure BDA0003900282370000082
A temperature spatial gradient;
in step 2, the calculating the temperature spatial gradient further comprises:
according to the actual situation, obtained from said calculation
Figure BDA0003900282370000083
K of the temperature space gradients are selected
Figure BDA0003900282370000084
Typical gradient composition measurement sample Δ T g (T), the measurement sample Δ T g (t) is expressed as:
Figure BDA0003900282370000091
recording inertial navigation starting time T =0 and corresponding to temperature gradient delta T g (T) smoothing to obtain a temperature gradient smoothing value Delta T gs Let the smoothing time be Δ t, expressed as:
Figure BDA0003900282370000092
the value Δ t =1s because the absolute temperature noise is small;
step 3, respectively carrying out n 1 Group inertial navigation cold environment starting test and n 2 Taking the group inertial navigation hot-state starting test as a training data set, taking two environments of a cold state and a hot state as two universes, and recording the universes as G 1 And G 2 (ii) a Calculating estimates of the mean of two populations of temperature gradients
Figure BDA0003900282370000093
And covariance estimation
Figure BDA0003900282370000094
The estimation of the ensemble mean of the temperature gradients of the two ensembles is calculated separately
Figure BDA0003900282370000095
And
Figure BDA0003900282370000096
the method comprises the following steps:
smoothing the value Δ T according to the test results and the temperature gradient gs Under the two population conditions, the temperature gradient elements belonging to the two populations are respectively expressed as
Figure BDA0003900282370000097
Figure BDA0003900282370000098
Estimation of the ensemble mean of the temperature gradients of two ensembles
Figure BDA0003900282370000099
And
Figure BDA00039002823700000910
calculated by the following formula:
Figure BDA00039002823700000911
Figure BDA00039002823700000912
the covariance estimate
Figure BDA0003900282370000101
Calculated by the following formula:
Figure BDA0003900282370000102
Figure BDA0003900282370000103
Figure BDA0003900282370000104
step 4, judging cold state and hot state starting based on a judgment threshold TSH according to the Mahalanobis distance overall judgment criterion;
the discrimination threshold TSH is expressed as
Figure BDA0003900282370000105
Wherein the content of the first and second substances,
Figure BDA0003900282370000106
in step 4, the cold state and the hot state start are determined based on the determination threshold TSH, and the method includes:
real-time temperature space gradient delta T obtained by using the discrimination threshold value to inertial navigation starting time gs Performing on-line discrimination, if the function omega (delta T) is discriminated gs ):
ω(ΔT gs )=a T ·ΔT gs ≥TSH
Judging that the inertial navigation is in a cold starting state; otherwise, judging that the inertial navigation is in a hot starting state;
step 5, if the hot start is judged through the step 4, directly selecting a hot start temperature error compensation model Mod 2 Compensating gyro output, mod 2 The error compensation amount calculated in real time for the pre-established thermal state temperature error model can be expressed as
Figure BDA0003900282370000107
Wherein T is g (t)=[T g1 (t)T g2 (t)...T gn (t)] T Is the real-time temperature collected by the gyro temperature sensor,
Figure BDA0003900282370000108
in order to establish a thermal state temperature error compensation quantity function according to the correlation between the error and the temperature, a polynomial can be specifically adopted for modeling: for example, if a single temperature sensor is used,
Figure BDA0003900282370000109
can be expressed as
Figure BDA0003900282370000111
(a 0 ,a 1 ,a 2 As model coefficients);
if the cold start is judged through the step 4, the inertial navigation working time t is less than t G Compensating the gyro output by using a cold start temperature error compensation model internally, t G The applicable duration of the model is started in a cold state; at t G ≤t≤t T Arranging the transition process in, t T For the arranged model transition duration, smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using a transition function; when t is>t T And then, compensating the output of the gyroscope by adopting a thermal state starting temperature error compensation model.
The compensating the gyro output by using the cold-start temperature error compensation model comprises the following steps: using cold start temperature error compensation model Mod 1 Compensating gyro output, mod 1 The error compensation amount calculated in real time for the pre-established cold state temperature error model can be expressed as
Figure BDA0003900282370000112
Figure BDA0003900282370000113
The function of the cold-state temperature error compensation quantity established according to the correlation between the error and the temperature is shown, and a polynomial can be used for modeling: if a single temperature measuring sensor is used,
Figure BDA0003900282370000114
can be expressed as
Figure BDA0003900282370000115
(b 0 ,b 1 ,b 2 As model coefficients).
The step of smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using the transition function comprises the following steps:
at t G ≤t≤t T Simultaneous cold state model Mod calculation in time period 1 With thermal state model Mod 2 Real-time compensation amount Δ M 1 (t) and Δ M 2 (t); using transitionsThe function F (t) smoothly and gradually switches the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model, and is represented as:
Figure BDA0003900282370000116
ΔM(t)=[1-F(t)]ΔM 1 (t)+F(t)ΔM 2 (t)t G ≤t≤t T
f (t) is a transition function, and delta M (t) is a transition process temperature error compensation quantity fusing a cold state model and a hot state model;
the adoption hot-state starts temperature error compensation model and compensates gyro output includes: adopt hot-state start-up temperature error compensation model Mod 2 And compensating the output of the gyroscope.
In the concrete implementation process, t G And t T The selection needs to be carried out according to the duration of the modeling data set of the cold-state starting model of the inertial navigation system and the temperature characteristic of the inertial navigation.
By adopting the scheme, the method can accurately judge the cold state and the hot state starting of the inertial navigation in a very short time when the inertial navigation is started, select a correct error model to compensate the output of the gyroscope, can complete the natural smooth switching of the compensation model from the cold state starting to the hot state stable working stage of the gyroscope, and is compatible with different working states of the cold state starting and the hot state starting.
The functions described above in this application may be performed at least in part by one or more hardware logic components. For example, without limitation, exemplary types of hardware logic components that may be used include: a Field Programmable Gate Array (FPGA), an Application Specific Integrated Circuit (ASIC), an Application Specific Standard Product (ASSP), a system on a chip (SOC), a load programmable logic device (CPLD), and the like.
Further, while operations are depicted in a particular order, this should be understood as requiring that such operations be performed in the particular order shown or in sequential order, or that all illustrated operations be performed, to achieve desirable results. Under certain circumstances, multitasking and parallel processing may be advantageous. Likewise, while several specific implementation details are included in the above discussion, these should not be construed as limitations on the scope of the disclosure. Certain features that are described in the context of separate embodiments can also be implemented in combination in a single implementation. Conversely, various features that are described in the context of a single implementation can also be implemented in multiple implementations separately or in any suitable subcombination.
Although the subject matter has been described in language specific to structural features and/or logical acts of devices, it is to be understood that the subject matter defined in the appended claims is not necessarily limited to the specific features or acts described above. Rather, the specific features and acts described above are disclosed as example forms of implementing the claims.

Claims (10)

1. A method for selecting and switching an optical gyro cold and hot start temperature error compensation model is characterized by comprising the following steps:
step 1, distributing a plurality of temperature sensors at different positions of different gyros of inertial navigation, and setting the total distribution number to be n;
step 2, calculating a temperature space gradient by using absolute temperature differences measured by temperature sensors at different positions;
step 3, respectively carrying out n 1 Group inertial navigation cold environment starting test and n 2 Taking the group inertial navigation hot-state starting test as a training data set, taking two environments of a cold state and a hot state as two universes, and recording the universes as G 1 And G 2 (ii) a Calculating estimates of the mean of two populations of temperature gradients
Figure FDA0003900282360000011
And covariance estimation
Figure FDA0003900282360000012
Step 4, judging cold state and hot state starting based on a judgment threshold TSH according to the Mahalanobis distance overall judgment criterion;
step 5, if the hot start is judged through the step 4, directly selecting a hot start temperature error compensation model Mod 2 Compensating gyro output, mod 2 The error compensation amount calculated in real time for the pre-established thermal state temperature error model can be expressed as
Figure FDA0003900282360000013
Wherein T is g (t)=[T g1 (t) T g2 (t) ... T gn (t)] T Is the real-time temperature collected by the gyro temperature sensor,
Figure FDA0003900282360000014
a thermal state temperature error compensation quantity function established according to the correlation between the error and the temperature;
if the cold start is judged through the step 4, the inertial navigation working time t is less than t G Internal-use cold-start temperature error compensation model Mod 1 Compensating gyro output, mod 1 The error compensation obtained in real time for the pre-established cold state temperature error model can be expressed as
Figure FDA0003900282360000015
Figure FDA0003900282360000016
Representing a cold temperature error compensation function, t, based on the error-temperature dependence G The applicable duration of the model is started in a cold state; at t G ≤t≤t T Arranging the transition process in, t T For the arranged model transition duration, smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using a transition function; when t is>t T Then, a hot-state starting temperature error compensation model Mod is adopted 2 And compensating the output of the gyroscope.
2. The method of claim 1, wherein in step 1, the disposing a plurality of temperature sensors at different positions of different gyros of inertial navigation comprises:
1 temperature sensor is respectively arranged on the gyroscope in X, Y, Z three directions, namely n =3 at the moment; or
For the laser gyroscope, a temperature sensor is respectively arranged on each gyroscope anode and the shell, namely n =6; or alternatively
For the fiber-optic gyroscope, a temperature sensor is respectively arranged at the position of each gyroscope close to the optical fiber and close to the surface of the shell, namely n =6.
3. The method of claim 1, wherein in step 2, said calculating a temperature spatial gradient comprises:
let the real-time temperature measurement values of n temperature sensors be T gi (t), i =1, …, n, the temperature gradient between the two temperature sensors is expressed as:
Figure FDA0003900282360000021
through the above calculation formula, the maximum value is obtained by calculation
Figure FDA0003900282360000022
A spatial gradient of temperature.
4. The method of claim 3, wherein in step 2, the calculating the temperature spatial gradient further comprises:
according to the actual situation, the calculation results
Figure FDA0003900282360000023
K of the temperature space gradients are selected
Figure FDA0003900282360000024
Typical gradient composition measurement sample Δ T g (T), the measurement sample Δ T g (t) is expressed as:
Figure FDA0003900282360000025
recording inertial navigation starting time T =0 and corresponding to temperature gradient delta T g (T) smoothing to obtain a temperature gradient smoothing value Delta T gs Let the smoothing time be Δ t, expressed as:
Figure FDA0003900282360000031
the value Δ t =1s since the absolute temperature noise is small.
5. The method of claim 1, wherein the separately calculating estimates of the ensemble mean of the temperature gradients of the two ensembles
Figure FDA0003900282360000032
And
Figure FDA0003900282360000033
the method comprises the following steps:
smoothing the value Δ T according to the test results and the temperature gradient gs Under the two population conditions, the temperature gradient elements belonging to the two populations are respectively expressed as
Figure FDA0003900282360000034
Figure FDA0003900282360000035
Estimation of the ensemble mean of the temperature gradients of two ensembles
Figure FDA0003900282360000036
And
Figure FDA0003900282360000037
calculated by the following formula:
Figure FDA0003900282360000038
Figure FDA0003900282360000039
6. the method of claim 5, wherein the covariance estimate
Figure FDA00039002823600000310
Calculated by the following formula:
Figure FDA0003900282360000041
Figure FDA0003900282360000042
Figure FDA0003900282360000043
7. the method according to claim 1, characterized in that said discrimination threshold TSH is expressed as
Figure FDA0003900282360000044
Wherein the content of the first and second substances,
Figure FDA0003900282360000045
8. the method according to one of claims 1 or 7, wherein determining cold and hot starts based on a discrimination threshold TSH in step 4 comprises:
real-time temperature space gradient delta T obtained by using the discrimination threshold value to inertial navigation starting time gs Performing on-line discrimination, if the function omega (delta T) is discriminated gs ):
ω(ΔT gs )=a T ·ΔT gs ≥TSH
Judging that the inertial navigation is in a cold starting state; otherwise, judging that the inertial navigation is in a hot starting state.
9. The method of claim 1, wherein smoothly and gradually switching the compensation amount output by the cold start temperature error compensation model to the compensation amount output by the hot start model by using a transition function comprises:
at t G ≤t≤t T Method for simultaneously calculating cold state model Mod in time period 1 And thermal state model Mod 2 Real-time compensation amount Δ M 1 (t) and; smoothly and gradually switching the compensation quantity output by the cold-state starting temperature error compensation model to the compensation quantity output by the hot-state starting model by using a transition function F (t), wherein the compensation quantity is represented as follows:
Figure FDA0003900282360000046
ΔM(t)=[1-F(t)]ΔM 1 (t)+F(t)ΔM 2 (t) t G ≤t≤t T
f (t) is a transition function, and delta M (t) is a transition process temperature error compensation quantity fusing a cold state model and a hot state model.
10. Method according to claim 9, characterized in that in step 5, t is G And t T And selecting according to the duration of a cold-state starting model modeling data set of the inertial navigation system and the temperature characteristic of the inertial navigation.
CN202211288342.9A 2022-10-20 2022-10-20 Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model Pending CN115824198A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211288342.9A CN115824198A (en) 2022-10-20 2022-10-20 Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211288342.9A CN115824198A (en) 2022-10-20 2022-10-20 Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model

Publications (1)

Publication Number Publication Date
CN115824198A true CN115824198A (en) 2023-03-21

Family

ID=85525135

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211288342.9A Pending CN115824198A (en) 2022-10-20 2022-10-20 Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model

Country Status (1)

Country Link
CN (1) CN115824198A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116026328A (en) * 2023-03-28 2023-04-28 航天星云(北京)科技有限公司 Construction method and compensation method of zero-bias hysteresis effect compensation model of micro inertial navigation

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116026328A (en) * 2023-03-28 2023-04-28 航天星云(北京)科技有限公司 Construction method and compensation method of zero-bias hysteresis effect compensation model of micro inertial navigation

Similar Documents

Publication Publication Date Title
CN107330149B (en) MIMU gyro random drift prediction method based on ARMA and BPNN combined model
US8947029B2 (en) Method and device for ascertaining the rotor temperature of a permanent-magnet synchronous machine
CN101115979B (en) Electronic thermometer
CN106370173B (en) A kind of modeling of lasergyro model of temperature compensation and verification method
CN109188026B (en) Automatic calibration deep learning method suitable for MEMS accelerometer
CN105866504B (en) A kind of optical fiber current mutual inductor temperature-compensation method based on Kalman filtering
CN115824198A (en) Method for selecting and switching optical gyroscope cold and hot start temperature error compensation model
CN101949710A (en) Rapid online dynamic calibration method for zero offset of GNSS (Global Navigation Satellite System) auxiliary MEMS (Micro Electro Mechanical Systems) inertial sensor
Abdelghani et al. Comparison study of subspace identification methods applied to flexible structures
CN110659755B (en) Modeling method, apparatus and storage medium for predicting motor temperature
KR20080081027A (en) Temperature estimation method and device
CN114061559B (en) Compensation method, system and computer storage medium for zero offset drift of fiber optic gyroscope
JP4693516B2 (en) Weight measurement method
CN106017511A (en) Method for testing and compensating temperature coefficient of fiber-optic gyroscope
CN109580033A (en) A kind of concrete dam distributed optical fiber temperature measurement data error compensation method
CN108225374A (en) A kind of Allan methods of analysis of variance of blending inheritance algorithm
Yuan et al. Statistical modeling of random walk errors for triaxial rate gyros
Xiong et al. Seamless global positioning system/inertial navigation system navigation method based on square-root cubature Kalman filter and random forest regression
CN110726852A (en) MEMS accelerometer temperature compensation method
CN113432612A (en) Navigation method, device and system of flying object
CN116026328A (en) Construction method and compensation method of zero-bias hysteresis effect compensation model of micro inertial navigation
KR101181755B1 (en) Apparatus and method for calculating error of inertial measurement unit
CN112525219B (en) Method and system for judging starting time of fiber-optic gyroscope
CN114001759A (en) Array type MEMS sensor control method and system
De Cock et al. Subspace system identification for mechanical engineering

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination