Next Article in Journal
Three-Dimensional Microscopic Image Reconstruction Based on Structured Light Illumination
Next Article in Special Issue
Doppler Measurement of Modulated Light for High Speed Vehicles
Previous Article in Journal
Automatic Inside Point Localization with Deep Reinforcement Learning for Interactive Object Segmentation
Previous Article in Special Issue
Optical Fibre-Based Sensors for Oil and Gas Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation and Error Analysis for Optomechanical Inertial Sensors

Department of Aerospace Engineering, Texas A&M University, 3141 TAMU, College Station, TX 77843-3141, USA
*
Author to whom correspondence should be addressed.
Sensors 2021, 21(18), 6101; https://doi.org/10.3390/s21186101
Submission received: 4 August 2021 / Revised: 31 August 2021 / Accepted: 4 September 2021 / Published: 11 September 2021

Abstract

:
A sensor model and methodology to estimate the forcing accelerations measured using a novel optomechanical inertial sensor with the inclusion of stochastic bias and measurement noise processes is presented. A Kalman filter for the estimation of instantaneous sensor bias is developed; the outputs from this calibration step are then employed in two different approaches for the estimation of external accelerations applied to the sensor. The performance of the system is demonstrated using simulated measurements and representative values corresponding to a bench-tested 3.76 Hz oscillator. It is shown that the developed methods produce accurate estimates of the bias over a short calibration step. This information enables precise estimates of acceleration over an extended operation period. These results establish the feasibility of reliably precise acceleration estimates using the presented methods in conjunction with state of the art optomechanical sensing technology.

1. Introduction

Accelerometers are used in countless applications including: satellite, aircraft, robotic, and automotive inertial navigation. Particularly, in the realm of satellite design these devices are used both as navigation equipment [1,2,3,4] and as scientific instrumentation [5,6,7,8,9,10,11]. In both cases the accuracy of the estimated acceleration directly impacts the performance of the system and its ability to meet critical mission requirements. There is presently a vast array of commercially available inertial sensors, developed using a diverse set of physical principles, providing a wide range of precision and sensitivity [12]. Notably, advances in optomechanics have facilitated the development of ultra-sensitive acceleration sensors with a comparatively small form factor [13]. The quality factor, a unitless parameter indicating the sharpness of the resonance peak of the oscillator [14], is of particular interest in quantifying the performance of these sensors. Previous studies indicate that such systems are capable of attaining an exceptionally high quality factor over a range of high and low frequencies, providing a broad dynamic range for acceleration sensing [15,16,17]. The use of highly accurate 3D manufacturing solutions has resulted in sensors with precise operating parameters that are uniquely suited to traditional estimation techniques without reliance on complex error models.
Inertial sensors do not directly observe forcing accelerations or torques, but rather the forces are estimated by assimilating the measurements of the dynamic response of a mechanical system strapped down to the system whose accelerations are to be measured. The fidelity of the sensor model directly impacts the quality of the estimate, and in turn the performance of the navigation system [18,19,20]. Sensor models can never perfectly characterize the behavior of the system and always possess some degree of uncertainty in the established parameters; for this reason these models are intrinsically stochastic. Since the mechanical system’s dynamical response to imposed forces is correlated with its own dynamics, which can be non-linear, the observed imposed accelerations are colored by the mechanical system dynamics. Stochastic processes used in the sensor modeling aim at decorrelating the sensor dynamics and associated modeling errors to estimate the imposed acceleration. Characterization of the sensor model’s response to these stochastic terms and their impact on performance is an essential component of practical sensor design. This analysis is routinely conducted for angular rate sensors by incorporating external measurements from attitude sensors into a Kalman filter to estimate the gyro biases and their associated steady-state statistics [21,22,23,24]. This analysis has been demonstrated to significantly improve performance in the area of spacecraft attitude determination. Similar procedures are required for estimation of accelerometer biases; namely some source of external information must be provided intermittently to conduct a calibration. This can be done using absolute position estimates from sources, such as GPS [25,26,27], or attitude estimates [28,29,30].
In this paper, a sensor model and estimators for use with state of the art optomechanical inertial sensors [16,17] are presented. Firstly, a method for sensor calibration using a generalized input acceleration is presented using a discrete time Kalman filter to estimate the instantaneous sensor bias. The information from this procedure is then fed forward to the operational phase where estimates of external accelerations are obtained. Two estimation approaches are presented for the task of estimating acceleration, each of which is compatible with the information obtained in the calibration step. These two phases constitute the operational flow of the sensor’s estimation program. Results are presented for each of these filters using representative oscillator parameters and simulated measurements in the presence of a sinusoidal external acceleration. It is shown that the calibration filter is able to provide accurate estimates of the instantaneous bias in the presence of an imperfect generalized force over a short window of time. Additionally, the operational filter (for estimates of acceleration) provides precise estimates of external acceleration whose errors follow the bias over extended periods of time. The bias contribution to the acceleration estimate error is well expressed by the associated estimate error variances.

2. Mathematical Formulation and Estimation Methods

2.1. Sensor Model

