CN109632974B

CN109632974B - Echo signal separation method for ultrasonic flaw detection

Info

Publication number: CN109632974B
Application number: CN201811597258.9A
Authority: CN
Inventors: 齐爱玲
Original assignee: Xian University of Science and Technology
Current assignee: Xian University of Science and Technology
Priority date: 2018-12-26
Filing date: 2018-12-26
Publication date: 2021-07-09
Anticipated expiration: 2038-12-26
Also published as: CN109632974A

Abstract

The invention discloses an echo signal separation method for ultrasonic flaw detection, which comprises the following steps: firstly, ultrasonic echo signal acquisition and synchronous uploading: ultrasonic detection is carried out on the object to be detected by adopting an ultrasonic flaw detection device, an ultrasonic echo signal F (t) of the object to be detected is obtained, and the obtained ultrasonic echo signal F (t) is synchronously transmitted to data processing equipment; the ultrasonic echo signals F (t) comprise ultrasonic echo signals at the positions of K defects in the tested object; secondly, signal sparse decomposition is carried out; thirdly, echo signal separation: and separating the ultrasonic echo signals F (t) by adopting data processing equipment according to the K optimal matching atoms in the optimal atom set by iterative decomposition at the moment to obtain K separation signals. The method has the advantages of simple steps, reasonable design, convenient implementation and good use effect, and realizes simple, convenient and quick separation of the superposed signals at the positions of the defects in the same ultrasonic echo signal by utilizing signal sparse decomposition.

Description

Echo signal separation method for ultrasonic flaw detection

Technical Field

The invention belongs to the technical field of ultrasonic flaw detection, and particularly relates to an echo signal separation method for ultrasonic flaw detection.

Background

Ultrasonic inspection (also called ultrasonic inspection or ultrasonic detection) is a nondestructive inspection method for inspecting internal defects of a material by using the difference of acoustic properties of the material and the defects thereof to the energy change of the reflection condition and the penetration time of an ultrasonic propagation waveform. The principle of ultrasonic flaw detection is that transmission loss of ultrasonic waves in a solid is small, the detection depth is large, the ultrasonic waves can reflect and refract on a heterogeneous interface, particularly, the ultrasonic waves can not pass through a gas-solid interface, if a metal has defects (gas in the defects) such as pores, cracks, layering and the like or impurities, the ultrasonic waves (also called ultrasonic signals or ultrasonic signals) can be totally or partially reflected when being transmitted to the interface between the metal and the defects, the reflected ultrasonic signals (also called ultrasonic echo signals, referred to as echo signals for short) are received by a probe, and the depth, the position and the shape of the defects in a workpiece can be judged according to the waveform change characteristics of the received ultrasonic echo signals.

Due to the influence of factors such as the characteristics of metal materials, inherent defects of machining processes and the like, internal defects of different degrees inevitably exist in the production process of industrial mechanical equipment, and common defects comprise cracks, air holes, shrinkage cavities, inclusions, sand inclusion cold shut and the like. For example, the coal mining machinery such as a coal mining machine box body inevitably has the defects in the production process, and as the manufacturing process of the coal mining machinery is complex, and the coal mining machinery works under the conditions of heavy load and alternating load for a long time, the working environment is severe, the stress condition of the equipment is poor, so that important parts of the equipment are gradually damaged from small to large in accidents in the using process, various safety hazards are generated, and particularly, the production is stopped due to the generation of the faults, and the yield of coal and the economic benefit of coal mines are directly influenced.

Nowadays, an ultrasonic detection method is widely used in the defect detection of a box body of a coal mining machine. Ultrasonic inspection is an important nondestructive inspection method, an ultrasonic signal is a broadband pulse signal modulated by the center frequency of a probe, the echo signal of the ultrasonic signal contains a large amount of information related to defects, but the ultrasonic signal (i.e. the ultrasonic echo signal) is often polluted by random noise and related noise of a measurement system and a test workpiece, particularly grain noise in coarse-grained materials, and the noise can make defect identification of the ultrasonic signal difficult, and seriously limits the precision and reliability of defect detection. Therefore, denoising the ultrasonic detection echo signal is very important to ensure the authenticity of the obtained defect signal. The method has important significance for improving the product defect detection rate, ensuring the product quality and prolonging the service life of the product for enterprises.

From the above, in the process of ultrasonic detection of defects of coal mining machines and other coal mining machines, the reliability and quality of the detection result are seriously affected by noise. The extraction of defect signals from signals mixed with various interference noises (i.e. ultrasonic detection signals) is the key to ensure the accuracy of echo signals, and when the defects in the material are detected by using ultrasonic waves, the defect information is represented by the amplitude, frequency or phase of the received ultrasonic echo signals. The defect signal refers to a detected ultrasonic echo signal, and the ultrasonic echo signal contains defect information. However, due to the influence of instrument noise and test environment factors, various interference noises such as various electrical noises, structural noises and pulse noises are accompanied with detection signals, and particularly when the microstructure noise of defect signals is very large or the crystal grains of the material are coarse, the defect and the signal intensity of the noise are weak, and the extraction of the defect signals under the background of strong noise is a difficult problem in the research of the field of ultrasonic signal processing. How to extract the desired information from these signals is a difficult and important issue. Therefore, denoising the ultrasonic detection echo signal is very important to ensure the authenticity of the obtained defect signal.

At present, a plurality of ultrasonic signal extraction methods, such as nonlinear filtering, fourier transform, wavelet transform and the like, are available, and the methods have a good effect of improving the signal-to-noise ratio of general ultrasonic signals, but have limitations on the extraction of defects under the background of small defects or strong noise, inaccurate detection results and low reliability. The sparse decomposition is a new signal analysis theory, a proper expansion function can be selected in a self-adaptive mode according to the characteristics of a signal to be extracted, the basic characteristics of the extracted signal can be represented by few functions, weak signals can be extracted better under the condition of low signal-to-noise ratio, and the original signal is approached to the maximum. The sparse decomposition algorithm was first proposed by Mal lat, which is a well-known matching pursuit algorithm. Therefore, the method is gradually popularized in the detection of the internal defects of the products. However, the algorithm has two defects, namely, the calculation amount of the sparse decomposition algorithm is large, the calculation time is very large under the existing calculation condition, and real-time detection cannot be carried out; secondly, the sparse decomposition algorithm is an optimal solution obtained under a continuous condition, and the detection precision of weak and small defects is still limited.

In ultrasonic nondestructive testing, an ultrasonic sensor in an ultrasonic flaw detection device sends pulses to pass through a discontinuous interface of a tested object, and received reflected echoes (namely ultrasonic echo signals) contain position information and flaw size information related to flaw characteristics, so that accurate detection of the position and size of a flaw and accurate estimation of the flaw are important contents of ultrasonic nondestructive evaluation. The ultrasonic echo signal is represented as a superimposed combination of defect waves (also referred to as defect signals) reflected from different interfaces at different depth positions along the time axis. In practical applications, a situation that one ultrasonic echo signal contains a plurality of defects often occurs, the positions of the plurality of defects are overlapped, and accordingly, the defect signals corresponding to the plurality of defects in the ultrasonic echo signal are also overlapped, and at this time, the defect signals corresponding to the plurality of defects in the ultrasonic echo signal need to be separated, otherwise, the defect signals cannot be extracted, and the defect classification and identification cannot be performed. After a plurality of defect signals in the ultrasonic echo signal are separated, respectively carrying out signal extraction, feature extraction and defect classification identification on each separated defect signal.

In addition, when a sparse decomposition method is adopted for signal drying (also called signal extraction) at present, the essence of the method is that a matching pursuit algorithm is adopted for signal extraction, and in the matching pursuit algorithm, firstly, the selection of Gabor atoms has no physical significance, and secondly, if an error Gabor atom is selected in a certain iteration, the error has to be corrected in the subsequent iteration process, so that the error tracking is caused, and the algorithm performance is poor. Moreover, according to the sparse signal representation theory, when the size of a dictionary atom is smaller than a threshold value, the sparse representation is unstable and unreliable due to the constraint of accumulative consistency. Thus, the existing sparse decomposition method has many defects.

Disclosure of Invention

The technical problem to be solved by the invention is to provide an echo signal separation method for ultrasonic flaw detection aiming at the defects in the prior art, the method has the advantages of simple steps, reasonable design, convenient implementation and good use effect, and the simple and fast separation of superposed signals at the positions of a plurality of defects in the same ultrasonic echo signal is realized by utilizing signal sparse decomposition.

In order to solve the technical problems, the invention adopts the technical scheme that: an echo signal separation method for ultrasonic flaw detection, characterized by comprising the steps of:

step one, ultrasonic echo signal acquisition and synchronous uploading: ultrasonic detection is carried out on the object to be detected by adopting an ultrasonic flaw detection device, an ultrasonic echo signal F (t) of the object to be detected is obtained, and the obtained ultrasonic echo signal F (t) is synchronously transmitted to data processing equipment;

