CN111628750A - Nonlinear filtering method for matching stochastic resonance in trap - Google Patents
Nonlinear filtering method for matching stochastic resonance in trap Download PDFInfo
- Publication number
- CN111628750A CN111628750A CN202010416149.3A CN202010416149A CN111628750A CN 111628750 A CN111628750 A CN 111628750A CN 202010416149 A CN202010416149 A CN 202010416149A CN 111628750 A CN111628750 A CN 111628750A
- Authority
- CN
- China
- Prior art keywords
- signal
- noise
- stochastic resonance
- potential well
- nonlinear
- 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
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
- H03H17/0213—Frequency domain filters using Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/141—Discrete Fourier transforms
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0248—Filters characterised by a particular frequency response or filtering method
- H03H17/0255—Filters based on statistics
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0248—Filters characterised by a particular frequency response or filtering method
- H03H17/0261—Non linear filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H2017/0072—Theoretical filter design
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Biophysics (AREA)
- General Engineering & Computer Science (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Computing Systems (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Discrete Mathematics (AREA)
- Biomedical Technology (AREA)
- Nonlinear Science (AREA)
- Databases & Information Systems (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Artificial Intelligence (AREA)
- Probability & Statistics with Applications (AREA)
- Computational Linguistics (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
- Complex Calculations (AREA)
Abstract
The invention provides a nonlinear filtering method for in-trap matching stochastic resonance, which is characterized by establishing a potential well constrained generalized second-order nonlinear model, deducing the matching relation of in-trap nonlinear parameters, establishing a parameter matching optimization method and improving the dynamic range of a filter frequency domain. Compared with the conventional stochastic resonance method, the method can improve the dynamic range of the frequency domain of the filter, breaks through the high sampling limit of the stochastic resonance, can improve the anti-noise performance, and greatly improves the flexibility and feasibility of the stochastic resonance in engineering application, especially for an embedded system with limited operational capability.
Description
Technical Field
The invention relates to the field of signal processing, in particular to a filtering and noise reduction method.
Background
Since the 21 st century, the world countries have increasingly competitive around the marine field in politics, economy and military, and have developed corresponding marine development strategies, and the protection, development and utilization of marine resources have become the focus of common world attention. Currently, the national core interests of China are mainly embodied in two aspects of economic development and safety interests, economic construction is a central task for reforming and opening the world, and ocean rights and interests are fundamentally guaranteed for realizing ocean power in a new period. Therefore, the research of the advanced weak signal processing method has great research value and practical significance for the detection and identification of the target in the long-distance water.
The filter is a general signal processing means for processing noise, and is widely applied in various fields, however, the filtering performance of the commonly used method is limited, and the processing capability for strong background noise is still insufficient. In recent years, a weak signal processing method of stochastic resonance has been drawing attention from domestic and foreign research institutes because of its enhanced characteristics for weak signals. The stochastic resonance is not denoised in a noise filtering way like the traditional weak signal processing method (high-order spectrum analysis, wavelet analysis, empirical mode decomposition analysis and the like), but noise is utilized, and a strong background noise signal is input into a special nonlinear system (resonance system), so that partial energy of the noise is converted into energy of the noise, the signal output is enhanced while the small noise energy is weakened, and the weak small signal detection can be effectively used.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a nonlinear filtering method for in-trap matching stochastic resonance, which is used for establishing a potential well constrained generalized second-order nonlinear model, deducing the matching relation of nonlinear parameters in a trap, establishing a parameter matching optimization method and improving the dynamic range of a filter frequency domain.
The technical scheme adopted by the invention for solving the technical problem comprises the following specific steps:
the first step is as follows: collecting sound signals in the sea by using a sonar, and recording the sound signals as g (t), namely input signals; the input signal simultaneously contains a mixture of a single-frequency line spectrum signal and a noise signal, namely:
g(t)=s(t)+n(t) (1)
wherein s (t) ═ Acos (2 pi f)0t), A is the input signal amplitude, f0Is the input signal frequency, n (t) is the ocean background noise signal;
the second step is that: and (3) noise intensity estimation, namely respectively estimating the noise variance of each reconstructed signal by adopting a maximum likelihood estimation method, wherein the calculation formula is as follows:
wherein the content of the first and second substances,is an estimated value of noise intensity D, N is the length of signal g (T), and T (x) is the test statistic;
the third step: constructing a second-order duffing nonlinear system:
wherein, x is the system output,is the second derivative of x and is,the first derivative of x, gamma is a damping factor, a and b are nonlinear potential parameters, and g (t) is an input signal;
the fourth step: constructing a generalized second-order matching stochastic resonance optimization model of potential well constraint:
wherein K is a potential well confinement factor, amatchAnd bmatchFor matching potential parameters, pi is the circumferential ratio, e is a natural constant,andrespectively representing the optimal potential well constraint factor and damping factor, wherein the SNRI is the output signal-to-noise ratio gain;
the fifth step: initializing and setting a potential well constraint factor and a damping factor, and setting a parameter search range; wherein the region of the potential well confinement factor isThe search interval for the damping factor is
And a sixth step: and (3) performing numerical solution on the formula (3) by using a four-order Rungestota method, wherein the initial value is determined to be (0,0), and the step length h is 1/fs,fsObtaining an output sequence x for sampling frequency, calculating a signal-to-noise ratio gain value, optimizing two parameters through a genetic algorithm until the signal-to-noise ratio gain is unchanged or the iteration times are maximum to obtain an optimal valueAnd
the seventh step: obtaining the optimal potential well constraint factor and damping factor according to the sixth stepAndcomputing optimal nonlinear filter output x using a four-order Runge Kutta methodopt。
In the fourth step, the calculation of the output signal-to-noise ratio gain is as follows:
the signal-to-noise ratio SNR defines: performing N-point Discrete Fourier Transform (DFT) on input g (t) and output x to obtain power S corresponding to each frequencyiAnd i represents an arbitrary value between 1 and N, the signal-to-noise ratio is calculated as follows:
respectively calculating the SNR of the input signal and the system output signal by using the formula (5)inAnd SNRoutThe SNR gain is then represented by the formula SNRI ═ SNRout-SNRinAnd (4) calculating.
The nonlinear filtering method for the in-trap matching stochastic resonance has the advantages that compared with a conventional stochastic resonance method, the dynamic range of a filter frequency domain can be enlarged, the high sampling limit of the stochastic resonance is broken through, the anti-noise performance can be improved, the flexibility and feasibility of the stochastic resonance in engineering application are greatly improved, and particularly, the method is suitable for an embedded system with limited operational capability.
Drawings
Figure 1 is a frequency response curve of the nonlinear filtering of the present invention.
FIG. 2 is a graph of the anti-noise response of the non-linear filtering of the present invention.
Detailed Description
The invention is further illustrated with reference to the following figures and examples.
In order to overcome the defects of the prior art, the invention provides a nonlinear filtering method for in-trap matching stochastic resonance, which is used for establishing a potential well constrained generalized second-order nonlinear model, deducing the matching relation of nonlinear parameters in a trap, establishing a parameter matching optimization method and improving the dynamic range of a filter frequency domain.
The technical scheme adopted by the invention for solving the technical problems is implemented in the following specific steps:
the first step is as follows: collecting sound signals in the sea by using a sonar, and recording the sound signals as g (t), namely input signals; the input signal containing both a single-frequency line-spectrum signal and a noise signal, i.e. mixing
g(t)=s(t)+n(t) (1)
Wherein s (t) ═ Acos (2 pi f)0t), A is the input signal amplitude, f0For the input signal frequency, n (t) is the ocean background noise signal.
The second step is that: and (3) noise intensity estimation, namely respectively estimating the noise variance of each reconstructed signal by adopting a maximum likelihood estimation method, wherein the calculation formula is as follows:
wherein the content of the first and second substances,is an estimated value of noise intensity D, N is the length of signal g (T), and T (x) is the test statistic;
the third step: constructing a second-order duffing nonlinear system:
wherein, x is the system output,is the second derivative of x and is,is the first derivative of x, gamma is the damping factor, a, b are nonlinear potential parameters, and g (t) is the input signal.
The fourth step: constructing a generalized second-order matching stochastic resonance optimization model of potential well constraint:
wherein K is a potential well confinement factor, amatchAnd bmatchTo match the potential parameters, π is the circumferential ratio, e is a natural constant, with a value of about 2.71828,andrespectively representing the optimal potential well constraint factor and the optimal damping factor, wherein the SNRI is the output signal-to-noise ratio gain, and the output signal-to-noise ratio gain is calculated by the following method:
the signal-to-noise ratio SNR defines: an N-point DFT (discrete Fourier transform) is performed on the input g (t) and the output x to obtain power S corresponding to each frequencyiAnd i represents an arbitrary value between 1 and N, the signal-to-noise ratio is calculated as follows:
Respectively calculating the SNR of the input signal and the output signal of the systeminAnd SNRoutThe SNR gain can be given by the following equation SNRI SNRout-SNRinAnd (4) calculating.
The fifth step: initializing and setting a potential well constraint factor and a damping factor, setting a parameter search range, wherein an empirical interval of the potential well constraint factor isThe search interval for the damping factor is
And a sixth step: the four-step Runge Kutta method is used for solving the numerical value of the formula (3), the initial value is determined to be (0,0), and the step length is determinedh=1/fs,fsObtaining an output sequence x for sampling frequency, calculating a signal-to-noise ratio gain value, optimizing two parameters through a genetic algorithm, and obtaining an optimal value when the signal-to-noise ratio gain is unchanged or the iteration times are maximumAndthe main parameter settings for genetic algorithm optimization are proposed as follows: population: 2. 100 population individuals, 10 maximum iterations, 0.95 cross probability and 0.01 mutation probability.
The seventh step: according to the optimal potential well constraint factor and damping factorAndcalculating the optimal nonlinear filter output by using a four-order Runge Kutta method, and using xoptAnd (4) showing.
The nonlinear filtering method for matching stochastic resonance in the trap can improve the filtering performance and the anti-noise performance of signals, the frequency response curve of nonlinear filtering is shown in figure 1, compared with the conventional stochastic resonance method, the dynamic range of the frequency domain of a filter can be improved, the high sampling limit of stochastic resonance is broken through, and the noise response curve of nonlinear filtering is shown in figure 2, compared with the conventional stochastic resonance method, the processing gain is improved under different noise intensities. The conventional Stochastic Resonance methods referenced in the figures are Haitao Dong, Haiyan Wang, Xiaohong Shen, et al. effects of second-ordered stored Resonance for Weak Signal Detection [ J ]. IEEEAccess,2018,6:46505-46515.
The method breaks through the high sampling limit of stochastic resonance, can improve the anti-noise performance, and greatly increases the flexibility and feasibility of the stochastic resonance in engineering application, especially for an embedded system with limited computing capability.
Claims (2)
1. A method of nonlinear filtering of in-trap matched stochastic resonance, comprising the steps of:
the first step is as follows: collecting sound signals in the sea by using a sonar, and recording the sound signals as g (t), namely input signals; the input signal simultaneously contains a mixture of a single-frequency line spectrum signal and a noise signal, namely:
g(t)=s(t)+n(t) (1)
wherein s (t) ═ Acos (2 pi f)0t), A is the input signal amplitude, f0Is the input signal frequency, n (t) is the ocean background noise signal;
the second step is that: and (3) noise intensity estimation, namely respectively estimating the noise variance of each reconstructed signal by adopting a maximum likelihood estimation method, wherein the calculation formula is as follows:
wherein the content of the first and second substances,is an estimated value of noise intensity D, N is the length of signal g (T), and T (x) is the test statistic;
the third step: constructing a second-order duffing nonlinear system:
wherein, x is the system output,is the second derivative of x and is,the first derivative of x, gamma is a damping factor, a and b are nonlinear potential parameters, and g (t) is an input signal;
the fourth step: constructing a generalized second-order matching stochastic resonance optimization model of potential well constraint:
wherein K is a potential well confinement factor, amatchAnd bmatchFor matching potential parameters, pi is the circumferential ratio, e is a natural constant,andrespectively representing the optimal potential well constraint factor and damping factor, wherein the SNRI is the output signal-to-noise ratio gain;
the fifth step: initializing and setting a potential well constraint factor and a damping factor, and setting a parameter search range; wherein the region of the potential well confinement factor isThe search interval for the damping factor is
And a sixth step: and (3) performing numerical solution on the formula (3) by using a four-order Rungestota method, wherein the initial value is determined to be (0,0), and the step length h is 1/fs,fsObtaining an output sequence x for sampling frequency, calculating a signal-to-noise ratio gain value, optimizing two parameters through a genetic algorithm until the signal-to-noise ratio gain is unchanged or the iteration times are maximum to obtain an optimal valueAnd
2. The method of claim 1, wherein the method comprises:
in the fourth step, the calculation of the output signal-to-noise ratio gain is as follows:
the signal-to-noise ratio SNR defines: performing N-point Discrete Fourier Transform (DFT) on input g (t) and output x to obtain power S corresponding to each frequencyiAnd i represents an arbitrary value between 1 and N, the signal-to-noise ratio is calculated as follows:
wherein Sf0Is the frequency f of the signal s (t)0The corresponding power value;
respectively calculating the SNR of the input signal and the system output signal by using the formula (5)inAnd SNRoutThe SNR gain is then represented by the formula SNRI ═ SNRout-SNRinAnd (4) calculating.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010416149.3A CN111628750B (en) | 2020-05-17 | 2020-05-17 | Nonlinear filtering method for random resonance matching in trap |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010416149.3A CN111628750B (en) | 2020-05-17 | 2020-05-17 | Nonlinear filtering method for random resonance matching in trap |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111628750A true CN111628750A (en) | 2020-09-04 |
CN111628750B CN111628750B (en) | 2023-04-18 |
Family
ID=72259809
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010416149.3A Active CN111628750B (en) | 2020-05-17 | 2020-05-17 | Nonlinear filtering method for random resonance matching in trap |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111628750B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502076A (en) * | 2023-06-30 | 2023-07-28 | 浙江大学 | Weak signal identification method and device and readable storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549105A (en) * | 2018-07-02 | 2018-09-18 | 西北工业大学 | A kind of seismic wave characteristic extracting method of self-adapting random resonant |
CN108875685A (en) * | 2018-07-02 | 2018-11-23 | 西北工业大学 | A kind of underwater AUV detection method of Adaptive matching accidental resonance |
CN110376575A (en) * | 2019-08-19 | 2019-10-25 | 西北工业大学 | A kind of low frequency spectrum lines detection method based on damping parameter matching accidental resonance |
-
2020
- 2020-05-17 CN CN202010416149.3A patent/CN111628750B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108549105A (en) * | 2018-07-02 | 2018-09-18 | 西北工业大学 | A kind of seismic wave characteristic extracting method of self-adapting random resonant |
CN108875685A (en) * | 2018-07-02 | 2018-11-23 | 西北工业大学 | A kind of underwater AUV detection method of Adaptive matching accidental resonance |
CN110376575A (en) * | 2019-08-19 | 2019-10-25 | 西北工业大学 | A kind of low frequency spectrum lines detection method based on damping parameter matching accidental resonance |
Non-Patent Citations (2)
Title |
---|
HAITAO DONG,HAIYAN WANG: "《Effects of Second-Order Matched Stochastic Resonance for Weak Signal Detection》", 《IEEE ACCESS》, 22 August 2018 (2018-08-22) * |
朱维娜,林敏: "《基于人工鱼群算法的轴承故障随机共振自适应检测方法》", 《振动与冲击》, 26 June 2014 (2014-06-26) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116502076A (en) * | 2023-06-30 | 2023-07-28 | 浙江大学 | Weak signal identification method and device and readable storage medium |
CN116502076B (en) * | 2023-06-30 | 2023-10-17 | 浙江大学 | Weak signal identification method and device and readable storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN111628750B (en) | 2023-04-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110376575B (en) | Low-frequency line spectrum detection method based on damping parameter matching stochastic resonance | |
WO2021056727A1 (en) | Joint noise reduction method based on variational mode decomposition and permutation entropy | |
CN104881567A (en) | Statistical model based bridge health monitoring data wavelet denoising method | |
CN106328120B (en) | Method for extracting abnormal sound features of public places | |
CN103995950A (en) | Wavelet coefficient partial discharge signal noise elimination method based on related space domain correction threshold values | |
CN102353952A (en) | Line spectrum detection method by coherent accumulation of frequency domains | |
CN111628750B (en) | Nonlinear filtering method for random resonance matching in trap | |
CN114992033B (en) | NLM-CEEMDAN-based hydroelectric generating set signal denoising method | |
Wang et al. | Joint noise and mask aware training for DNN-based speech enhancement with sub-band features | |
CN111638501B (en) | Spectral line enhancement method for self-adaptive matching stochastic resonance | |
CN114038470A (en) | Method, device, equipment and medium for extracting voiceprint features of transformer signals | |
CN110515063A (en) | Underwater acoustic signal processing method and apparatus based on the steady wavelet transform of iteration | |
CN115859054B (en) | MIC-CEEMDAN-based hydroelectric generating set draft tube pressure pulsation signal denoising method | |
Lei et al. | High voltage shunt reactor acoustic signal denoising based on the combination of VMD parameters optimized by coati optimization algorithm and wavelet threshold | |
CN117609702A (en) | Pipeline leakage acoustic emission signal denoising method, system, equipment and medium | |
CN117111155A (en) | Microseism data denoising method based on integrated framework | |
CN112539887A (en) | WT-LCD-WD-based pipeline leakage signal denoising method | |
Zehtabian et al. | A novel speech enhancement approach based on singular value decomposition and genetic algorithm | |
CN112926504A (en) | Acoustic emission signal denoising method based on noise reduction self-encoder | |
CN113065387B (en) | Wavelet denoising method and system for pump system | |
CN110703089B (en) | Wavelet threshold denoising method for low-frequency oscillation Prony analysis | |
Yang et al. | Iterative Morlet wavelet with SOSO boosting strategy for impulsive feature extraction | |
Yang et al. | Noise Modeling and Deep Learning Noise Suppression of Mud Signal | |
CN109871784B (en) | Full-wave nuclear magnetic resonance signal noise filtering method for optimizing matching pursuit by genetic algorithm | |
CN103854658A (en) | Steel plate corrosion acoustic emission signal de-noising method based on short-time fractal dimension enhancing method |
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 |