The sensor is composed of two monolithic fused silica resonators, each of which are made up of a dual-flexure supported proof mass. The dynamics governing the motion of the proof mass undergoing small deflections can be modeled as a perturbed linear harmonic oscillator. The effects of neglected vibrational modes, as well as thermomechanical and optical noise sources are then captured by a Gaussian white noise process and a bias term. The resultant equation of motion is given as:
x ¨ + 2 ω ζ x ˙ + ω 2 x = b ( t ) + n v ( t ) + g ( t )
where x is the proof mass’s displacement from the equilibrium position, g ( t ) is the forcing acceleration, ω is the oscillators natural frequency, and ζ is the damping ratio, which is inversely proportional to the quality factor ( Q = 1 / 2 ζ ). The bias term b ( t ) is modeled as a Wiener process:
b ˙ ( t ) = n u ( t )
where { n v ( t ) , n u ( t ) } are uncorrelated, zero-mean Gaussian random processes with autocorrelation functions:
E [ n v ( t i ) n v ( t j ) ] = σ v 2 δ ( t i t j ) E [ n u ( t i ) n u ( t j ) ] = σ u 2 δ ( t i t j )
where δ ( t ) is the Dirac delta function. The contributing process noise sources, including: thermomechanical losses and cavity drift have been investigated in references [16,17], respectively. Highly precise manufacturing processes coupled with pressure and temperature control allow for these noise sources to be modeled as Gaussian. The effect of radiation pressure from the optical measurement device is also conservatively introduced to the system using the Wiener process model for the bias. It is worth noting that strictly speaking the contribution of radiation pressure is not Gaussian, however it is commonly linearized for reduced order modeling [31,32]. For this single degree of freedom model the effect of off-axis accelerations are not considered, but could be augmented to a three-axis accelerometer configuration using this model in a straightforward manner.
The equations of motion can be expressed in state-space representation as [33]:
X ˙ ( t ) = 0 1 0 ω 2 2 ω ζ 1 0 0 0 X ( t ) + 0 g ( t ) 0 + 0 n v ( t ) n u ( t ) .
In this formulation the state is defined as: X ( t ) = x ( t ) x ˙ ( t ) b ( t ) T . The continuous time dynamics for the linear time invariant system can be discretized in a straightforward manner [33,34]:
X k + 1 = Φ ( t k + 1 , t k ) X k + Γ ( t k + 1 , t k ) g k + w k
Ψ ( t k + 1 , t k ) = t k t k + 1 Φ ( t k + 1 , τ ) d τ Γ ( t k + 1 , t k ) = Ψ 12 ( t k + 1 , t k ) Ψ 22 ( t k + 1 , t k ) 0 T
where X k = X ( t k ) , and Φ ( t k + 1 , t k ) = e A ( t k + 1 t k ) is the state transition matrix, which is derived explicitly in Appendix A following the developments of Skelton [35]. It is assumed that g ( t ) is constant over a sufficiently small interval ( t k + 1 t k ) . The discrete random process w k encompasses the effects of both n v ( t ) and n u ( t ) , and is expressed as:
w k = t k t k + 1 Φ ( t k + 1 , τ ) f ( τ ) d τ f ( t ) = 0 n v ( t ) n u ( t ) T .
The process noise can be sampled using the covariance ( Q = E [ w k w k T ] ) provided in Appendix B. In conjunction with Equation (5), this allows for the state to be determined at discrete time intervals.
Observations of the oscillator are obtained using a highly sensitive laser displacement measuring interferometer which operates at a fixed sampling frequency. These measurements are expressed as a linear function of the state:
y ˜ k = H X k + ν k = x k + ν k
where H = 1 0 0 is the measurement mapping matrix. The measurement noise ( ν k ) made up of shot noise, electric readout noise, and optical frequency noise are taken to be zero-mean Gaussian with variance σ m 2 , and ν k is Gaussian white noise with variance σ m 2 . It is important to note that g ( t ) is made up of all frequencies of interest, while the observable frequencies of the sensor are limited by the oscillator frequency, the sampling frequency, and the relevant model uncertainties. Therefore, only a component of the forcing acceleration is sensed by the accelerometer.

2.2. Calibration Phase: Estimating the Sensor Bias

The accelerometer bias must be monitored to prevent the accumulated build-up of errors in the estimate of acceleration. To do this a calibration step must be conducted on a reasonable basis. One such method is to carefully apply an input to the sensor over a specified time and observe the response of the system. Knowing that the acceleration cannot be perfectly delivered, we define the provided input as:
g ( t ) = g ¯ ( t ) + n g ( t )
where g ¯ ( t ) is the desired acceleration at time t, and n g ( t ) is a zero-mean Gaussian process with variance σ g 2 . The discrete sensor dynamics given in Equation (5) can now be augmented with the imperfect input acceleration:
X k + 1 = Φ ( t k + 1 , t k ) X k + Γ ( t k + 1 , t k ) g ¯ k + w k
where w k is defined in the same way as Equation (7) with the addition of n g ( t ) :
f ( t ) = 0 n g ( t ) + n v ( t ) n u ( t ) T .
Because the input noise n g ( t ) and the process noise terms { n v ( t ) , n u ( t ) } are independent, the Q matrix derived in Appendix B would be of the same form, with the sum of σ g 2 and σ v 2 taking the place of σ v 2 . With the dynamics in this form, the discrete-time Kalman filter can be used to sequentially estimate the state using displacement measurements defined in Equation (8) [34,36,37]. The update step, which is conducted when a new measurement is incorporated is given as:
X ^ k + = X ^ k + K k y ˜ k H X ^ k P k + = I K k H P k . K k = P k H T H P k H T + R 1
where the caret corresponds to the estimate of the indicated quantity and P is the estimate error covariance matrix. After a measurement is incorporated it is propagated over the sampling interval ( t k + 1 t k ) . The propagation equations are given as:
X ^ k + 1 = Φ ( t k + 1 , t k ) X ^ k + + Γ ( t k + 1 , t k ) g ¯ k P k + 1 = Φ ( t k + 1 , t k ) P k + Φ T ( t k + 1 , t k ) + Q .
Using this method, an estimate of the sensor bias and the variance of the bias estimate error can be obtained sequentially in the calibration step. Provided a reasonably steady input of an admissible magnitude is applied, the filter will converge to steady-state behavior quite rapidly. With this in mind, the convergent value of P can be found by solving the discrete Riccati equation [38,39,40]:
P = Φ ( t k + 1 , t k ) I K H P Φ T ( t k + 1 , t k ) + Q .
In doing so, a closed-form solution for the expected filter performance can be obtained based on oscillator parameters, sampling frequency, and process noise parameters. This steady-state analysis is important for sensor design and component selection. Accuracies of the displacement-sensing interferometers, process models of the mechanical elements, and the signal processing specifications can be designed using the mathematics and methods of the steady-state analysis.

2.3. Operation Phase: Estimating the Forcing Acceleration