wherein f (t) ═ f (t)₁),f(t₂),...,f(t_N')]^TT represents a time parameter, t_iThe ith sampling time of the ultrasonic flaw detector, f (t)_i) Signal value sampled at ith sampling time of ultrasonic flaw detection deviceI is a positive integer and i ═ 1, 2, 3, …, N' is a positive integer and it is the signal length of the ultrasound echo signal f (t);

the ultrasonic echo signals F (t) comprise ultrasonic echo signals at the positions of K defects in the measured object;

step two, signal sparse decomposition: adopting data processing equipment and calling a sparse decomposition module to carry out iterative decomposition processing on the ultrasonic echo signal F (t) in the step one, and converting the ultrasonic echo signal F (t) into the ultrasonic echo signal F (t)

And obtaining the optimal atom set of iterative decomposition at the moment; the iterative decomposition of the best set of atoms at this time contains K best matching atoms,

decomposing for the iteration the z-th best matching atom in the set of best atoms;

in the formula R^K(t) is a residual quantity of the ultrasonic echo signal f (t) after K iterative decompositions, where K is a preset total iterative decomposition number and K is a positive integer, z is a positive integer and z is 1, 2, …, and K; a is_zThe expansion coefficient of the best matching atom after the z-th iterative decomposition and the residual error after the last iterative decomposition is obtained;

for the best matching atom found at the z-th iterative decomposition,

is a Gabor atom and

wherein the function ψ (t) is a Gaussian window function and

r_zis composed of

Of time-frequency parameters r_z＝(s_z,u_z,v_z,w_z)，s_zAs a scale parameter, u_zAs a displacement parameter, v_zAs a frequency parameter, w_zIs a phase parameter;

step three, echo signal separation: separating the ultrasonic echo signals F (t) in the first step by adopting data processing equipment according to the K optimal matching atoms in the iterative decomposition optimal atom set at the moment to obtain K separation signals;

the K' th of the K separation signals is denoted as F_k'(t)，

Wherein the content of the first and second substances,

to this end, the iterative decomposition of the K ' th best matching atom in the best set of atoms, K ' being a positive integer and K ' ═ 1, 2, …, K;

represents R^k'^-1(t) and

inner product of (d); r^k'^-1(t) is a residual quantity R after calling a matching tracking algorithm module and performing k '-1 iterative decomposition on the ultrasonic echo signal F (t) by utilizing the first k' -1 best matching atoms in the iterative decomposition best atom set at the moment⁰(t)＝F(t)。

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: in the first step, the ultrasonic flaw detection device is connected with the data processing equipment in a wired or wireless mode.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: in the fourth step, each separation signal is an ultrasonic echo signal at the position of a defect in the detected object;

overlapping ultrasonic echo signals at positions where K defects in the measured object are located in the first step;

in the first step, the ultrasonic echo signals f (t) are signals obtained by superimposing the ultrasonic echo signals at the positions of the K defects in the object to be tested, and each of the separated signals in the fourth step is a superimposed signal.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: before signal sparse decomposition is carried out in the second step, determining the number K of the superposed signals in the ultrasonic echo signal F (t);

when the number K of the superimposed signals in the ultrasonic echo signals f (t) is determined, a data processing device is used for determining, and the process is as follows:

step B1, determining wave crests and wave troughs: respectively determining all wave crests and all wave troughs of the ultrasonic echo signal F (t) in the step one by adopting data processing equipment, and respectively synchronously recording the sampling time and the signal value of each determined wave crest and each determined wave trough;

in this step, each determined peak and each determined valley are an extreme point of the ultrasonic echo signal f (t);

step B2, removing extreme points: adopting data processing equipment and calling a time domain extreme point eliminating module or a frequency domain extreme point eliminating module to eliminate extreme points to obtain M 'eliminated extreme points, and arranging the M' extreme points from front to back according to the sampling time sequence of the extreme points; wherein M' is a positive integer and is the total number of extreme points obtained after the extreme points are removed in the step;

when the data processing equipment is adopted and a time domain extreme point rejecting module is called to reject the extreme points, rejecting the extreme points of which the absolute values of the signal values in all the extreme points determined in the step B1 are smaller than beta 'to obtain M' number of rejected extreme points; wherein, β' is a constant and the value range thereof is 0.1-0.35, and max | f (t) | is the maximum absolute value of the signal value in the ultrasonic echo signal f (t);

when the data processing equipment is adopted and a frequency domain extreme point eliminating module is called to eliminate the extreme points, eliminating the extreme points of which the absolute values of the signal values are smaller than beta in all the extreme points determined in the step B1 to obtain M' eliminated extreme points; wherein β is a preset rejection threshold, β ═ α × max | y (f) |, α is a constant and has a value range of 0.25 to 0.35, y (f) is a frequency spectrum of the ultrasonic echo signal f (t), and max | y (f) | is a maximum absolute value of an amplitude value in the frequency spectrum of the ultrasonic echo signal f (t);

step B3, superposition judgment: respectively performing superposition judgment on two adjacent extreme points in the M 'extreme points in the step B2 from front to back by adopting data processing equipment according to the sampling time sequence, and determining the total number of superposed extreme point pairs in the M' extreme points according to a superposition judgment result;

the superposition judgment methods of two adjacent extreme points in the M' extreme points are the same; when the M ' th extreme point and the M ' +1 th extreme point in the M ' extreme points are subjected to superposition judgment, adopting data processing equipment and calling a wave crest judgment module or a wave trough judgment module to perform superposition judgment; wherein M 'is a positive integer and M' is 1, 2, …, M '-2, M' -1;

when data processing equipment is adopted and a wave crest judging module is called for superposition judgment, when the m 'th extreme point and the m' +1 th extreme point are both wave crests, the m 'th extreme point and the m' +1 th extreme point are judged to be one superposition extreme point pair;

when data processing equipment is adopted and a wave trough judging module is called to carry out superposition judgment, when the m 'th extreme point and the m' +1 th extreme point are both wave troughs, the m 'th extreme point and the m' +1 th extreme point are judged to be one superposition extreme point pair;

after all adjacent two extreme points in the M ' extreme points are subjected to superposition judgment, obtaining all superposed extreme point pairs in the M ' extreme points, obtaining the total number of the superposed extreme point pairs in the M ' extreme points, and marking the total number of the superposed extreme point pairs in the M ' extreme points as K '; and determining the number K of the superimposed signals in the ultrasonic echo signal f (t) according to the total number K ' of the superimposed extremum point pairs in the M ' extremum points and according to a formula K ═ K ' + 1.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: before the superposition judgment in the step B3, initializing the total number K ' of the superposed extremum point pairs in the M ' extremum points by using a data processing device, where K ' is 0;

when the M 'th extreme point and the M' +1 th extreme point in the M 'extreme points are subjected to superposition judgment, when the M' th extreme point and the M '+1 th extreme point are both peaks or both valleys, judging that the M' th extreme point and the M '+1 th extreme point are one superposition extreme point pair, and adding 1 to K'; otherwise, K' at the moment is kept unchanged;

and after all adjacent two extreme points in the M ' extreme points are subjected to superposition judgment, K ' at the moment is the total number of the superposed extreme point pairs in the M ' extreme points.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: when signal sparse decomposition is carried out in the second step, respectively searching K optimal matching atoms in the iterative decomposition optimal atom set from beginning to end by adopting data processing equipment;

said

To find the best matching atom with the data processing apparatus,

time-frequency parameter r of_zIs the best time frequency parameter;

adopt data processing equipment and pair

When searching, pair

Time-frequency parameter r of_zThe search is performed as follows:

c1, optimizing time frequency parameters: calling an optimization algorithm module, finding out the optimal time frequency parameter with the maximum fitness value according to the preset value ranges of the scale parameter, the displacement parameter, the frequency parameter and the phase parameter in the time frequency parameters, and finding out the fitness value greater than thrd_zA plurality of preferred time-frequency parameters;

in this step, the number of the preferred time-frequency parameters found out is denoted as k_z，k_zIs an integer and k_z≥0；

Wherein the content of the first and second substances,

c' is a constant and is more than or equal to 2 and less than or equal to 3; r_z-1Residual amounts, R, for the first z-1 of said best matching atoms that have been found by the data processing apparatus at this time_z-1＝F(t)-ψ_z-1(t)，ψ_z-1(t) for this purpose use is made of the sum of the first z-1 of said best matching atoms, | | R, which has been determined by the data processing device_z-1||₂Is R_z-12-norm of (d); when z is 1, R₀＝F(t)；

The optimal time frequency parameter is the time frequency parameter with the maximum found fitness value, and the better time frequency parameter is the found fitness value larger than thrd_zThe fitness value of each time-frequency parameter is the Gabor atom and R corresponding to the time-frequency parameter_z-1Inner product of (d);

step C2, determining the optimal time frequency parameter: firstly, the number k of the better time-frequency parameters in the step C1_zAnd (4) judging: when k is_z0 or k_zWhen the time frequency is 1, the optimal time frequency parameter found in the step C1 is used as the optimal time frequency parameter r_z(ii) a When k is_zWhen k is not less than 2, k found in the step C1_zCalculating the sparsity of the better time-frequency parameters respectively, and taking the better time-frequency parameter with the minimum sparsity as the optimal time-frequency parameter r_z；

To k is paired_zThe rarity of the d-th better time-frequency parameter in the better time-frequency parametersSparsity XS_dWhen doing the calculation, XS_d＝||R^d||_ξ，R^dIs psi_dResidual amount of (t), R^d＝F(t)-ψ_d(t)，ψ_d(t) is k_zA Gabor atom corresponding to the d-th better time frequency parameter in the better time frequency parameters; r | |^d||_ξRepresents R^dXi-norm of which xi is a constant and is more than or equal to 0 and less than or equal to 1; d is a positive integer and 1, 2, …, k_z。

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: before echo signal separation in the third step, updating judgment of the optimal atom set of iterative decomposition is needed;

when the best atom set updating judgment of iterative decomposition is actually carried out, judging R^K(t)||₂Whether less than epsilon: when | | | R^K(t)||₂If the number is less than epsilon, judging that the optimal atom set of iterative decomposition in the step two does not need to be updated, and entering a step three; otherwise, when R^K(t)||₂When the number of the atoms is more than or equal to epsilon, judging that the optimal atom set of the iterative decomposition in the step two needs to be updated, and entering a step D;

wherein, | | R^K(t)||₂Is the R in the step two^K(t) 2-norm, epsilon is a preset residual quantity judgment threshold;

step D, updating the iterative decomposition optimal atom set: updating the iterative decomposition optimal atomic set in the second step by adopting data processing equipment, wherein the method comprises the following steps:

step D1, optimizing the optimal matching atoms, wherein the process is as follows:

step D11, atom random selection: randomly taking out an optimal matching atom from the iterative decomposition optimal atom set at the moment by adopting data processing equipment as an atom to be optimized, wherein the atom to be optimized is marked as

Wherein j1 is a positive integer and j1 is not less than 1 and not more than K;

k-1 best matching atoms except the atom to be optimized in the iterative decomposition optimal atom set are all atoms to be processed, and K-1 atoms to be processed form the atom set to be processed;

step D12, finding the best matching atom: the best matching atom found is recorded as

Is recorded as a time-frequency parameter r_j1', time-frequency parameter r_j1'＝(s_j1',u_j1',v_j1',w_j1')；

For the best matching atom

When searching, according to the preset s_j1'、u_j1'、v_j1' and w_j1' the value range is obtained by adopting data processing equipment and calling the optimizing algorithm module to find out the fitness value fitness (r)_j1') the largest optimal time-frequency parameter, the found optimal time-frequency parameter being the time-frequency parameter r_j1'; according to the formula

Solving for the best matching atom

Wherein the content of the first and second substances,

to represent

And

inner product of (d);

ψ'₀(t) is the sum of K-1 of said atoms to be treated in step D11;

d13, atom replacement judgment, atom replacement and iterative decomposition optimal atom set updating: adopting data processing equipment and calling a residual value judging module, a fitness value judging module or a sparsity judging module to judge whether the atoms to be optimized in the step D11 need to be replaced or not, and replacing the atoms to be optimized according to the judgment result;

adopting data processing equipment and calling a residual value judging module to judge whether the atoms to be optimized in the step D11 need to be replaced or not, and according to the residual value R after replacement_j1'^K(t)||_ξWhether it is less than the residue before replacement | | R_j1 ^K(t)||_ξAnd (4) judging: when | | | R_j1'^K(t)||_ξ＜||R_j1 ^K(t)||_ξIf so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

Obtaining the updated iterative decomposition optimal atom set; otherwise, judging that the atoms to be optimized in the step D11 do not need to be replaced, and entering the step D15;

wherein R is_j1'^K(t)＝F(t)-ψ_j1'(t)，

R_j1 ^K(t)＝F(t)-ψ_j1(t)，ψ_j1(t) is the sum of the K best matching atoms in the iterative decomposition best atom set before atom replacement judgment in this step; r | |_j1'^K(t)||_ξRepresents R_j1'^KXi-norm of (t, | | R_j1 ^K(t)||_ξRepresents R_j1 ^KXi-norm of (t), xi is constant and xi is more than or equal to 0 and less than or equal to 1;

d11, when the data processing equipment is adopted and the Fitness value judging module is called to judge whether the atoms to be optimized in the step D11 need to be replaced, the Fitness value Fitness (r) after replacement is carried out_j1') whether greater than pre-replacement fitnessValue Fitness (r)_j1) And (4) judging: when Fitness (r)_j1')＞Fitness(r_j1) If so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

wherein the content of the first and second substances,

represents R^j1-1(t) and

inner product of (A), R^j1-1(t)＝F(t)-ψ_j1-1(t)，ψ_j1-1(t) is the sum of the first j1-1 best matching atoms in the set of best atoms for this time of the iterative decomposition;

represents R^j1-1(t) and

inner product of (d);

adopting data processing equipment and calling a sparsity judging module to judge whether to replace the atoms to be optimized in the step D11 according to | | R^j1'||_ξWhether or not less than R^j1||_ξAnd (4) judging: when | | | R^j1'||_ξ＜||R^j1||_ξIf so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

Obtaining the updated iterative decomposition optimal atom set; otherwise, it is determined that the wait in step D11 is not neededOptimizing atoms for replacement and entering step D15;

wherein R is^j1' is

Amount of residual error of

R^j1Is composed of

Amount of residual error of

||R^j1'||_ξRepresents R^j1The xi-norm, | | R of `^j1||_ξRepresents R^j1ξ -norm of;

in this step, after completing the atom replacement judgment, the atom replacement, and the updating of the optimal atom set of iterative decomposition, the optimization process of the optimal matching atom selected in step D11 is completed;

step D14, residual error quantity judgment: and D13, judging the optimized residual quantity of the best matching atoms: when | | R'_j1 ^K(t)||₂If the number is less than epsilon, finishing the updating process of the iterative decomposition optimal atom set, and entering the step three; otherwise, when | | R'_j1 ^K(t)||₂When the value is more than or equal to epsilon, entering the step D15;

wherein, | R'_j1 ^K(t)||₂Is R'_j1 ^K(t) 2-norm; r'_j1 ^K(t) is the residual quantity after K iterative decompositions are performed on F (t) according to K best matching atoms in the iterative decomposition best atom set at the moment;

step D15, next best matching atom optimization: optimizing one of said best matching atoms in said iteratively decomposed set of best atoms that was not optimized at that time, according to the method described in steps D11 through D13;

step D16, residual error quantity judgment: optimizing the best matching atom in the step D15The residual error quantity is judged: when | | | R "_j1 ^K(t)||₂If the number is less than epsilon, finishing the updating process of the iterative decomposition optimal atom set, and entering the step three; otherwise, when | | R "_j1 ^K(t)||₂When the value is more than or equal to epsilon, returning to the step D15;

wherein, | | R "_j1 ^K(t)||₂Is R'_j1 ^K(t) 2-norm; r'_j1 ^K(t) is the residual quantity after K iterative decompositions are performed on F (t) according to K best matching atom pairs in the iterative decomposition best atom set at the moment.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: after signal sparse decomposition is carried out in the step D1, synchronously storing the iterative decomposition optimal atom set into a data memory by adopting data processing equipment, wherein the data memory is connected with the data processing equipment;

after atom replacement judgment and atom replacement are carried out in the step D13, the updated iterative decomposition optimal atom set is synchronously stored by adopting data processing equipment;

after signal sparse decomposition is performed in step D1, when the optimal atom set for iterative decomposition is synchronously stored in a data memory by using data processing equipment, storing K optimal matching atoms in the optimal atom set for iterative decomposition according to an iterative decomposition order; wherein the content of the first and second substances,

the best matching atom found when the z-th iterative decomposition is performed on the ultrasonic echo signal f (t) in step D1.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: when optimizing the best matching atoms in the step D1, optimizing the best matching atoms in the iterative decomposition best atom set according to the storage sequence;

when the best matching atom in the step D1 is optimized, the best matching atom which is optimized first is the 1 st best matching atom in the iterative decomposition best atom set in the step D1.

The echo signal separation method for ultrasonic flaw detection is characterized by comprising the following steps: after the echo signals are separated in the third step, respectively extracting the K separated signals; the signal extraction methods of the K separation signals are the same;

when any one of the K separation signals is subjected to signal extraction, signal extraction is carried out on the separation signal by adopting data processing equipment, and the separation signal is a signal to be processed and is marked as a signal f (t);

when the data processing equipment is adopted to extract the signals f (t), the process is as follows:

step A1, sparse decomposition based on an optimization algorithm: adopting data processing equipment and calling a sparse decomposition module to carry out iterative decomposition processing on the signal f (t), and converting the signal f (t) into the signal f (t)

And obtaining the optimal atom set of iterative decomposition at the moment; the iterative decomposition of the best set of atoms at this time contains m best matching atoms,

decomposing the nth best matching atom in the best atom set for the iteration;

in the formula R^m(t) is the residual error of the signal f (t) after m iterative decompositions, where m is the preset total iterative decompositions times and m is a positive integer, n is a positive integer and n is 1, 2, …, m; a is_nThe expansion coefficient of the best matching atom after the nth iterative decomposition and the residual error after the last iterative decomposition is obtained;

adopting data processing equipment and calling an optimization algorithm module to find out the best matching atom for the nth iterative decomposition;

is a Gabor atom and

wherein the function ψ (t) is a Gaussian window function and

r_nis composed of

Of time-frequency parameters r_n＝(s_n,u_n,v_n,w_n)，s_nAs a scale parameter, u_nAs a displacement parameter, v_nAs a frequency parameter, w_nIs a phase parameter;

in this step, the best matching atom is found

According to a preset s_n、u_n、v_nAnd w_nThe value range of the adaptive value Fitness (r) is found out by adopting data processing equipment and calling an optimization algorithm module_n) The maximum optimal time frequency parameter, the found optimal time frequency parameter is the time frequency parameter r_n；

Wherein, Fitness (r)_n) Is a time-frequency parameter r_nThe value of the fitness value of (a) is,

represents Rⁿ-¹(t) and

inner product of (d); r^n-1(t) is the residual quantity of the signal f (t) after n-1 iterative decompositions, R⁰(t)＝f(t)；

Step A2, residual error amount judgment: judgment | | | R^m(t)||₂Whether less than epsilon: when | | | R^m(t)||₂If the value is less than epsilon, the step A4 is entered; otherwise, when R^m(t)||₂When the value is more than or equal to epsilon, entering the step A3;

wherein, | | R^m(t)||₂Is described in step A1R^m(t) 2-norm, epsilon is a preset residual quantity judgment threshold;

step A3, optimizing the best matching atom, the process is as follows:

step A31, atom random selection: randomly taking out an optimal matching atom from the iterative decomposition optimal atom set at the moment by adopting data processing equipment as an atom to be optimized, wherein the atom to be optimized is marked as

Wherein j is a positive integer and j is more than or equal to 1 and less than or equal to m;

at the moment, m-1 best matching atoms except the atom to be optimized in the iterative decomposition optimal atom set are all atoms to be processed, and m-1 atoms to be processed form the atom set to be processed at the moment;

step A32, finding the best matching atom: the best matching atom found is recorded as

Is recorded as a time-frequency parameter r_j', time-frequency parameter r_j'＝(s_j',u_j',v_j',w_j')；

For the best matching atom

When searching, according to the preset s_j'、u_j'、v_j' and w_j' the value range is obtained by adopting data processing equipment and calling the optimizing algorithm module to find out the fitness value fitness (r)_j') the largest optimal time-frequency parameter, the found optimal time-frequency parameter being the time-frequency parameter r_j'; according to the formula

Solving for the best matching atom

Wherein the content of the first and second substances,

to represent

And

inner product of (d);

ψ₀(t) is the sum of m-1 of said atoms to be treated in step A31;

step A33, atom replacement judgment and atom replacement: adopting data processing equipment and calling a residual value judging module, a fitness value judging module or a sparsity judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and replacing the atoms to be optimized according to the judgment result;

when the data processing equipment is adopted and the residual value judging module is called to judge whether the atoms to be optimized in the step A31 need to be replaced, the residual value R after the replacement is carried out_j'^m(t)||_ξWhether it is less than the residue before replacement | | R_j ^m(t)||_ξAnd (4) judging: when | | | R_j'^m(t)||_ξ＜||R_j ^m(t)||_ξWhen the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

Obtaining the updated iterative decomposition optimal atom set; otherwise, judging that the atoms to be optimized in the step A31 do not need to be replaced, and entering the step A35;

wherein R is_j'^m(t)＝f(t)-ψ_j'(t)，

R_j ^m(t)＝f(t)-ψ_j(t)，ψ_j(t) is the sum of m best matching atoms in the iterative decomposition best atom set before atom replacement judgment in the step; r | |_j'^m(t)||_ξRepresents R_j'^mXi-norm of (t, | | R_j ^m(t)||_ξRepresents R_j ^mXi-norm of (t), xi is constant and xi is more than or equal to 0 and less than or equal to 1;

adopting data processing equipment and calling a Fitness value judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and according to the replaced Fitness value Fitness (r)_j') whether it is greater than the pre-replacement Fitness value Fitness (r)_j) And (4) judging: when Fitness (r)_j')＞Fitness(r_j) When the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein the content of the first and second substances,

represents R^j-1(t) and

inner product of (A), R^j-1(t)＝f(t)-ψ_j-1(t)，ψ_j-1(t) is the sum of the first j-1 best matching atoms in the set of best atoms for this time of the iterative decomposition;

represents R^j-1(t) and

inner product of (d);

adopting data processing equipment and calling a sparsity judging module to judge whether the sparsity judging module needs to be used for judging whether the step A31 is neededWhen the atoms to be optimized are replaced, according to | | | R^j'||_ξWhether or not less than R^j||_ξAnd (4) judging: when | | | R^j'||_ξ＜||R^j||_ξWhen the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein R is^j' is

Amount of residual error of

R^jIs composed of

Amount of residual error of

R^j'||_ξRepresents R^jThe xi-norm, | | R of `^j||_ξRepresents R^jξ -norm of;

in this step, after the atom replacement judgment and the atom replacement are completed, the optimization process of the best matching atom selected in step a31 is completed;

step A34, residual error amount judgment: and B, judging the optimized residual quantity of the best matching atom in the step A33: when | | R'_j ^m(t)||₂If the value is less than epsilon, the step A4 is entered; otherwise, when | | R'_j ^m(t)||₂When the value is more than or equal to epsilon, entering the step A35;

