CN114400986A - Cascade filtering method and device for portable relative gravimeter - Google Patents

Cascade filtering method and device for portable relative gravimeter Download PDF

Info

Publication number
CN114400986A
CN114400986A CN202111486896.5A CN202111486896A CN114400986A CN 114400986 A CN114400986 A CN 114400986A CN 202111486896 A CN202111486896 A CN 202111486896A CN 114400986 A CN114400986 A CN 114400986A
Authority
CN
China
Prior art keywords
gravimeter
sampling
gravity
filtering
output data
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
CN202111486896.5A
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.)
717th Research Institute of CSIC
Original Assignee
717th Research Institute of CSIC
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 717th Research Institute of CSIC filed Critical 717th Research Institute of CSIC
Priority to CN202111486896.5A priority Critical patent/CN114400986A/en
Publication of CN114400986A publication Critical patent/CN114400986A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0283Filters characterised by the filter structure
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V7/00Measuring gravitational fields or waves; Gravimetric prospecting or detecting
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/06Non-recursive filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H2017/0072Theoretical filter design
    • H03H2017/0081Theoretical filter design of FIR filters

Abstract

The invention relates to a cascade filtering method and a device of a portable relative gravimeter, and the scheme comprises two-stage filtering, wherein the output data of the gravimeter is subjected to down-sampling, and then the down-sampled output data of the gravimeter is sent to an FIR filter for filtering to obtain the final output result of the gravimeter. Through the technical scheme, the output precision of the gravimeter is high, the down-sampling treatment of the cascade filter can inhibit high-frequency noise in advance, and the measurement precision is improved. Meanwhile, the scheme saves computing resources, and the order of the filter can be reduced through the down-sampling processing of the cascade filter, so that the operation pressure is reduced.

Description