The results obtained from the calibration step can now be incorporated to develop procedure for the estimation of the forcing acceleration using a series of displacement measurements defined in Equation (8). Upon completion of the calibration step, it is assumed that the bias estimate ( b ^ ) is itself unbiased. Under this assertion, the sensor bias is given as:
b ( t 0 ) = b ^ + n b 0 b ( t k ) = b ( t 0 ) + t 0 t k n u ( τ ) d τ
where t 0 is the time of calibration and n b 0 is a zero-mean Gaussian random variable with variance σ b 0 2 . Between calibrations the bias estimate is not updated, so b ^ does not change. Because n b 0 and n u ( t ) are both zero mean, the error in the bias estimate between calibrations will remain zero mean; however the variance in the estimate will grow [41]. The statistical quantities governing this growth are given by:
E [ b k b ^ ] = 0 E [ ( b k b ^ ) ( b k b ^ ) T ] = σ b 0 2 + ( t k t 0 ) σ u 2 .
This growth in bias uncertainty between calibrations will be shown to degrade the accuracy of the subsequent acceleration estimate as the time since calibration increases. Using this instantaneous definition for the bias term, the discrete dynamics in Equation (5) can be rewritten as:
x k + 1 x ˙ k + 1 = Φ 11 Φ 12 Φ 21 Φ 22 x k x ˙ k + Ψ 12 Ψ 22 g k + Φ 13 Φ 23 b k + t k t k + 1 n v ( τ ) Φ 12 ( t k + 1 , τ ) n v ( τ ) Φ 22 ( t k + 1 , τ ) d τ .
Using this reduced form of the discrete dynamics, two estimators will now be presented.

2.3.1. Least Squares Formulation

The first estimator considered is a moving least squares method using a batch of N + 1 measurements obtained from the laser interferometer, which are defined in Equation (8). A measurement obtained at n { 0 , 1 , 2 , , N | N 2 } steps ahead of t k is then expressed as the following:
y k + n = Φ 11 ( n ) Φ 12 ( n ) Ψ 12 ( n ) x k x ˙ k g k + Φ 13 ( n ) b ^ + ν ˜ k + n ν ˜ k + n = ν k + n + Φ 13 ( n ) ( b k b ^ ) + t k t k + n Φ 13 ( t ( k + n ) , τ ) n v ( τ ) d τ ,
where Φ i j ( n ) and Ψ i j ( n ) are used as shorthand for Φ i j ( t ( k + n ) , t k ) and Ψ i j ( t ( k + n ) , t k ) , respectively. Using this, a set of measurements can then be expressed as a function of the instantaneous oscillator state and forcing acceleration:
y k y k + 1 y k + N = 1 0 0 Φ 11 ( 1 ) Φ 12 ( 1 ) Ψ 12 ( 1 ) Φ 11 ( n ) Φ 12 ( n ) Ψ 12 ( n ) x k x ˙ k g k + 0 Φ 13 ( 1 ) Φ 13 ( n ) b ^ + ν ˜ k ν ˜ k + 1 ν ˜ k + N .
Now, the minimum variance estimate of x k x ˙ k g k T can be determined from the set of measurements:
x ^ k x ˙ ^ k g ^ k = [ H ˜ k T R 1 H ˜ k ] 1 H ˜ k T R 1 y ˜ η k b ^ P = [ H ˜ k T R 1 H ˜ k ] 1
H ˜ k = 1 0 0 Φ 11 ( 1 ) Φ 12 ( 1 ) Ψ 12 ( 1 ) Φ 11 ( n ) Φ 12 ( n ) Ψ 12 ( n ) y ˜ = y k y k + 1 y k + N η k = 0 Φ 13 ( 1 ) Φ 13 ( n )
where the caret corresponds to the estimate of the indicated term and R is the measurement noise covariance. The elements of R are given in indicial notation as:
R i j = σ m 2 δ i j + Φ 13 ( i ) Φ 13 ( j ) σ b 0 2 + σ u 2 ( t k t 0 ) + σ v 2 t k t k + i Φ 12 2 ( t k + i , τ ) d τ
where i j and δ i j is the Kronecker delta. A closed form solution for the integral in this expression is derived in Appendix B. The measurement noise has an explicit dependence on time, and as a result will increase with the time since last calibration. As a result the certainty in the estimate of acceleration will intuitively degrade during operation, necessitating periodic calibration. Provided the user has some operational requirements on the acceleration estimates, the calibration cycle can be determined explicitly in this manner.

2.3.2. Kalman Filter Formulation

An estimator for the determination of the oscillator state and forcing acceleration can also be developed using a Kalman filter approach. Returning to Equation (17), the dynamics can be expressed as:
x k + 1 x ˙ k + 1 = Φ 11 Φ 12 Φ 21 Φ 22 x k x ˙ k + Ψ 12 Ψ 22 ( g k + b ^ ) + w ˜ k w ˜ k = Φ 13 Φ 23 ( b k b ^ ) + t k t k + 1 n v ( τ ) Φ 12 ( t k + 1 , τ ) n v ( τ ) Φ 22 ( t k + 1 , τ ) d τ .
The new process noise covariance Q ˜ , which has been augmented with the uncertainty in the instantaneous value of the bias, is given by the following:
Q ˜ = Φ 13 2 Φ 13 Φ 23 Φ 23 Φ 13 Φ 23 2 σ b 0 2 + σ u 2 ( t k t 0 ) + σ v 2 t k t k + 1 Φ 12 2 ( t k + 1 , τ ) d τ t k t k + 1 Φ 12 ( t k + 1 , τ ) Φ 22 ( t k + 1 , τ ) d τ t k t k + 1 Φ 12 ( t k + 1 , τ ) Φ 22 ( t k + 1 , τ ) d τ t k t k + 1 Φ 22 2 ( t k + 1 , τ ) d τ
where the integrals contained in the second matrix are derived in Appendix B. The discrete dynamics can then be rearranged to include the instantaneous forcing acceleration as an additional state. In the absence of a model for g ( t ) , the estimate must be driven by the incorporation of new measurements. This can be done by tuning the process noise associated with the channel to an appropriately scaled value ( α ).
x k + 1 x ˙ k + 1 g k + 1 = Φ 11 Φ 12 Ψ 12 Φ 21 Φ 22 Ψ 22 0 0 1 x k x ˙ k g k + Ψ 12 Ψ 22 0 b ^ + w ˜ k 0
Using these discrete dynamics, the Kalman filter can be run with slight modifications to the update and propagation equations given in Equations (12) and (13), respectively.
X ^ k + = X ^ k + K k y ˜ k H X ^ k P k + = I K k H P k . K k = P k H T H P k H T + R 1 ,
X ^ k + 1 = Φ ( t k + 1 , t k ) X ^ k + + Γ ( t k + 1 , t k ) b ^ P k + 1 = Φ ( t k + 1 , t k ) P k + Φ T ( t k + 1 , t k ) + Q ˜ 0 0 α ,
Note here that Φ ( t k + 1 , t k ) provided in Appendix A is identical to matrix in Equation (25), which allows for its use in the propagation equations above. Like the calibration step, the filter is run using displacement measurements of the proof mass given in Equation (8), which each possess measurement noise covariance R = σ m 2 . Figure 1 summarizes the developments for the calibration and operation filter and provides the framework for intended use.