wherein, | R'_j ^m(t)||₂Is R'_j ^m(t) 2-norm; r'_j ^m(t) is according toThe residual quantity of m best matching atom pairs f (t) in the best atom set of iterative decomposition after m iterative decomposition is carried out at the moment;

step A35, next best matching atom optimization: optimizing one of said best matching atoms in said iteratively decomposed set of best atoms that was not optimized at that time, according to the method described in steps a31 through a 33;

step A36, residual error amount judgment: and B, judging the optimized residual quantity of the best matching atom in the step A35: when | | | R "_j ^m(t)||₂If the value is less than epsilon, the step A4 is entered; otherwise, when | | R "_j ^m(t)||₂When the value is more than or equal to epsilon, returning to the step A35;

wherein, | | R "_j ^m(t)||₂Is R'_j ^m(t) 2-norm; r'_j ^m(t) is the residual quantity after m iterative decompositions are performed on f (t) according to m best matching atoms in the iterative decomposition best atom set at the moment;

step A4, signal reconstruction: obtaining an approximate signal f' (t) of a signal f (t) by adopting data processing equipment according to the iterative decomposition optimal atom set at the moment; wherein the approximation signal f "(t) is a signal extracted from the signal f (t),

wherein

To this end the iteration decomposes the nth 'best matching atom in the best set of atoms, n' being a positive integer and n ═ 1, 2, …, m; a is_n' is

And f (t) is subjected to n '-1 times of iterative decomposition according to the first n' -1 best matching atoms in the iterative decomposition best atom set at the moment, and then the residual quantity is expanded.

Compared with the prior art, the invention has the following advantages:

1. the method has the advantages of simple steps, reasonable design, convenient implementation and low input cost.

2. The signal separation speed is high, the data processor is adopted to automatically complete the signal separation process, and the signal separation process can be completed in several minutes or even shorter time, so that the real-time signal separation is realized.

3. The signal sparse decomposition method based on the optimization algorithm is simple, reasonable in design, convenient to implement and good in using effect, improves the signal extraction speed, can effectively improve the quality and performance indexes of original signals after signal extraction, and particularly plays an important role in ultrasonic nondestructive inspection. Meanwhile, the value range of the frequency parameter v is limited to

The unit of fo is MHz, so that on one hand, the calculated amount of a sparse decomposition algorithm can be effectively reduced, and real-time detection is realized; on the other hand, the performance of the MP algorithm (namely, the matching pursuit algorithm) is effectively improved, so that the sparsely represented signal can effectively meet the detection precision of weak and small defects, and the purpose of more simply and accurately obtaining effective information contained in the signal is achieved. By limiting the range of values of the frequency parameter v to

Effective information contained in the signals can be further highlighted, the sparsely represented signals can express the effective information more heavily, redundant information is weakened, signal intrinsic characteristics can be expressed more accurately, and the precision of signal separation can be effectively guaranteed.

4. When the optimal matching atom is searched based on the optimization algorithm, besides the index of the fitness value, the signal sparsity is also considered so as to select the atom optimally matched with the signal F (t), so that the signal (or the extracted signal) after drying can more accurately express the intrinsic characteristics of the signal, and the signal extraction precision can be effectively ensured.

5. The method is simple and easy to implement, after sparse decomposition of the signals, whether the optimal atom set of iterative decomposition at the moment meets the preset signal extraction precision requirement or not needs to be judged through residual quantity judgment, and optimal matching atom optimization is carried out according to the judgment result, so that the accuracy of signal extraction can be further improved, the signal extraction precision is further improved, the extracted signals further approach the original signals, the optimal matching with the original signals is realized, the extraction accuracy and the extraction speed of the signals are improved, and the signal separation effect is ensured.

6. The adopted optimal matching atom optimization method is reasonable in design, convenient to implement and good in using effect, one optimal matching atom is randomly selected from the iterative decomposition optimal atom set for optimization, after the optimization is completed, whether the iterative decomposition optimal atom set meets the signal extraction precision requirement or not is judged through the residual quantity, and whether the optimization needs to be continuously carried out on the rest optimal matching atoms or not is determined according to the judgment result. Therefore, the method is simple and convenient to realize, can realize the combination of quick optimization and real-time judgment of an optimization result, can effectively simplify the optimization process of the optimal matching atom, can quickly achieve the aim of optimizing the optimal matching atom, and further effectively improves the signal extraction precision. Meanwhile, the adopted atom replacement judgment method is reasonable in design, simple and convenient to realize and good in use effect, any one of the methods of residual value judgment, fitness value judgment or sparsity judgment is adopted for atom replacement judgment, any one of the methods can be selected for atom replacement judgment, the use mode is flexible, and each atom replacement judgment method can realize effective atom replacement judgment.

7. The improved artificial bee colony algorithm is adopted for optimizing to realize the best matching atom search, all atoms in an over-complete dictionary do not need to be generated before signal sparse decomposition, only the positions of honey sources need to be generated to replace Gabor atoms in an atom library, and the storage space is greatly saved. In addition, the artificial bee colony algorithm searches the optimal matching atoms in a continuous space, and the matching tracking algorithm searches the atoms in a discrete search space, so that the search range of the artificial bee colony algorithm is wider, the extracted atoms can better reflect the characteristics of original signals, the calculation speed is improved, and the accuracy of parameter extraction is improved because the atoms are optimized in the continuous solution space range. Compared with a discrete space range, the method can more accurately extract the optimal matching atoms from signal matching, thereby improving the precision of signal extraction and effectively extracting useful signals under a strong noise background.

8. The method for determining the number K of the superposed signals in the ultrasonic echo signals F (t) is reasonable in design, simple and convenient to implement and good in using effect, the traveling wave peak and the wave trough are determined firstly, and all extreme points in the ultrasonic echo signals F (t) are obtained correspondingly; and then, effective extreme point elimination is carried out according to a specific extreme point elimination method, so that the calculation amount of the superposition judgment in the step B3 is effectively reduced, the superposed signals can be more highlighted, the total number K ' of the superposed extreme point pairs in the M ' extreme points after the extreme point elimination is judged and identified through superposition in the step B3, and the number K of the superposed signals in the ultrasonic echo signal F (t) can be determined according to a formula K ═ K ' + 1. Only if the effective extreme point elimination in the step B2 is effectively combined with the superposition judgment in the step B3, the superposition extreme point pairs in the M' extreme points and the superposition positions of the two superposed signals are in one-to-one correspondence, so that the determination of the number K of the superposed signals in the ultrasonic echo signal f (t) is realized, and the number K of the superposed signals in the ultrasonic echo signal f (t) is neither too large nor too small. After the effective extreme points are removed in the step B2, the extreme points of a superimposed signal that have no practical analysis significance can be removed, and only the extreme points with real value in the superimposed signal are retained, so that the condition that the number K of the superimposed signals is determined to be too large can be effectively avoided; meanwhile, in the step B3, by performing the overlapping judgment on two adjacent extreme points in the M' extreme points respectively, any overlapping position is not missed, so that the situation that the number K of the determined overlapping signals is too small can be effectively avoided, and the accuracy of the number K of the determined overlapping signals is very high.

9. After signal separation, each separated signal needs to be extracted, and the adopted signal extraction method has simple steps, reasonable design, convenient realization and lower investment costThe method has the advantages of low signal extraction speed, automatic signal extraction process completion by adopting the data processor, signal extraction process completion in several minutes or even shorter time, and signal real-time extraction. When the signal is actually extracted, the signal sparse decomposition method based on the optimization algorithm is simple, reasonable in design, convenient to implement and good in using effect, the signal extraction speed is improved, the quality and performance indexes of original signals can be effectively improved after the signal is extracted, and the method has an important effect particularly in ultrasonic nondestructive flaw detection. Meanwhile, the value range of the frequency parameter v is limited to

And f_oThe unit of (A) is MHz, on one hand, the calculated amount of a sparse decomposition algorithm can be effectively reduced, and real-time detection is realized; on the other hand, the performance of the MP algorithm (namely, the matching pursuit algorithm) is effectively improved, so that the sparsely represented signal can effectively meet the detection precision of weak and small defects, and the purpose of more simply and accurately obtaining effective information contained in the signal is achieved. By limiting the range of values of the frequency parameter v to

Effective information contained in the signals can be further highlighted, the sparsely represented signals can express the effective information more heavily, redundant information is weakened, signal intrinsic characteristics can be expressed more accurately, and signal extraction precision can be effectively guaranteed.

Meanwhile, a signal extraction precision verification step is added, the verification method is simple and easy to implement, after signal sparse decomposition, whether the current iteration decomposition optimal atomic set meets the preset signal extraction precision requirement or not needs to be judged through residual quantity judgment, and optimal matching atomic optimization is carried out according to the judgment result, so that the accuracy of signal extraction can be further improved, the signal extraction precision is further improved, the extracted signal further approaches to the original signal, the optimal matching with the original signal is realized, and the signal extraction accuracy and the signal extraction speed are improved.

The method can adopt an improved artificial bee colony algorithm to carry out optimization to realize extraction of the separated signals, particularly provides a reliable basis for accurate detection of defects for ultrasonic nondestructive inspection, simultaneously provides a theoretical basis for qualitative and quantitative analysis of the defects, can effectively solve the problem of difficult extraction of weak defects under a strong noise background, can solve the problems of extraction speed and precision of the weak defects, can accurately extract defect information under the strong noise background, improves the speed of ultrasonic signal extraction, and provides technical support for real-time automatic detection, thereby effectively solving the problems of high algorithm complexity, over-matching and the like of the existing matching tracking algorithm. Therefore, the invention utilizes the artificial bee colony algorithm to select the atoms which are optimally matched with the ultrasonic signals from the continuous dictionary library, thereby recovering the signals to be processed.

10. The method has the advantages of good use effect, high practical value, simple steps, reasonable design, convenient implementation and good use effect, and realizes simple, convenient and quick separation of superposed signals at positions of a plurality of defects in the same ultrasonic echo signal by utilizing signal sparse decomposition.

The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.

Drawings

FIG. 1 is a block diagram of the process flow of the present invention.

Fig. 2 is a schematic block diagram of a signal separation processing system according to the present invention.

Description of reference numerals:

1-ultrasonic flaw detection device; 2-a data processing device; 3-data memory.

Detailed Description

An echo signal separation method for ultrasonic testing shown in fig. 1 includes the following steps:

step one, ultrasonic echo signal acquisition and synchronous uploading: ultrasonic detection is carried out on the object to be detected by adopting an ultrasonic flaw detection device 1, an ultrasonic echo signal F (t) of the object to be detected is obtained, and the obtained ultrasonic echo signal F (t) is synchronously transmitted to a data processing device 2;

wherein f (t) ═ f (t)₁),f(t₂),...,f(t_N')]^TT represents a time parameter, t_iThe ith sampling time, f (t) of the ultrasonic flaw detector 1_i) A signal value sampled at the ith sampling time of the ultrasonic flaw detection apparatus 1, where i is a positive integer, i is 1, 2, 3, …, N' is a positive integer, and is a signal length of the ultrasonic echo signal f (t);

wherein, the [ alpha ], [ beta ]]^TRepresents a transpose of a matrix; [ f (t)₁),f(t₂),...,f(t_N')]^TIs a matrix [ f (t)₁),f(t₂),...,f(t_N')]Transposing;

step two, signal sparse decomposition: adopting a data processing device 2 and calling a sparse decomposition module to carry out iterative decomposition processing on the ultrasonic echo signal F (t) in the step one, and converting the ultrasonic echo signal F (t) into the ultrasonic echo signal F (t)

for the best matching atom found at the z-th iterative decomposition,

is a Gabor atom and

wherein the function ψ (t) is a Gaussian window function and

r_zis composed of

step three, echo signal separation: separating the ultrasonic echo signals F (t) in the first step by adopting data processing equipment 2 according to the K optimal matching atoms in the iterative decomposition optimal atom set at the moment to obtain K separation signals;

the K' th of the K separation signals is denoted as F_k'(t)，

Wherein the content of the first and second substances,

represents R^k'^-1(t) and

The sparse decomposition algorithm (also called MP algorithm) is an approximation of atoms by multiple iterations, and each iteration adds a single atom to match with the ultrasound echo signal, so that the inner product of the selected atom and the residual of the signal is maximized, and the energy of the residual is minimized. Thus, the MP algorithm recovers the signal by matching the signal features with the minimum subset of atoms, which is not only computationally intensive, but more importantly, lacks of physical interpretation. If an accurate physical meaning of the selected atoms can be given, the atoms correspond to the ultrasonic echo signals one by one, dominant atoms are sequentially selected from an over-complete atom library and added into a signal approximation representation to recover a plurality of defect signals, and the selection of the atoms with physical explanation enables the extraction of the plurality of defect signals to have a reliable basis.

In the invention, the K optimal matching atoms in the iterative decomposition optimal atom set are respectively in one-to-one correspondence with K separation signals (namely, the ultrasonic echo signals at the positions of K defects in the measured object, also called as K defect signals) contained in the ultrasonic echo signals F (t), and the separation of the K defect signals is realized through signal sparse decomposition so as to respectively and accurately identify the K defect signals.

And in the fourth step, each separation signal is an ultrasonic echo signal at the position of one defect in the measured object.

And step one, overlapping the ultrasonic echo signals at the positions of the K defects in the measured object.

In the present embodiment, in the first step, the ultrasonic testing apparatus 1 and the data processing device 2 are connected by a wired or wireless method.

Step two described

In practical use, when the number K of the superimposed signals included in the ultrasonic echo signal f (t) is known, the method described in the step one to the step three can simply, conveniently and quickly complete the separation process of the ultrasonic echo signals (i.e., K superimposed signals) at the positions of K defects in the ultrasonic echo signal f (t); however, when the number K of the superimposed signals included in the ultrasound echo signal f (t) is unknown, the number K of the superimposed signals in the ultrasound echo signal f (t) needs to be determined.

In this embodiment, before performing signal sparse decomposition in the second step, the number K of the superimposed signals in the ultrasonic echo signal f (t) is determined;

when the number K of the superimposed signals in the ultrasonic echo signals f (t) is determined, the data processing device 2 is used for determining, and the process is as follows:

step B1, determining wave crests and wave troughs: respectively determining all wave crests and all wave troughs of the ultrasonic echo signal F (t) in the step one by adopting a data processing device 2, and respectively synchronously recording the sampling time and the signal value of each determined wave crest and each determined wave trough;

step B2, removing extreme points: adopting data processing equipment 2 and calling a time domain extreme point eliminating module or a frequency domain extreme point eliminating module to eliminate extreme points to obtain M 'eliminated extreme points, and arranging the M' extreme points from front to back according to the sampling time sequence of the extreme points; wherein M' is a positive integer and is the total number of extreme points obtained after the extreme points are removed in the step;

when the data processing device 2 is adopted and the time domain extreme point eliminating module is called to eliminate the extreme points, eliminating the extreme points of which the absolute values of the signal values in all the extreme points determined in the step B1 are smaller than beta 'to obtain M' number of eliminated extreme points; wherein, β' is a constant and the value range thereof is 0.1-0.35, and max | f (t) | is the maximum absolute value of the signal value in the ultrasonic echo signal f (t);

when the data processing device 2 is adopted and the frequency domain extreme point eliminating module is called to eliminate the extreme points, eliminating the extreme points of which the absolute values of the signal values are smaller than beta in all the extreme points determined in the step B1 to obtain M' eliminated extreme points; wherein β is a preset rejection threshold, β ═ α × max | y (f) |, α is a constant and has a value range of 0.25 to 0.35, y (f) is a frequency spectrum of the ultrasonic echo signal f (t), and max | y (f) | is a maximum absolute value of an amplitude value in the frequency spectrum of the ultrasonic echo signal f (t);

step B3, superposition judgment: respectively performing superposition judgment on two adjacent extreme points in the M 'extreme points in the step B2 from front to back by adopting the data processing equipment 2 according to the sampling time sequence, and determining the total number of superposed extreme point pairs in the M' extreme points according to a superposition judgment result;

the superposition judgment methods of two adjacent extreme points in the M' extreme points are the same; when the M ' th extreme point and the M ' +1 th extreme point in the M ' extreme points are subjected to superposition judgment, the data processing equipment 2 is adopted and a wave crest judgment module or a wave trough judgment module is called for carrying out superposition judgment; wherein M 'is a positive integer and M' is 1, 2, …, M '-2, M' -1;