Cascade filtering method and device for portable relative gravimeter
Technical Field
The invention relates to the technical field of surveying and mapping, in particular to a cascade filtering method and device of a portable relative gravimeter.
Background
The portable relative gravimeter is used as a measuring instrument capable of sensing environmental gravity information, can directly reflect the spatial distribution characteristics of mineral resources so as to provide convenience for resource exploration and utilization, and can also provide original reference gravity information for military striking. However, in the actual gravity measurement process, the external motion disturbance is often coupled with the gravity field and is simultaneously measured by the gravity sensor, thereby causing measurement errors. Therefore, an effective method capable of separating the motion disturbance from the actual gravity field needs to be adopted to improve the gravity measurement precision.
Disclosure of Invention
Aiming at the technical problems in the prior art, the invention provides a cascade filtering method and a cascade filtering device for a portable relative gravimeter, which can solve the problem that the measurement precision of the portable relative gravimeter is influenced by motion disturbance.
The technical scheme for solving the technical problems is as follows:
in a first aspect, the present invention provides a method for cascading filtering of a portable relative gravimeter, including:
down-sampling the output data of the gravimeter;
and sending the output data of the gravity meter after the down sampling into an FIR filter for filtering to obtain a final output result of the gravity meter.
Further, the gravimeter is digitally leveled prior to the gravimeter measurement.
Further, the digital leveling of the gravimeter comprises:
according to the Euler angle definition, the output value of the gravity meter and the pitching and rolling states of the equipment have the following relations:
Figure BDA0003397827890000021
wherein [ a ] isx ay az]The output value of the gravimeter is g, the reference gravity is g, the pitch to be solved of the equipment is alpha, and the roll to be solved of the equipment is beta;
and the relation between the output value of the gravity meter and the pitching and rolling states of the equipment is utilized to complete pitching and rolling calculation and realize digital leveling.
Further, down-sampling the gravimeter output data includes: down-sampling the gravimeter output data according to the following formula:
Figure BDA0003397827890000022
in the formula (f)smoothIs a down-sampling result; f. ofiThe gravity signal is an original gravity signal in the day direction at the moment i; n is the down-sampling time span.
Further, the down-sampling time span N is taken to be 1 s.
Further, the cut-off frequency f of the FIR filtersCalculated by the following formula:
d=0.5×v×(1/fs)
in the formula: d is the spatial resolution of the gravimeter and v is the average horizontal movement velocity of the gravimeter carrier.
In a second aspect, the present invention provides a cascade filter device for a portable relative gravimeter, comprising:
the down-sampling module is used for down-sampling the data output by the gravimeter;
and the FIR filter module is used for filtering the gravity meter output data after the down sampling to obtain the final output result of the gravity meter.
Further, the device also comprises a digital leveling module which is used for carrying out digital leveling on the gravimeter before the gravimeter measures.
In a third aspect, the present invention provides an electronic device comprising:
a memory for storing a computer software program;
and the processor is used for reading and executing the computer software program so as to realize the cascade filtering method of the portable relative gravimeter according to the first aspect of the invention.
In a fourth aspect, the present invention provides a non-transitory computer readable storage medium having stored therein a computer software program for implementing the method of cascaded filtering for a portable relative gravimeter according to the first aspect of the present invention.
The invention has the beneficial effects that: the invention comprises two stages of low-pass filters, wherein the first stage of low-pass filter is used for realizing data down-sampling, and the second stage of low-pass filter is used for finishing gravity measurement filtering. The technology can be used for gravity measurement filtering of a portable relative gravimeter and has certain engineering value.
Considering that the gravity information is in a low frequency band, namely the mean value is down sampled and the effective gravity information is not attenuated, the data processing flow adopts a flow sequential processing mode. After original gravity information is measured, mean value down-sampling is carried out on high-frequency sampling data in advance, and disturbance influence is initially inhibited; and then, the low-frequency gravity data after down sampling is processed by an FIR low-pass filter to complete gravity calculation.
Through the technical scheme, the output precision of the gravimeter is high, the down-sampling treatment of the cascade filter can inhibit high-frequency noise in advance, and the measurement precision is improved. Meanwhile, the scheme saves computing resources, and the order of the filter can be reduced through the down-sampling processing of the cascade filter, so that the operation pressure is reduced.
Drawings
Fig. 1 is a schematic structural diagram of a cascade filter device of a portable relative gravimeter according to an embodiment of the present invention;
fig. 2 is a schematic flow chart of a cascade filtering method of a portable relative gravimeter according to an embodiment of the present invention.
Fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of a non-transitory computer-readable storage medium according to an embodiment of the present invention.
Detailed Description
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth by way of illustration only and are not intended to limit the scope of the invention.
As shown in fig. 1, an embodiment of the present invention provides a cascade filter device for a portable relative gravimeter. The device comprises a digital leveling module, a down-sampling module and an FIR filter module.
The digital leveling module is used for performing digital leveling on the gravity meter before measurement of the gravity meter;
the down-sampling module is used for down-sampling the output data of the gravimeter;
and the FIR filter module is used for filtering the gravity meter output data after the down sampling to obtain the final output result of the gravity meter.
Based on the device, the embodiment of the invention also provides a cascade filtering method of the portable relative gravimeter, and the flow chart of the method is shown in fig. 2. The method comprises the following steps:
carrying out digital leveling on the gravimeter;
down-sampling the output data of the gravimeter;
and sending the output data of the gravity meter after the down sampling into an FIR filter for filtering to obtain a final output result of the gravity meter.
The portable relative gravimeter adopts three accelerometers as gravity sensors, and digital leveling is required to be completed according to the pitching and rolling states of the equipment in the actual measurement process, so that the output value of the three-axis accelerometer can be correctly projected to the sky direction. According to the Euler angle definition, when the equipment has pitching and rolling, the measured value of the three-axis accelerometer and the pitching and rolling states of the equipment have the following relations:
Figure BDA0003397827890000041
wherein [ a ] isx ay az]The output value of the gravimeter is g, the reference gravity is g, the pitch to be solved of the equipment is alpha, and the roll to be solved of the equipment is beta;
and the relation between the output value of the gravity meter and the pitching and rolling states of the equipment is utilized to complete pitching and rolling calculation and realize digital leveling.
In the gravity measurement process, a raw sensor signal generally has a high-frequency sampling frequency of 50-500 HZ, and the signal contains high-frequency noise with large amplitude and external disturbance noise. The direct low-pass filtering calculation by using the high-frequency signal not only consumes the calculation resources of the processor seriously, but also causes the phenomenon that the low-frequency gravity information cannot be reflected correctly, and is not favorable for attenuation suppression of high-frequency disturbance. Therefore, it is necessary to adopt a down-sampling technique to process the raw data of the sensor:
Figure BDA0003397827890000051
in the formula (f)smoothIs a down-sampling result; f. ofiThe gravity signal is an original gravity signal in the day direction at the moment i; n is the down-sampling time span.
To fit satellite navigation measurements, which are typically at 1Hz sampling frequency, N is typically chosen to be 1 s. After the down-sampling technology is adopted, the suppression of high-frequency noise can be directly realized, and the subsequent gravity measurement filtering processing can be facilitated.
In order to further increase the accuracy of the gravity measurement result, it is necessary to design a low-pass filter for gravity measurement based on the down-sampling result based on the frequency band characteristics and measurement requirements of the gravity information. Currently, scholars mainly use finite impulse response filters (FIR filters) to implement gravity measurement, and the FIR filters can be equivalent to a differential equation form convenient for encoding implementation:
Figure BDA0003397827890000052
in the formula: a isk,bkIs the filter coefficient; m and N jointly form the order of the filter; x () is the input sequence; y () is the output sequence. In the process of actually applying the FIR filter, the spatial resolution requirement (or referred to as half-wavelength resolution) of the gravity measurement result needs to be combined, and the spatial resolution of the gravimeter is as follows:
d=0.5×v×(1/fs)
in the formula: d is the spatial resolution of the gravimeter, v is the average horizontal movement velocity of the gravimeter carrier, fsThe frequency is cut for a filter corresponding to an attenuation gain of-3 dB.
Spatial resolution defined by the cut-off frequency fsDetermine, and fsThe cut-off frequency f of the filter can be calculated according to the average motion speed v when the carrier measures gravity with the required spatial resolution d as inputs
For motion disturbance suppression, a mainstream solution is to separate motion acceleration from gravity based on signal frequency domain characteristics starting from a low-pass digital filter. Based on the technical scheme, the portable cascade filtering technology for the relative gravimeter comprises two stages of low-pass filters, wherein the first stage of low-pass filter is used for achieving data down-sampling, and the second stage of low-pass filter is used for achieving gravity measurement filtering. The technology can be used for gravity measurement filtering of a portable relative gravimeter and has certain engineering value.
Considering that the gravity information is in a low frequency band, namely the mean value is down sampled and the effective gravity information is not attenuated, the data processing flow adopts a flow sequential processing mode. After original gravity information is measured, mean value down-sampling is carried out on high-frequency sampling data in advance, and disturbance influence is initially inhibited; and then, the low-frequency gravity data after down sampling is processed by an FIR low-pass filter to complete gravity calculation.
Through the technical scheme, the output precision of the gravimeter is high, the down-sampling treatment of the cascade filter can inhibit high-frequency noise in advance, and the measurement precision is improved. Meanwhile, the scheme saves computing resources, and the order of the filter can be reduced through the down-sampling processing of the cascade filter, so that the operation pressure is reduced.
Referring to fig. 3, fig. 3 is a schematic diagram of an embodiment of an electronic device according to an embodiment of the invention. As shown in fig. 3, an embodiment of the present invention provides an electronic device 500, which includes a memory 510, a processor 520, and a computer program 511 stored in the memory 520 and executable on the processor 520, wherein the processor 520 executes the computer program 511 to implement the following steps:
down-sampling the output data of the gravimeter;
and sending the output data of the gravity meter after the down sampling into an FIR filter for filtering to obtain a final output result of the gravity meter.
Referring to fig. 4, fig. 4 is a schematic diagram illustrating an embodiment of a computer-readable storage medium according to the present invention. As shown in fig. 4, the present embodiment provides a computer-readable storage medium 600 having a computer program 611 stored thereon, the computer program 611, when executed by a processor, implementing the steps of:
down-sampling the output data of the gravimeter;
and sending the output data of the gravity meter after the down sampling into an FIR filter for filtering to obtain a final output result of the gravity meter.
It should be noted that, in the foregoing embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to relevant descriptions of other embodiments for parts that are not described in detail in a certain embodiment.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
While preferred embodiments of the present invention have been described, additional variations and modifications in those embodiments may occur to those skilled in the art once they learn of the basic inventive concepts. Therefore, it is intended that the appended claims be interpreted as including preferred embodiments and all such alterations and modifications as fall within the scope of the invention.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.