3. Numerical Simulation Results

The combined calibration and acceleration estimation procedure is implemented and tested using a simulated dataset generated using the sensor model developed in Section 2.1. The oscillator parameters that define the sensor response are provided in Table 1. These parameters correspond to experimental results obtained from benchtesting of a prototype sensor [16,17]. The parameters that define the stochastic processes driving the system { σ v , σ u } are provided in Table 2. Using this information a truth model has been generated, from which the displacement time history has been extracted and corrupted by noise to create a simulated measurement set. These noisy measurements are then provided to the calibration and acceleration estimation programs.
Firstly, the calibration step was conducted over a span of 10 min using a constant input acceleration of g ¯ = 1 × 10 6 m/s 2 . The input noise density corresponds to 10 % error in delivered acceleration, and is provided in Table 2. The initial estimate of the bias is significantly offset from the true value (1 m/s 2 ). The true and estimated bias during the calibration are plotted with the bias estimate error and corresponding 3 σ bounds in Figure 2. The filter is able to converge in the span of a few seconds, producing an unbiased estimate with root mean square error of 1.9816 ng. The converged variance of the bias estimate error (provided in Table 2) is then fed forward to the acceleration estimation procedure in the next step.
After the calibration phase is implemented, the acceleration estimation step is carried out over an extended operation interval of 8 h. The initial true bias was set to zero, while b ^ was set as a random drawing from N ( 0 , σ b 0 2 ) . The external acceleration was defined as a sinusoid of amplitude 1 × 10 5 m/s 2 and frequency 1/100 Hz. The results for the error in estimated acceleration and the 3 σ bounds are provided in Figure 3 for the moving least squares procedure and in Figure 4 for the Kalman filter. The estimate errors for both filters start as zero mean; however, as time increases the influence of the bias on the estimate error becomes apparent in both cases. This is expected, as the filter is initialized with a constant b ^ for the initial time, and beyond that point the uncertainty in b ( t ) is reflected by the growth in the covariance. In both cases the estimate error is contained within the plotted 3 σ bounds. The acceleration estimates can be deemed reliable until the 3 σ bounds become sufficiently large to necessitate calibration, at which point the new bias estimate and σ b 0 can be fed forward to continue operation.
To quantify the filter performance in a manner that is independent of the specific realization of the bias random process, the procedure has been conducted for both estimation approaches in a series of 10,000, 1 h trials to generate statistical measures of performance. The results from this Monte Carlo analysis are presented in Table 3.

4. Discussion

The sensor model developed in Section 2 has been employed to simulate noisy sensor dynamics and measurements. Using these measurements, both the calibration and operation phase estimation procedures have been conducted in series. Figure 2 shows the results obtained in the calibration phase. It is shown that the filter is able to converge to steady-state operation in less than 5 s, in spite of poor a priori knowledge of the bias state. It is shown that the estimate tracks the true value well over the extended 10 min operating window.
Figure 3 shows the results from the operation stage using the moving least squares formulation over a period of 8 h without calibration. During this time, the covariance grows significantly as the certainty in the bias estimate provided by the calibration is diminished. The resultant estimate error is shown to be dominated by the effect of the bias, which is plotted in yellow. This same procedure was conducted over a 1 h window for 10,000 Monte Carlo iterations, 100 of which are plotted with the 3 σ bounds in Figure 5. The figure shows that the covariance is consistent with the expected results in that the errors are bounded by red curves. An important observation is that the user can determine the required calibration cycle based on required level of accuracy of the acceleration estimates. Figure 4 shows the results from the operation stage using the Kalman filter formulation. The results here mirror those of the least squares formulation. The estimate error is primarily driven by the bias term over the 8 h period, and the estimate error is contained within the 3 σ bounds which are plotted in red. The main distinction between the two is the tuning of the parameter α for the process noise covariance. Figure 6 shows the results from 100 Monte Carlo iterations using the Kalman filter. As is shown, the covariance scaling is consistent with the acceleration estimate error results obtained.

5. Conclusions

Inertial navigation is a crucial component to the future of autonomy across several disciplines. It is required that accurate and precise measurements of acceleration and angular rates be acquired in spite of various sources of error to conduct this procedure effectively. In this paper a stochastic model for novel optomechanical inertial sensors is presented and utilized to derive procedures for the estimation of sensor bias during a calibration step, and forcing accelerations during operation. These procedures are developed to be fully compatible, as outputs from the calibration filter are used directly in the operation phase, thus constituting a full sequence for accurate acceleration measurement. This procedure is demonstrated using representative sensor parameters and simulated process and measurement noise. The simulation results verify that the Kalman filter used for calibration is able to produce a highly accurate estimate of the instantaneous bias in the presence of an imperfect input, and that both operational filters produce precise estimates of acceleration. The acceleration estimation errors induced by the bias as it evolves over the operational window are shown to be fully captured by the covariance; this result is proposed as a quantitative method to determine the required calibration cycle to maintain a desired degree of accuracy. In summary, this work addresses the critical concerns of: how to estimate the forcing acceleration in the presence of a stochastic bias term, how to estimate this bias term during a calibration phase, and how regularly must calibration phases be conducted to retain the desired degree of precision.

Author Contributions

Conceptualization, P.K., M.M. and F.G.; methodology, P.K. and M.M.; software, P.K.; validation, P.K.; formal analysis, P.K. and M.M.; investigation, P.K.; writing—original draft preparation, P.K.; writing—review and editing, P.K., M.M. and F.G.; visualization, P.K.; supervision, M.M.; funding acquisition, M.M. and F.G. All authors have read and agreed to the published version of the manuscript.

Funding

Support from National Geospatial Intelligence Agency (NGA) through Grant No. HM0476-19-1-2015 is gratefully acknowledged. Dr. Scott A. True is thanked for serving as a technical monitor and research advocate for this work. Office of Naval Research Grant Number N00014-19-1-2435 is acknowledged for their support. The authors are thankful to Dr. Brian Holm-Hansen and Dr. David Gonzales for their encouragement and support. The Jet Propulsion Laboratory (JPL) subcontract 1659175 Texas Intelligent Space Systems (TISS) initiative and JPL’s Strategic University Research Partnership (SURP) initiative are gratefully acknowledged.