when the data processing device 2 is adopted and a wave crest judging module is called to carry out superposition judgment, when the m 'th extreme point and the m' +1 th extreme point are both wave crests, the m 'th extreme point and the m' +1 th extreme point are judged to be one superposition extreme point pair;

when the data processing device 2 is adopted and a wave trough judging module is called to carry out superposition judgment, when the m 'th extreme point and the m' +1 th extreme point are both wave troughs, the m 'th extreme point and the m' +1 th extreme point are judged to be one superposition extreme point pair;

In this embodiment, in step B3, the data processing device 2 is adopted and a peak determining module is invoked to perform superposition determination, and when the m 'th extremum point and the m' +1 th extremum point are both peaks, it is determined that the m 'th extremum point and the m' +1 th extremum point are one superimposed extremum point pair.

In practical use, in step B3, the data processing device 2 may also be adopted and a valley judging module is invoked to perform superposition judgment, and when the m 'th extreme point and the m' +1 th extreme point are both valleys, it is judged that the m 'th extreme point and the m' +1 th extreme point are one superimposed extreme point pair.

In this embodiment, in step B2, the data processing device 2 is adopted and the time domain extreme point eliminating module is called to eliminate the extreme point, and the value size of α can be adjusted accordingly according to actual needs.

In practical use, the data processing device 2 may be adopted in step B2, and the frequency domain extreme point removing module is called to remove the extreme points, and the value of α' may be adjusted accordingly according to practical needs.

In the signal sampling process, due to the influence of various factors such as environmental factors and sampling systems, a plurality of non-real extreme points (namely interference extreme points) inevitably exist in the ultrasonic echo signal f (t), and the extreme points are not the extreme points of defect signals in the ultrasonic echo signal f (t), so that the non-real extreme points need to be removed, thus not only effectively reducing the calculated amount, but also removing the interference extreme points so as to accurately determine the number K of the existing superposed signals.

The two eliminating methods of time domain eliminating and frequency domain eliminating can simply, conveniently, quickly and effectively eliminate the interference extreme points, and the two eliminating methods can be automatically completed by adopting the data processing equipment 2, so the method is convenient to realize and high in processing speed, and the eliminating result can be effectively ensured.

When the time domain eliminating method or the frequency domain eliminating method is adopted to eliminate the extreme points, the adopted eliminating threshold value beta' and the eliminating threshold value beta are set reasonably. The determination of the rejection threshold β' and the rejection threshold β are both directly related to the ultrasonic echo signal f (t), and the determination is performed by the currently processed ultrasonic echo signal f (t) instead of a fixed value for different signals to be processed, so that certain adaptability is achieved, the rejection effect can be effectively ensured, and the real extreme points cannot be rejected too much. And in addition, the actual extreme point removing effect and the ultrasonic defect signal identification purpose are combined, and the actual using effect of the extreme point removing method is very good.

Wherein β ' ═ α ' × max | f (t) |, β ' is determined by the maximum absolute value max | f (t)) | of the signal values in the ultrasound echo signal f (t), the signal values of the superimposed signals included in the currently processed ultrasound echo signal f (t) can be known according to max | f (t)) |, and the extreme points whose absolute values are smaller than β ' are substantially not associated with the superimposed signals included in the ultrasound echo signal f (t), and the meaning of the actual analysis is substantially absent, so that the extreme points whose absolute values are smaller than β ' are eliminated.

Accordingly, β ═ α × max | y (f) |, β is determined by the maximum absolute value of the amplitude in the ultrasound echo signal f (t) max | y (f) |, and the amplitude of the superimposed signal included in the currently processed ultrasound echo signal f (t) can be known according to max | y (f) |, and the extreme points where the absolute value of the signal value is smaller than β are substantially unrelated to the superimposed signal included in the ultrasound echo signal f (t), and the meaning of the actual analysis does not exist, so the extreme points where the absolute value of the signal value is smaller than β are eliminated.

When the superposition judgment is performed in the step B3, the adopted superposition judgment method is reasonable in design, convenient to implement, and good in use effect, the total number K ' of the superposition extreme point pairs in the M ' extreme points is identified by respectively performing superposition judgment on the eliminated M ' extreme points, and the number K of the superposition signals in the ultrasonic echo signal f (t) can be directly determined according to the total number K ' of the superposition extreme point pairs in the M ' extreme points.

The determination method of the number K of the superimposed signals in the ultrasonic echo signals f (t) adopted in the steps B1 to B3 is reasonable in design, simple and convenient to implement, and good in use effect, and the traveling wave peak and the wave trough are determined first, and all extreme points in the ultrasonic echo signals f (t) are obtained correspondingly; and then, effective extreme point elimination is carried out according to a specific extreme point elimination method, so that the calculation amount of the superposition judgment in the step B3 is effectively reduced, the superposed signals can be more highlighted, the total number K ' of the superposed extreme point pairs in the M ' extreme points after the extreme point elimination is judged and identified through superposition in the step B3, and the number K of the superposed signals in the ultrasonic echo signal F (t) can be determined according to a formula K ═ K ' + 1. Only if the effective extreme point elimination in the step B2 is effectively combined with the superposition judgment in the step B3, the superposition extreme point pairs in the M' extreme points and the superposition positions of the two superposed signals are in one-to-one correspondence, so that the determination of the number K of the superposed signals in the ultrasonic echo signal f (t) is realized, and the number K of the superposed signals in the ultrasonic echo signal f (t) is neither too large nor too small. After the effective extreme points are removed in the step B2, the extreme points of a superimposed signal that have no practical analysis significance can be removed, and only the extreme points with real value in the superimposed signal are retained, so that the condition that the number K of the superimposed signals is determined to be too large can be effectively avoided; meanwhile, in the step B3, by performing the overlapping judgment on two adjacent extreme points in the M' extreme points respectively, any overlapping position is not missed, so that the situation that the number K of the determined overlapping signals is too small can be effectively avoided, and the accuracy of the number K of the determined overlapping signals is very high.

In this embodiment, before the overlap determination in step B3, a data processing device 2 is first used to initialize the total number K ' of the overlap extremum point pairs in the M ' extremum points, where K ' is 0;

In this embodiment, when performing signal sparse decomposition in step two, the data processing device 2 is adopted to search K best matching atoms in the iterative decomposition best atom set from first to last;

said

To use the best matching atoms found by the data processing apparatus 2,

time-frequency parameter r of_zIs the best time frequency parameter;

using a data processing device 2 and

when searching, pair

Time-frequency parameter r of_zThe search is performed as follows:

Wherein the content of the first and second substances,

c' is a constant and is more than or equal to 2 and less than or equal to 3; r_z-1 residual amounts, R, of the first z-1 of said best matching atoms that have been found for this time with the data processing device 2_z-1＝F(t)-ψ_z-1(t)，ψ_z-1(t) for this purpose the sum of the first z-1 of said best matching atoms, | | R, which has been determined by the data processing device 2 is used_z-1||₂Is R_z-12-norm of (d); when z is 1, R₀＝F(t)；

To k is paired_zSparsity XS of the d-th one of the better time-frequency parameters_dWhen doing the calculation, XS_d＝||R^d||_ξ，R^dIs psi_dResidual amount of (t), R^d＝F(t)-ψ_d(t)，ψ_d(t) is k_zA Gabor atom corresponding to the d-th better time frequency parameter in the better time frequency parameters; r | |^d||_ξRepresents R^dXi-norm of which xi is a constant and is more than or equal to 0 and less than or equal to 1; d is a positive integer and 1, 2, …, k_z。

The 2-norm is the sum of the squares of the elements of the vector and then the square root (i.e., the L2 norm) is taken, as is well known in the art.

R in step C1_z-1Is a vector of dimension N' x 1, | | R_z-1||₂Is R_z-1The power of 1/2 of the sum of the 2 nd power of the absolute values of the N' elements.

R in step C2^dIs a vector of dimension N' x 1, | | R^d||_ξIs R^dThe absolute values of the N' elements are 1/xi to the power of the xi power sum.

In this embodiment, the ultrasonic flaw detector 1 is an a-type digital ultrasonic flaw detector. In practice, other types of ultrasonic flaw detection equipment may be used.

Before signal sparse decomposition is carried out in the second step, the value ranges of the time-frequency parameter mesoscale parameter, the displacement parameter, the frequency parameter and the phase parameter are respectively determined according to the value range determination method of the time-frequency parameter mesoscale parameter, the displacement parameter, the frequency parameter and the phase parameter when the conventional signal sparse decomposition is carried out. Step two described

Is the best matching atom when the signal F (t) to be processed is subjected to the z-th iterative decomposition.

Each Gabor atom corresponds to its time-frequency parameter, and each Gabor atom corresponds to its time-frequency parameter one to one.

The article "implementing MP-based sparse signal decomposition using FFT" (author: infusory) published in the journal of electronics and information (vol.28, No. 4) at 4.2006 discloses: "…, r ═ s, u, v, w), the time-frequency parameters can be discretized as follows: r ═ alpha^j,pα^jΔu,kα^-jΔ v, i Δ w), where α ═ 2, Δ u ═ 1/2, Δ v ═ pi, Δ w ═ pi/6, 0 < j < log₂N,0≤p≤N2^-j+1，0≤k≤N2^j+1I is more than or equal to 0 and less than or equal to 12. The above description gives a specific over-complete atom library ". From the above, the frequency parameter v is based on k α^-jΔ v is discretized by k is more than or equal to 0 and less than or equal to N2^j+1、0＜j＜log₂N, α ═ 2, and Δ v ═ pi, in which case the frequency parameter v has a very large value range, the minimum value of the frequency parameter v being 0 and the maximum value thereof being

The frequency parameter v thus has a value range of

Even with discretization, the range of values of the frequency parameter v is still very large.

In this embodiment, the scale parameter in step C1 has a value range of [1, N']The value range of the displacement parameter is [0, N']The value range of the frequency parameter is

The value range of the phase parameter is [0,2 pi]. Wherein f is_oSampling frequency f of the ultrasonic flaw detector 1_oIn MHz.

According to the common knowledge in the field, the sparse decomposition algorithm (also called MP algorithm) has two defects, one is that the calculation amount of the sparse decomposition algorithm is large, the calculation time is huge under the existing calculation condition, and real-time detection cannot be carried out; secondly, the sparse decomposition algorithm is an optimal solution obtained under a continuous condition, and the detection precision of weak and small defects is still limited.

The purpose of signal sparse representation is to represent a signal by using as few atoms as possible in a given overcomplete dictionary, and a more concise representation mode of the signal can be obtained, so that information contained in the signal can be more easily obtained, and the signal can be more conveniently processed, such as compression, encoding and the like. The research focus of the signal sparse representation direction mainly focuses on the aspects of sparse decomposition algorithm, overcomplete atom dictionary (also called atom library, Gabor dictionary), application of sparse representation and the like. Two major tasks of signal sparse representation are dictionary generation and sparse decomposition of signals. However, existing research has demonstrated that searching for atoms in scale and frequency from a coarse scale to a fine scale can significantly improve the performance of MP algorithms (i.e., matching pursuit algorithms) without increasing the size of the atom pool. Thus, the range of the frequency parameter v

Further reduction of energyThe performance of the MP algorithm (namely the matching pursuit algorithm) is effectively improved. Especially for the frequency parameter, the value range has a larger influence on the performance of the MP algorithm (i.e. the matching pursuit algorithm).

Because the value range of the frequency parameter (i.e., the frequency parameter v) is related to the actual sampling frequency of the signal, on the basis of the research experience of years of sparse decomposition, the value range of the time-frequency parameter is used for fully and long-term research and verification on the influence of improving the performance of the MP algorithm (i.e., the matching pursuit algorithm), and then the value range of the frequency parameter and the actual sampling frequency of the processed signal (i.e., the sampling frequency f of the ultrasonic flaw detection device 1) are obtained_o) Closely related and not completely corresponding one to one, and the value range of the frequency parameter is limited to be in the comprehensive view of simplifying the calculated amount of the sparse decomposition algorithm and refining the value range of the time-frequency parameter and improving the performance of the matching tracking algorithm

And f_oThe unit of (A) is MHz, on one hand, the calculated amount of a sparse decomposition algorithm can be effectively reduced, and real-time detection is realized; on the other hand, the performance of the MP algorithm (namely, the matching pursuit algorithm) is effectively improved, so that the sparsely represented signal can effectively meet the detection precision of weak and small defects, and the purpose of more simply and accurately obtaining effective information contained in the signal is achieved. By limiting the value range of the frequency parameter to

According to the common general knowledge in the art, and combining the article "implementing MP-based signal sparse decomposition by FFT" (author: yiloy) published in the journal of electronics and information in 4.2006 (vol.28, No. 4), it is known that, before sparse decomposition is currently performed, four parameters of a time-frequency parameter are usually discretized respectively, and an overcomplete atom library is generated, but the number of atoms in the overcomplete atom library is usually very large, the occupied storage space is very large, the calculation amount is large, the calculation engineering is complex, and all atoms in the overcomplete atom library need to be analyzed and judged respectively, and the best matching atom is found; meanwhile, the parameter value range and the discretization method also have great influence on the generated overcomplete atom library, which inevitably causes poor accuracy of the generated overcomplete atom library (also called an overcomplete dictionary, Gabor dictionary), so that the intrinsic characteristics of the signals cannot be accurately expressed, and the signal extraction precision cannot be guaranteed.

Before signal sparse decomposition is carried out in the second step, all atoms in the over-complete dictionary do not need to be generated, and the data processing equipment 2 is only needed to be adopted and the optimization algorithm module is called to carry out optimization, so that the best matching atoms can be simply, conveniently and quickly found out one by one, and the storage space is greatly saved. In addition, the optimizing algorithm module searches for the best matching atoms in the value range of each parameter (specifically in a continuous space), and searches for the best matching atoms in a discrete search space (namely an over-complete dictionary or an over-complete atom library obtained through discretization) with the traditional matching tracking algorithm, so that the searching range of the optimizing algorithm module is wider, and the searched best matching atoms can better reflect the characteristics of an original signal (namely the ultrasonic echo signal F (t)), so that the accuracy of signal extraction can be further ensured, and the separation effect of K superposed signals can be correspondingly ensured.

The optimizing algorithm module in the step C1 is a genetic algorithm module, an artificial fish swarm algorithm module or an artificial bee swarm algorithm module. In practice, the optimization algorithm module may be other types of optimization algorithm modules. When a genetic algorithm module is called for optimizing, a conventional genetic algorithm is adopted; when an artificial fish school algorithm module is called for optimizing, a conventional artificial fish school algorithm is adopted; and when the genetic algorithm module and the artificial bee colony algorithm module are called for optimizing, the conventional artificial bee colony algorithm is adopted.

The method for determining the best matching atom by adopting the data processing equipment 2 and calling the optimizing algorithm module for optimizing has the following advantages: firstly, the defect that the traditional methods such as Fourier transform, wavelet transform and the like can only carry out decomposition on orthogonal basis is overcome, and the intrinsic characteristics of signals can be more accurately expressed, so that the accuracy of signal extraction is improved; secondly, the generation of local optimal values can be effectively avoided, optimization searching in a continuous space can be carried out, and compared with the optimization searching in a discrete space carried out by the original matching tracking algorithm, the searching range is expanded, so that the accuracy of signal extraction is further effectively improved; thirdly, the optimal matching atoms are found out through optimization of the optimization algorithm module, the method is simple and convenient to realize and high in extraction speed, the problem of high complexity of an original matching algorithm can be effectively solved, the convergence speed of noise reduction processing and the signal extraction speed are greatly improved, and the real-time performance of signal extraction is improved; fourthly, the accuracy of signal extraction can be effectively improved, and the problems of signal extraction under the background of strong noise and extraction of weak and small signals are solved; fifthly, the use effect is good, the detection problems of weak and small defects and the like in the field of ultrasonic nondestructive inspection can be solved, the product quality of production enterprises is improved, and potential safety hazards are avoided; sixth, application scope is wide, can effectively be applicable to the extraction process of multiple signal, especially can effectively extract to the sound signal of difficult inspection of nonstationary. Therefore, the method for determining the best matching atoms by calling the optimization algorithm module for optimization has the advantages of reasonable design, good effect and high practical value, not only improves the speed of signal extraction, but also can effectively improve the quality and performance indexes of original signals after signal extraction, and has an important effect particularly in ultrasonic nondestructive inspection.

