Home Analog Devices Feedback Subscribe Archives 简体中文 日本語 | |

Download this article in PDF format. (737 KB) Analyzing Frequency Response of Inertial MEMS in Stabilization Systems
Since many stabilization systems require more than one axis of active correction, An accelerometer responds to both Another type of sensor is the A third type of sensor is the 3-axis The generic block diagram of Figure 2 shows how gyroscope and accelerometer measurements can be employed in a manner that uses their basic strengths but minimizes the impact of their weaknesses. The pole locations of the low-pass accelerometer and high-pass gyroscope filters are typically application-dependent, with accuracy goals, phase delay, vibration, and “normal” motion expectations, all contributing to these decisions. System-dependent behaviors will also affect the weighting factors, which also have an impact on how these two measurements are combined. The extended Kalman filter is one example of an algorithm that combines the filtering and weighting functions to calculate dynamic angle estimates.
A strategy for evaluating IMU bandwidth is determining what information is available in product documentation, analyzing the impact of this information on the system’s response to inertial motion, and stabilizing the system’s response. This analysis, and any corrective actions it entails, will become the basis for preliminary testing. Frequency response is often represented as “bandwidth” in specification tables for IMUs and gyroscopes. As a performance parameter, it represents the frequency at which the output magnitude drops to about 70% (–3 dB) of the actual magnitude of motion that the sensor is experiencing. In some cases, bandwidth may also be defined by the frequency at which the output response lags the actual motion by 90 degrees (for a 2-pole system). Both of these metrics can directly impact an important stability criterion for a control loop: Analyzing frequency response in a system starts with a high-level, “black box” view, which describes the system’s response to inputs over the entire frequency range of interest. In electronic circuits, where the input and output are defined in common terms, such as signal level (volts), this typically involves developing a The analysis process starts with understanding all of the components associated with a sensor signal chain. Figure 3 offers an overall diagram of the typical functions. The signal chain starts with a core sensor element, which translates the inertial motion into a representative electrical signal. If the bandwidth is not limited in the sensor element, it is often limited by filters in the signal-conditioning circuit preceding the ADC. After the signals are digitized, a processor typically applies correction (calibration) formulas and digital filtering. The secondary digital filters reduce the bandwidth and sample rates that the feedback systems use in their control routines. All of these stages can influence the gain and phase of the sensor signal, with respect to frequency. Figure 3 provides an example of an IMU that has multiple filters in a mixed-signal processing system. This system will serve as an example for illustrating some useful analysis techniques. Figure 3. ADIS16488 sensor in a signal chain for frequency analysis.
The accelerometer’s single-pole ( These formulas provide the basis for numerical analysis in programs that can manage the complex numbers associated with the Fmax = 9840/2; % one-half of the sample rate For a quick assessment of the time delay associated with these filters, notice that the phase delay of a single-pole filter is equal to 45° at its −3-dB frequency, or 1/8 of the corner frequency’s period. In this case, the time delay of the accelerometer’s filter is approximately equal to 0.38 ms. For the gyroscope, the delay is equal to the sum of the time delays of the two stages, or about 0.47 ms.
When using MATLAB for analysis, use a sample rate ( N), along with the same frequency array (f) used to analyze the analog filter. Using a common frequency array will make it easier to combine the results of each stage. Use the following code to analyze this first stage:Fmax = 9840/2; % one-half of the sample rate Analyzing the second averaging/decimating filter will require prior knowledge of the control system’s sample rate but will use the same relationships. For example, if a control loop requires a sample rate that is close to 400 SPS, the second filter’s average and decimation rate would be equal to six (for a sample rate of 410 SPS and four samples, 9840/[410 × 4] = 6). Use the same
Fortunately, many modern programs contain specific tools or commands for this type of analysis, based on these basic relationships. It is still useful to understand them when verifying results of the automatic assessment tools and in developing an intuitive feel for when to question the outputs of a FIR design tool. The MATLAB “fdatool” command launches its filter analysis and design package, which helps design and analyze the system FIR filter implementation.
For early analysis validation without a rate table, measuring the spectral noise over the frequency band of interest can provide useful insights. This simplified approach does not require sophisticated test equipment but only a secure mechanical connection to a stable platform and data collection instrumentation. However, it does rely on the mechanical noise having a “flat” noise magnitude with respect to frequency. Figure 6 illustrates two examples that both use the same 2-pole, low-pass filter. The first example (ADIS16375) uses a gyroscope that has a flat response over its usable frequency range. The second example (ADIS16488) uses a gyroscope that has a modest amount of peaking at 1.2 kHz, which actually extends the –3-dB frequency to approximately 380 Hz. Recognizing this resonant behavior can be valuable for those in the process of modeling and simulating a control loop. Identifying this behavior in a simple test can also help explain noise levels that are higher than expected when performing a more thorough system characterization. When understood and identified early in a project, these behaviors can normally be managed with adjustments to the filter poles. When measuring noise density, make sure that the sample rate is at least twice the highest frequency of interest to meet the Nyquist criterion. Also, take enough data samples to reduce the uncertainty of the measurements. The plots in Figure 6 were derived from FFT analysis of a time record with a length of 256k samples at a maximum rate of 2.46 kSPS.
Another approach uses a gyroscope’s self-test function. The self-test function provides an opportunity to stimulate the sensor’s mechanical structure, using an electrical signal, without requiring the device to be subjected to external inertial motion. The self-test function forces a change in the sensor core that simulates its response to actual motion, producing a corresponding change in the electrical output. Not all products provide real-time access to this, but it can be a useful tool when available, or if the manufacturer can provide data from this type of frequency-response test. In the simplest approach, the self-test, which simulates response to a step, is compared with the analytical expectation. Repeating the self-test assertion at specific frequencies provides a direct method for studying the magnitude of the sensor response at each frequency. Consider the two different responses in Figure 7. At the lower frequency, the gyroscope output looks like a square wave, with the exception of the transient response at each transition. The transient response follows the expectation of a “step response” for the filter network in the sensor signal chain. In the second example, where the frequency of the self-test is high enough to prevent full settling, a decrease in magnitude occurs. Notice the difference in magnitudes between the blue and black-dotted responses, on the bottom signal in this figure. There are a number of methods for estimating the magnitude of these time records. A discrete Fourier transform (DFT) separates the primary frequency content (self-test frequency) from the harmonic content, which can contribute errors to the magnitude/frequency response.
Figure 7. Self test.
Developing an understanding of an IMU’s bandwidth and its role in system stability should employ analysis, modeling, test data, and iteration of these factors. Start by quantifying the information available, make assumptions to close any gaps, and then develop a plan to refine these assumptions. I invite you to comment on inertial MEMS in stabilization systems in the Analog Dialogue Community on EngineerZone.
WEBCAST: Using MEMS Sensors for Industrial Platform Stabilization Systems The Five Motion Senses: Using MEMS Inertial Sensing to Transform Applications Improving Industrial Control with Integrated MEMS Inertial Sensors
Copyright 1995- Analog Devices, Inc. All rights reserved. |