Acknowledgments

The authors would like to thank John L. Junkins for his support on the project and Adam Hines for correspondence regarding his work on optomechanical inertial sensors. Additional thanks to Anup B. Katake, Brian T. Young, Tim McElrath, Fred Y. Hadaegh and Sarah U. Stevens of Jet Propulsion Laboratory for their technical inputs and kind help.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. State Transition Matrix

The state transition matrix Φ ( t , t 0 ) provides the fundamental solution to a homogeneous linear state-space system, mapping the state at some reference time ( t 0 ) into the solution space at a future time ( t ) . For a state-space defined as:
x ˙ ( t ) = A x ( t ) + B u ( t )
where { A , B } are constant matrices, the state transition matrix is determined by the matrix exponential:
Φ ( t , t 0 ) = e A ( t t 0 ) .
The resultant solution flow of the non-homogeneous system given in Equation (A1) is then determined with the STM:
x ( t ) = Φ ( t , t 0 ) x ( t 0 ) + B t 0 t Φ ( t , τ ) u ( τ ) d τ
as is seen in the formulation of Equation (5). For the continuous LTI system defined in Equation (4) the matrix exponential provides the following analytic STM:
Φ = Φ 11 Φ 12 Φ 13 Φ 21 Φ 22 Φ 23 0 0 1
p 1 = ω 1 ζ 2 p 2 = ω ζ
Φ 11 = 1 p 1 e p 2 ( t t 0 ) p 1 cos p 1 ( t t 0 ) + p 2 sin p 1 ( t t 0 ) Φ 12 = 1 p 1 e p 2 ( t t 0 ) sin p 1 ( t t 0 ) Φ 13 = 1 ω 2 1 1 p 1 e p 2 ( t t 0 ) p 1 cos p 1 ( t t 0 ) + p 2 sin p 1 ( t t 0 )
Φ 21 = ω 2 p 1 e p 2 ( t t 0 ) sin p 1 ( t t 0 ) Φ 22 = 1 p 1 e p 2 ( t t 0 ) p 1 cos p 1 ( t t 0 ) p 2 sin p 1 ( t t 0 ) Φ 23 = 1 p 1 e p 2 ( t t 0 ) sin p 1 ( t t 0 ) .
Additionally, to account for the effect of the forcing function, the integral of Φ ( t , t 0 ) is required, particularly the elements that make up Γ ( t , t 0 ) which account for the effect of the forcing acceleration and bias terms.
Ψ ( t , t 0 ) = t 0 t Φ ( t , τ ) d τ
Ψ 11 = 1 p 1 ω 2 2 p 1 p 2 + e p 2 d t ( p 1 2 p 2 2 ) sin p 1 ( t t 0 ) 2 p 1 p 2 cos p 1 ( t t 0 ) Ψ 12 = 1 ω 2 1 1 p 1 e p 2 ( t t 0 ) p 1 cos p 1 ( t t 0 ) + p 2 sin p 1 ( t t 0 ) Ψ 13 = 1 ω 2 1 p 1 ω 4 2 p 1 p 2 + e p 2 d t ( p 1 2 p 2 2 ) sin p 1 ( t t 0 ) 2 p 1 p 2 cos p 1 ( t t 0 ) Ψ 21 = 1 1 p 1 e p 2 ( t t 0 ) p 1 cos p 1 ( t t 0 ) + p 2 sin p 1 ( t t 0 ) Ψ 22 = 1 p 2 e p 2 ( t t 0 ) sin p 1 ( t t 0 ) Ψ 23 = 1 ω 2 1 1 p 1 e p 2 ( t t 0 ) p 1 cos p 1 ( t t 0 ) + p 2 sin p 1 ( t t 0 ) Ψ 31 = 0 Ψ 32 = 0 Ψ 33 = ( t t 0 )

Appendix B. Process Noise Covariance

