CN106841402A - A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm - Google Patents

A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm Download PDF

Info

Publication number
CN106841402A
CN106841402A CN201710005147.3A CN201710005147A CN106841402A CN 106841402 A CN106841402 A CN 106841402A CN 201710005147 A CN201710005147 A CN 201710005147A CN 106841402 A CN106841402 A CN 106841402A
Authority
CN
China
Prior art keywords
signal
algorithm
optimal
reconstruction
ultrasonic phase
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201710005147.3A
Other languages
Chinese (zh)
Inventor
白志亮
陈世利
肖启阳
贾乐成
曾周末
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201710005147.3A priority Critical patent/CN106841402A/en
Publication of CN106841402A publication Critical patent/CN106841402A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/44Processing the detected response signal, e.g. electronic circuits specially adapted therefor

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

The invention discloses a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, the signal reconstruction optimization method is comprised the following steps:Ultrasonic phase array defect detecting system is built, the ultrasonic echo reflected via the defective locations of test specimen is obtained, and is extracted A and sweep signal;Signal is swept to A using orthogonal basis carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication;According to optimal sparse base, using match tracing, orthogonal matching pursuit, the orthogonal matching pursuit of regularization, progressively orthogonal matching pursuit, compression sampling match tracing are reconstructed respectively to ultrasonic phase array signal;The ultrasonic phase array signal reconstruction errors that above-mentioned five kinds greedy restructing algorithms are used under different compression ratios are calculated, and optimal algorithm is selected according to result.The present invention is also contrasted come the quality reconstruction and applicability of verification algorithm by with what Traditional Wavelet compressed.

Description

A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm
Technical field
Field, more particularly to a kind of phased array based on greedy algorithm are perceived the present invention relates to ultrasonic phase array Signal Compression Ultrasonic signal reconstruction and optimization method.
Background technology
Ultrasonic phase array is a technology of field of ultrasonic nondestructive detection newly-developed, in recent years with its detection speed it is fast, Imaging precision is high, adapt to complex topography component the advantages of detection and get growing concern for.The technology is by controlling each hair Penetrate the actuation duration of array element to change the phase that sound wave reaches target area, so as to be deflected and focused beam.Use ultrasound Phased-array technique can obtain the defect distribution situation in whole detected region and can be real-time in the case of not mobile probe Display testing result, for industrial nondestructive testing provides convenience.
With continuing to develop for the technology, it is a burst of that phased array probe has developed into Two-Dimensional Moment by traditional one-dimensional linear array Row, annular array, fan-shaped array, flexible array etc. greatly improve detection suitable for the array format of more high detection demand Efficiency and imaging precision.But, the increase of matrix number and improving for shape also bring a very big problem simultaneously --- The surge of data volume, increased the difficulty of data transmission and processing and the pressure of front-end collection sensor.
The content of the invention
The invention provides a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, the present invention uses pressure Contracting perception algorithm is compressed to signal and reconstructs and determine optimal greedy algorithm by experimental result, the present invention also by with tradition The contrast of wavelet compression carrys out the quality reconstruction and applicability of verification algorithm, described below:
A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, the signal reconstruction optimization method includes Following steps:
Ultrasonic phase array defect detecting system is built, the ultrasonic echo reflected via the defective locations of test specimen is obtained, And extract A and sweep signal;
Signal is swept to A using orthogonal basis carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication;
According to optimal sparse base, using match tracing, orthogonal matching pursuit, the orthogonal matching pursuit of regularization, progressively just Match tracing, compression sampling match tracing is handed over to be reconstructed ultrasonic phase array signal respectively;
Calculate and the above-mentioned five kinds ultrasonic phase array signal reconstruction errors of greedy restructing algorithm and root are used under different compression ratios Optimal algorithm is selected according to result;
The result of optimal algorithm is compared with traditional wavelet compression result, the applicability of verification algorithm.
The ultrasonic phase array defect detecting system includes:The host computer that is sequentially connected electrically, ultrasonic phase array detector, with And ultrasonic phase array probe.
The use orthogonal basis sweeps signal to A and carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication The step of be specially:
Using the degree of rarefication of each sparse transformation of formula quantitative description between L1 norms and L2 norms;
Signal is swept to A carry out discrete Fourier transform and obtain X (k), carry out discrete cosine transform and obtain D (k), and calculate phase The degree of rarefication answered;Wavelet transform is carried out to x (n) using the sym4 wavelet basis of four layers of decomposition and obtains WTf(m, n), and calculate Its degree of rarefication;
Choose common db, totally 54 kinds of wavelet basis are swept signal and decomposed to A for sym, bior, rbio and coif family, point The solution number of plies is set as 2 to 6 layers, and optimal sparse base is chosen according to degree of rarefication result of calculation.
It is described to be specially according to optimal sparse base, the step of be reconstructed to ultrasonic phase array signal using match tracing:
Obtain the coefficient n of corresponding column vector in matrix Ak, column vector a is chosen from matrix An, make it that there is highest with residual error Correlation;
Update and rebuild echo signal;Residual error is updated, is iteratively repeated, update residual values;Return to reconstruction signal s1
It is described according to optimal sparse base, it is specific the step of be reconstructed to ultrasonic phase array signal using orthogonal matching pursuit For:
Column vector a is chosen from matrix An, make it that there is highest correlation, record coefficient correlation n with residual errork;Calculate and work as Optimal approximation coefficient under preceding column vector;
It is iteratively repeated:Residual values are updated, reconstruction signal s is returned2
It is described according to optimal sparse base, ultrasonic phase array signal is reconstructed using the orthogonal matching pursuit of regularization Step is specially:
Coefficient correlation u is calculated, the index position corresponding to K maximum is found from u, be deposited into indexed set J;K It is sparse angle value;
Postsearch screening, the maximum corresponding atom index value of one group of coefficient correlation of selection energy are carried out using regularization method It is stored in J0In, J0It is first element of set J;
Make Ω=Ω ∪ J0, update supported collection AΩ;Signal approximation is carried out using least square method and residual error updates;Ω is rope Draw collection;
| Ω | >=2K is steps be repeated alternatively until, reconstruction signal s is returned3
It is described according to optimal sparse base, the step of being reconstructed to ultrasonic phase array signal using progressively orthogonal matching pursuit Specially:
Set up soft-threshold and calculate the atom set J under the threshold valuek={ j:|Ck(j) | > tkσk, wherein2≤t of threshold parameterk≤3;CkJ () represents interim index set CkIn j-th element;rk-1It is residual error Value;M is the line number of calculation matrix Φ;
Coefficients of Approximation vector is calculated, residual values are updated, reconstruction signal s is returned4
It is described according to optimal sparse base, the step of being reconstructed to ultrasonic phase array signal using compression sampling match tracing Specially:
α K column vectors are chosen from matrix A, makes it that there is highest correlation, recording indexes Ω with residual error;α is tuning Parameter, K is sparse angle value;
Supported collection merges;Calculate the optimal approximation coefficient under current column vector;Preserve maximum preceding K coefficient, iteration weight It is multiple, residual values are updated, return to reconstruction signal s5
Described calculating use the above-mentioned five kinds ultrasonic phase array signal reconstruction errors of greedy restructing algorithm under different compression ratios And be specially according to the step of result selection optimal algorithm:
It the ratio between is the signal length compressed with original signal strength to define compression ratio, chooses percentage mean square error and quantifies Evaluate the precision of restructing algorithm;Part primary signal is removed by random, it is 20%~80% to set compression ratio scope, every 5% 1 grade, calculate the A under different compression ratios and sweep signal reconstruction error and optimal algorithm is selected according to result.
The result by optimal algorithm compares with traditional wavelet compression result, the step of the applicability of verification algorithm Specially:
By the compression ratio for setting different threshold values to adjust wavelet compression, all defect signal is entered using wavelet compression Row is decomposed and reconstituted, and the result of compressed sensing algorithm and wavelet compression, the usability of verification algorithm are compared under identical compression ratio.
The beneficial effect of technical scheme that the present invention is provided is:
1st, the compressed sensing algorithm that the present invention will be invented in recent years is applied to ultrasonic phase array Signal Compression field, achieves Good effect;
2nd, the present invention is reconstructed using five kinds of greedy algorithms to ultrasonic phase array flaw echo, is compared in different compressions Reconstruction accuracy under rate, have selected the optimal restructing algorithm (OMP) of suitable such signal;
3rd, the present invention is tested using multigroup flaw indication, and by the compression effectiveness of new algorithm and classical wavelet compression Contrasted.Result shows, when compression ratio reaches 70%, the average percent reconstructed error of compressed sensing algorithm is only 3.0054%, the reconstruction accuracy with wavelet compression is suitable, fully meets industrial detection demand.
Brief description of the drawings
Fig. 1 is a kind of flow chart of the phased array supersonic signal reconstruction optimization method based on greedy algorithm;
Fig. 2 is the structural representation of ultrasonic phase array defect detecting system;
Fig. 3 is the schematic diagram that defect sets and numbers;
Fig. 4 (a) is the schematic diagram that No. 8 defect A sweep signal;
Fig. 4 (b) is No. 8 schematic diagrames of defect discrete Fourier transform result;
Fig. 4 (c) is No. 8 schematic diagrames of defect discrete cosine transform result;
Fig. 4 (d) is No. 8 schematic diagrames of defect wavelet transform result;
Fig. 5 is the reconstructed error result figure that different restructing algorithms are used under different compression ratios;
Fig. 6 for when compression ratio is 70% all defect use the reconstructed error comparison diagram of OMP algorithms and wavelet compression respectively.
Specific embodiment
To make the object, technical solutions and advantages of the present invention clearer, further is made to embodiment of the present invention below Ground is described in detail.
In recent years, a kind of emerging compressed sensing algorithm receives many such as signal transacting, computer science, applied mathematics The concern of area research person.In classical signal process field, in order to avoid distorted signals, sample frequency necessarily be greater than signal highest More than the twice of frequency, and compressed sensing can keep the prototype structure of signal and by most by non-self-adapting linear projection Optimization method accurate reconstruction signal, such that it is able to break through the limitation of aromatic sampling thheorem.In terms of ultrasonic phase array data compression, What researcher before used is all the conventional methods such as wavelet compression, huffman coding compression, although these methods can also be obtained Good compression effectiveness and reconstruction accuracy, but aromatic sampling thheorem is still followed, sampled signal is not reduced fundamentally.This hair It is bright that compressed sensing algorithm is innovatively applied to ultrasonic phase array Signal Compression field, probe into its effect and feasibility.
The core theory of compressed sensing includes three aspects, is respectively the reconstruct of sparse signal representation, encoding measurement and algorithm. Focal point is placed on algorithm reconstruct aspect by the embodiment of the present invention.The reconstruct of signal is vital one during compressed sensing Step, and quality reconstruction and reconstruction accuracy are to weigh the most important index of compressed sensing application effect.Signal reconstruction is typically to solve for One optimization problem of minimum L0 norms, but the presenter Donoho of compressed sensing etc. points out that the problem is NP-hard problems And cannot solve.In consideration of it, researcher proposes a series of alternate algorithms, mainly include minimal L1 norm method, greedy algorithm, most The small full calculus of variations, iteration method etc..
In restructing algorithm mentioned above, greedy algorithm because fast operation, computation complexity is low, reconstruction accuracy is high and It is widely adopted.Therefore, the embodiment of the present invention have selected a series of greedy algorithms and ultrasonic phase array signal be reconstructed, including Matching pursuit algorithm (MP), orthogonal matching pursuit algorithm (OMP), the orthogonal matching pursuit algorithm (ROMP) of regularization, progressively just Matching pursuit algorithm (StOMP) and compression sampling matching pursuit algorithm (CoSaMP) are handed over, by comparing under different compression ratios Reconstruct mean square error to select optimal algorithm, there is provided a kind of ultrasonic phase array signal reconstruction optimization side based on greedy algorithm Method.Finally, the embodiment of the present invention is also contrasted come the quality reconstruction and applicability of verification algorithm by with what Traditional Wavelet compressed.
Embodiment 1
A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, referring to Fig. 1, signal reconstruction optimization bag Include following steps:
101:Ultrasonic phase array defect detecting system is built, the ultrasound reflected via the defective locations of test specimen is obtained and is returned Ripple, and extract A and sweep signal;
102:Signal is swept to A using orthogonal basis carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication;
103:According to optimal sparse base, using match tracing, orthogonal matching pursuit, the orthogonal matching pursuit of regularization, by Step orthogonal matching pursuit, compression sampling match tracing are reconstructed to ultrasonic phase array signal;
Wherein, the embodiment of the present invention is not limited to the consequence for reconstruction of above-mentioned five kinds of algorithms.
104:The ultrasonic phase array signal reconstruction errors that above-mentioned five kinds greedy restructing algorithms are used under different compression ratios are calculated, And optimal algorithm is selected according to result;
105:The result of optimal algorithm is compared with traditional wavelet compression result, the applicability of verification algorithm.
Referring to Fig. 2, the ultrasonic phase array defect detecting system includes:Host computer 1, the ultrasonic phase array inspection being sequentially connected electrically Survey instrument 2 and ultrasonic phase array probe 3.
Wherein, the use orthogonal basis in step 102 sweeps signal to A and carries out sparse transformation, and is selected by calculating degree of rarefication The step of taking optimal sparse base is specially:
Using the degree of rarefication of each sparse transformation of formula quantitative description between L1 norms and L2 norms;
Signal is swept to A carry out discrete Fourier transform and obtain X (k), carry out discrete cosine transform and obtain D (k), and calculate phase The degree of rarefication answered;Wavelet transform is carried out to x (n) using the sym4 wavelet basis of four layers of decomposition and obtains WTf(m, n), and calculate Its degree of rarefication;
Choose common db, totally 54 kinds of wavelet basis are swept signal and decomposed to A for sym, bior, rbio and coif family, point The solution number of plies is set as 2 to 6 layers, and optimal sparse base is chosen according to degree of rarefication result of calculation.
Wherein, in step 103 according to optimal sparse base, ultrasonic phase array signal is reconstructed using match tracing Step is specially:
Obtain the coefficient n of corresponding column vector in matrix Ak, column vector a is chosen from matrix An, make it that there is highest with residual error Correlation;
Update and rebuild echo signal;Residual error is updated, is iteratively repeated, update residual values;Return to reconstruction signal s1
Wherein, in step 103 according to optimal sparse base, weight is carried out to ultrasonic phase array signal using orthogonal matching pursuit The step of structure, is specially:
Column vector a is chosen from matrix An, make it that there is highest correlation, record coefficient correlation n with residual errork;Calculate and work as Optimal approximation coefficient under preceding column vector;
It is iteratively repeated:Residual values are updated, reconstruction signal s is returned2
Wherein, in step 103 according to optimal sparse base, the orthogonal matching pursuit using regularization is believed ultrasonic phase array The step of number being reconstructed is specially:
Coefficient correlation u is calculated, the index position corresponding to K maximum is found from u, be deposited into indexed set J;K It is sparse angle value;
Postsearch screening, the maximum corresponding atom index value of one group of coefficient correlation of selection energy are carried out using regularization method It is stored in J0In, J0It is first element of set J;
Make Ω=Ω ∪ J0, update supported collection AΩ;Signal approximation is carried out using least square method and residual error updates;Ω is rope Draw collection;
| Ω | >=2K is steps be repeated alternatively until, reconstruction signal s is returned3
Wherein, in step 103 according to optimal sparse base, ultrasonic phase array signal is entered using progressively orthogonal matching pursuit The step of line reconstruction, is specially:
Set up soft-threshold and calculate the atom set J under the threshold valuek={ j:|Ck(j) | > tkσk, wherein 2≤t of threshold parameterk≤3;CkJ () represents interim index set CkIn j-th element;rk-1It is residual values;M is calculation matrix The line number of Φ;
Coefficients of Approximation vector is calculated, residual values are updated, reconstruction signal s is returned4
Wherein, in step 103 according to optimal sparse base, ultrasonic phase array signal is entered using compression sampling match tracing The step of line reconstruction, is specially:,
α K column vectors are chosen from matrix A, makes it that there is highest correlation, recording indexes Ω with residual error;α is tuning Parameter, K is sparse angle value;
Supported collection merges;Calculate the optimal approximation coefficient under current column vector;Preserve maximum preceding K coefficient, iteration weight It is multiple, residual values are updated, return to reconstruction signal s5
Wherein, the above-mentioned five kinds ultrasonic phase arrays of greedy restructing algorithm are used under the different compression ratios of calculating in step 104 Signal reconstruction error is simultaneously specially according to the step of result selection optimal algorithm:
It the ratio between is the signal length compressed with original signal strength to define compression ratio, chooses percentage mean square error and quantifies Evaluate the precision of restructing algorithm;Part primary signal is removed by random, it is 20%~80% to set compression ratio scope, every 5% 1 grade, calculate the A under different compression ratios and sweep signal reconstruction error and optimal algorithm is selected according to result.
Wherein, the result by optimal algorithm in step 105 compares with traditional wavelet compression result, verification algorithm The step of applicability, is specially:
By the compression ratio for setting different threshold values to adjust wavelet compression, all defect signal is entered using wavelet compression Row is decomposed and reconstituted, and the result of compressed sensing algorithm and wavelet compression, the usability of verification algorithm are compared under identical compression ratio.
In sum, the invention provides a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, this Invention is compressed to signal using compressed sensing algorithm and reconstructs and determine optimal greedy algorithm by experimental result, and the present invention is also Contrasted come the quality reconstruction and applicability of verification algorithm by with what Traditional Wavelet compressed.
Embodiment 2
The scheme in embodiment 1 is described in detail with reference to specific computing formula, example, it is as detailed below to retouch State:
201:Ultrasonic phase array defect detecting system is built, the ultrasound reflected via the defective locations of test specimen 4 is obtained Echo, and extract A and sweep signal;
The detailed operation of the step is:
1) ultrasonic phase array defect detecting system is built, the system includes:It is host computer 1, ultrasonic phase array detector 2, super Sound phased array probe 3 and test specimen 4, detecting system are as shown in Fig. 2 defect sets and numbering is as shown in Figure 3.
2) first in the surface smear couplant of test specimen 4, all defect detection is finished, and extracts the A of each defective locations Sweeping signal carries out data processing.
In order to improve detection coverage rate, the embodiment of the present invention carries out defects detection, sample frequency using sector scan mode It is 100MHz.Test specimen 4 can be containing 9 aluminum test blocks of the different-diameter different depth flat-bottom hole of artificial.
The embodiment of the present invention is with 2,36 degree of resin glass of ultrasonic phase array detector of the MULTI2000 models of M2M companies Illustrated as a example by the ultrasonic phase array probe 3 of glass voussoir.When implementing, the embodiment of the present invention is to ultrasonic phase array detector 2nd, the model of ultrasonic phase array probe 3 is not limited, as long as the device of above-mentioned functions can be completed.
202:Signal is swept to A using orthogonal basis carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication;
To be that signal will have openness for important priori conditions in compression sensing algorithm, therefore this step takes common Orthogonal basis sweeps signal to A and carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication, the detailed operation of the step For:
1) the degree of rarefication computational methods of definition signal, embodiment of the present invention selection one kind is between L1 norms and L2 norms Formula carry out the degree of rarefication of each sparse transformation of quantitative description:
In the embodiment of the present invention, N represents the length of signal, θiIt is each coefficient after conversion.
2) signal x (n) is swept to selected A carry out discrete Fourier transform and obtain X (k), and calculate its degree of rarefication.
WN=e-j2π/N
3) discrete cosine transform is carried out to x (n) and obtains D (k), and calculate its degree of rarefication.
4) wavelet transform is carried out to x (n) using the sym4 wavelet basis of four layers of decomposition and obtains WTf(m, n), and calculate it Degree of rarefication.
Wherein, wavelet transform can be obtained by the scale parameter a and translation parameters b in discretization continuous wavelet. Takem,n∈Z,a0≠ 1, a0、b0It is initial coefficients, typically takes a0>1、b0>0;M, n are respectively scale factor And shift factor;Z is set of integers.
By wavelet basis functionObtain It is mother Small echo.
Corresponding discrete wavelet transformer is changed to:
Wherein, ψ * (a0 -mt-nb0) it is the wavelet function generated by morther wavelet;* it is conjugate of symbol.
5) result of wavelet transformation is different and different with wavelet basis and Decomposition order, chooses common db, sym, Totally 54 kinds of wavelet basis are swept signal and are decomposed to A in 2) for bior, rbio and coif family, and Decomposition order is set as 2 to 6 layers, root Optimal sparse base is chosen according to degree of rarefication result of calculation.
203:Ultrasonic phase array signal is reconstructed using match tracing, reconstruction signal is designated as s1
The embodiment of the present invention is calculation matrix using gaussian random matrix in subsequent reconstructIn selection State analysis and show that optimal sparse base is designated asUltrasonic phase array signal is designated as x, and its nonlinear measurement value is y=Ax, its Middle A=Φ Ψ.Ultrasonic phase array signal is reconstructed first by match tracing (MP) algorithm, reconstruction signal is designated as s1
The detailed operation of the step is:
1) initialize:Make s1=0, residual values r0=y, approximation x0=0, iteration count k=1.
2) the coefficient n of corresponding column vector in A is obtainedk.Column vector a is chosen from matrix An, make it that there is highest with residual error Correlation:
Wherein,<>It is inner product operation symbol;rk-1It is current residue.
3) update and rebuild echo signal:
Wherein, xkIt is current approximation value;xk-1It is previous step approximation;For the column vector that previous step is selected.
4) residual error is updated:
5) k=k+1, repeat the above steps 2)~4), until meeting end condition.
6) s is made1=xk, s1Reconstruction signal obtained by as.
204:Ultrasonic phase array signal is reconstructed using orthogonal matching pursuit (OMP) algorithm, reconstruction signal is designated as s2
The detailed operation of the step is:
1) initialize:Make indexed setResidual values r0=y, iteration count k=1;It is empty set.
2) search identification:Column vector a is chosen from matrix An, make it that there is highest correlation with residual error, record is related Coefficient nk
Ωkk-1∪{nk}
Wherein,<rk-1,an>It is current residue rk-1With column vector anInner product;Ωk-1It is indexed set before;ΩkFor new Indexed set.
3) parameter Estimation:Calculate the optimal approximation coefficient x under current column vectork
Wherein,For the column vector that previous step is selected.
4) it is iteratively repeated:Update residual values:
K=k+1, repeat the above steps 2)~4), until meeting end condition.
5) reconstruction signal s is returned2=xk, s2Reconstruction signal obtained by as.
205:Ultrasonic phase array signal is reconstructed using orthogonal matching pursuit (ROMP) algorithm of regularization, reconstruct letter Number it is designated as s3
The detailed operation of the step is:
1) initialize:Residual error r0=y, indexes value setIndex value set Represent empty set.Estimate letter Number degree of rarefication is K=68, iteration count k=1.
2) coefficient correlation u is calculated:U={ un|un=|<rk-1,an>|, n=1,2 ..., N }, wherein rk-1It is current residue; anThe column vector of representing matrix A;unIt is the element in u;<>It is inner product operation symbol.
3) index position corresponding to K maximum is found from u, is deposited into J.
4) carry out postsearch screening using regularization method, the coefficient correlation of the corresponding atom of index value in J is divided into it is some, So that | u (i) |≤2 | u (j) |, i, j ∈ J, then select the corresponding atom index value of one group of coefficient correlation of energy maximum to be stored in J0In, J0It is first element of set J.
5) Ω=Ω ∪ J are made0
6) carry out Signal approximation using least square method and residual error updates:
xk=arg minx||y-AΩx||2
rk=y-AΩxk
Wherein, AΩFor the column vector that previous step is selected;xkIt is current approximation value.
7) k=k+1, repeat the above steps 2)~6), until | Ω | >=2K.
8) s is made3=xk, s3Reconstruction signal obtained by as.
206:Ultrasonic phase array signal is reconstructed using progressively orthogonal matching pursuit (StOMP) algorithm, reconstruction signal It is designated as s4
The detailed operation of the step is:
1) initialize:Residual error r0=y, indexes value setIndex value set Represent empty set.Estimate letter Number degree of rarefication is K=68, iteration count k=1.
2) C is calculatedk=<A,rk-1>, wherein CkIt is interim indexed set, rk-1It is current residue;<>It is inner product operation symbol.
3) soft-threshold h=t is set upkσk, calculate the atom set J under the threshold valuek={ j:|Ck(j) | > tkσk, wherein2≤t of threshold parameterk≤3;CkJ () represents set CkIn j-th element;rk-1It is current residue;M is The line number of calculation matrix Φ.
4) Ω=Ω ∪ J are made0, calculate Coefficients of Approximation xk
xk=(AΩ TAΩ)-1AΩ Ty
Wherein AΩIt is the column vector selected;AΩ TRepresent AΩTransposition.
5) updating residual values is:rk=y-Axk
6) k=k+1, repeat the above steps 2)~5), until meeting end condition.
7) s is made4=xk, s4Reconstruction signal obtained by as.
207:Ultrasonic phase array signal is reconstructed using compression sampling match tracing (CoSaMP) algorithm, reconstruction signal It is designated as s5
The detailed operation of the step is:
1) initialize:Make residual values r0=y, indexed setApproximation x0=0, estimate that signal degree of rarefication is K=68, Iteration count k=1, tuning parameter is α.
2) search identification:α K column vectors a is chosen from matrix An, make it that there is highest correlation with residual error, record Index Ω:
Ω∈arg min|P|≤αKn∈P|<rk-1,an>|
Wherein, P is indexed set;rk-1It is current residue;<>It is inner product operation symbol.
3) supported collection merges:P=supp (xk-1)∪Ω;Supp is supported collection symbol, xk-1Represent previous step approximation collection Close.
4) parameter Estimation:Calculate the optimal approximation coefficient under current column vector:
uk=arg minu||rk-1-APu||2
Wherein, APFor the column vector that previous step determines.
5) screening retains:Preserve maximum preceding K coefficient, xk=[u]K
6) it is iteratively repeated:Update residual values, rk=y-Axk
7) k=k+1, repeat the above steps 2)~6), until meeting end condition.
8) s is made5=xk, s5Reconstruction signal obtained by as.
208:Calculate the ultrasonic phase array signal reconstruction errors that above-mentioned five kinds greedy restructing algorithms are used under different compression ratios And optimal algorithm is selected according to result;
The detailed operation of the step is:
1) it the ratio between is the signal length compressed with original signal strength to define compression ratio (CR):
Wherein, M is the length of measured value.
2) embodiment of the present invention chooses the precision of percentage mean square error (PRD) quantitative assessment restructing algorithm:
It follows that PRD values are smaller, represent that reconstruction accuracy is higher.
3) part primary signal is removed by random, it is 20%~80% to set compression ratio scope, every 5% 1 grade.Make The A under different compression ratios is calculated with foregoing five kinds of greedy algorithms sweep signal reconstruction error.
4) according to result, selection optimal algorithm is considered.
209:All defect signal is compressed using traditional wavelet algorithm and calculates reconstruction accuracy, by optimal algorithm The result of reconstruct compared with wavelet compression result, the applicability of verification algorithm.
The detailed operation of the step is:
1) compression ratio is adjusted in wavelet compression by setting different threshold values, soft-threshold mode is used in the present invention.
2) all defect signal is carried out using wavelet compression it is decomposed and reconstituted,
3) result of compressed sensing algorithm and wavelet compression is compared under identical compression ratio.
In sum, the invention provides a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, this Invention is compressed to signal using compressed sensing algorithm and reconstructs and determine optimal greedy algorithm by experimental result, and the present invention is also Contrasted come the quality reconstruction and applicability of verification algorithm by with what Traditional Wavelet compressed.
Embodiment 3
Feasibility checking is carried out to the scheme in embodiment 1 and 2 with reference to specific test data, it is as detailed below to retouch State:
Device parameters used in the present embodiment are:Centre frequency is popped one's head in for the 64 array element ultrasonic phase arrays of 5MHz, battle array First centre-to-centre spacing is 0.6mm, and material for test to be measured is aluminium, and AD sample frequencys are 100MHz.Host computer CPU is AMD Athlon X4 Four cores, 4GB internal memories, operating system is the 64bit of Windows 7.
1) subsequent analysis are carried out as a example by No. 8 flat-bottom holes of selection first, a diameter of 5mm of the flat-bottom hole, depth is 2mm.Defect Echo A is swept shown in signal such as Fig. 4 (a), signal length N=1024.
2) discrete Fourier transform, discrete cosine transform and wavelet transform are carried out respectively to above-mentioned signal, change is got in return The result for arriving is respectively as shown in Fig. 4 (b), (c), (d).The degree of rarefication calculated using formula described in step 202 is respectively 0.8140th, 0.7559 and 0.8438.By comparing, it is evident that the degree of rarefication of wavelet transformation is higher.
3) carry out Its Sparse Decomposition using other wavelet basis described in step 202 and calculate degree of rarefication, the result for obtaining shows 4 The bior3.1 small echos that layer is decomposed have best degree of rarefication, are 0.8973.
4) to other 8 flat-bottom hole flaw indications also using it is above-mentioned 2)~3) step be analyzed, wherein 5 at 4 layers There is best degree of rarefication under bior3.1 wavelet decompositions, 3 degree of rarefications under the decomposition are also highly desirable in addition.For unified point Analysis, selects 4 layers of bior3.1 small echos of decomposition as the sparse base of whole ultrasonic phase array signals in embodiments of the present invention.
5) using step 203 to five kinds of restructing algorithms described in step 207 to signal described in 1) under all compression ratios It is reconstructed, obtains result as shown in Figure 5.Special instruction, due to the randomness of Gauss measurement matrix, each run knot Fruit has small difference, and the data point in Fig. 5 is 100 average values of computing.
6) by Fig. 5 it will be evident that OMP algorithms all show very low reconstructed error in each compression ratio level, say Understand superiority of the algorithm in the perception of ultrasonic phase array Signal Compression.Particularly when compression ratio reaches 70%, reconstruct is missed Difference is only 3.0901%.
7) using OMP algorithms and wavelet compression all defect signal is carried out respectively it is decomposed and reconstituted, from greedy algorithm 4 layers of bior3.1 small echos of identical are used as wavelet basis.When compression ratio is 70%, OMP algorithms and the weight of wavelet compression algorithm are used Application condition result is as shown in Figure 6.
8) contrast understands that the reconstruction accuracy that OMP algorithms reach can compare favourably with classical wavelet compression, further demonstrates Applicability of the algorithm in ultrasonic phase array Signal Compression.It is noted that the compressed sensing algorithm based on OMP can be by Sampling and compression are carried out simultaneously, and wavelet compression still follows traditional mode.
The foregoing is only preferable preferable example of the invention, be not intended to limit the invention, it is all it is of the invention spirit and Within principle, any modification, equivalent substitution and improvements made etc. should be included within the scope of the present invention.
To the model of each device in addition to specified otherwise is done, the model of other devices is not limited the embodiment of the present invention, As long as the device of above-mentioned functions can be completed.
It will be appreciated by those skilled in the art that accompanying drawing is a schematic diagram for preferred embodiment, the embodiments of the present invention Sequence number is for illustration only, and the quality of embodiment is not represented.
The foregoing is only presently preferred embodiments of the present invention, be not intended to limit the invention, it is all it is of the invention spirit and Within principle, any modification, equivalent substitution and improvements made etc. should be included within the scope of the present invention.

Claims (10)

1. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm, it is characterised in that the signal reconstruction is excellent Change method is comprised the following steps:
Ultrasonic phase array defect detecting system is built, the ultrasonic echo reflected via the defective locations of test specimen is obtained, and carry Take A and sweep signal;
Signal is swept to A using orthogonal basis carries out sparse transformation, and optimal sparse base is chosen by calculating degree of rarefication;
According to optimal sparse base, using match tracing, orthogonal matching pursuit, the orthogonal matching pursuit of regularization, progressively orthogonal Ultrasonic phase array signal is reconstructed respectively with tracking, compression sampling match tracing;
Calculate and the above-mentioned five kinds ultrasonic phase array signal reconstruction errors of greedy restructing algorithm are used under different compression ratios and according to knot Fruit selection optimal algorithm;
The result of optimal algorithm is compared with traditional wavelet compression result, the applicability of verification algorithm.
2. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is that the ultrasonic phase array defect detecting system includes:The host computer that is sequentially connected electrically, ultrasonic phase array detector and Ultrasonic phase array is popped one's head in.
3. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is that the use orthogonal basis sweeps signal to A and carries out sparse transformation, and the step of optimal sparse base is chosen by calculating degree of rarefication It is rapid to be specially:
Using the degree of rarefication of each sparse transformation of formula quantitative description between L1 norms and L2 norms;
Signal is swept to A carry out discrete Fourier transform and obtain X (k), carry out discrete cosine transform and obtain D (k), and calculate corresponding Degree of rarefication;Wavelet transform is carried out to x (n) using the sym4 wavelet basis of four layers of decomposition and obtains WTf(m, n), and it is dilute to calculate its Dredge degree;
Choose common db, totally 54 kinds of wavelet basis are swept signal and decomposed to A for sym, bior, rbio and coif family, decomposition layer Number is set as 2 to 6 layers, and optimal sparse base is chosen according to degree of rarefication result of calculation.
4. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is, it is described according to optimal sparse base, it is specially the step of be reconstructed to ultrasonic phase array signal using match tracing:
Obtain the coefficient n of corresponding column vector in matrix Ak, column vector a is chosen from matrix An, make it that there is highest phase with residual error Guan Xing;
Update and rebuild echo signal;Residual error is updated, is iteratively repeated, update residual values;Return to reconstruction signal s1
5. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is, it is described according to optimal sparse base, it is specially the step of be reconstructed to ultrasonic phase array signal using orthogonal matching pursuit:
Column vector a is chosen from matrix An, make it that there is highest correlation, record coefficient correlation n with residual errork;Prostatitis is worked as in calculating Optimal approximation coefficient under vector;
It is iteratively repeated, updates residual values, returns to reconstruction signal s2
6. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is, it is described according to optimal sparse base, the step being reconstructed to ultrasonic phase array signal using the orthogonal matching pursuit of regularization It is rapid to be specially:
Coefficient correlation u is calculated, the index position corresponding to K maximum is found from u, be deposited into indexed set J;K is dilute Dredge angle value;
Postsearch screening is carried out using regularization method, the maximum corresponding atom index value of one group of coefficient correlation of selection energy is stored in J0In, J0It is first element of set J;
Make Ω=Ω ∪ J0, update supported collection AΩ;Signal approximation is carried out using least square method and residual error updates;Ω is indexed set;
| Ω | >=2K is steps be repeated alternatively until, reconstruction signal s is returned3
7. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is, it is described according to optimal sparse base, have the step of be reconstructed to ultrasonic phase array signal using progressively orthogonal matching pursuit Body is:
Set up soft-threshold and calculate the atom set J under the threshold valuek={ j:|Ck(j) | > tkσk, wherein 2≤t of threshold parameterk≤3;CkJ () represents interim index set CkIn j-th element;rk-1It is residual values;M is calculation matrix The line number of Φ;
Coefficients of Approximation vector is calculated, residual values are updated, reconstruction signal s is returned4
8. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is, it is described according to optimal sparse base, have the step of be reconstructed to ultrasonic phase array signal using compression sampling match tracing Body is:
α K column vectors are chosen from matrix A, makes it that there is highest correlation, recording indexes Ω with residual error;α joins for tuning Number, K is sparse angle value;
Supported collection merges;Calculate the optimal approximation coefficient under current column vector;Maximum preceding K coefficient is preserved, is iteratively repeated, more New residual values, return to reconstruction signal s5
9. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is that described calculating use the above-mentioned five kinds ultrasonic phase array signal reconstruction errors of greedy restructing algorithm and root under different compression ratios The step of selecting optimal algorithm according to result is specially:
It the ratio between is the signal length compressed with original signal strength to define compression ratio, chooses percentage mean square error quantitative assessment The precision of restructing algorithm;Part primary signal is removed by random, it is 20%~80% to set compression ratio scope, every 5% 1 Shelves, calculate the A under different compression ratios and sweep signal reconstruction error and select optimal algorithm according to result.
10. a kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm according to claim 1, its feature It is that the result by optimal algorithm compares with traditional wavelet compression result, has the step of the applicability of verification algorithm Body is:
By the compression ratio for setting different threshold values to adjust wavelet compression, all defect signal is divided using wavelet compression Solution reconstruct, compares the result of compressed sensing algorithm and wavelet compression, the usability of verification algorithm under identical compression ratio.
CN201710005147.3A 2017-01-04 2017-01-04 A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm Pending CN106841402A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710005147.3A CN106841402A (en) 2017-01-04 2017-01-04 A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710005147.3A CN106841402A (en) 2017-01-04 2017-01-04 A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm

Publications (1)

Publication Number Publication Date
CN106841402A true CN106841402A (en) 2017-06-13

Family

ID=59118575

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710005147.3A Pending CN106841402A (en) 2017-01-04 2017-01-04 A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm

Country Status (1)

Country Link
CN (1) CN106841402A (en)

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108416296A (en) * 2018-03-09 2018-08-17 东北电力大学 A kind of ultrasound signal signatures acquisition methods based on improvement sparse decomposition
CN108426945A (en) * 2018-02-10 2018-08-21 天津大学 A kind of phased array supersonic image reconstructing method based on time domain data compression perception
CN108445083A (en) * 2018-02-11 2018-08-24 天津大学 A kind of phased array supersonic image reconstruction optimization method based on frequency domain compressed sensing
CN109632974A (en) * 2018-12-26 2019-04-16 西安科技大学 A kind of ultrasonic examination echo-signal separation method
CN110133628A (en) * 2019-06-13 2019-08-16 哈尔滨工业大学 A kind of restoration methods of LFM Signal Embeded in Strong Interference in short-term
CN110261472A (en) * 2019-03-08 2019-09-20 李言钦 The method that wavelet basis function is rebuild for the measurement of furnace profile flow field sonic method
GB2582545A (en) * 2019-03-14 2020-09-30 Darkvision Tech Compressing ultrasound data in a downhole tool
CN112187404A (en) * 2020-09-02 2021-01-05 北京电子工程总体研究所 Signal reconstruction method and system based on channelized orthogonal matching pursuit
CN112653472A (en) * 2020-12-15 2021-04-13 哈尔滨工程大学 Dolphin whistle call signal reconstruction method based on block sparse compressed sensing
CN112697269A (en) * 2020-12-07 2021-04-23 福州大学 Sound intensity sparse measurement high-resolution imaging method based on simulation sound field prior information
CN113252779A (en) * 2021-04-20 2021-08-13 桂林电子科技大学 Method for improving defect inversion imaging quality
CN113395181A (en) * 2021-06-11 2021-09-14 中国人民解放军陆军勤务学院 Signal measurement method and device, and state monitoring method and device of Internet of things network

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
吴迪等: "分段正则化正交匹配追踪算法", 《光学精密工程》 *
唐朝伟等: "一种稀疏度自适应分段正交匹配追踪算法", 《中南大学学报(自然科学版)》 *
姜剑: "相控阵结构健康监测的数据压缩采样与损伤识别研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
宋丽娟: "一种改进的压缩采样匹配追踪算法的研究", 《电脑知识与技术》 *
杨晓霞: "超声相控阵汽车发动机内腔腐蚀检测关键技术研究", 《中国博士学位论文全文数据库工程科技Ⅱ辑》 *
郭亮等: "基于压缩感知理论的时变信号压缩技术", 《西南交通大学学报》 *

Cited By (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108426945A (en) * 2018-02-10 2018-08-21 天津大学 A kind of phased array supersonic image reconstructing method based on time domain data compression perception
CN108445083A (en) * 2018-02-11 2018-08-24 天津大学 A kind of phased array supersonic image reconstruction optimization method based on frequency domain compressed sensing
CN108416296A (en) * 2018-03-09 2018-08-17 东北电力大学 A kind of ultrasound signal signatures acquisition methods based on improvement sparse decomposition
CN108416296B (en) * 2018-03-09 2021-12-14 东北电力大学 Ultrasonic signal feature acquisition method based on improved sparse decomposition
CN109632974B (en) * 2018-12-26 2021-07-09 西安科技大学 Echo signal separation method for ultrasonic flaw detection
CN109632974A (en) * 2018-12-26 2019-04-16 西安科技大学 A kind of ultrasonic examination echo-signal separation method
CN110261472B (en) * 2019-03-08 2021-12-24 李言钦 Method for measuring and reconstructing hearth section flow field by using wavelet basis function through acoustic wave method
CN110261472A (en) * 2019-03-08 2019-09-20 李言钦 The method that wavelet basis function is rebuild for the measurement of furnace profile flow field sonic method
US11795816B2 (en) 2019-03-14 2023-10-24 Darkvision Technologies Inc. Compressing ultrasound data in a downhole tool
GB2582545B (en) * 2019-03-14 2021-04-14 Darkvision Tech Compressing ultrasound data in a downhole tool
GB2582545A (en) * 2019-03-14 2020-09-30 Darkvision Tech Compressing ultrasound data in a downhole tool
CN110133628A (en) * 2019-06-13 2019-08-16 哈尔滨工业大学 A kind of restoration methods of LFM Signal Embeded in Strong Interference in short-term
CN112187404B (en) * 2020-09-02 2023-08-01 北京电子工程总体研究所 Signal reconstruction method and system based on channelized orthogonal matching pursuit
CN112187404A (en) * 2020-09-02 2021-01-05 北京电子工程总体研究所 Signal reconstruction method and system based on channelized orthogonal matching pursuit
CN112697269B (en) * 2020-12-07 2022-05-31 福州大学 Sound intensity sparse measurement high-resolution imaging method based on simulation sound field prior information
CN112697269A (en) * 2020-12-07 2021-04-23 福州大学 Sound intensity sparse measurement high-resolution imaging method based on simulation sound field prior information
CN112653472A (en) * 2020-12-15 2021-04-13 哈尔滨工程大学 Dolphin whistle call signal reconstruction method based on block sparse compressed sensing
CN113252779A (en) * 2021-04-20 2021-08-13 桂林电子科技大学 Method for improving defect inversion imaging quality
CN113252779B (en) * 2021-04-20 2023-08-01 桂林电子科技大学 Method for improving defect inversion imaging quality
CN113395181A (en) * 2021-06-11 2021-09-14 中国人民解放军陆军勤务学院 Signal measurement method and device, and state monitoring method and device of Internet of things network

Similar Documents

Publication Publication Date Title
CN106841402A (en) A kind of phased array supersonic signal reconstruction optimization method based on greedy algorithm
CN106841401A (en) A kind of phased array supersonic signal reconstruction optimization method based on sensing matrix
Lorintiu et al. Compressed sensing reconstruction of 3D ultrasound data using dictionary learning and line-wise subsampling
CN106814141A (en) A kind of phased array supersonic compression method based on orthogonal matching pursuit
CN105676168B (en) A kind of acoustic vector sensor array direction estimation method
CN106770691A (en) A kind of turbine rotor ultrasonic phase array compression method based on compressed sensing
Lavarello et al. A study on the reconstruction of moderate contrast targets using the distorted Born iterative method
CN109712077B (en) Depth dictionary learning-based HARDI (hybrid automatic repeat-based) compressed sensing super-resolution reconstruction method
CN111956180B (en) Method for reconstructing photoacoustic endoscopic tomographic image
CN109875606B (en) Ultrasonic CT sound velocity imaging method based on prior reflection imaging
CN112754527B (en) Data processing method for low-frequency ultrasonic thoracic imaging
Bai et al. Ultrasonic defect characterization using the scattering matrix: A performance comparison study of Bayesian inversion and machine learning schemas
Jin et al. Sparse deconvolution method for ultrasound images based on automatic estimation of reference signals
Zhang et al. Unrolled convolutional neural network for full-wave inverse scattering
CN114820847A (en) Amplitude extraction method for transmission attenuation ultrasonic tomography
CN115236206A (en) High-frequency ultrasonic denoising method and system based on multipath matching pursuit algorithm improvement
CN112465924B (en) Rapid medical image reconstruction method based on multi-feature fusion
CN113449683A (en) High-frequency ultrasonic sparse denoising method and system based on K-SVD training local dictionary
CN110706298B (en) Regularization weighted least square transmission-reflection dual-mode ultrasonic imaging reconstruction method
CN117147694A (en) Inverse problem-based ultrasonic full-focusing imaging sparse regularization reconstruction method and equipment
Hua et al. Sparse‐View Ultrasound Diffraction Tomography Using Compressed Sensing with Nonuniform FFT
CN111462077A (en) Method for characterizing biological tissues by utilizing nonlinear information entropy
CN108426945A (en) A kind of phased array supersonic image reconstructing method based on time domain data compression perception
Hughes et al. Characterization of digital waveforms using thermodynamic analogs: Applications to detection of materials defects
CN113974560A (en) Sparse array element optimization selection and compressive sensing imaging method for annular photoacoustic tomography system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170613

WD01 Invention patent application deemed withdrawn after publication