In the actual use process, the best matching atom is determined whether sparse decomposition is carried out after an over-complete atom library is established in a conventional matching and tracking algorithm or when signal sparse decomposition is completed by using an optimization algorithm module to find the best matching atom, and the fitness value of the atom (namely the inner product of the Gabor atom and the corresponding residual quantity) is used as the unique matching judgment standard, so that the method has one-sidedness. Correspondingly, the determined best matching atoms are only Gabor atoms with the best fitness value, and when the signal is extracted, in order to select the atoms which are optimally matched with the signal F (t) to be processed, the signal sparsity needs to be considered besides the index of the fitness value, so that the intrinsic characteristics of the signal can be more accurately expressed, and the signal extraction precision can be effectively ensured.

In the second step of the invention, the data processing equipment 2 pairs are adopted

Time-frequency parameter r of_zWhen searching is carried out, firstly, the time frequency parameter is optimized according to the step C1, the optimal time frequency parameter with the maximum adaptability value is found out, and the adaptability value is more than thrd_zA plurality of preferred time-frequency parameters, thrd_zDetermining a threshold value for the fitness value; and determining the optimal time-frequency parameters by combining the sparsity of the optimal time-frequency parameters according to the step C2. Therefore, when determining the best matching atom, not only the matching judgment standard of the fitness value is considered, but also the sparsity (also called robust support) judgment of each atom is considered, specifically according to the formula XS_d＝||R^d||_ξCalculating to obtain the sparsity XS of the d-th better time frequency parameter_dAnd according to the calculation result (namely the robust support of each atom), selecting the optimal time-frequency parameter with the minimum sparsity as the optimal time-frequency parameter r_z. Therefore, the determined best matching atom is more accurate, and the intrinsic characteristics of the signal can be more accurately expressed.

To further ensure the accuracy of signal extraction and the accuracy of each of the separated superimposed signals, the best matching atoms need to be selected more accurately. In this embodiment, before the echo signal separation in step three, an iterative decomposition optimal atom set updating judgment needs to be performed, so as to select a more accurate optimal matching atom.

When the best atom set updating judgment of iterative decomposition is actually carried out, judging R^K(t)||₂Whether less than epsilon: when | | | R^K(t)||₂If the number is less than epsilon, judging that the optimal atom set of iterative decomposition in the step two does not need to be updated, and entering a step three; otherwise, when R^K(t)||₂When the number of atoms is more than or equal to epsilon, judging that the optimal atom set of the iterative decomposition in the step two needs to be updated, and enteringStep D;

step D, updating the iterative decomposition optimal atom set: updating the iterative decomposition optimal atomic set in the step two by adopting the data processing equipment 2, wherein the method comprises the following steps:

step D11, atom random selection: randomly taking out an optimal matching atom from the iterative decomposition optimal atom set at the moment by adopting the data processing equipment 2 as an atom to be optimized, wherein the atom to be optimized is marked as

Wherein j1 is a positive integer and j1 is not less than 1 and not more than K;

For the best matching atom

When searching, according to the preset s_j1'、u_j1'、v_j1' and w_j1' the value range is found out by adopting the data processing equipment 2 and calling the optimizing algorithm module to enable the fitness value to be fitness (r)_j1') the largest optimal time-frequency parameter, the found optimal time-frequency parameter being the time-frequency parameter r_j1'; according to the formula

Solving for the best matching atom

Wherein the content of the first and second substances,

to represent

And

inner product of (d);

ψ'₀(t) is the sum of K-1 of said atoms to be treated in step D11;

d13, atom replacement judgment, atom replacement and iterative decomposition optimal atom set updating: adopting a data processing device 2 and calling a residual value judging module, a fitness value judging module or a sparsity judging module to judge whether the atoms to be optimized in the step D11 need to be replaced or not, and replacing the atoms to be optimized according to the judgment result;

when the data processing device 2 is adopted and the residual value judging module is called to judge whether the atoms to be optimized in the step D11 need to be replaced, the residual value R after the replacement is carried out_j1'^K(t)||_ξWhether it is less than the residue before replacement | | R_j1 ^K(t)||_ξAnd (4) judging: when | | | R_j1'^K(t)||_ξ＜||R_j1 ^K(t)||_ξIf so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

Obtaining the updated iterative decomposition optimal atom set; otherwise, it is judged as noneReplacing the atoms to be optimized in the step D11, and entering a step D15;

wherein R is_j1'^K(t)＝F(t)-ψ_j1'(t)，

adopting the data processing equipment 2 and calling a Fitness value judging module to judge whether the atoms to be optimized in the step D11 need to be replaced, and according to the replaced Fitness value Fitness (r)_j1') whether it is greater than the pre-replacement Fitness value Fitness (r)_j1) And (4) judging: when Fitness (r)_j1')＞Fitness(r_j1) If so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

wherein the content of the first and second substances,

represents R^j1-1(t) and

represents R^j1-1(t) and

inner product of (d);

when the data processing equipment 2 is adopted and the sparsity judgment module is called to judge whether the atoms to be optimized in the step D11 need to be replaced, according to | | R^j1'||_ξWhether or not less than R^j1||_ξAnd (4) judging: when | | | R^j1'||_ξ＜||R^j1||_ξIf so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

wherein R is^j1' is

Amount of residual error of

R^j1Is composed of

Amount of residual error of

step D16, residual error quantity judgment: and D15, judging the optimized residual quantity of the best matching atoms: when | | | R "_j1 ^K(t)||₂If the number is less than epsilon, finishing the updating process of the iterative decomposition optimal atom set, and entering the step three; otherwise, when | | R "_j1 ^K(t)||₂When the value is more than or equal to epsilon, returning to the step D15;

Wherein, R is^K(t) is an N' x 1-dimensional vector, | | R^K(t)||₂Is R^KThe 1/2 th power of the 2 nd power sum of the absolute values of the N' elements in (t).

R 'in step D14'_j1 ^K(t) is a vector of N 'x 1 dimension, | R'_j1 ^K(t)||₂Is R'_j1 ^K(t)R'_j ^mThe 1/2 th power of the 2 nd power sum of the absolute values of the N elements in (t).

R in step D16 "_j1 ^K(t) is an N' x 1-dimensional vector, | | R "_j1 ^K(t)||₂Is R'_j1 ^KThe 1/2 th power of the 2 nd power sum of the absolute values of the N' elements in (t).

R in the step D13_j1'^K(t) is an N' x 1-dimensional vector, | | R_j1'^K(t)_ξIs R_j1'^KAnd (t) the power of 1/xi of the xi power sum of the absolute values of the N' elements.

Said R_j1 ^K(t) is an N' x 1-dimensional vector, | | R_j1 ^K(t)_ξIs R_j1 ^KAnd (t) the power of 1/xi of the xi power sum of the absolute values of the N' elements.

Said R^j1Is a vector of N' x 1 dimension, | | R^j1'||_ξIs R^j1The absolute values of the N elements in the' are 1/xi to the power of the xi power sum. Said R^j1Is a vector of dimension N' x 1, | | R^j1||_ξIs R^j1The absolute values of the N' elements are 1/xi to the power of the xi power sum.

Before searching for the best matching atom in step D12, s is first searched_j1'、u_j1'、v_j1' and w_j1' the value ranges are set respectively, and the set s_j1' the value range is the same as the value range of the scale parameter set in the step C1, and the set u is_j1' the value range is the same as the value range of the displacement parameter set in the step C1, and v is set_j1' the value range is the same as the value range of the frequency parameter set in the step C1, and w is set_j1The range of values of' is the same as the range of values of the phase parameter set in step C1. Thus, s_j1'is selected from the range of [1, N']And s_j1'∈[1,N']，u_j1'is selected from the range of [0, N']And u is_j1'∈[0,N']，v_j1' has a value range of

And is

w_j1' has a value range of[0,2π]And w_j1'∈[0,2π]。

In this example, s_zIs in the value range of [1, N']And s_z∈[1,N']，u_zIs [0, N']And u is_z∈[0,N']，v_zHas a value range of

And is

w_zHas a value range of [0,2 pi]And w_z∈[0,2π]。

In the actual use process, no matter the conventional matching pursuit algorithm carries out sparse decomposition after establishing an over-complete atom library, or the optimizing algorithm module is utilized to optimize and find out the optimal matching atom to complete signal sparse decomposition, the method has certain limitations, and the optimal matching atom is obtained under certain limiting conditions, so that when the two methods are adopted to carry out signal extraction, the accuracy of signal extraction is only relatively high. When the overcomplete atom library is adopted for sparse decomposition, the value range of each parameter in the time-frequency parameters and the discretization method both have great influence on the generated overcomplete atom library, and the finally determined overcomplete atom library cannot include all atoms and inevitably omits one or more optimal matching atoms, so that the accuracy of signal extraction is influenced. When the optimization algorithm module is used for optimizing and finding out the best matching atoms, although the signal extraction speed can be improved, the search in a continuous interval can be realized, the found best matching atoms are only the best matching atoms to a certain extent or within a certain range under the influence of the advantages and the disadvantages of the algorithm in the optimization algorithm module, such as the search step length, the search strategy, the search termination condition and the like, so that the accuracy of signal extraction can be influenced to a certain extent, and the accuracy of the best matching atoms can be correspondingly influenced, thereby influencing the accuracy of signal separation.

As can be seen from the above, after the signal sparse decomposition is completed in the second step, residual quantity judgment is further performed to judge whether the current iterative decomposition optimal atom set meets the preset signal extraction precision requirement, and if not, the step D1 is performed to perform optimal matching atom optimization, so as to further improve the accuracy of the optimal matching atoms. Therefore, after the signal sparse decomposition is completed in the step two, whether the optimal atom set subjected to iterative decomposition after the signal sparse decomposition in the step two meets the preset requirement is judged according to the residual error quantity judgment result, and the verification link of the signal extraction precision (also called atom selection precision) is added, so that the signal extraction precision can be further improved, and the accuracy of signal separation is ensured.

When the optimal matching atom optimization is performed in the step D1, the optimal matching atom optimization method is reasonable in design, convenient to implement, and good in use effect, and a best matching atom is randomly selected from the iterative decomposition optimal atom set for optimization, and after the optimization is completed by a best matching atom, whether the iterative decomposition optimal atom set satisfies the signal extraction precision requirement at that time is judged through the residual error amount, and whether the optimization of the rest of the best matching atoms is required to be continued is determined according to the judgment result. Therefore, the method is simple and convenient to realize, can randomly select an optimal matching atom for optimization, has no limit on the atom optimization sequence, can judge the residual error amount once when the optimization process of the optimal matching atom is finished, can realize the combination of quick optimization and real-time judgment of the optimization result, can effectively simplify the optimization process of the optimal matching atom, can quickly achieve the aim of optimizing the optimal matching atom, and effectively improves the signal separation precision.

When the atom to be optimized is optimized, the method for searching the best matching atom corresponding to the atom to be optimized (i.e., the method for searching the best matching atom in step D12) is reasonably designed, and the best matching atom better than the atom to be optimized can be simply, conveniently and quickly found out.

Found time frequency parameter r_j1' to make the fitness value fitness (r)_j1') the largest best time-frequency parameter;

due to the fact that

And psi'₀(t) is the sum of K-1 of the atoms to be treated in step D11, thus

Subtracting the residual error of K-1 atoms to be processed except the atom to be optimized from the ultrasonic echo signal F (t) (namely the original signal), thereby obtaining the ultrasonic echo signal

Is a residual signal directly related to the atom to be optimized, thus making use of

Finding the time-frequency parameter r as an evaluation_j1' the index is more targeted, except for the residual signal of K-1 atoms to be processed in the optimal atom set of iterative decomposition at the moment

Directly related to the atom to be optimized, and finding out the time-frequency parameter r by using an optimization algorithm module_j1The process of' is not affected by other atoms (i.e., K-1 of the atoms to be processed) and the probability of finding the best matching atom better than the atom to be optimized is higher, while the best matching atom obtained

The K-1 atoms to be processed in the optimal atom set subjected to iterative decomposition at the moment are not affected, signal sparse decomposition is not needed again, the atom replacement of the atoms to be optimized is completed only according to the method in the step D13, the using effect is very good, the signal extraction precision can be simply, conveniently and quickly improved, and the accuracy of the selected atoms is ensured.

When the atom replacement judgment and the atom replacement are performed in the step D13, any one of the methods of residual value judgment, fitness value judgment and sparsity judgment is used for atom replacement judgment, any one of the methods can be selected for atom replacement judgment, the use mode is flexible, and each atom replacement judgment method can realize effective atom replacement judgment.

When the residual value judgment module is called to judge whether the atoms to be optimized in the step D11 need to be replaced, the residual value R after the replacement is carried out_j1'^K(t)||_ξWhether it is less than the residue before replacement | | R_j1 ^K(t)||_ξThe judgment result judges whether to replace the atoms to be optimized, and the atoms with smaller residual values are selected to ensure that the residual quantity of the signals is smaller, thereby effectively improving the signal extraction precision and ensuring that the extracted signals further approach the original signals.

Calling a Fitness value judging module to judge whether the atoms to be optimized in the step D11 need to be replaced or not, and according to the Fitness value Fitness (r) after replacement_j1') whether it is greater than the pre-replacement Fitness value Fitness (r)_j1) And judging whether to replace the atoms to be optimized, and selecting the atoms with larger fitness value to reduce the residual quantity of the signal, thereby effectively improving the signal extraction precision and enabling the extracted signal to further approach the original signal.

And D11, when the sparsity judging module is called to judge whether the atoms to be optimized need to be replaced, judging whether the atoms need to be replaced according to the minimum robust support, and selecting the atoms with lower robust support, so that the signal characteristics can be better matched, the representation of the signal is sparser, the purpose of effectively improving the signal extraction precision is achieved, and the extracted signal is further close to the original signal.

Wherein the content of the first and second substances,

R^j1'(t_i) Is R^j1' the signal value at the ith sampling instant, i.e. R^j1The ith signal value of.

In this embodiment, after performing signal sparse decomposition in step D1, the data processing device 2 is used to synchronously store the iterative decomposition optimal atom set into the data storage 3, and the data storage 3 is connected to the data processing device 2. The ultrasonic flaw detection device 1, the data processing equipment 2 and the data storage 3 form a signal separation processing system, which is shown in detail in fig. 2.

After the atom replacement judgment and the atom replacement are performed in step D13, the updated iterative decomposition optimal atom set is synchronously stored by using the data processing device 2.

After signal sparse decomposition is performed in step D1, when the data processing device 2 is used to synchronously store the iterative decomposition optimal atom set into the data memory 3, the K optimal matching atoms in the iterative decomposition optimal atom set are respectively stored according to the iterative decomposition order; wherein the content of the first and second substances,

In this embodiment, when the best matching atom in step D1 is optimized, the best matching atom in the iterative decomposition best atom set is optimized according to the storage sequence;

In actual use, when the best matching atom in step D1 is optimized, the best matching atom in the iterative decomposition best atom set may not be optimized in the storage order.

The epsilon is a preset residual quantity judgment threshold, and the value of the epsilon can be limited according to specific requirements during actual use.

In this embodiment, the epsilon ═ e is described^-5。

In actual use, the value of epsilon can be correspondingly adjusted according to specific requirements.

In this embodiment, ξ ═ 1 in step C2.

When the device is actually used, the value of xi can be correspondingly adjusted according to specific requirements.

When signal sparse decomposition is performed in the second step, finding out K optimal matching atoms in the iterative decomposition optimal atom set from first to last by adopting the data processing equipment 2 according to the methods from the step C1 to the step C2;

the optimizing algorithm module in the step C1 is an artificial bee colony algorithm module.

In practical use, the optimization algorithm module can also be other optimization algorithm modules, such as a genetic algorithm module, an artificial fish swarm algorithm module and the like.

In this embodiment, step two is as described

The best matching atom is found when the z-th iterative decomposition is carried out on the ultrasonic echo signal F (t).

The optimizing algorithm module in the step C1 is an artificial bee colony algorithm module;

when the time-frequency parameter optimization is performed in the step C1, the process is as follows:

step 2011, parameter initialization: setting the maximum iteration times MC, the number SN of honey sources, the number of employed bees, the number of observation bees and the maximum exploitation times limit of the honey sources of the artificial bee colony algorithm module by adopting the data processing equipment 2; meanwhile, SN different honey sources are randomly generated by adopting the data processing equipment 2, the SN honey sources are all honey sources to be mined, and the pth honey source in the generated SN honey sources is recorded as a 4-dimensional vector X_p＝(X_1p,X_2p,X_3p,X_4p) Each honey source is a time-frequency parameter; the number of the employed bees and the number of the observation bees are SN, and each generated honey source is distributed to one employed bee;

wherein p is a positive integer and p ═ 1, 2, …, SN; x_1pAnd s preset in step 201_nHave the same value range of X_2pAnd u preset in step 201_nHave the same value range of X_3pAnd v preset in step 201_nHave the same value range of X_4pAnd w preset in step 201_nThe value ranges of (A) are the same;

the initialization of the parameters in this step is completedAfter the completion, the data processing equipment 2 is adopted to respectively calculate the fitness values of the SN honey sources, and all the fitness values in the SN honey sources are larger than thrd_zAnd each fitness value > thrd is recorded synchronously by means of the data processing device 2_zThe honey sources are marked as the better time frequency parameters;

step 2012, hiring bee neighborhood search: each hiring bee carries out neighborhood search on the allocated honey source, if the fitness value of the searched new honey source is larger than that of the original honey source, the new honey source is used as the honey source to be exploited, which is searched by the hiring bee, and the exploited frequency is set to be 0; otherwise, adding 1 to the mined times of the original honey source;

in this step, in the process of neighborhood searching for the allocated honey source by each hiring bee, all the searched fitness values > thrd are obtained by the data processing equipment 2_zAnd each recorded fitness value > thrd is recorded synchronously with the data processing device 2_zThe honey sources are marked as the better time frequency parameters;

step 2013, searching of bee neighborhood observation: calculating the selection probability of each honey source searched by the hiring bee according to the fitness values of all the honey sources searched by the hiring bee in the step 2012; the observation bees select honey sources for honey collection from all the honey sources searched by the employment bees as new honey sources according to the calculated selection probability of each honey source;

the observation bee carries out neighborhood search on the selected honey source, if the fitness value of the searched new honey source is larger than that of the original honey source, the observation bee is changed into a employment bee, the new honey source is used as the searched honey source to be exploited, and the exploited frequency is set to be 0; otherwise, if the honey source and the employed bees are not changed, adding 1 to the mined times of the original honey source;

in this step, in the process of neighborhood searching for the selected honey source by the observation bees, all the searched fitness values > thrd are obtained by the data processing equipment 2_zAnd each recorded fitness value > thrd is recorded synchronously with the data processing device 2_zAll honey sources are marked as the better time frequencyA parameter;

step 2014, recording the optimal honey source in real time: after the search of the employed bee neighborhood and the search of the observation bee neighborhood are finished, obtaining the optimal honey source at the moment and synchronously recording, wherein the iteration times of the artificial bee colony algorithm module is added with 1;

in the process of employing bee neighborhood searching and observing bee neighborhood searching, if the mined times of the honey source reach the maximum mined times l imit of the honey source, the observing bee is converted into a detecting bee, a new honey source is generated through the detecting bee, and the mined times are set to be 0;

step 2015, repeating steps 2012 to 2014 for multiple times until the iteration number of the artificial bee colony algorithm module reaches the maximum iteration number MC, and the obtained optimal honey source is the optimal time-frequency parameter in step C1;

when the fitness values of the SN honey sources are respectively calculated in the step 2011, when the employment bee neighborhood search is performed in the step 2012 and the observation bee neighborhood search is performed in the step 2013, the fitness values of any honey source are the Gabor atom and the R atom corresponding to the honey source^z-1(t) inner product;

after repeating steps 2012 to 2014 a plurality of times in step 2015, k described in step C1 is obtained_zThe better time frequency parameters are selected;

in step 2011, step 2012, step 2013 and step 2015, the total number of all honey sources marked as the better time-frequency parameter by the data processing device 2 is k in step C1_zEach honey source marked as the preferred time-frequency parameter by the data processing device 2 is one of the preferred time-frequency parameters.

Wherein R is^z-1And (t) is the residual quantity of the ultrasonic echo signal F (t) after z-1 times of iterative decomposition in the step two.

The time frequency parameter r_zThe corresponding Gabor atom is

The optimal honey source obtained in step 2014 is the optimal honey source obtained in one iteration process, and the optimal honey source obtained in step 2015 is the optimal honey source with the maximum fitness value in the optimal honey sources obtained in the MC iteration process.

In this embodiment, the fitness value of the original honey source

Represents R^z-1(t) and

inner product of (d);

in this step, the number of the honey sources to be mined, which are searched by the hiring bees, is multiple, and all the honey sources to be mined, which are searched by the hiring bees, are the honey sources searched by the hiring bees.

The fitness value of any searched new honey source is the Gabor atom and R corresponding to the honey source^z-1(t) inner product.

In this embodiment, the best matching atom in step D12

When searching, the data processing equipment 2 is adopted and the optimizing algorithm module is called to search

Time-frequency parameter r of_j1'The optimizing algorithm module is an artificial bee colony algorithm module, and the process is as follows:

step D21, parameter initialization: setting the maximum iteration times MC ', the number SN ' of honey sources, the number of employed bees, the number of observation bees and the maximum exploitation times limit ' of the honey sources of the artificial bee colony algorithm module by adopting the data processing equipment 2; meanwhile, SN 'different honey sources are randomly generated by the data processing equipment 2, the SN' honey sources are all honey sources to be mined, and the p 'th honey source in the generated SN' honey sources is recorded as a 4-dimensional vector X_p'＝(X_1p',X_2p',X_3p',X_4p') Each honey source is a time-frequency parameter(ii) a The number of the employed bees and the number of the observation bees are SN', and each generated honey source is distributed to one employed bee;

wherein p ' is a positive integer and p ' ═ 1, 2, …, SN '; x_1p'And s preset in the step two_zHave the same value range of X_2p'And u preset in the step two_zHave the same value range of X_3p'And v preset in the step two_zHave the same value range of X_4p'And w preset in the step two_zThe value ranges of (A) are the same;

step D22, hiring bee neighborhood search: each hiring bee carries out neighborhood search on the allocated honey source, if the fitness value of the searched new honey source is larger than that of the original honey source, the new honey source is used as the honey source to be exploited, which is searched by the hiring bee, and the exploited frequency is set to be 0; otherwise, adding 1 to the mined times of the original honey source;

step D23, observation bee neighborhood searching: calculating the selection probability of each honey source searched by the employed bees according to the fitness values of all the honey sources searched by the employed bees in the step D22; the observation bees select honey sources for honey collection from all the honey sources searched by the employment bees as new honey sources according to the calculated selection probability of each honey source;

step D24, recording the optimal honey source in real time: after the search of the employed bee neighborhood and the search of the observation bee neighborhood are finished, obtaining the optimal honey source at the moment and synchronously recording, wherein the iteration times of the artificial bee colony algorithm module is added with 1;

d25, repeating the steps D22-D23 for multiple times until the iteration number of the artificial bee colony algorithm module reaches the maximum iteration number MC, and obtaining the optimal honey source as

Time-frequency parameter r of_j1'，r_j1'＝(s_j1',u_j1',v_j1',w_j1')；

When the employed bee neighborhood search is performed in the step D22 and the observation bee neighborhood search is performed in the step D23, the fitness value of any honey source is the Gabor atom and R corresponding to the honey source^z-1(t) inner product.

In this example, X_1p'Is in the value range of [1, N']And X_1p∈[1,N']，X_2p'Is [0, N']And X_2p∈[0,N']，X_3p'Has a value range of

And is

X_4pHas a value range of [0,2 pi]And X_4p'∈[0,2π]。

The time frequency parameter r_j1'The corresponding Gabor atom is

The optimal honey source obtained in the step D24 is the optimal honey source obtained in the one-time iterative process, and the optimal honey source obtained in the step D25 is the optimal honey source with the maximum fitness value among the optimal honey sources obtained in the MC' time iterative process.

In addition, the method adopts border-crossing retracing processing, adopts the hiring bees and the observation bees to carry out neighborhood search, carries out boundary detection on a new honey source after the new honey source is generated, and carries out border-crossing retracing operation on the new honey source if the new honey source exceeds the upper and lower bounds. And when the boundary-crossing retracing operation is carried out on the new honey source, the boundary-crossing retracing operation is respectively carried out on the 4 elements of the new honey source according to the maximum value and the minimum value of the four elements of the honey source. Carrying out boundary detection on the new honey source, and respectively carrying out-of-range judgment on 4 elements of the new honey source according to the maximum value and the minimum value of the four elements of the honey source; and performing border-crossing retracing operation on 4 elements of the new honey source according to the border-crossing judgment result, and obtaining the honey source after the border-crossing retracing operation, thereby avoiding the phenomenon of error search.

Wherein the new honey source

For new honey source

Q element of (2)

When the out-of-range judgment is performed, when

When it is determined that

Does not exceed the boundary, does not need to be aligned

Performing boundary-crossing retracing operation; when in use

When it is determined that

Beyond the lower bound, according to the formula

After obtaining an out-of-range retracing operation

When in use

When it is determined that

Beyond the upper bound, according to the formula

After obtaining an out-of-range retracing operation

In step 2013, when the selection probability of each honey source searched by the employed bees is calculated according to the fitness values of all the honey sources searched by the employed bees in step 2012, the selection probability of each honey source is calculated according to the roulette mode. The selected probability of any honey source is the ratio of the fitness value of the honey source to the sum of the fitness values of all honey sources searched by the employment bees. And 2013, selecting the honey source with the highest selected probability as a new honey source when the honey source for honey collection is selected from all the honey sources searched by the employed bees according to the calculated selected probability of each honey source by the observation bees.

Correspondingly, when the selection probability of each honey source searched by the employed bees is calculated according to the fitness values of all the honey sources searched by the employed bees in the step D22 in the step D23, the selection probability of each honey source is calculated according to a roulette mode. The selected probability of any honey source is the ratio of the fitness value of the honey source to the sum of the fitness values of all honey sources searched by the employment bees. And D23, selecting the honey source with the maximum selected probability as a new honey source when the honey source for honey collection is selected from all the honey sources searched by the employed bees according to the calculated selected probability of each honey source by the observation bees.

When the search of the neighborhood of the observer bees is performed in the step 2013 and the search of the neighborhood of the observer bees is performed in the step D23, in order to accelerate the search speed, the search mode is changed from the random search to the following search mode: judging whether the fitness value of the honey source searched next randomly is larger than the fitness value of the honey source at the center position of the bee colony at the moment, and taking the honey source searched next randomly as a new honey source when the fitness value of the honey source searched next randomly is larger than the fitness value of the honey source at the center position of the bee colony at the moment; otherwise, the honey source at the center position of the bee colony at the moment is used as a new honey source to improve the searching speed of the algorithm. And the honey source at the central position of the bee colony at the moment is the average value of the sum of all the honey sources searched at the moment.

Because the distance from the optimal atom is closer and closer along with the increase of the search times of the bee colony, in order to accelerate the optimization speed and avoid trapping in local optimization, the honey source concentration (namely fitness value) of the next search position and the central position of the bee is compared when the bee is observed for searching, and a new honey source is determined according to the comparison result, so that the search step length is increased, and the speed of the bee moving towards the optimal atom direction is accelerated.

In this embodiment, when the parameters are initialized in step 2011 and step D21, the initial bee colony is generated by using a uniform distribution method.

The randomness of initial swarm distribution in the original artificial swarm algorithm can cause uncertainty of a search space, and if the initial swarm search space does not contain a global optimal solution and cannot cover a region of the global optimal solution in limited searches, premature convergence can be caused. The initial bee colony is generated by adopting a uniform distribution method, so that the premature convergence phenomenon can be effectively avoided.

In this embodiment, in the process of initializing parameters in step 2011, when SN honey sources are generated, the SN honey sources are generated according to a formula

Calculating the qth element X of the pth honey source in SN honey sources_qpWherein q is a positive integer and q is 1, 2, 3 or 4; x_qupIs the maximum value of the qth element of the honey source, X_qlowIs the minimum value of the qth element of the honey source.

Wherein the 1 st element of the honey source has a maximum value of N' and a minimum value of 1, thus X_1upN' and X_1low＝1。The maximum value of the 2 nd element of the honey source is N' and its minimum value is 0, thus X_2upN' and X_2low0. The maximum value of the 3 rd element of the honey source is

And its minimum value is 0, thus

And X_3low0. The maximum value of the 4 th element of the honey source is 2 pi and the minimum value thereof is 0, thus X _4up2 pi and X_4low＝0。

In this embodiment, in step 2012, when the neighboring search of the employer bee is performed, the neighboring search is performed by the employer bee near the current honey source location to generate a new honey source, and the new honey source location is determined according to formula X_p*＝X_p+φ_p(X_p-X_l) Making a determination wherein X_pFor the currently searched source of raw honey, phi_pIs [ -1,1 [ ]]A random number in the range, X_lIs a random honey source, X_p*Is a new source of honey, passes through phi_pThe range of new honey sources is limited.

In this embodiment, in the process of initializing parameters in step D21, when SN '(i.e., SN') honey sources are generated, the SN honey sources are generated according to a formula

Calculating the qth element X of the pth honey source in SN honey sources_qp'。

In this embodiment, in the step D22, when the neighboring search of the employed bee is performed, the neighboring search is performed by the employed bee near the current honey source location to generate a new honey source, and the new honey source location is determined according to the formula X_p'*＝X_p'+φ_p(X_p'*-X_l) Making a determination wherein X_p'For the currently searched source of raw honey, phi_pIs [ -1,1 [ ]]A random number in the range, X_lIs a random honey source, X_p'*Is a new source of honey, passes through phi_pThe range of new honey sources is limited.

In the present embodiment, the first and second electrodes are,r 'in step D14'_j1 ^K(t) according to the formula

Performing a calculation in which

To this end the iteration decomposes the n1' th best matching atom in the best set of atoms, n1' is a positive integer and n1' is 1, 2, …, K; a is_n1'Is composed of

Expansion coefficients of residual quantities after n1'-1 times of iterative decomposition of F (t) according to the first n1' -1 best matching atom pairs in the iterative decomposition best atom set at this time;

r' described in step D16 "_j1 ^K(t) according to the formula

Performing a calculation in which

To this end the iteration decomposes the n2' th best matching atom in the best set of atoms, n2' is a positive integer and n2' is 1, 2, …, K; a is_n2'Is composed of

And (c) an expansion coefficient of a residual quantity after n2'-1 times of iterative decomposition of the first n2' -1 best matching atom pairs F (t) in the optimal atom set according to the iterative decomposition at the moment.

In this embodiment, R in step D13^j1-1(t) is the residual quantity after j1-1 times of iterative decomposition is carried out on the first j1-1 best matching atom pairs F (t) in the iterative decomposition best atom set before atom replacement judgment in the step.

For R in step D13^j1-1(t) when calculating, iteratively decomposing the optimal set of atoms according to the values of the set of atoms before performing the atom replacement judgment in this stepAnd

performing a calculation wherein k ' is a positive integer and k ' ═ 1, 2, …, j1-1, k ' < j 1;

for the k' th best matching atom in the iterative decomposition best atom set before the atom replacement judgment in the step, a_k'Is composed of

And f (t) is subjected to k '-1 times of iterative decomposition according to the expansion coefficient of residual quantity after the first k' -1 best matching atoms in the iterative decomposition best atom set are subjected to atom replacement judgment in the step.

In this embodiment, after the echo signals are separated in step three, signal extraction needs to be performed on the K separated signals respectively; the signal extraction methods of the K separation signals are the same.

When any one of the K separation signals is subjected to signal extraction, the data processing equipment 2 is adopted to extract the separation signal, and the separation signal is a signal to be processed and is marked as a signal f (t);

when the data processing device 2 is used for extracting the signal f (t), the process is as follows:

step A1, sparse decomposition based on an optimization algorithm: adopting a data processing device 2 and calling a sparse decomposition module to carry out iterative decomposition processing on the signal f (t), and converting the signal f (t) into the signal f (t)

decomposing the nth best matching atom in the best atom set for the iteration;

adopting data processing equipment 2 and calling an optimization algorithm module to find out the best matching atom for the nth iterative decomposition;

is a Gabor atom and

wherein the function ψ (t) is a Gaussian window function and

r_nis composed of

in this step, the best matching atom is found

According to a preset s_n、u_n、v_nAnd w_nThe data processing device 2 is adopted and an optimization algorithm module is called to find out the Fitness value Fitness (r)_n) The maximum optimal time frequency parameter, the found optimal time frequency parameter is the time frequency parameter r_n；