The process noise vector w k defined in Equation (5) is composed of a Gaussian white noise process, and a random walk. The instantaneous value of this term is given by:
w k = t k t k + 1 Φ 11 ( t k + 1 , τ ) Φ 12 ( t k + 1 , τ ) Φ 13 ( t k + 1 , τ ) Φ 21 ( t k + 1 , τ ) Φ 22 ( t k + 1 , τ ) Φ 23 ( t k + 1 , τ ) Φ 31 ( t k + 1 , τ ) Φ 32 ( t k + 1 , τ ) Φ 33 ( t k + 1 , τ ) 0 n v ( τ ) n u ( τ ) d τ = t k t k + 1 Φ 12 ( t k + 1 , τ ) n v ( τ ) + Φ 13 ( t k + 1 , τ ) n u ( τ ) Φ 22 ( t k + 1 , τ ) n v ( τ ) + Φ 23 ( t k + 1 , τ ) n u ( τ ) n u ( τ ) d τ
For any admissible value of k, w k is a zero-mean Gaussian random vector with covariance Q = E [ w k w k T ] . As such, the elements of this covariance matrix are required to define the multivariate distribution from which drawings of w k are taken. By applying the definitions in Equation (3), the elements of Q are given as:
Q 11 = σ v 2 t k t k + 1 Φ 12 2 ( t k + 1 , τ ) d τ + σ u 2 t k t k + 1 Φ 13 2 ( t k + 1 , τ ) d τ Q 12 = Q 21 = σ v 2 t k t k + 1 Φ 12 ( t k + 1 , τ ) Φ 22 ( t k + 1 , τ ) d τ + σ u 2 t k t k + 1 Φ 13 ( t k + 1 , τ ) Φ 23 ( t k + 1 , τ ) d τ Q 13 = Q 31 = σ u 2 t k t k + 1 Φ 13 ( t k + 1 , τ ) d τ Q 22 = σ v 2 t k t k + 1 Φ 22 2 ( t k + 1 , τ ) d τ + σ u 2 t k t k + 1 Φ 23 2 ( t k + 1 , τ ) d τ Q 23 = Q 32 = σ u 2 t k t k + 1 Φ 23 ( t k + 1 , τ ) d τ Q 33 = σ u 2 t k t k + 1 1 d τ .
Note that the units of { σ v , σ u } must agree with those provided in Table 2 to ensure that these products are dimensionally consistent. Neglecting the contribution of these integrals will result in misrepresentation of the random process parameters. To proceed, the integrals in the above expressions must be evaluated. These have each been evaluated analytically, giving the following:
t k t k + 1 Φ 12 2 ( t k + 1 , τ ) d τ = 1 4 p 1 2 1 p 2 1 e 2 p 2 d t 1 ω 2 p 2 e 2 p 2 d t p 2 cos ( 2 p 1 d t ) p 1 sin ( 2 p 1 d t ) t k t k + 1 Φ 13 2 ( t k + 1 , τ ) d τ = 1 4 p 2 p 1 2 ω 10 [ e 2 p 2 d t p 2 ( p 2 3 3 p 1 2 p 2 ) cos 2 p 1 d t + p 2 ( p 1 3 3 p 1 p 2 2 ) sin 2 p 1 d t ω 4 + e p 2 d t 16 p 1 2 p 2 2 cos p 1 d t + 8 p 1 p 2 ( p 1 + p 2 ) ( p 2 p 1 ) sin p 1 d t + p 1 2 p 2 2 4 p 2 d t 11 + p 1 2 1 + 4 p 2 d t ]
t k t k + 1 Φ 22 2 ( t k + 1 , τ ) d τ = 1 4 p 2 1 4 p 1 2 p 2 e 2 p 2 d t ω 2 p 2 2 cos ( 2 p 1 d t ) p 1 p 2 sin ( 2 p 1 d t ) t k t k + 1 Φ 23 2 ( t k + 1 , τ ) d τ = 1 4 p 1 2 1 p 2 1 e 2 p 2 d t 1 ω 2 p 2 e 2 p 2 d t p 2 cos ( 2 p 1 d t ) p 1 sin ( 2 p 1 d t )
t k t k + 1 Φ 12 ( t k + 1 , τ ) Φ 22 ( t k + 1 , τ ) d τ = 1 4 p 1 2 e 2 p 2 d t 1 cos ( 2 p 1 d t ) t k t k + 1 Φ 13 ( t k + 1 , τ ) Φ 23 ( t k + 1 , τ ) d τ = e 2 p 2 d t 2 p 1 2 ω 6 p 1 cos p 1 d t + p 2 sin p 1 d t p 1 e p 2 d t 2
t k t k + 1 Φ 13 ( t k + 1 , τ ) d τ = 1 ω 4 ω 2 d t 2 P 2 + e P 2 d t 2 P 2 cos ( P 1 d t ) + P 2 2 P 1 P 1 sin ( P 1 d t ) t k t k + 1 Φ 23 ( t k + 1 , τ ) d τ = 1 P 1 ω 2 P 1 e P 2 d t P 2 sin ( P 1 d t ) + P 1 cos ( P 1 d t )
where d t = t k + 1 t k . Alternatively, elegant solutions for the computation of integrals involving the matrix exponential have been developed by Van Loan [42]. Using the system matrix in Equation (4), the covariance can be computed simply as:
Q = 0 d t e A τ Q c e A T τ d τ = F 3 T ( d t ) G 2 ( d t )
Q c = 0 0 0 0 σ v 2 0 0 0 σ u 2
exp A Q c 0 A T d t = F 2 ( d t ) G 2 ( d t ) 0 F 3 ( d t )
Using either of these closed form solutions for the covariance matrix of the oscillator state and bias, the stochastic dynamics can be realized using drawings from a zero-mean Gaussian random vector of unit variance:
X k + 1 = Φ ( t k + 1 , t k ) X k + Γ ( t k + 1 , t k ) g ¯ k + Q N 3 × 1 ( 0 | 1 )