Claims (10)

1. A cascade filtering method of a portable relative gravimeter is characterized by comprising the following steps:
down-sampling the output data of the gravimeter;
and sending the output data of the gravity meter after the down sampling into an FIR filter for filtering to obtain a final output result of the gravity meter.
2. The method of claim 1, wherein the gravimeter is digitally leveled prior to the gravimeter measurement.
3. The method of claim 2, wherein the digitally leveling the gravimeter comprises:
according to the Euler angle definition, the output value of the gravity meter and the pitching and rolling states of the equipment have the following relations:
Figure FDA0003397827880000011
wherein [ a ] isx ay az]The output value of the gravimeter is g, the reference gravity is g, the pitch to be solved of the equipment is alpha, and the roll to be solved of the equipment is beta;
and the relation between the output value of the gravity meter and the pitching and rolling states of the equipment is utilized to complete pitching and rolling calculation and realize digital leveling.
4. The method of claim 1, wherein down-sampling the gravimeter output data comprises: down-sampling the gravimeter output data according to the following formula:
Figure FDA0003397827880000012
in the formula (f)smoothIs a down-sampling result; f. ofiThe gravity signal is an original gravity signal in the day direction at the moment i; n is the down-sampling time span.
5. The method of claim 4, wherein the downsampled time span N is taken to be 1 s.
6. Method according to claim 1, characterized in that the cut-off frequency f of the FIR filtersCalculated by the following formula:
d=0.5×v×(1/fs)
in the formula: d is the spatial resolution of the gravimeter and v is the average horizontal movement velocity of the gravimeter carrier.
7. A cascaded filter arrangement for a portable relative gravimeter, comprising:
the down-sampling module is used for down-sampling the data output by the gravimeter;
and the FIR filter module is used for filtering the gravity meter output data after the down sampling to obtain the final output result of the gravity meter.
8. The apparatus of claim 7, further comprising a digital leveling module for digitally leveling the gravimeter prior to the gravimeter measurement.
9. An electronic device, comprising:
a memory for storing a computer software program;
a processor for reading and executing the computer software program to implement the cascade filtering method of a portable relative gravimeter according to any of claims 1 to 6.
10. A non-transitory computer readable storage medium having stored therein a computer software program for implementing the method of cascaded filtering for a portable relative gravimeter according to any of claims 1-6.
CN202111486896.5A 2021-12-07 2021-12-07 Cascade filtering method and device for portable relative gravimeter Pending CN114400986A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111486896.5A CN114400986A (en) 2021-12-07 2021-12-07 Cascade filtering method and device for portable relative gravimeter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111486896.5A CN114400986A (en) 2021-12-07 2021-12-07 Cascade filtering method and device for portable relative gravimeter

Publications (1)

Publication Number Publication Date
CN114400986A true CN114400986A (en) 2022-04-26

Family

ID=81227498

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111486896.5A Pending CN114400986A (en) 2021-12-07 2021-12-07 Cascade filtering method and device for portable relative gravimeter

Country Status (1)

Country Link
CN (1) CN114400986A (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101083461A (en) * 2007-07-13 2007-12-05 北京创毅视讯科技有限公司 Drop sampling filtering process and drop sampling filter
CN102033163A (en) * 2009-10-06 2011-04-27 精工爱普生株式会社 Frequency measurement method, frequency measurement device and apparatus equipped with frequency measurement device
CN107359868A (en) * 2016-05-10 2017-11-17 深圳市中兴微电子技术有限公司 Pulse density modulated change-over circuit and method
US10784881B1 (en) * 2019-11-15 2020-09-22 Liquid Instruments Pty Ltd. Blended anti-aliasing analog-to-digital conversion for digital test and measurement devices
CN111900953A (en) * 2020-07-16 2020-11-06 上海富芮坤微电子有限公司 System and method for realizing sampling conversion and filtering of digital microphone
CN212726967U (en) * 2020-09-04 2021-03-16 成都泰格微波技术股份有限公司 Digital decimation filter for analog-to-digital converter

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101083461A (en) * 2007-07-13 2007-12-05 北京创毅视讯科技有限公司 Drop sampling filtering process and drop sampling filter
CN102033163A (en) * 2009-10-06 2011-04-27 精工爱普生株式会社 Frequency measurement method, frequency measurement device and apparatus equipped with frequency measurement device
CN107359868A (en) * 2016-05-10 2017-11-17 深圳市中兴微电子技术有限公司 Pulse density modulated change-over circuit and method
US10784881B1 (en) * 2019-11-15 2020-09-22 Liquid Instruments Pty Ltd. Blended anti-aliasing analog-to-digital conversion for digital test and measurement devices
CN111900953A (en) * 2020-07-16 2020-11-06 上海富芮坤微电子有限公司 System and method for realizing sampling conversion and filtering of digital microphone
CN212726967U (en) * 2020-09-04 2021-03-16 成都泰格微波技术股份有限公司 Digital decimation filter for analog-to-digital converter

Similar Documents

Publication Publication Date Title
Sun et al. Mooring alignment for marine SINS using the digital filter
El-Sheimy et al. Wavelet de-noising for IMU alignment
CN101840529B (en) Optic fiber gyroscope random drift modeling method based on locally variable integrated neural network
CN108759845A (en) A kind of optimization method based on inexpensive multi-sensor combined navigation
CN104596514A (en) Real-time noise reduction system and real-time noise reduction method of accelerometer and gyroscope
CN105300381A (en) Rapid convergence method based on improved complementary filter for attitude of self-balance mobile robot
CN103414451B (en) A kind of EKF method being applied to attitude of flight vehicle and estimating
CN104280047A (en) Gyroscope shift filtering system and method integrating multiple sensors
CN106802143A (en) A kind of hull deformation angle measuring method based on inertial instruments and Iterative-Filtering Scheme
CN110703205A (en) Ultrashort baseline positioning method based on adaptive unscented Kalman filtering
CN108007457B (en) Monitoring navigation system asynchronous data fusion method based on subdivided time slices
CN103940447B (en) Mooring state initial aligning method based on self-adaptive digital filter
CN109186571A (en) A kind of gyroscope filtering and noise reduction method
CN114400986A (en) Cascade filtering method and device for portable relative gravimeter
CN111121820B (en) MEMS inertial sensor array fusion method based on Kalman filtering
CN110533005B (en) Method for measuring deformation of ship body under complex sea condition
CN112611382A (en) Strapdown inertial navigation system heave measurement method with phase compensation
Yongpan et al. Real-time zero phase filtering for heave measurement
CN116608859A (en) Navigation method, storage medium and device of self-adaptive unscented Kalman filtering based on threshold processing
CN112051595A (en) Backward differential filtering method for solving motion acceleration of carrier by utilizing DGPS (differential global positioning system) position information
CN112629540B (en) Heave measurement method based on carrier attitude information
CN107664499B (en) On-line noise reduction method for accelerometer of ship strapdown inertial navigation system
Sun et al. Novel approach to GPS/SINS integration for IMU alignment
CN206479651U (en) A kind of air-ground TEM reception systems with posture writing function
RU2643201C2 (en) Strap down inertial attitude-and-heading reference

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
CB02 Change of applicant information
CB02 Change of applicant information

Country or region after: China

Address after: 430000 981 Xiongchu street, Hongshan District, Wuhan City, Hubei Province

Applicant after: Huazhong Optoelectronic Technology Research Institute (717 Research Institute of China Shipbuilding Corp.)

Address before: 430000 981 Xiongchu street, Hongshan District, Wuhan City, Hubei Province

Applicant before: HUAZHONG OPTOELECTRONIC TECHNOLOGY Research Institute (THE 717TH RESEARCH INSTITUTE OF CSIC)

Country or region before: China