represents Rⁿ-¹(t) and

wherein | Rm (t) | non-woven cells₂For R in the step A1^m(t) 2-norm, epsilon is a preset residual quantity judgment threshold;

step A3, optimizing the best matching atom, the process is as follows:

step A31, atom random selection: randomly taking out an optimal matching atom from the iterative decomposition optimal atom set at the moment by adopting the data processing equipment 2 as an atom to be optimized, wherein the atom to be optimized is marked as

For the best matching atom

When searching, according to the preset s_j'、u_j'、v_j' and w_j' the value range is found out by adopting the data processing equipment 2 and calling the optimizing algorithm module to enable the fitness value to be fitness (r)_j') the largest optimal time-frequency parameter, the found optimal time-frequency parameter being the time-frequency parameter r_j'; according to the formula

Solving for the best matching atom

Wherein the content of the first and second substances,

to represent

And

inner product of (d);

ψ₀(t) is the sum of m-1 of said atoms to be treated in step A31;

step A33, atom replacement judgment and atom replacement: adopting a data processing device 2 and calling a residual value judging module, a fitness value judging module or a sparsity judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and replacing the atoms to be optimized according to the judgment result;

when the data processing device 2 is adopted and the residual value judging module is called to judge whether the atoms to be optimized in the step A31 need to be replaced, the residual value R after the replacement is carried out_j'^m(t)||_ξWhether it is less than the residue before replacement | | R_j ^m(t)||_ξAnd (4) judging: when | | | R_j'^m(t)||_ξ＜||R_j ^m(t)||_ξWhen the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein R is_j'^m(t)＝f(t)-ψ_j'(t)，

adopting the data processing equipment 2 and calling a Fitness value judging module to judge whether the atoms to be optimized in the step A31 need to be replaced, and according to the replaced Fitness value Fitness (r)_j') whether it is greater than the pre-replacement Fitness value Fitness (r)_j) And (4) judging: when Fitness (r)_j')＞Fitness(r_j) When the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein the content of the first and second substances,

represents R^j-1(t) and

represents R^j-1(t) and

inner product of (d);

when the data processing device 2 is adopted and the sparsity judgment module is called to judge whether the atoms to be optimized in the step A31 need to be replaced, according to | | R^j'||_ξWhether or not less than R^j||_ξAnd (4) judging: when | | | R^j'||_ξ＜||R^j||_ξWhen the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein R is^j' is

Amount of residual error of

R^jIs composed of

Amount of residual error of

R^j'||_ξRepresents R^jThe xi-norm, | | R of `^j||_ξRepresents R^jξ -norm of;

wherein, | R'_j ^m(t)||₂Is R'_j ^m(t) 2-norm; r'_j ^m(t) is the residual quantity after m iterative decompositions are performed on f (t) according to m best matching atoms in the iterative decomposition best atom set at the moment;

step A4, signal reconstruction: obtaining an approximate signal f' (t) of the signal f (t) by adopting the data processing device 2 according to the iterative decomposition optimal atom set at the moment; wherein the approximation signal f "(t) is a signal extracted from the signal f (t),

wherein

R in step A2^m(t) is a vector of dimension Nx 1, | | R^m(t)||₂Is R^mThe 1/2 th power of the 2 nd power sum of the absolute values of the N elements in (t). Where N is a positive integer and is the signal length of the signal f (t). Wherein N is N'.

R 'in step A34'_j ^m(t) is an Nx 1-dimensional vector, | R'_j ^m(t)||₂Is R'_j ^mThe 1/2 th power of the 2 nd power sum of the absolute values of the N elements in (t).

R in step A36 "_j ^m(t) is an Nx 1-dimensional vector, | | R "_j ^m(t)||₂Is R'_j ^mThe 1/2 th power of the 2 nd power sum of the absolute values of the N elements in (t).

R in step A33_j'^m(t) is a vector of dimension Nx 1, | | R_j'^m(t)||_ξIs R_j'^mAnd (t) the power of 1/xi of the xi power sum of the absolute values of the N elements.

Said R_j ^m(t) is a vector of dimension Nx 1, | | R_j ^m(t)||_ξIs R_j ^mAnd (t) the power of 1/xi of the xi power sum of the absolute values of the N elements.

Said R^j' is a vector of dimension Nx 1, | | R^j'||_ξIs R^jThe absolute values of the N elements in the' are 1/xi to the power of the xi power sum. Said R^jIs in the Nx 1 dimension directionAmount, | | R^j||_ξIs R^jThe absolute values of the N elements are 1/xi to the power of the xi power sum.

Wherein, in step A1

Described in step A1

Is the best matching atom when the signal f (t) to be processed is decomposed by the nth iteration.

Before searching for the best matching atom in step A32, s is first searched_j'、u_j'、v_j' and w_j' the value ranges are set respectively, and the set s_j' the value range is equal to s set in the step A1_nAre in the same value range, and the set u_j' the value range is equal to u set in the step A1_nHave the same value range, and set v_j' the value range is equal to v set in the step A1_nHave the same value range, set w_j' the value range is equal to w set in the step A1_nThe value ranges of (A) are the same.

In this embodiment, s in step A32_jThe value range of' and s_nHave the same value range of u_j' range of values and u_nHave the same value range, v_j' value range and v_nHave the same value range, w_j' range of values and w_nThe value ranges of (A) are the same. Thus, s_j' has a value range of [1, N]And s_j'∈[1,N]，u_j' has a value range of [0, N]And u is_j'∈[0,N]，v_j' has a value range of

And is

w_j' has a value range of [0,2 π]And w_j'∈[0,2π]。

After the signal sparse decomposition is completed in the step a1, the step a2 is further required to perform residual quantity judgment, and whether the iterative decomposition optimal atom set at this time meets the preset signal extraction precision requirement is judged, if not, the step A3 is required to perform optimal matching atom optimization, so as to further improve the accuracy of signal extraction. Therefore, after the signal sparse decomposition is completed in the step a1, according to the residual quantity judgment result in the step a2, whether the optimal atom set subjected to iterative decomposition after the signal sparse decomposition in the step a1 meets the preset signal extraction precision requirement is judged, and the verification link of the signal extraction precision is added, so that the signal extraction precision can be further improved, and the extracted signal further approaches to the original signal.

When the optimal matching atom optimization is performed in the step a3, the adopted optimal matching atom optimization method is reasonable in design, convenient to implement and good in use effect, a best matching atom is randomly selected from the iterative decomposition optimal atom set for optimization, after the optimization is completed by one best matching atom, whether the iterative decomposition optimal atom set meets the signal extraction precision requirement or not is judged through the residual quantity, and whether the optimization needs to be continuously performed on the rest of the best matching atoms or not is determined according to the judgment result. Therefore, the method is simple and convenient to realize, can randomly select an optimal matching atom for optimization, has no limit on the atom optimization sequence, can judge the residual error amount once when the optimization process of the optimal matching atom is finished, can realize the combination of quick optimization and real-time judgment of the optimization result, can effectively simplify the optimization process of the optimal matching atom, can quickly achieve the aim of optimizing the optimal matching atom, and effectively improves the signal extraction precision.

When the atom to be optimized is optimized, the method for searching the best matching atom corresponding to the atom to be optimized (i.e., the method for searching the best matching atom in step a 32) is reasonably designed, and the best matching atom better than the atom to be optimized can be simply, conveniently and quickly found out.

Found time frequency parameter r_j' to make the fitness value fitness (r)_j') the largest best time-frequency parameter;

due to the fact that

And psi₀(t) is the sum of m-1 of said atoms to be treated in step A31, thus

Subtracting the residual error of m-1 atoms to be processed except the atom to be optimized from the signal f (t) (i.e. the original signal) to be processed, thereby obtaining the final product

Finding the time-frequency parameter r as an evaluation_j' the index is more targeted, except for the residual signal of m-1 atoms to be processed in the best atom set in iterative decomposition at the moment

Directly related to the atom to be optimized, and finding out the time-frequency parameter r by using an optimization algorithm module_jThe process of' is not affected by other atoms (i.e., m-1 of the atoms to be processed) and the probability of finding the best matching atom better than the atom to be optimized is higher, while the best matching atom obtained

The m-1 atoms to be processed in the optimal atom set subjected to iterative decomposition at the moment are not affected, signal sparse decomposition is not needed again, the atom replacement of the atoms to be optimized is completed only according to the method in the step A33, and finally the step A4 is directly performed to reconstruct the signals, so that the method has a very good using effect and can simply, conveniently and quickly improve the signal extraction precision.

When the atom replacement judgment and the atom replacement are performed in the step a33, any method of residual value judgment, fitness value judgment or sparsity judgment is adopted for atom replacement judgment, any method can be selected for atom replacement judgment, the use mode is flexible, and each atom replacement judgment method can realize effective atom replacement judgment.

When the residual value judgment module is called to judge whether the atoms to be optimized in the step A31 need to be replaced, the residual value R after the replacement is carried out_j'^m(t)||_ξWhether it is less than the residue before replacement | | R_j ^m(t)||_ξThe judgment result judges whether to replace the atoms to be optimized, and the atoms with smaller residual values are selected to ensure that the residual quantity of the signals is smaller, thereby effectively improving the signal extraction precision and ensuring that the extracted signals further approach the original signals.

Calling a Fitness value judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and according to the Fitness value Fitness (r) after replacement_j') whether it is greater than the pre-replacement Fitness value Fitness (r)_j) And judging whether to replace the atoms to be optimized, and selecting the atoms with larger fitness value to reduce the residual quantity of the signal, thereby effectively improving the signal extraction precision and enabling the extracted signal to further approach the original signal.

And calling a sparsity judging module to judge whether to replace the atoms to be optimized in the step A31, judging whether to replace the atoms according to the minimum robust support, and selecting the atoms with lower robust support, so that the signal characteristics can be better matched, the signal representation is sparser, the purpose of effectively improving the signal extraction precision is achieved, and the extracted signal is further close to the original signal.

Wherein the content of the first and second substances,

R^j'(t_i) Is R^j' the signal value at the ith sampling instant, i.e. R^jThe ith signal value of.

In this embodiment, after signal sparse decomposition is performed in step a1, the optimal atom set for iterative decomposition is synchronously stored in a data storage 3 by using a data processing device 2, and the data storage 3 is connected to the data processing device 2;

after the atom replacement judgment and the atom replacement are performed in step a33, the updated iterative decomposition optimal atom set is synchronously stored by using the data processing device 2.

The signal sampling system 1, the data processing device 2 and the data memory 3 form a signal sampling and extracting system.

The best matching atom optimized in step A35 is one of the best matching atoms in the set of iteratively decomposed best atoms in step A1. The best matching atom that has completed optimization cannot be optimized again.

In this embodiment, after the optimization process of the best matching atom is completed in step a33, the best matching atom selected in step a31 is labeled as an optimized atom. Thus, the best matching atom optimized in step a35 is one of the best matching atoms in the best set of atoms other than the optimized atom for this time. Wherein, at this time, one of the best matching atoms in the set of iteratively decomposed best atoms that is not optimized is one of the best matching atoms in the set of iteratively decomposed best atoms other than the optimized atom at this time.

In this embodiment, after signal sparse decomposition is performed in step a1, when the data processing device 2 is used to synchronously store the iterative decomposition optimal atom set into the data memory 3, m optimal matching atoms in the iterative decomposition optimal atom set are respectively stored according to the iterative decomposition order; wherein the content of the first and second substances,

the best matching atom is found when the nth iterative decomposition is performed on the signal f (t) to be processed in step a 1.

In this embodiment, when the best matching atom in step a3 is optimized, the best matching atom in the iterative decomposition best atom set is optimized according to the storage sequence;

when the best matching atom in the step A3 is optimized, the best matching atom which is optimized first is the 1 st best matching atom in the iterative decomposition best atom set in the step a 1.

In actual use, when the best matching atom in step a3 is optimized, the best matching atom in the iterative decomposition best atom set may not be optimized in the order of storage.

The epsilon described in the step a2 is a preset residual quantity judgment threshold, and the value of epsilon can be limited according to specific needs when in actual use.

In this embodiment, epsilon ═ e in step a2^-5。

In this embodiment, ξ ═ 1 in step a 33.

In this embodiment, step A1

Finding the best matching atom when the signal f (t) to be processed is subjected to the nth iterative decomposition in the step A1;

when signal sparse decomposition is carried out in the step A1, finding m best matching atoms in the iterative decomposition best atom set in the step A1 from first to last by adopting data processing equipment (2);

the optimizing algorithm module in the step A1 is an artificial bee colony algorithm module. In practical use, the optimization algorithm module can also be other optimization algorithm modules, such as a genetic algorithm module, an artificial fish swarm algorithm module and the like.

In this example, R 'in step A34'_j ^m(t) according to the formula

To carry outIs calculated, wherein

To this end the iterative decomposition of the n1 th best matching atom in the best set of atoms, n1 being a positive integer and n1 being 1, 2, …, m; a is_n1Is composed of

Expansion coefficients of residual quantities after n1-1 iterative decompositions on f (t) according to the first n1-1 best matching atom pairs in the iterative decomposition best atom set at the moment;

r' described in step A36 "_j ^m(t) according to the formula

Performing a calculation in which

To this end the iterative decomposition of the n2 th best matching atom in the best set of atoms, n2 being a positive integer and n2 being 1, 2, …, m; a is_n2Is composed of

And (c) expansion coefficients of residual quantities after n2-1 times of iterative decomposition according to the first n2-1 best matching atom pairs f (t) in the iterative decomposition best atom set at the moment.

In this embodiment, R in step A33^j-1(t) is the residual quantity after j-1 iterative decompositions are performed on the first j-1 best matching atom pairs f (t) in the iterative decomposition best atom set according to the judgment of atom replacement in the step.

For R in step A33^j-1(t) when calculating, iteratively decomposing the optimal set of atoms and

performing a calculation, wherein k is a positive integer and k is 1, 2, …, j-1, k < j;

a is the k-th best matching atom in the iterative decomposition best atom set before the atom replacement judgment in the step_kIs composed of

And f (t) is subjected to k-1 times of iterative decomposition according to the expansion coefficient of residual quantity after the first k-1 best matching atoms in the iterative decomposition best atom set are subjected to atom replacement judgment in the step.

The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and all simple modifications, changes and equivalent structural changes made to the above embodiment according to the technical spirit of the present invention still fall within the protection scope of the technical solution of the present invention.

Claims

1. An echo signal separation method for ultrasonic flaw detection, characterized by comprising the steps of:

step one, ultrasonic echo signal acquisition and synchronous uploading: ultrasonic detection is carried out on a detected object by adopting an ultrasonic flaw detection device (1), an ultrasonic echo signal F (t) of the detected object is obtained, and the obtained ultrasonic echo signal F (t) is synchronously transmitted to a data processing device (2);

wherein f (t) ═ f (t)₁),f(t₂),...,f(t_N')]^TT represents a time parameter, t_iIs the ith sampling time f (t) of the ultrasonic flaw detection device (1)_i) The signal value is a signal value sampled at the ith sampling time of the ultrasonic flaw detection device (1), wherein i is a positive integer, i is 1, 2, 3, …, N' is a positive integer, and the signal value is the signal length of an ultrasonic echo signal F (t);

step two, signal sparse decomposition: adopting a data processing device (2) and calling a sparse decomposition module to carry out pair on the ultrasonic echo signals F (t)Performing iterative decomposition processing to convert the ultrasonic echo signal F (t)

for the best matching atom found at the z-th iterative decomposition,

is a Gabor atom and

wherein the function ψ (t) is a Gaussian window function and

r_zis composed of

step three, echo signal separation: separating the ultrasonic echo signals F (t) in the step I by adopting data processing equipment (2) according to K optimal matching atoms in the iterative decomposition optimal atom set at the moment to obtain K separation signals;

the K' th of the K separation signals is denoted as F_k'(t)，

Wherein the content of the first and second substances,

represents R^k'^-1(t) and

inner product of (d); r^k'^-1(t) is a residual quantity R after calling a matching tracking algorithm module and performing k '-1 iterative decomposition on the ultrasonic echo signal F (t) by utilizing the first k' -1 best matching atoms in the iterative decomposition best atom set at the moment⁰(t)＝F(t)；

In the third step, each separation signal is an ultrasonic echo signal at the position of a defect in the detected object;

in the first step, the ultrasonic echo signals F (t) are signals obtained by superposing the ultrasonic echo signals at the positions of K defects in the tested object, and each separated signal in the third step is a superposed signal;

before signal sparse decomposition is carried out in the second step, determining the number K of the superposed signals in the ultrasonic echo signal F (t);

when the number K of the superposed signals in the ultrasonic echo signals F (t) is determined, a data processing device (2) is adopted for determining, and the process is as follows:

step B1, determining wave crests and wave troughs: respectively determining all wave crests and all wave troughs of the ultrasonic echo signal F (t) in the step one by adopting a data processing device (2), and respectively synchronously recording the sampling time and the signal value of each wave crest and each wave trough;

step B2, removing extreme points: adopting data processing equipment (2) and calling a time domain extreme point eliminating module or a frequency domain extreme point eliminating module to eliminate extreme points to obtain M 'eliminated extreme points, and arranging the M' extreme points from front to back according to the sampling time sequence of each extreme point; wherein M' is a positive integer and is the total number of extreme points obtained after the extreme points are removed in the step;

when the data processing equipment (2) is adopted and a time domain extreme point rejecting module is called to reject the extreme points, rejecting the extreme points of which the absolute values of the signal values in all the extreme points determined in the step B1 are smaller than beta 'to obtain M' number of rejected extreme points; wherein, β' is a constant and the value range thereof is 0.1-0.35, and max | f (t) | is the maximum absolute value of the signal value in the ultrasonic echo signal f (t);

when the data processing equipment (2) is adopted and a frequency domain extreme point eliminating module is called to eliminate the extreme points, eliminating the extreme points of which the absolute values of the signal values in all the extreme points determined in the step B1 are smaller than beta to obtain M' eliminated extreme points; wherein β is a preset rejection threshold, β ═ α × max | y (f) |, α is a constant and has a value range of 0.25 to 0.35, y (f) is a frequency spectrum of the ultrasonic echo signal f (t), and max | y (f) | is a maximum absolute value of an amplitude value in the frequency spectrum of the ultrasonic echo signal f (t);

step B3, superposition judgment: respectively carrying out superposition judgment on two adjacent extreme points in the M 'extreme points in the step B2 from front to back by adopting data processing equipment (2) according to the sampling time sequence, and determining the total number of superposed extreme point pairs in the M' extreme points according to a superposition judgment result;

the superposition judgment methods of two adjacent extreme points in the M' extreme points are the same; when the M ' th extreme point and the M ' +1 th extreme point in the M ' extreme points are subjected to superposition judgment, a data processing device (2) is adopted and a peak judgment module or a trough judgment module is called for carrying out superposition judgment; wherein M 'is a positive integer and M' is 1, 2, …, M '-2, M' -1;

when data processing equipment (2) is adopted and a wave crest judging module is called to carry out superposition judgment, when the m 'th extreme point and the m' +1 th extreme point are both wave crests, the m 'th extreme point and the m' +1 th extreme point are judged to be one superposition extreme point pair;

when data processing equipment (2) is adopted and a wave trough judging module is called to carry out superposition judgment, when the m 'th extreme point and the m' +1 th extreme point are both wave troughs, the m 'th extreme point and the m' +1 th extreme point are judged to be one superposition extreme point pair;

2. An echo signal separating method for ultrasonic testing according to claim 1, characterized in that: in the first step, the ultrasonic flaw detection device (1) is connected with the data processing equipment (2) in a wired or wireless mode.

3. An echo signal separating method for ultrasonic testing according to claim 1, characterized in that: before the superposition judgment in the step B3, initializing the total number K ' of the superposition extremum point pairs in the M ' extremum points by using a data processing device (2), where K ' is 0;

4. An echo signal separating method for ultrasonic testing according to claim 1 or 2, characterized in that: when signal sparse decomposition is carried out in the step two, respectively searching K optimal matching atoms in the iterative decomposition optimal atom set from beginning to end by adopting data processing equipment (2);

said

For the best matching atom found with the data processing device (2),

time-frequency parameter r of_zIs the best time frequency parameter;

using a data processing device (2) and

when searching, pair

Time-frequency parameter r of_zThe search is performed as follows:

step C1, time frequency parameter optimization: calling an optimization algorithm module, finding out the optimal time frequency parameter with the maximum fitness value according to the preset value ranges of the scale parameter, the displacement parameter, the frequency parameter and the phase parameter in the time frequency parameters, and finding out the fitness value greater than thrd_zThe better time-frequency parameter;

Wherein the content of the first and second substances,

c' is a constant and is more than or equal to 2 and less than or equal to 3; r_z-1For the residual quantity, R, of the first z-1 of said best matching atoms which have been found by the data processing device (2) at that time_z-1＝F(t)-ψ_z-1(t)，ψ_z-1(t) for this purpose the sum of the first z-1 of said best matching atoms, | | R, determined by the data processing device (2) is used_z-1||₂Is R_z-12-norm of (d); when z is 1, R₀＝F(t)；

To k is paired_zSparsity XS of the d-th one of the better time-frequency parameters_dWhen doing the calculation, XS_d＝||R^d||_ξ，R^dIs psi_dResidual error amount of (t)，R^d＝F(t)-ψ_d(t)，ψ_d(t) is k_zA Gabor atom corresponding to the d-th better time frequency parameter in the better time frequency parameters; r | |^d||_ξRepresents R^dXi-norm of which xi is a constant and is more than or equal to 0 and less than or equal to 1; d is a positive integer and 1, 2, …, k_z。

5. An echo signal separating method for ultrasonic testing according to claim 1 or 2, characterized in that: before echo signal separation in the third step, updating judgment of the optimal atom set of iterative decomposition is needed;

step D, updating the iterative decomposition optimal atom set: updating the iterative decomposition optimal atomic set in the second step by adopting a data processing device (2), wherein the method comprises the following steps:

step D11, atom random selection: randomly taking out a best matching atom from the iterative decomposition best atom set at the moment by adopting a data processing device (2) as an atom to be optimized, wherein the atom to be optimized is marked as

Wherein j1 is a positive integer and j1 is not less than 1 and not more than K;

Is recorded as a time-frequency parameter r_j1'Time-frequency parameter r_j1'＝(s_j1',u_j1',v_j1',w_j1')；

For the best matching atom

When searching, according to the preset s_j1'、u_j1'、v_j1'And w_j1'The value range of the adaptive value is found out by adopting the data processing equipment (2) and calling an optimization algorithm module_j1') The maximum optimal time frequency parameter, the found optimal time frequency parameter is the time frequency parameter r_j1'(ii) a According to the formula

Solving for the best matching atom

Wherein the content of the first and second substances,

to represent

And

inner product of (d);

ψ'₀(t) is of K-1 of said atoms to be treated in step D11And;

d13, atom replacement judgment, atom replacement and iterative decomposition optimal atom set updating: adopting data processing equipment (2) and calling a residual value judging module, an adaptability value judging module or a sparsity judging module to judge whether the atoms to be optimized in the step D11 need to be replaced or not, and replacing the atoms to be optimized according to the judging result;

adopting a data processing device (2) and calling a residual value judging module to judge whether the atoms to be optimized in the step D11 need to be replaced or not, and according to the replaced residual value R_j1' ^K(t)||_ξWhether it is less than the residue before replacement | | R_j1 ^K(t)||_ξAnd (4) judging: when | | | R_j1' ^K(t)||_ξ＜||R_j1 ^K(t)||_ξIf so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

wherein R is_j1' ^K(t)＝F(t)-ψ_j1'(t)，ψ_j1'(t)＝ψ₀(t)+ψ_rj1'(t)；R_j1 ^K(t)＝F(t)-ψ_j1(t)，ψ_j1(t) is the sum of the K best matching atoms in the iterative decomposition best atom set before atom replacement judgment in this step; r | |_j1' ^K(t)||_ξRepresents R_j1' ^KXi-norm of (t, | | R_j1 ^K(t)||_ξRepresents R_j1 ^KXi-norm of (t), xi is constant and xi is more than or equal to 0 and less than or equal to 1;

adopting a data processing device (2) and calling a Fitness value judging module to judge whether the atoms to be optimized in the step D11 need to be replaced, and according to the replaced Fitness value Fitness (r)_j1') Whether greater than pre-replacement fitness value Fitness(r_j1) And (4) judging: when Fitness (r)_j1')＞Fitness(r_j1) If so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

wherein the content of the first and second substances,

represents R^j1-1(t) and

represents R^j1-1(t) and

inner product of (d);

adopting data processing equipment (2) and calling a sparsity judging module to judge whether to replace the atoms to be optimized in the step D11 according to | | R^j1'||_ξWhether or not less than R^j1||_ξAnd (4) judging: when | | | R^j1'||_ξ＜||R^j1||_ξIf so, judging that the atom to be optimized in the step D11 needs to be replaced, and replacing the atom to be optimized in the step D11 with the best matching atom in the step D12

wherein R is^j1'Is composed of

Amount of residual error of

R^j1Is composed of

Amount of residual error of

||R^j1'||_ξRepresents R^j1'Xi-norm, | | R^j1||_ξRepresents R^j1ξ -norm of;

6. An echo signal separating method for ultrasonic testing according to claim 5, characterized in that: after signal sparse decomposition is carried out in the step D1, synchronously storing the iterative decomposition optimal atom set into a data storage (3) by adopting data processing equipment (2), wherein the data storage (3) is connected with the data processing equipment (2);

after atom replacement judgment and atom replacement are carried out in the step D13, the updated iterative decomposition optimal atom set is synchronously stored by adopting data processing equipment (2);

after signal sparse decomposition is carried out in the step D1, when the optimal atom set of iterative decomposition is synchronously stored in a data storage (3) by adopting data processing equipment (2), respectively storing K optimal matching atoms in the optimal atom set of iterative decomposition according to the iterative decomposition sequence; wherein the content of the first and second substances,

7. An echo signal separating method for ultrasonic testing according to claim 6, characterized in that: when optimizing the best matching atoms in the step D1, optimizing the best matching atoms in the iterative decomposition best atom set according to the storage sequence;

8. An echo signal separating method for ultrasonic testing according to claim 1 or 2, characterized in that: after the echo signals are separated in the third step, respectively extracting the K separated signals; the signal extraction methods of the K separation signals are the same;

when any one of the K separation signals is subjected to signal extraction, a data processing device (2) is adopted to extract the separation signal, and the separation signal is a signal to be processed and is marked as a signal f (t);

when the data processing device (2) is adopted to extract the signal f (t), the process is as follows:

step A1, sparse decomposition based on an optimization algorithm: adopting a data processing device (2) and calling a sparse decomposition module to carry out iterative decomposition processing on the signal f (t), and converting the signal f (t) into the signal f (t)

decomposing the nth best matching atom in the best atom set for the iteration;

in the formula R^m(t) is the residual error of the signal f (t) after m iterative decompositions, where m is the preset total iterative decompositions times and m is a positive integer, n is a positive integer and n is 1, 2, …, m;a_nthe expansion coefficient of the best matching atom after the nth iterative decomposition and the residual error after the last iterative decomposition is obtained;

adopting data processing equipment (2) for nth iteration decomposition and calling an optimization algorithm module to find out the best matching atom;

is a Gabor atom and

wherein the function ψ (t) is a Gaussian window function and

r_nis composed of

in this step, the best matching atom is found

According to a preset s_n、u_n、v_nAnd w_nThe value range of the adaptive value Fitness (r) is found out by adopting the data processing equipment (2) and calling an optimization algorithm module_n) The maximum optimal time frequency parameter, the found optimal time frequency parameter is the time frequency parameter r_n；

represents R^n-1(t) and

wherein, | | R^m(t)||₂For R in the step A1^m(t) 2-norm, epsilon is a preset residual quantity judgment threshold;

step A3, optimizing the best matching atom, the process is as follows:

step A31, atom random selection: randomly taking out a best matching atom from the iterative decomposition best atom set at the moment by adopting a data processing device (2) as an atom to be optimized, wherein the atom to be optimized is marked as

Is recorded as a time-frequency parameter r_j'Time-frequency parameter r_j'＝(s_j',u_j',v_j',w_j')；

For the best matching atom

When searching, according to the preset s_j'、u_j'、v_j'And w_j'The data processing equipment (2) is adopted and the optimizing algorithm module is called to find out the fitness value fitness (r)_j') The maximum optimal time frequency parameter, the found optimal time frequency parameter is the time frequency parameter r_j'(ii) a According to the formula

Solving for the best matching atom

Wherein the content of the first and second substances,

to represent

And

inner product of (d);

ψ₀(t) is the sum of m-1 of said atoms to be treated in step A31;

step A33, atom replacement judgment and atom replacement: adopting data processing equipment (2) and calling a residual value judging module, an adaptability value judging module or a sparsity judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and replacing the atoms to be optimized according to the judging result;

adopting a data processing device (2) and calling a residual value judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and according to the replaced residual value R_j' ^m(t)||_ξWhether is less thanResidual value R before replacement_j ^m(t)||_ξAnd (4) judging: when | | | R_j' ^m(t)||_ξ＜||R_j ^m(t)||_ξWhen the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein R is_j' ^m(t)＝f(t)-ψ_j'(t)，

R_j ^m(t)＝f(t)-ψ_j(t)，ψ_j(t) is the sum of m best matching atoms in the iterative decomposition best atom set before atom replacement judgment in the step; r | |_j' ^m(t)||_ξRepresents R_j' ^mXi-norm of (t, | | R_j ^m(t)||_ξRepresents R_j ^mXi-norm of (t), xi is constant and xi is more than or equal to 0 and less than or equal to 1;

adopting a data processing device (2) and calling a Fitness value judging module to judge whether the atoms to be optimized in the step A31 need to be replaced or not, and according to the replaced Fitness value Fitness (r)_j') Whether greater than the pre-replacement Fitness value Fitness (r)_j) And (4) judging: when Fitness (r)_j')＞Fitness(r_j) When the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein the content of the first and second substances,

represents R^j-1(t) and

represents R^j-1(t) and

inner product of (d);

adopting data processing equipment (2) and calling a sparsity judging module to judge whether to replace the atoms to be optimized in the step A31 according to | | R^j'||_ξWhether or not less than R^j||_ξAnd (4) judging: when | | | R^j'||_ξ＜||R^j||_ξWhen the atom to be optimized in the step A31 needs to be replaced, the atom to be optimized in the step A31 is replaced by the best matching atom in the step A32

wherein R is^j'Is composed of

Amount of residual error of

R^jIs composed of

Amount of residual error of

||R^j'||_ξRepresents R^j'Xi-norm, | | R^j||_ξRepresents R^jξ -norm of;

wherein, | | R "_j ^m(t)||₂Is R'_j ^m(t) 2-norm; r'_j ^m(t) is according to the iteration at this timeDecomposing residual quantity after m times of iterative decomposition of m best matching atom pairs f (t) in a best atom set;

step A4, signal reconstruction: according to the iterative decomposition optimal atom set at the moment, obtaining an approximate signal f' (t) of a signal f (t) by adopting a data processing device (2); wherein the approximation signal f "(t) is a signal extracted from the signal f (t),

wherein

To this end the iteration decomposes the nth 'best matching atom in the best set of atoms, n' being a positive integer and n ═ 1, 2, …, m; a is_n'Is composed of