References

  1. Zanetti, R.; Holt, G.; Gay, R.; D’Souza, C.; Sud, J.; Mamich, H.; Begley, M.; King, E.; Clark, F.D. Absolute Navigation Performance of the Orion Exploration Flight Test 1. J. Guid. Control Dyn. 2017, 40, 1106–1116. [Google Scholar] [CrossRef]
  2. Verras, A.; Eapen, R.T.; Simon, A.B.; Majji, M.; Bhaskara, R.R.; Restrepo, C.I.; Lovelace, R. Vision and Inertial Sensor Fusion for Terrain Relative Navigation. In Proceedings of the AIAA Scitech 2021 Forum, 19–21 January 2021. [Google Scholar]
  3. Demcak, S.; Young, B.; Graat, E.; Beswick, R.; Criddle, K.; Ionasescu, R.; Hughes, R.S.; Lee, J.; Haggard, M.; Sealy, N.; et al. Navigation Design and Operations of MAVEN Aerobraking. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020; p. 0471. [Google Scholar]
  4. Jah, M.K.; Lisano, M.E.; Born, G.H.; Axelrad, P. Mars Aerobraking Spacecraft State Estimation by Processing Inertial Measurement Unit Data. J. Guid. Control Dyn. 2008, 31, 1802–1812. [Google Scholar] [CrossRef] [Green Version]
  5. Zurek, R.; Tolson, R.; Bougher, S.; Lugo, R.; Baird, D.; Bell, J.; Jakosky, B. Mars Thermosphere as Seen in MAVEN Accelerometer Data. J. Geophys. Res. Space Phys. 2017, 122, 3798–3814. [Google Scholar] [CrossRef]
  6. Christophe, B.; Foulon, B.; Liorzou, F.; Lebat, V.; Boulanger, D.; Huynh, P.A.; Zahzam, N.; Bidel, Y.; Bresson, A. Status of Development of the Future Accelerometers for Next Generation Gravity Missions. In International Symposium on Advancing Geodesy in a Changing World; Springer: Cham, Switzerland, 2019; pp. 85–89. [Google Scholar]
  7. Christophe, B.; Boulanger, D.; Foulon, B.; Huynh, P.A.; Lebat, V.; Liorzou, F.; Perrot, E. A new generation of ultra-sensitive electrostatic accelerometers for GRACE follow-on and towards the next generation gravity missions. Acta Astronaut. 2015, 117, 1–7. [Google Scholar] [CrossRef]
  8. Peron, R.; Lorenzini, E.C. METRIC: A Dedicated Earth-Orbiting Spacecraft for Investigating Gravitational Physics and the Space Environment. Aerospace 2017, 4, 38. [Google Scholar] [CrossRef] [Green Version]
  9. Antonucci, F.; Armano, M.; Audley, H.; Auger, G.; Benedetti, M.; Binetruy, P.; Boatella, C.; Bogenstahl, J.; Bortoluzzi, D.; Bosetti, P.; et al. LISA Pathfinder: Mission and status. Class. Quantum Gravity 2011, 28, 094001. [Google Scholar] [CrossRef]
  10. Middlemiss, R.; Samarelli, A.; Paul, D.; Hough, J.; Rowan, S.; Hammond, G. Measurement of the Earth tides with a MEMS gravimeter. Nature 2016, 531, 614–617. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Bai, Y.; Li, Z.; Hu, M.; Liu, L.; Qu, S.; Tan, D.; Tu, H.; Wu, S.; Yin, H.; Li, H.; et al. Research and Development of Electrostatic Accelerometers for Space Science Missions at HUST. Sensors 2017, 17, 1943. [Google Scholar] [CrossRef] [PubMed]
  12. El-Sheimy, N.; Youssef, A. Inertial sensors technologies for navigation applications: State of the art and future trends. Satell. Navig. 2020, 1, 1–21. [Google Scholar] [CrossRef] [Green Version]
  13. Guzmán, F.; Kumanchik, L.M.; Taylor, J.M.; Pratt, J.R. Optomechanial Gravimeter. U.S. Patent 10,545,259B2, 28 December 2017. [Google Scholar]
  14. Taylor, J.R. Classical Mechanics; University Science Books: Sausalito, CA, USA, 2005. [Google Scholar]
  15. Guzmán Cervantes, F.; Kumanchik, L.; Pratt, J.; Taylor, J.M. High sensitivity optomechanical reference accelerometer over 10 kHz. Appl. Phys. Lett. 2014, 104, 221111. [Google Scholar] [CrossRef] [Green Version]
  16. Hines, A.; Richardson, L.; Wisniewski, H.; Guzmán, F. Optomechanical inertial sensors. Appl. Opt. 2020, 59, G167–G174. [Google Scholar] [CrossRef]
  17. Wisniewski, H.; Richardson, L.; Hines, A.; Laurain, A.; Guzmán, F. Optomechanical lasers for inertial sensing. J. Opt. Soc. Am. 2020, 37, B87–B92. [Google Scholar] [CrossRef]
  18. DeMars, K.J.; Ward, K.C. Impact of Considering and Neglecting States on Descent-to-Landing Navigation. In Proceedings of the AIAA Scitech 2020 Forum, Orlando, FL, USA, 6–10 January 2020. [Google Scholar]
  19. Kratzer, K.M.; Helmuth, J.C.; Ward, K.C.; DeMars, K.J. Impact of Sensor Model Fidelity and Scheduling on Navigation Performance. In Proceedings of the 2018 AIAA Guidance Navigation, and Control Conference, Kissimmee, FL, USA, 8–12 January 2018. [Google Scholar]
  20. Christian, J.A.; Lightsey, E.G. High-Fidelity Measurement Models for Optical Spacecraft Navigation. In Proceedings of the 22nd International Technical Meeting of the Satellite Division of The Institute of Navigation (ION GNSS 2009), Savannah, GA, USA, 22–25 September 2009; pp. 1486–1503. [Google Scholar]
  21. Farrenkopf, R. Analytic steady-state accuracy solutions for two common spacecraft attitude estimators. J. Guid. Control 1978, 1, 282–284. [Google Scholar] [CrossRef]
  22. Markley, F.L.; Reynolds, R. Analytic steady-state accuracy of a spacecraft attitude estimator. J. Guid. Control Dyn. 2000, 23, 1065–1067. [Google Scholar] [CrossRef] [Green Version]
  23. Markley, F.L. Analytic Steady-State Accuracy of a Three-Axis Spacecraft Attitude Estimator. J. Guid. Control Dyn. 2017, 40, 2393–2400. [Google Scholar] [CrossRef]
  24. Dianetti, A.D.; Crassidis, J.L. Extension of Farrenkopf Steady-State Solutions with Estimated Angular Rate. In Proceedings of the 2018 AIAA Guidance Navigation, and Control Conference, Kissimmee, FL, USA, 8–12 January 2018; p. 2095. [Google Scholar]
  25. Crassidis, J. Sigma-Point Kalman Filtering for Integrated GPS and Inertial Navigation. IEEE Trans. Aerosp. Electron. Syst. 2006, 42, 750–756. [Google Scholar] [CrossRef]
  26. Van Helleputte, T.; Doornbos, E.; Visser, P. CHAMP and GRACE Accelerometer Calibration by GPS-Based Orbit Determination. Adv. Space Res. 2009, 43, 1890–1896. [Google Scholar] [CrossRef]
  27. Chen, K.; Shen, F.; Zhou, J.; Wu, X. Simulation Platform for SINS/GPS Integrated Navigation System of Hypersonic Vehicles Based on Flight Mechanics. Sensors 2020, 20, 5418. [Google Scholar] [CrossRef]
  28. Yu, Z.; Crassidis, J.L. Accelerometer Bias Calibration Using Attitude and Angular Velocity Information. J. Guid. Control Dyn. 2016, 39, 741–753. [Google Scholar] [CrossRef] [Green Version]
  29. Olsson, F.; Kok, M.; Halvorsen, K.; Schön, T.B. Accelerometer Calibration Using Sensor Fusion with a Gyroscope. In Proceedings of the 2016 IEEE Statistical Signal Processing Workshop (SSP), Palma de Mallorca, Spain, 26–29 June 2016. [Google Scholar]
  30. Edwan, E.; Knedlik, S.; Loffeld, O. Angular Motion Estimation Using Dynamic Models in a Gyro-Free Inertial Measurement Unit. Sensors 2012, 12, 5310–5327. [Google Scholar] [CrossRef]
  31. Brawley, G.; Vanner, M.; Larsen, P.E.; Schmid, S.; Boisen, A.; Bowen, W. Nonlinear optomechanical measurement of mechanical motion. Nat. Commun. 2016, 7, 10988. [Google Scholar] [CrossRef] [PubMed]
  32. Xiong, H.; Si, L.G.; Wu, Y. Precision measurement of electrical charges in an optomechanical system beyond linearized dynamics. Appl. Phys. Lett. 2017, 110, 171102. [Google Scholar] [CrossRef]
  33. Kailath, T. Linear Systems; Prentice-Hall: Englewood Cliffs, NJ, USA, 1980. [Google Scholar]
  34. Crassidis, J.L.; Junkins, J.L. Optimal Estimation of Dynamic Systems; CRC Press: Boca Raton, FL, USA, 2011. [Google Scholar]
  35. Skelton, R.E. Dynamic Systems Control: Linear Systems Analysis and Synthesis; John Wiley & Sons: New York, NY, USA, 1988. [Google Scholar]
  36. Kalman, R.E. A New Approach to Linear Filtering and Prediction Problems. J. Basic Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  37. Kalman, R.E.; Bucy, R.S. New Results in Linear Filtering and Prediction Theory. J. Basic Eng. 1961, 83, 95–108. [Google Scholar] [CrossRef]
  38. Vaughan, D. A Nonrecursive Algebraic Solution for the Discrete Riccati Equation. IEEE Trans. Autom. Control 1970, 15, 597–599. [Google Scholar] [CrossRef]
  39. Laub, A. A Schur Method for Solving Algebraic Riccati Equations. IEEE Trans. Autom. Control 1979, 24, 913–921. [Google Scholar] [CrossRef] [Green Version]
  40. Arnold, W.; Laub, A. Generalized Eigenproblem Algorithms and Software for Algebraic Riccati Equations. Proc. IEEE 1984, 72, 1746–1754. [Google Scholar] [CrossRef]
  41. Carpenter, J.R.; D’Souza, C.N. Navigation Filter Best Practices; Technical Report; NASA/TP–2018–219822; NASA Langley Research Center: Hampton, VA, USA, 2018. [Google Scholar]
  42. Van Loan, C. Computing Integrals Involving the Matrix Exponential. IEEE Trans. Autom. Control 1978, 23, 395–404. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Flowchart of the calibration and operation Kalman filters. As shown, the outputs of each filter are provided as inputs to the other. Equation (20) can easily be substituted into the calibration step to use the lerast squares solution in place of the Kalman filter.
