CN117928528A - Ship heave measurement method based on self-adaptive time-delay-free complementary band-pass filter - Google Patents
Ship heave measurement method based on self-adaptive time-delay-free complementary band-pass filter Download PDFInfo
- Publication number
- CN117928528A CN117928528A CN202410331990.0A CN202410331990A CN117928528A CN 117928528 A CN117928528 A CN 117928528A CN 202410331990 A CN202410331990 A CN 202410331990A CN 117928528 A CN117928528 A CN 117928528A
- Authority
- CN
- China
- Prior art keywords
- band
- frequency
- heave
- obtaining
- delay
- 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.)
- Granted
Links
- 230000000295 complement effect Effects 0.000 title claims abstract description 38
- 238000000691 measurement method Methods 0.000 title claims abstract description 14
- 238000000034 method Methods 0.000 claims abstract description 87
- 230000001133 acceleration Effects 0.000 claims abstract description 41
- 230000033001 locomotion Effects 0.000 claims abstract description 26
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 238000001514 detection method Methods 0.000 claims abstract description 11
- 238000007619 statistical method Methods 0.000 claims abstract description 7
- 230000003595 spectral effect Effects 0.000 claims description 39
- 230000003044 adaptive effect Effects 0.000 claims description 8
- 238000005311 autocorrelation function Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 230000009466 transformation Effects 0.000 claims description 5
- 230000003111 delayed effect Effects 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 10
- 230000010355 oscillation Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/10—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
- G01C21/12—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
- G01C21/16—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
- G01C21/165—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/20—Instruments for performing navigational calculations
- G01C21/203—Specially adapted for sailing ships
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C23/00—Combined instruments indicating more than one navigational value, e.g. for aircraft; Combined measuring devices for measuring two or more variables of movement, e.g. distance, speed or acceleration
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Navigation (AREA)
Abstract
The invention discloses a ship heave measurement method based on a self-adaptive time-delay-free complementary band-pass filter, which belongs to the technical field of inertial navigation measurement and is used for ship heave measurement, and comprises the steps of obtaining a power spectrum density flattening value and a frequency corresponding to the power spectrum density flattening value according to a statistical method and a nearest neighbor method so as to obtain a stop band bandwidth; calculating according to the center cut-off frequency and the stop band bandwidth to obtain a lower limit cut-off frequency and an upper limit cut-off frequency, and taking the low-pass filter as a medium to obtain a transfer function of the band-stop filter; according to the complementary principle, a transfer function of a non-delay complementary band-pass filter is obtained, and further, the heave acceleration is filtered and integrated twice to obtain the heave amount of the ship. The method provided by the invention can not only effectively reduce the dependence on priori knowledge or experience values and peak detection range setting, but also realize the filter design by adaptively acquiring the filter parameters, and remarkably reduce the influence of delay errors, thereby realizing high-precision measurement of ship heave motion.
Description
Technical Field
The invention discloses a ship heave measurement method based on a self-adaptive time-delay-free complementary band-pass filter, and belongs to the technical field of inertial navigation measurement.
Background
The existing inertial navigation-based heave measurement method is mainly divided into two types, wherein the first type is an inertial navigation-based combined measurement method, such as combined measurement of heave by using inertial navigation and satellite, although high-precision 3-dimensional positioning of the ship is easy to realize, and m-level absolute positioning precision can be achieved; but short-term relative measurement accuracy on the order of cm is often required in heave applications, which is difficult to achieve with conventional integrated navigation methods. Therefore, the ship heave motion measurement also mainly depends on inertial navigation, and the second type of method is inertial measurement and a digital filter, wherein the method utilizes the comprehensive acceleration output of the inertial navigation in the heave direction, and the heave motion is measured through integration and filtering by means of the digital filter. However, the parameter acquisition of the conventional filter is often obtained by means of priori knowledge or empirical values, and the same filter parameters are adopted in the whole period, so that the conventional filter is not suitable for the situation of different sea conditions, and the conventional filter has the influence of time delay errors.
Disclosure of Invention
The invention aims to provide a ship heave measurement method based on a self-adaptive time-delay-free complementary band-pass filter, so as to solve the problem of low ship heave measurement precision in the prior art.
A ship heave measurement method based on a self-adaptive time-delay-free complementary band-pass filter comprises the following steps:
S1, mechanically arranging according to a strapdown inertial navigation system to obtain roll angle, pitch angle and course angle information, and constructing a direction cosine matrix converted from a carrier coordinate system b system to a navigation coordinate system n system Vertical acceleration information output by strapdown inertial navigation system/>And local gravitational acceleration/>Obtaining heave acceleration/>, of a vessel;
S2, carrying out power spectral density analysis on heave acceleration in different windows by adopting a sliding window, and obtaining the stop band center cut-off frequency by using a maximum value detection methodObtaining a value of which the power spectral density tends to be flat and the corresponding frequency according to two methods of statistical analysis and nearest neighbor so as to obtain the stop band bandwidth/>;
S3 according toAnd/>Calculating to obtain the lower limit cut-off angle frequency/>, of the band-stop filterAnd upper cutoff frequency/>Taking the low-pass filter as a medium to obtain a transfer function of the band-stop filter;
and S4, according to the transfer function of the band-stop filter, obtaining a transfer function of a time delay-free complementary band-pass filter based on complementarity, and filtering and twice integrating the ship heave acceleration to obtain a ship heave motion value.
S1 includes the group consisting of,The method comprises the following steps:
;
;/>;
;
In the method, in the process of the invention, Heading angle, roll angle and pitch angle respectively,/>For vertical acceleration in navigational coordinate system,/>Representing the local latitude.
S2 comprises that the length of the heave acceleration data isSliding window length is/>Number of windows/>The method comprises the following steps:
;
If it is For integers, heave acceleration data is divided into/>Windows, carrying out power spectrum density analysis on the window data;
If it is For floating point number or decimal, performing power spectral density analysis on window data of integer part, and performing power spectral density analysis on/>Obtaining window number/>, by rounding downThe method comprises the following steps:
;
In the method, in the process of the invention, The method comprises the steps of representing downward rounding, and taking the rounded data of the insufficient integer as a window to perform power spectrum density analysis.
The power spectral density analysis is that of random signalsIs the autocorrelation function/>The method comprises the following steps:
;
In the method, in the process of the invention, For time delay,/>For observing time,/>For the current time,/>For subtracting the time delayed random signal;
fourier transforming the two sides of the autocorrelation function to obtain the power spectral density of the signal :
;
In the method, in the process of the invention,Representing the sampling frequency,/>Representing a complex number of units.
Obtaining stop band center cut-off frequency using maximum detection methodThe method comprises the steps of carrying out power spectrum density analysis according to ship heave acceleration in different windows, and obtaining frequency/>, corresponding to the maximum value, through maximum value detectionWill/>Center cut-off frequency/>, as a band stop filter:
;
Obtaining a value of which the power spectral density tends to be flat and the corresponding frequency according to two methods of statistical analysis and nearest neighbor so as to obtain the stop band bandwidthThe method comprises obtaining average value of power spectral density value range corresponding to maximum signal count by histogram statistics, and obtaining frequency/>, closest to center cut-off frequency, corresponding to power spectral density of interference signal by nearest neighbor method according to relationship between power spectral density and frequencyStop band bandwidth/>, of a band stop filterThe method comprises the following steps:
。
s3, obtaining lower and upper limit cut-off angular frequencies of the analog band-stop filter according to the center cut-off frequency and the stop band bandwidth, wherein the lower and upper limit cut-off angular frequencies are as follows:
;
In the method, in the process of the invention, The cut-off angle frequency is corresponding to the center cut-off frequency;
Obtaining the stop band cut-off angular frequency of the simulated prototype low-pass filter according to the lower limit cut-off angular frequency and the upper limit cut-off angular frequency of the simulated band-stop filter And passband angular frequency/>;
According to、/>And filter order and type, obtaining a digital prototype low-pass filter transfer function/>, based on bilinear transformationAnd obtain transfer function/>, of the digital band-stop filter according to the spectral transformation:
;
In the method, in the process of the invention,、/>、/>、/>、/>、/>Is a function parameter.
S4 comprises that the input signal isTransfer function/>, of a non-time-delay complementary bandpass filterThe method comprises the following steps:
;
Heave motion of ship The method comprises the following steps:
。
Compared with the prior art, the invention has the following beneficial effects: the method provided by the invention can not only effectively reduce the dependence on priori knowledge or experience values and peak detection range setting, but also realize the filter design by adaptively acquiring the filter parameters, and remarkably reduce the influence of delay errors, thereby realizing high-precision measurement of ship heave motion.
Drawings
FIG. 1 is a flow chart of transfer function acquisition for a non-delay complementary bandpass filter;
FIG. 2 is a graph of power spectral density analysis of vertical acceleration of a vessel;
FIG. 3 is an output heave motion diagram of the vessel;
FIG. 4 is a diagram of a conventional bandpass filter obtaining heave motion of a vessel;
FIG. 5 is a heave motion diagram of a vessel obtained with a fixed parameter time delay-free complementary bandpass filter;
FIG. 6 is a heave motion of a vessel obtained by a slip adaptive delay-free complementary bandpass filter;
FIG. 7 is a graph of error between heave of a vessel and a reference true value obtained by a conventional bandpass filter, and the RMS value;
FIG. 8 is a plot of error and RMS values between heave of the vessel and a reference true value obtained by a fixed parameter non-delayed complementary bandpass filter;
FIG. 9 is a plot of error versus RMS value between heave of the vessel and a reference true value obtained by sliding an adaptive delay-free complementary bandpass filter.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the present invention more apparent, the technical solutions in the present invention will be clearly and completely described below, and it is apparent that the described embodiments are some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
A ship heave measurement method based on a self-adaptive time-delay-free complementary band-pass filter comprises the following steps:
S1, mechanically arranging according to a strapdown inertial navigation system to obtain roll angle, pitch angle and course angle information, and constructing a direction cosine matrix converted from a carrier coordinate system b system to a navigation coordinate system n system Vertical acceleration information output by strapdown inertial navigation system/>And local gravitational acceleration/>Obtaining heave acceleration/>, of a vessel;
S2, carrying out power spectral density analysis on heave acceleration in different windows by adopting a sliding window, and obtaining the stop band center cut-off frequency by using a maximum value detection methodObtaining a value of which the power spectral density tends to be flat and the corresponding frequency according to two methods of statistical analysis and nearest neighbor so as to obtain the stop band bandwidth/>;
S3 according toAnd/>Calculating to obtain the lower limit cut-off angle frequency/>, of the band-stop filterAnd upper cutoff frequency/>Taking the low-pass filter as a medium to obtain a transfer function of the band-stop filter;
and S4, according to the transfer function of the band-stop filter, obtaining a transfer function of a time delay-free complementary band-pass filter based on complementarity, and filtering and twice integrating the ship heave acceleration to obtain a ship heave motion value.
S1 includes the group consisting of,The method comprises the following steps:
;
;/>;
;
In the method, in the process of the invention, Heading angle, roll angle and pitch angle respectively,/>For vertical acceleration in navigational coordinate system,/>Representing the local latitude.
S2 comprises that the length of the heave acceleration data isSliding window length is/>Number of windows/>The method comprises the following steps:
;
If it is For integers, heave acceleration data is divided into/>Windows, carrying out power spectrum density analysis on the window data;
If it is For floating point number or decimal, performing power spectral density analysis on window data of integer part, and performing power spectral density analysis on/>Obtaining window number/>, by rounding downThe method comprises the following steps:
;
In the method, in the process of the invention, The method comprises the steps of representing downward rounding, and taking the rounded data of the insufficient integer as a window to perform power spectrum density analysis.
The power spectral density analysis is that of random signalsIs the autocorrelation function/>The method comprises the following steps:
;
In the method, in the process of the invention, For time delay,/>For observing time,/>For the current time,/>For subtracting the time delayed random signal;
fourier transforming the two sides of the autocorrelation function to obtain the power spectral density of the signal :
;
In the method, in the process of the invention,Representing the sampling frequency,/>Representing a complex number of units.
Obtaining stop band center cut-off frequency using maximum detection methodThe method comprises the steps of carrying out power spectrum density analysis according to ship heave acceleration in different windows, and obtaining frequency/>, corresponding to the maximum value, through maximum value detectionWill/>Center cut-off frequency/>, as a band stop filter:
;
Obtaining a value of which the power spectral density tends to be flat and the corresponding frequency according to two methods of statistical analysis and nearest neighbor so as to obtain the stop band bandwidthThe method comprises obtaining average value of power spectral density value range corresponding to maximum signal count by histogram statistics, and obtaining frequency/>, closest to center cut-off frequency, corresponding to power spectral density of interference signal by nearest neighbor method according to relationship between power spectral density and frequencyStop band bandwidth/>, of a band stop filterThe method comprises the following steps:
。
s3, obtaining lower and upper limit cut-off angular frequencies of the analog band-stop filter according to the center cut-off frequency and the stop band bandwidth, wherein the lower and upper limit cut-off angular frequencies are as follows:
;
In the method, in the process of the invention, The cut-off angle frequency is corresponding to the center cut-off frequency;
Obtaining the stop band cut-off angular frequency of the simulated prototype low-pass filter according to the lower limit cut-off angular frequency and the upper limit cut-off angular frequency of the simulated band-stop filter And passband angular frequency/>;
According to、/>And filter order and type, obtaining a digital prototype low-pass filter transfer function/>, based on bilinear transformationAnd obtain transfer function/>, of the digital band-stop filter according to the spectral transformation:
;
In the method, in the process of the invention,、/>、/>、/>、/>、/>Is a function parameter.
S4 comprises that the input signal isTransfer function/>, of a non-time-delay complementary bandpass filterThe method comprises the following steps:
;
Heave motion of ship The method comprises the following steps:
。
The definition of heave motion of a vessel is: periodic movement of the vessel relative to the local sea level. Heave motions are defined herein as motions similar to the above, i.e. motions of the vessel up and down relative to the local sea level, which motions are along the vertical axis and have periodicity, the frequency band of the motions being called heave frequency band. The acceleration along heave direction of inertial navigation output inevitably contains oscillation type system errors such as a schler oscillation period term and an earth oscillation period term, and compared with heave motion with a period of about ten seconds, the system error period is long and the frequency is low. If the integral operation is directly carried out on the vertical acceleration, the system error is seriously accumulated, and the twice integral natural displacement mainly shows that part of the Schlemer oscillation is in a divergent state in a short time.
The transfer function acquisition flow of the non-delay complementary band-pass filter is shown in figure 1, according to the heave acceleration power spectral density, the frequency corresponding to the maximum power spectral density obtained by a histogram statistics and nearest neighbor method is used as the center cut-off frequency and the stop band bandwidth, and according to the relational expression of the lower and upper limit cut-off angular frequencies, the center cut-off frequency and the stop band bandwidth of the analog band-stop filter, the lower and upper limit cut-off angular frequencies of the analog band-stop filter are obtained through the pre-distortion treatment; frequency conversion is carried out on the filter to obtain stop band cut-off angular frequency and pass band angular frequency of the analog low-pass filter, so that the transfer function of the analog low-pass filter is obtained; according to bilinear transformation, a digital low-pass filter transfer function is obtained, and spectrum transformation is carried out on the digital low-pass filter transfer function to obtain a band-stop filter transfer function; and finally, according to the complementarity, obtaining the transfer function of the non-delay complementary band-pass filter.
In order to verify the effectiveness of the algorithm, the inertial navigation vertical acceleration data are processed, and four methods of POS MV, a traditional band-pass filter, a fixed parameter non-delay complementary band-pass filter and a sliding self-adaptive non-delay complementary band-pass filter are compared. And taking the result of the POS MV as a reference true value, and obtaining the RMS value of the error between the three methods and the reference true value.
Carrying out power spectral density analysis on the vertical acceleration of the ship as shown in fig. 2, carrying out difference between the vertical acceleration output by strapdown inertial navigation and the local gravity acceleration to obtain the heave acceleration of the ship, and solving the power spectral density of the ship, wherein the horizontal axis represents frequency and the vertical axis represents power spectral density; the ship heave output by the POS MV is adopted, the horizontal axis represents time, and the vertical axis represents ship heave motion, as shown in figure 3; the ship heave acceleration is filtered and integrated twice by adopting a traditional band-pass filter to obtain the ship heave as shown in fig. 4, wherein the horizontal axis represents time and the vertical axis represents ship heave motion; filtering and twice integrating the ship heave acceleration by adopting a fixed parameter time delay-free complementary band-pass filter to obtain ship heave as shown in fig. 5, wherein the horizontal axis represents time and the vertical axis represents ship heave motion; filtering and twice integrating the ship heave acceleration by adopting a sliding self-adaptive non-delay complementary band-pass filter to obtain ship heave, wherein the horizontal axis represents time and the vertical axis represents ship heave motion as shown in fig. 6; the difference between the ship heave obtained by POS MV and the ship heave obtained by the traditional band-pass filter is shown in figure 7, a ship heave error map of the traditional band-pass filter is obtained, and the RMS value of the ship heave error map is obtained to be 0.044741m; the difference between the ship heave obtained by the POS MV and the ship heave obtained by the fixed parameter non-delay complementary band-pass filter is shown as figure 8, a ship heave error map of the fixed parameter non-delay complementary band-pass filter is obtained, and the RMS value of the ship heave error map is 0.029791m; the difference between the ship heave obtained by the POS MV and the ship heave obtained by the sliding self-adaptive non-delay complementary band-pass filter is shown in figure 9, a ship heave error map of the sliding self-adaptive non-delay complementary band-pass filter is obtained, and the RMS value of the ship heave error map is obtained to be 0.026692m. Compared with the traditional band-pass filter and the fixed-parameter time-delay-free complementary band-pass filter, the method provided by the invention has the advantage that the heave precision of the ship is improved.
The above embodiments are only for illustrating the technical aspects of the present invention, not for limiting the same, and although the present invention has been described in detail with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may be modified or some or all of the technical features may be replaced with other technical solutions, which do not depart from the scope of the technical solutions of the embodiments of the present invention.
Claims (7)
1. The ship heave measurement method based on the self-adaptive time-delay-free complementary band-pass filter is characterized by comprising the following steps of:
S1, mechanically arranging according to a strapdown inertial navigation system to obtain roll angle, pitch angle and course angle information, and constructing a direction cosine matrix converted from a carrier coordinate system b system to a navigation coordinate system n system Vertical acceleration information output by strapdown inertial navigation system/>And local gravitational acceleration/>Obtaining heave acceleration/>, of a vessel;
S2, carrying out power spectral density analysis on heave acceleration in different windows by adopting a sliding window, and obtaining the stop band center cut-off frequency by using a maximum value detection methodObtaining a value of which the power spectral density tends to be flat and the corresponding frequency according to two methods of statistical analysis and nearest neighbor so as to obtain the stop band bandwidth/>;
S3 according toAnd/>Calculating to obtain the lower limit cut-off angle frequency/>, of the band-stop filterAnd upper cutoff frequency/>Taking the low-pass filter as a medium to obtain a transfer function of the band-stop filter;
and S4, according to the transfer function of the band-stop filter, obtaining a transfer function of a time delay-free complementary band-pass filter based on complementarity, and filtering and twice integrating the ship heave acceleration to obtain a ship heave motion value.
2. The method for measuring heave of a vessel based on an adaptive time-lapse complementary bandpass filter according to claim 1, wherein S1 comprises,The method comprises the following steps:
;
;/>;
;
In the method, in the process of the invention, Heading angle, roll angle and pitch angle respectively,/>For vertical acceleration in navigational coordinate system,/>Representing the local latitude.
3. The method for measuring heave of a vessel based on the adaptive time-lapse-free complementary bandpass filter according to claim 2, wherein S2 comprises that the length of the heave acceleration data isSliding window length is/>Number of windows/>The method comprises the following steps:
;
If it is For integers, heave acceleration data is divided into/>Windows, carrying out power spectrum density analysis on the window data;
If it is For floating point number or decimal, performing power spectral density analysis on window data of integer part, and performing power spectral density analysis on/>Obtaining window number/>, by rounding downThe method comprises the following steps:
;
In the method, in the process of the invention, The method comprises the steps of representing downward rounding, and taking the rounded data of the insufficient integer as a window to perform power spectrum density analysis.
4. A ship heave measurement method based on an adaptive time-delay-free complementary bandpass filter according to claim 3, characterized in that the power spectral density analysis is a random signalIs the autocorrelation function/>The method comprises the following steps:
;
In the method, in the process of the invention, For time delay,/>For observing time,/>For the current time,/>For subtracting the time delayed random signal;
fourier transforming the two sides of the autocorrelation function to obtain the power spectral density of the signal :
;
In the method, in the process of the invention,Representing the sampling frequency,/>Representing a complex number of units.
5. The ship heave measurement method based on the adaptive time-delay-free complementary band-pass filter according to claim 4, wherein the stop band center cut-off frequency is obtained by using a maximum value detection methodThe method comprises the steps of carrying out power spectrum density analysis according to ship heave acceleration in different windows, and obtaining frequency/>, corresponding to the maximum value, through maximum value detectionWill/>Center cut-off frequency/>, as a band stop filter:
;
Obtaining a value of which the power spectral density tends to be flat and the corresponding frequency according to two methods of statistical analysis and nearest neighbor so as to obtain the stop band bandwidthThe method comprises obtaining average value of power spectral density value range corresponding to maximum signal count by histogram statistics, and obtaining frequency/>, closest to center cut-off frequency, corresponding to power spectral density of interference signal by nearest neighbor method according to relationship between power spectral density and frequencyStop band bandwidth/>, of a band stop filterThe method comprises the following steps:
。
6. The ship heave measurement method based on the adaptive time-delay-free complementary band-pass filter according to claim 5, wherein S3 comprises obtaining lower and upper limit cut-off angular frequencies of the analog band-stop filter according to the center cut-off frequency and the stop band bandwidth as follows:
;
In the method, in the process of the invention, The cut-off angle frequency is corresponding to the center cut-off frequency;
Obtaining the stop band cut-off angular frequency of the simulated prototype low-pass filter according to the lower limit cut-off angular frequency and the upper limit cut-off angular frequency of the simulated band-stop filter And passband angular frequency/>;
According to、/>And filter order and type, obtaining a digital prototype low-pass filter transfer function/>, based on bilinear transformationAnd obtain transfer function/>, of the digital band-stop filter according to the spectral transformation:
;
In the method, in the process of the invention,、/>、/>、/>、/>、/>Is a function parameter.
7. The adaptive time-delay-free complementary bandpass filter-based ship heave measurement method according to claim 6, wherein S4 comprises that the input signal isTransfer function/>, of a non-time-delay complementary bandpass filterThe method comprises the following steps:
;
Heave motion of ship The method comprises the following steps:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410331990.0A CN117928528B (en) | 2024-03-22 | 2024-03-22 | Ship heave measurement method based on self-adaptive time-delay-free complementary band-pass filter |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410331990.0A CN117928528B (en) | 2024-03-22 | 2024-03-22 | Ship heave measurement method based on self-adaptive time-delay-free complementary band-pass filter |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117928528A true CN117928528A (en) | 2024-04-26 |
CN117928528B CN117928528B (en) | 2024-05-31 |
Family
ID=90767006
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410331990.0A Active CN117928528B (en) | 2024-03-22 | 2024-03-22 | Ship heave measurement method based on self-adaptive time-delay-free complementary band-pass filter |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117928528B (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030225491A1 (en) * | 2002-05-30 | 2003-12-04 | Sowada Delroy J. | Methods and systems for determining heave and heave rate of vessels |
CN107101631A (en) * | 2017-05-10 | 2017-08-29 | 哈尔滨工程大学 | A kind of ship heave measuring method based on auto-adaptive filtering technique |
CN110319838A (en) * | 2019-07-09 | 2019-10-11 | 哈尔滨工程大学 | A kind of adaptive athletic posture frame of reference heave measurement method |
US20220326019A1 (en) * | 2021-04-21 | 2022-10-13 | Harbin Engineering University | Strapdown Inertial Navigation Heave Measurement Method Using Multiple Low-Pass Filter Units |
CN115876194A (en) * | 2022-12-25 | 2023-03-31 | 哈尔滨工程大学 | Self-adaptive ship heave measurement method based on strapdown inertial navigation |
CN115950423A (en) * | 2023-01-10 | 2023-04-11 | 江苏科技大学 | Ship heave motion measurement method based on adaptive filtering |
-
2024
- 2024-03-22 CN CN202410331990.0A patent/CN117928528B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030225491A1 (en) * | 2002-05-30 | 2003-12-04 | Sowada Delroy J. | Methods and systems for determining heave and heave rate of vessels |
CN107101631A (en) * | 2017-05-10 | 2017-08-29 | 哈尔滨工程大学 | A kind of ship heave measuring method based on auto-adaptive filtering technique |
CN110319838A (en) * | 2019-07-09 | 2019-10-11 | 哈尔滨工程大学 | A kind of adaptive athletic posture frame of reference heave measurement method |
US20220326019A1 (en) * | 2021-04-21 | 2022-10-13 | Harbin Engineering University | Strapdown Inertial Navigation Heave Measurement Method Using Multiple Low-Pass Filter Units |
CN115876194A (en) * | 2022-12-25 | 2023-03-31 | 哈尔滨工程大学 | Self-adaptive ship heave measurement method based on strapdown inertial navigation |
CN115950423A (en) * | 2023-01-10 | 2023-04-11 | 江苏科技大学 | Ship heave motion measurement method based on adaptive filtering |
Non-Patent Citations (3)
Title |
---|
Z GANR, CHI UMG, MANEAE: "Ship manoeuvr ability prediction using neural networks analysis", ADVANCED MATERIALS RESEARCH, 31 December 2014 (2014-12-31) * |
严恭敏;苏幸君;翁浚;秦永元;: "基于惯导和无时延滤波器的舰船升沉测量", 导航定位学报, no. 02, 30 June 2016 (2016-06-30) * |
刘锡祥;黄荣;王启明;汪宋兵;: "一种基于惯性测量与自适应滤波的舰船升沉计算方法", 中国惯性技术学报, no. 01, 28 February 2019 (2019-02-28) * |
Also Published As
Publication number | Publication date |
---|---|
CN117928528B (en) | 2024-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105222809B (en) | The method of the boat magnetic disturbance penalty coefficient estimation of a kind of ground magnetic gradient robust | |
CN106487359B (en) | The building method of Morphologic filters based on self-adapting multi-dimension AVG-Hat transformation | |
CN109883392B (en) | Strapdown inertial navigation heave measurement method based on phase compensation | |
CN114998419B (en) | Gully volume change slicing method based on terrain point cloud | |
CN113552565B (en) | Phase unwrapping method for SAR data high noise and large gradient change area | |
CN112347992B (en) | Remote sensing estimation method for time sequence AGB in desert area | |
CN108108692A (en) | A kind of sinusoidal signal frequency recognition methods | |
CN112269192B (en) | Quick self-adaptive dynamic Beidou monitoring real-time resolving denoising method | |
CN114662348B (en) | Method for establishing time-varying sea tide load displacement model based on GNSS technology | |
CN104933235B (en) | A kind of method for merging coastal waters multi-satellite sea level height abnormal data | |
CN117928528B (en) | Ship heave measurement method based on self-adaptive time-delay-free complementary band-pass filter | |
CN111369620A (en) | Horizontal oil-water two-phase flow velocity measurement method based on improved PIV | |
CN114964219A (en) | Hybrid EMD algorithm based on parameter optimization | |
CN112710470B (en) | Self-adaptive wavelet threshold MEMS gyroscope noise reduction method based on maneuver identification | |
CN110082787A (en) | A method of extracting half Sunday of Sunday tide signal from GNSS time sequence | |
CN113702666A (en) | Signal joint noise reduction method for fiber optic gyroscope inertial measurement unit | |
CN106972840B (en) | Sampling rate conversion method and device | |
CN102589551A (en) | Real-time filtering method based on wavelet transformation for ship optical-fiber gyro signals | |
CN115950423A (en) | Ship heave motion measurement method based on adaptive filtering | |
CN111983579B (en) | Method for eliminating radar speed ambiguity by using pulse repetition time of difference | |
CN112230296B (en) | Gravity-related time reciprocal determination method | |
CN103398799A (en) | Fabry-Perot interference ring image processing method | |
CN113310470A (en) | Wave measuring method based on inertia measuring system | |
CN113433596B (en) | Gravity gradient dynamic measurement filtering method based on spatial domain | |
CN112629850B (en) | Fault diagnosis method for variable-speed planetary gearbox |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |