CN111178325A - Signal extraction and processing method - Google Patents
Signal extraction and processing method Download PDFInfo
- Publication number
- CN111178325A CN111178325A CN202010028226.8A CN202010028226A CN111178325A CN 111178325 A CN111178325 A CN 111178325A CN 202010028226 A CN202010028226 A CN 202010028226A CN 111178325 A CN111178325 A CN 111178325A
- Authority
- CN
- China
- Prior art keywords
- signal
- norm
- input signal
- variable
- maximum
- 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.)
- Withdrawn
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/08—Feature extraction
Landscapes
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Artificial Intelligence (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Signal Processing For Digital Recording And Reproducing (AREA)
Abstract
The invention discloses a signal extraction and processing method, which relates to the technical field of signal extraction and processing, and specifically comprises the following steps: firstly, extracting upper and lower outlines of an input signal, averaging, carrying out Newton iteration to adjust norm to minimize the norm to obtain a variable-mode signal, then carrying out reduction operation on the variable-mode signal and the input signal for multiple times to obtain multiple groups of variable-mode signals, and carrying out Hilbert change on the multiple groups of variable-mode signals to obtain the amplitude and phase angle of all the variable-mode signals. The invention decomposes the signal to obtain the amplitude and the phase angle of the variable-mode signal, lays a solid foundation for the further processing of the subsequent signal, solves the defects of poor robustness, poor separation effect and easy error introduction of the traditional method, and can be widely applied to the existing filtering methods such as noise elimination, mode identification and the like.
Description
Technical Field
The invention relates to the technical field of signal processing, in particular to a method for extracting and processing signals.
Background
With the continuous development of the digital signal processing field, the requirements of people on signal quality are continuously improved, in order to ensure the accurate extraction of signals in a complex environment and avoid the loss of important information, at present, various methods are adopted, but how to remove interference signals and clutter signals in original signals to enable the signals to meet the follow-up requirements is a troublesome problem. The modulation mode of the signal is to superpose a high-frequency target signal on a selected low-frequency carrier signal, and during signal transmission, clutter signals influencing the accurate extraction and processing of the signal are used. The traditional denoising method generally performs Fourier transform on signals, and the processing mode adopting the Fourier transform is more suitable for stable signals and has poor processing effect on non-stable signals. The processing of non-stationary signals is mainly to remove the existing clutter signals directly before demodulation, and requires either a high signal-to-noise ratio or a clean original signal. However, in the actual design, the algorithm is vague and difficult to understand, is not easy to implement, has poor effect, and is easy to lose important information, and aliasing signals are processed according to the method, because the aliasing signals are mutually overlapped on a time frequency domain and a frequency domain, the signals cannot be separated by the traditional time domain or frequency domain filtering method, so that the method for separating or extracting the aliasing signals has important theoretical significance and practical application value.
Disclosure of Invention
In order to solve the above problems, the present invention provides a method for extracting and processing signals.
The specific technical scheme provided by the invention is as follows: a method of signal extraction and processing, comprising the steps of:
s1, obtaining an original signal x (t), and simultaneously setting a Z value, wherein the Z value is the number of the variable modulus signals to be obtained;
s2, inputting an original signal x (t) as an input signal u (t), and setting i = 1;
s3, extracting the upper and lower profiles of the input signal u (t) and averaging to obtain the characteristic signal of the input signal u (t);
S4, judging the characteristic signalWhether the difference between the number of the medium maximum value points and the number of the minimum value points is 1 or 0: if not, the characteristic signal is sentReturn to execution of step S3 as an input signal;
s5, setting a proportionality coefficient and aiming at the maximum characteristic signalFitting the maximum value point and the minimum value point to obtain a signalThen make a request forWith input signal u (t)A norm; using Newton iteration to fine-tune the norm to minimize the normAs a variable mode signal;
S6, judging the mode-changing signalIs equal to Z: if yes, sequentially outputting all the variable modulus signals; if not, the mode-changing signal is adjustedPerforming a subtraction operation with the input signal u (t) as a new input signal, returning to step S3, while letting i = i + 1;
s7, performing Hilbert transform on all the variable-mode signals output in the step S6, and then fusing to obtain the amplitude and phase angle of the variable-mode signals
2. Extraction of the signal according to claim 1And a processing method, wherein the step S3 extracts the upper and lower contours of the input signal u (t) and averages them to obtain a characteristic signal of the input signal u (t)The method comprises the following steps: collecting a maximum value point and a minimum value point of an input signal u (t), and recording the maximum value point as Xmax(1) To Xmax(n 1), the minimum value point is denoted as Xmin(1) To Xmin(n 2), n1 and n2 being the number of maximum and minimum points, respectively;
maximum value point X to be collectedmax(1) To Xmax(n 1) fitting to obtain the upper envelope x of the input signal u (t)1(t) collecting the minimum value point Xmin(1) To Xmin(n 2) fitting to obtain a lower envelope x of x (t)2(t) obtaining the upper envelope x by the formula (1)1(t) and the lower envelope x2(t) obtaining a characteristic signal from the mean value,
Further, the upper envelope x1(t) and the lower envelope x2(t) the fitting formula (1) is:
further, the setting of the scaling factor is performed for the maximum characteristic signalFitting the maximum value point and the minimum value point to obtain a fitted signalThe method comprises the following steps:
s11, extracting maximum characteristic signalIs recorded as the maximum value point ofToMinimum value point is marked asTo,Andrespectively representing the number of maximum points and minimum points;
s12, k1 toEach of the scale coefficients in the set of scale coefficients corresponds to a maximum point of the maximum characteristic signalToOr minimum pointToMultiplying, fitting by equation (2) to obtain a signalWherein k1 is set toThe initial values of (a) are all 1;
wherein the content of the first and second substances,and ts is the norm, and the termination time of the current data acquisition.
Further, the performing norm fine-tuning to minimize the norm in step S5 using newton iteration includes:
s13, iterating the norm through a formula (4) to obtain new proportionality coefficients k1 toStep S12 is executed to obtain a new oneAccording to newCalculating the new value again according to equation (3)With input signal u (t)Norm, and repeatedly executing S13 until the norm is minimum after iteration u times;
wherein m isiAnd (t) is a variable-mode signal, and t and w are respectively a time axis and a frequency axis of the current result.
The invention has the beneficial effects that:
the invention provides a signal extraction and processing method, which solves the defects that the traditional method is poor in robustness and easy to introduce errors, so that the detection result is closer to a true value, aliasing signals are decomposed efficiently, and the amplitude values and phase angles of all variable-mode signals of the original signals are obtained. The method has the advantages of simple algorithm and accurate separation, two broadband signal waveforms containing a plurality of frequency components can be separated by simple algorithm execution, and the separation of a plurality of signal components, namely the separation of multiple input and multiple output signals, can be completed by multiple times of algorithm execution.
In addition to the objects, features and advantages described above, other objects, features and advantages of the present invention are also provided. The present invention will be described in further detail below.
Description of the drawings:
fig. 1 is a flow chart of a signal extraction and processing method according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail with reference to the following embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
The signal extracting and processing method provided by any embodiment of the invention comprises the following steps:
s1, obtaining an original signal x (t), and simultaneously setting a Z value, wherein the Z value is the number of the variable modulus signals to be obtained;
a signal extracting and processing method according to any embodiment of the present invention, wherein the original signal x (t) may be any continuous mixed signal, for example, radar echo, bio-impedance signal;
the Z value may be set to any one of 1, 2, 3, 4, and 100.
S2, inputting an original signal x (t) as an input signal u (t), and setting i = 1;
s3, extracting the upper and lower profiles of the input signal u (t) and averaging to obtain the characteristic signal of the input signal u (t);
In a signal extracting and processing method according to any embodiment of the present invention, the upper and lower contours of the input signal u (t) are extracted and averaged to obtain the characteristic signal of the input signal u (t)The method comprises the following steps: collecting a maximum value point and a minimum value point of an input signal u (t), and recording the maximum value point as Xmax(1) To Xmax(n 1), the minimum value point is denoted as Xmin(1) To Xmin(n 2), n1 and n2 being the number of maximum and minimum points, respectively;
maximum value point X to be collectedmax(1) To Xmax(n 1) fitting by equation (2) to obtain the upper envelope of the input signal u (t)Line x1(t) collecting the minimum value point Xmin(1) To Xmin(n 2) fitting by equation (2) to obtain the lower envelope x of x (t)2(t);
Calculating the upper envelope line x by formula (1)1(t) and the lower envelope x2(t) obtaining a characteristic signal from the mean value,
S4, judging the characteristic signalWhether the difference between the number of the medium maximum value points and the number of the minimum value points is 1 or 0: if not, the characteristic signal is sentReturn to execution of step S3 as an input signal;
s5, setting a proportionality coefficient and aiming at the maximum characteristic signalFitting the maximum value point and the minimum value point to obtain a signalThen make a request forWith input signal u (t)A norm; using Newton iteration to fine-tune the norm to minimize the normAs a variable mode signal;
A signal extracting and processing method according to any of the embodiments of the present invention, wherein the scaling factor is set to k1 toTo ensure the continuity of the signal, k1 is set toAll are 1, and the maximum characteristic signal is obtainedIs recorded as the maximum value point ofToMinimum value point is marked asTo,Andeach representing a maximum point and a minimum pointCounting;
k1 toEach of the scale coefficients in the set of scale coefficients corresponds to a maximum point of the maximum characteristic signalToOr minimum valueToMultiplying, fitting by equation (2) to obtain a signal。
Exemplary, k1Andmultiplication, k2Andmultiplication,,,,,, kplus-1Andmultiplication, kplusAndmultiplication.
A method of extracting and processing a signal according to any of the embodiments of the present invention, wherein the method comprisesWith input signal u (t)The norm is given by equation (3):
wherein the content of the first and second substances,and ts is the norm, and the termination time of the current data acquisition.
In the method for extracting and processing a signal according to any embodiment of the present invention, a specific way of finely adjusting a norm by using newton iteration to minimize the norm is as follows:
the norm is iterated through equation (4) to obtain new scaling factors k1 toNew k1 toEach of the scale coefficients is in one-to-one correspondence withMaximum value point ofToOr minimum valueToMultiplying, fitting by equation (2) to obtain a new signalAccording to newCalculating the new value again according to equation (3)With input signal u (t)Performing norm, and repeatedly executing the step until the norm is minimum after iteration for u times;
s6, judging the mode-changing signalIs equal to Z: if yes, sequentially outputting all the variable modulus signals; if not, the mode-changing signal is adjustedPerforming a subtraction operation with the input signal u (t) as a new input signal, returning to step S3, while letting i = i + 1;
illustratively, when the value of Z is set to 4, when i = Z, all the modulus-varying signals output in sequence are、、、(ii) a When the value of Z is set to be 10, when i = Z, all the modulus-changing signals output in sequence are、、、、、、、、、。
S7, performing Hilbert transform on all the variable-mode signals output in the step S6, and then fusing to obtain the amplitude and phase angle of the variable-mode signals
A method for extracting and processing signals according to any embodiment of the present invention, wherein the amplitude and phase angle of the modulus-variable signal are calculatedEquation (5) of (a) is:
wherein m isiAnd (t) is a variable-mode signal, and t and w are respectively a time axis and a frequency axis of the current result.
The above description is only for the purpose of illustrating the preferred embodiments of the present invention and should not be taken as limiting the invention, but rather the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the present invention.
Claims (6)
1. A method for extracting and processing signals, comprising the steps of:
s1, obtaining an original signal x (t), and simultaneously setting a Z value, wherein the Z value is the number of the variable modulus signals to be obtained;
s2, inputting an original signal x (t) as an input signal u (t), and setting i = 1;
s3, extracting the upper and lower profiles of the input signal u (t) and averaging to obtain the characteristic signal of the input signal u (t);
S4, judging the characteristic signalWhether the difference between the number of the medium maximum value points and the number of the minimum value points is 1 or 0: if not, the characteristic signal is sentReturn to execution of step S3 as an input signal;
s5, setting a proportionality coefficient and aiming at the maximum characteristic signalFitting the maximum value point and the minimum value point to obtain a signalThen make a request forWith input signal u (t)A norm; fine tuning norm using newton iterationObtained when norm is minimizedAs a variable mode signal;
S6, judging the mode-changing signalIs equal to Z: if yes, sequentially outputting all the variable modulus signals; if not, the mode-changing signal is adjustedPerforming a subtraction operation with the input signal u (t) as a new input signal, returning to step S3, while letting i = i + 1;
2. The method for extracting and processing the signal according to claim 1, wherein the step S3 is to extract the upper and lower contours of the input signal u (t) and average the upper and lower contours to obtain the characteristic signal of the input signal u (t)The method comprises the following steps: collecting a maximum value point and a minimum value point of an input signal u (t), and recording the maximum value point as Xmax(1) To Xmax(n 1), the minimum value point is denoted as Xmin(1) To Xmin(n 2), n1 and n2 being the number of maximum and minimum points, respectively;
maximum value to be collectedPoint Xmax(1) To Xmax(n 1) fitting to obtain the upper envelope x of the input signal u (t)1(t) collecting the minimum value point Xmin(1) To Xmin(n 2) fitting to obtain a lower envelope x of x (t)2(t) obtaining the upper envelope x by the formula (1)1(t) and the lower envelope x2(t) obtaining a characteristic signal from the mean value,
4. a method for extracting and processing signals according to claim 3, wherein the scaling factor is set for the maximum feature signalFitting the maximum value point and the minimum value point to obtain a fitted signalThe method comprises the following steps:
s11, extracting maximum characteristic signalIs recorded as the maximum value point ofToMinimum value point is marked asTo,Andrespectively representing the number of maximum points and minimum points;
s12, k1 toEach of the scale coefficients in the set of scale coefficients corresponds to a maximum point of the maximum characteristic signalToOr minimum pointToMultiplying, fitting by equation (2) to obtain a signalWherein k1 is set toThe initial values of (a) are all 1;
5. The method for extracting and processing the signal according to claim 4, wherein the step S5 of performing norm fine-tuning using newton iteration to minimize the norm comprises:
s13, iterating the norm through a formula (4) to obtain new proportionality coefficients k1 toStep S12 is executed to obtain a new oneAccording to newCalculating the new value again according to equation (3)With input signal u (t)Norm, and repeatedly executing S13 until the norm is minimum after iteration u times;
6. method for extracting and processing signals according to claim 1, characterized in that the amplitude and phase angle of said modulus-varying signalEquation (5) for (t, w) is:
wherein m isiAnd (t) is a variable-mode signal, and t and w are respectively a time axis and a frequency axis of the current result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010028226.8A CN111178325A (en) | 2020-01-10 | 2020-01-10 | Signal extraction and processing method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010028226.8A CN111178325A (en) | 2020-01-10 | 2020-01-10 | Signal extraction and processing method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111178325A true CN111178325A (en) | 2020-05-19 |
Family
ID=70656281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010028226.8A Withdrawn CN111178325A (en) | 2020-01-10 | 2020-01-10 | Signal extraction and processing method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111178325A (en) |
-
2020
- 2020-01-10 CN CN202010028226.8A patent/CN111178325A/en not_active Withdrawn
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108680786B (en) | Adaptive filtering envelope extraction method for pulse signal frequency domain | |
CN108364659B (en) | Frequency domain convolution blind signal separation method based on multi-objective optimization | |
CN108469602B (en) | Pulse signal type automatic discrimination method based on spectral feature extraction | |
CN109285561B (en) | Ship propeller cavitation noise modulation spectrum feature fidelity enhancement method based on self-adaptive window length | |
CN109557429A (en) | Based on the GIS partial discharge fault detection method for improving wavelet threshold denoising | |
TW438981B (en) | Method for separating a desired signal from an interfering signal | |
CN114025379B (en) | Broadband multi-signal detection method, device and equipment | |
WO2014043833A1 (en) | Image denoising method based on phase consistency for magnetic resonance imaging | |
CN114325598A (en) | Signal denoising method and device | |
CN114895248A (en) | Sinusoidal frequency modulation signal parameter estimation method, system and medium | |
CN109087657B (en) | Voice enhancement method applied to ultra-short wave radio station | |
CN111178325A (en) | Signal extraction and processing method | |
CN115835192B (en) | Accompanying carrying type hidden communication method, system, medium, equipment and terminal | |
CN111611686A (en) | Detection method for communication signal time-frequency domain | |
CN109460614B (en) | Signal time-frequency decomposition method based on instantaneous bandwidth | |
CN113627398B (en) | Signal characteristic detection method based on adaptive reconstruction filtering | |
CN110210081A (en) | A kind of SS-OCT system k-clock delay correcting algorithm | |
CN113325401B (en) | Distortion towing linear array signal reconstruction method based on line spectrum phase difference deblurring | |
CN210347836U (en) | Wave velocity calculating device for cable oscillatory wave partial discharge detection system | |
CN110198196B (en) | Amplitude estimation method based on signal intensity in communication system | |
CN114674410A (en) | Component number time-varying underwater acoustic signal instantaneous frequency estimation method | |
CN113822363A (en) | Convolutional neural network time delay estimation method based on quadratic rational kernel function | |
CN109683142B (en) | Method for estimating parameters of triangular linear frequency modulation continuous signals based on differential envelope detection | |
CN110426711B (en) | Time delay estimation method and system based on polarity zero detection | |
CN114062992A (en) | Ultrasonic extraction method and device for primary side detection process of current transformer |
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 | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200519 |
|
WW01 | Invention patent application withdrawn after publication |