Figure 1. Flowchart of the calibration and operation Kalman filters. As shown, the outputs of each filter are provided as inputs to the other. Equation (20) can easily be substituted into the calibration step to use the lerast squares solution in place of the Kalman filter.
Sensors 21 06101 g001
Figure 2. Results for the bias estimate obtained from the simulated calibration phase. The top figure shows the true and estimated bias over the calibration time interval of 10 min. On the bottom, the error in the bias estimate is plotted in blue with the corresponding 3 σ bounds in red.
Figure 2. Results for the bias estimate obtained from the simulated calibration phase. The top figure shows the true and estimated bias over the calibration time interval of 10 min. On the bottom, the error in the bias estimate is plotted in blue with the corresponding 3 σ bounds in red.
Sensors 21 06101 g002
Figure 3. Results for the forcing acceleration estimate error from the moving least squares procedure. Estimate error is plotted in blue with the associated 3 σ bounds in red. This is accompanied by the true bias time history in yellow.
Figure 3. Results for the forcing acceleration estimate error from the moving least squares procedure. Estimate error is plotted in blue with the associated 3 σ bounds in red. This is accompanied by the true bias time history in yellow.
Sensors 21 06101 g003
Figure 4. Results for the forcing acceleration estimate error from the Kalman filter. Estimate error is plotted in blue with the associated 3 σ bounds in red. This is accompanied by the true bias time history in cyan.
Figure 4. Results for the forcing acceleration estimate error from the Kalman filter. Estimate error is plotted in blue with the associated 3 σ bounds in red. This is accompanied by the true bias time history in cyan.
Sensors 21 06101 g004
Figure 5. Results for the forcing acceleration estimate error from the moving least squares procedure for 100 Monte Carlo iterations. Estimate error is plotted in blue with the associated 3 σ bounds in red.
Figure 5. Results for the forcing acceleration estimate error from the moving least squares procedure for 100 Monte Carlo iterations. Estimate error is plotted in blue with the associated 3 σ bounds in red.
Sensors 21 06101 g005
Figure 6. Results for the forcing acceleration estimate error from the Kalman filter procedure for 100 Monte Carlo iterations. Estimate error is plotted in blue with the associated 3 σ bounds in red.
Figure 6. Results for the forcing acceleration estimate error from the Kalman filter procedure for 100 Monte Carlo iterations. Estimate error is plotted in blue with the associated 3 σ bounds in red.
Sensors 21 06101 g006
Table 1. Oscillator parameters for optomechanical inertial sensor.
Table 1. Oscillator parameters for optomechanical inertial sensor.
TermValueUnits
ω 3.76Hz
Q 1.14 × 10 5
f s 30.5Hz
Table 2. Random process parameters.
Table 2. Random process parameters.
TermValueUnits
σ v 1 × 10 9 m/s Hz
σ u 1 × 10 8 m/s 2 Hz
σ m 1 × 10 11 m
σ g 1.8107 × 10 8 m/s Hz
σ b 0 1.3506 × 10 8 m/s 2
Table 3. Key performance statistics from 10,000 Monte Carlo trials of both acceleration estimation procedures.
Table 3. Key performance statistics from 10,000 Monte Carlo trials of both acceleration estimation procedures.
TermValueUnits
Least Squares Error Mean 7.7990 × 10 5 μg
Least Squares Error Std Dev 3.6446 × 10 2 μg
Kalman Filter Error Mean 3.7493 × 10 5 μg
Kalman Filter Error Std Dev 3.5476 × 10 2 μg
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kelly, P.; Majji, M.; Guzmán, F. Estimation and Error Analysis for Optomechanical Inertial Sensors. Sensors 2021, 21, 6101. https://doi.org/10.3390/s21186101

AMA Style

Kelly P, Majji M, Guzmán F. Estimation and Error Analysis for Optomechanical Inertial Sensors. Sensors. 2021; 21(18):6101. https://doi.org/10.3390/s21186101

Chicago/Turabian Style

Kelly, Patrick, Manoranjan Majji, and Felipe Guzmán. 2021. "Estimation and Error Analysis for Optomechanical Inertial Sensors" Sensors 21, no. 18: 6101. https://doi.org/10.3390/s21186101

APA Style

Kelly, P., Majji, M., & Guzmán, F. (2021). Estimation and Error Analysis for Optomechanical Inertial Sensors. Sensors, 21(18), 6101. https://doi.org/10.3390/s21186101

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop