CN100507555C - Fractal elimination method for mass data in ultrasonic intelligent detection - Google Patents

Abstract

The related fractal elimination method for great much data in ultrasonic intelligent detection comprises: resolving the great much data storage into a duality recognition problem; with finite noise-contained sequence fractal R/S algorithm, extracting Hurst index from signal, presenting a fractal linear level index to distinguish signal. This invention stores only the defective ultrasonic echo signal, has well real-time feature, and high accuracy and stability.

Description

The fractal elimination method of mass data during ultrasonic intelligent detects

Technical field

The present invention relates to the fractal elimination method of mass data in a kind of ultrasonic intelligent detection, is a kind of long-time, long magnanimity detection processing method of data that produces apart from ultrasonic automatic detection, when supersonic array detects automatically that is specially adapted to.Can distinguish the kind of data accurately, fast, reject useless part in the data, keep useful data, can significantly reduce the storage data volume, solve the technical bottleneck of present data storage.

Background technology

Ultrasonic intelligent detects, detects automatically and the multisensor array detection has wide application market, and these test processs are long-time often, long distance (or large tracts of land) is carried out, require test macro record test data effectively, for ex-post analysis and assessment.Because the ultrasound wave that is used to detect often has very high frequency of operation, also very high to the signals sampling frequency requirement simultaneously in order to obtain information exactly, this data volume that makes test produce is very big, and the real-time storage of data is brought great challenge.Existing memory technology is difficult to competent this work.For example, if sensor array adopts 72 probes, the center sensor frequency is 5MHz, and sample frequency is 40MHz, and the A/D conversion accuracy is 16, and then the data volume of per second is the 5.625G byte, and data volume hourly is the 20250G byte.In fact, flaw indication part shared ratio in whole detection signal is less, the signal of the overwhelming majority is the noise signal that does not have defective, there is no need so write down all echoed signals, only defective echoed signal need be noted just passable.But, how to allow detection system determine independently, in real time that under the situation that nobody participates in it is a technical barrier that signal has zero defect.

At present, reduce the method for storing data volume and mainly contain data compression method and threshold value differentiation scalping method.Although data compression method has been obtained great success at communication field and multimedia technology field, from all being difficult to satisfy the requirement of big data quantity Ultrasonic Detection occasion on compression ratio still is real-time.Moreover, also there is no need all data are compressed.Traditional threshold value is differentiated scalping method and is generally compared according to the time domain amplitude of signal or energy and pre-set threshold.Signal less than threshold value is disallowable, and is retained greater than the signal of threshold value.This method is easier, real-time also can meet the demands, but the amplitude of signal and energy are subject to disturb, all can cause the variation of received signal amplitude and energy such as the subtle change of interference of noise, sensor water journey distance (sensor surface is to detecting surperficial distance) and irradiating angle etc., and these factors all are unpredictable.Moreover, may be in same magnitude with noise for its amplitude of little defective and energy, have in addition be lower than noise (as the situation of negative signal to noise ratio (S/N ratio)), these factors have all limited the discrimination of this method to defective.

Since the eighties, the applied research of fractal technology has entered into the signal Processing field, particularly the fractal research to one dimension finite time sequence has also obtained significant progress (1.R.Esteller, et al.A comparisonof waveform fractal dimension algorithms, IEEE TRANSACTIONS ON CIRCUITS ANDSYSTEM, Vol.48, No.2,2001, pp.177-183.2. Zhao Hui, Hou Jianrong executes the Bole, the Study on Similarity of nonstationary time series data at random, software journal, Vol.15, No.5,, pp.633-640 in 2004.3.P.Barat，Fractal?characterization?of?ultrasonic?signals?frompolycrystalline?materials，Chaos，Solitons?&?Fractals，1998，Vol.9，No.11，pp.1827-1834。4.Rafat?Weron，Estimating?long-range?dependence：finitesample?properties?and?confidence?intervals，PHYSICA?A，2002，pp.285-299)。All there is the report of application in field such as identification of targets in such as water wave data rule, the dividend yield over time diagnosis, water of identification, the speaker's of relation, the voice signal of the time domain waveform of time response, the heart, the electroencephalogram of rule, atmospheric parameter and disease identification, mechanical fault over time.But be used for fractal algorithm and application thereof that ultrasonic signal is handled, relevant report is also considerably less.

Summary of the invention

The objective of the invention is to solve when producing the substantive test data in the present Ultrasonic Detection, the choice problem of data is with the technical bottleneck of reply mass data storage.This is a kind of method based on the statistics fractal theory, is mainly used in the time series of band noise, it by the time domain space on the segmentation self-similarity nature of analytic signal distinguish signal.The amplitude and the energy that do not relate to signal are not subjected to The noise yet, and algorithm is easy, has good real-time, can program realize, can realize by hardware yet.By a large amount of actual measurement ultrasonic signals is verified, prove that this algorithm has very high discrimination, even under the situation of negative signal to noise ratio (S/N ratio), preferable performance is arranged also.

The present invention includes following concrete steps:

1) data to be analyzed is rejected the question resolves itself into binary identification problem, judge at first that promptly detection signal is defective echoed signal, does not still have the noise signal of flaw echo;

2) with data sementation amplitude to be analyzed normalization, utilize segmentation Hurst exponentiation algorithm based on (R/S), extract the Hurst index of ultrasonic signal; R is the discrete-time series cumulative departure, and S is the standard deviation of discrete-time series;

3) extract 2) in the fractal linearity index of respective signal, promptly based on cumulative departure based on the straight line y=kx+b of the segmentation Hurst exponentiation algorithm match of (R/S);

Set up data fractal characteristic staqtistical data base:

4) according to the practical study object, with above-mentioned 1)-3) method different types of known target body is carried out the fractal characteristic statistical study, try to achieve its statistics fractal characteristic, be stored in the property data base use when differentiating;

Measured data is judged classification:

5) to measured data by 1)-3) and method try to achieve its fractal property, utilize 4 again) in the fractal statistical feature data are classified; The zero defect signal, promptly noise signal is disallowable, and the defectiveness echoed signal is retained, and store.

Description of drawings

Fig. 1 is a FB(flow block) of the present invention

Fig. 2 among the present invention as the ultrasonic signal of special case

Fig. 3 among the present invention as the ultrasonic signal fractal characteristic distributed area figure of special case

Fig. 4 among the present invention as the fractal Hurst index map of ultrasonic signal (zero defect echo) of special case

Fig. 5 among the present invention as the fractal Hurst index map of ultrasonic signal (defectiveness echo) of special case

Embodiment

Below in conjunction with drawings and Examples technical scheme of the present invention is further described.

1) problem description

If the signal that any one sensor is received is:

x _i(t)＝s _i(t)+n _i(t)，i＝1，2，…，M (1)

Wherein, s _i(t) be i the flaw echoes that sensor is received, n _i(t) be i the noise signal that sensor is received, M is the total number of array upper sensor.Work as s _iWhen (t) being zero, x _i(t) be the zero defect signal; Work as s _iWhen (t) non-vanishing, x _i(t) be the defectiveness signal.And noise n _i(t) exist always.

x _i(t) after the A/D conversion, obtain discrete series a: X _i=[x _I1, x _I2..., x _IN], i=1,2 ..., N, N are that the signal intercepting is counted.

Utilize this N discrete data point to judge the existence that has or not flaw indication in real time, for the data in real time storage provides safeguard.If from N dimension time series, extract the d dimensional feature vector, can be expressed as: Y _i=[Y _I1, Y _I2..., Y _Id].

If this classification problem has two classifications, be designated as H respectively ₀And H ₁Wherein, H ₀Expression does not have defective, H ₁There is defective in expression.If can find d dimensional feature vector Y _I0And Y _I1, then this binary classification problems can be described as: $\left\{\begin{array}{cc}{X}_i\&Element;{\mathrm{<figure-callout\; id="0"\; label="H"\; filenames="C200610040544E00101.png,C200610040544E00111.png"\; state="\{\{state\}\}">H</figure-callout>}}_0,& Y_i\&Subset;Y_{i0}\\ X_i\&Element;{\mathrm{<figure-callout\; id="1"\; label="H"\; filenames="C200610040544E00092.png,C200610040544E00111.png"\; state="\{\{state\}\}">H</figure-callout>}}_1,& Y_i\&Subset;Y_{i1}\end{array}\right.---\left(2\right)$

Differentiate problem for the ultrasonic signal of being studied, (shown in Fig. 2 (a)) classifies as H1 when having flaw echo and noise simultaneously in the signal.(as Fig. 2 (b)) classifies as H when having only noise in the signal ₀

2) fractal algorithm

If discrete-time series X _i, its dimension is N, average is

(N), standard deviation is S (N), cumulative departure X (n N), then has following relational expression:

$\tilde{x}\left(N\right)=\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{x}_i---\left(3\right)$

$S\left(N\right)={\left[\frac{1}{N}\underset{i=1}{\overset{N}{\mathrm{\&Sigma;}}}{(x_i-\tilde{x}\left(N\right))}^2\right]}^{1/2}---\left(4\right)$

$X(n,N)=\underset{i=1}{\overset{n}{\mathrm{\&Sigma;}}}[x_i-\tilde{x}\left(N\right)],0<n\&le;N---\left(5\right)$

N is tieed up time series (length is N) is divided into d the segment that length is n, get R (N)=max{X (n, N) }-min{X (n, N) }, then

${\left(\frac{R}{S}\right)}_n=\frac{R\left(N\right)}{S\left(N\right)}---\left(6\right)$

Following formula is the ratio of discrete-time series cumulative departure and standard deviation.For fractal sequences, (R/S) _nFollowing statistical property is arranged:

(R/S) _n～cn ^H (7)

Taking the logarithm in both sides, can get

log(R/S) _n＝logc+Hlogn (8)

Utilize the least square fitting straight line, the slope of this straight line is exactly this seasonal effect in time series Hurst index, represents with H.Specifically can adopt following step:

(a) on whole collection, try to achieve (R/S) ₁:

(b) with whole set pair half point, try to achieve a R/S on each collection, two just are averaged have obtained (R/S) ₂

(c) can divide down always, as N/4, N/8 ..., obtain (R/S) ₃, (R/S) ₄

(d) make log-log figure, use the least square fitting straight line, its slope is exactly the estimated value of Hurst index.

For the finite length sequence that has fractal property, its fitting a straight line has good linearity, and its Hurst index can be represented this fractal property well; For the finite length sequence that does not have fractal property, its fitting a straight line linear very poor, its Hurst index is general scatter bigger.

3) fractal linearity index and algorithm

Find that under study for action although said method can find fractal phenomenon in time series, the FRACTAL DIMENSION numerical value of trying to achieve in this way is not exclusively reliable.Because this method adopts the least square fitting straight line, the slope by fitting a straight line provides the Hurst exponential quantity.And in the reality, do not have the point set of trends of straight line can simulate straight line equally, but such point set does not possess fractal property, and the Hurst index that provides does not have practical significance, and this is the deficiency of bringing with least square method.

Adopt a kind of new method to describe the feasibility of this fractal property in the invention first.If feasible, the FRACTAL DIMENSION numerical value that then calculates is of practical significance, otherwise, then nonsensical.Specifically describe as follows:

For the R/S algorithm, establish N dimension time series { x _i, i=1,2 ..., N is divided into 2 ^j(j=0,1,2 ... d) interval of section non-overlapping copies, the length n=N/2 of each section ^j, then can try to achieve d+1 to log (R/S), logn is y=kx+b by the straight line of top method match, then the cumulative departure of fitting a straight line is:

$\mathrm{\&Delta;}=\frac{1}{d+1}\underset{j=1}{\overset{d+1}{\mathrm{\&Sigma;}}}|y_i-\log {(R/S)}_j|---\left(9\right)$

Claim in this research that Δ is fractal linearity index.It can be used as an index weighing the data fractal characteristic.

4) instance analysis and conclusion

When utilizing the R/S method that ultrasonic signal is analyzed,, in the test signal has been carried out normalized in order to reduce the influence of signal amplitude dimension.For the fractal property to these signals has one to get more information about, in Fig. 4 and Fig. 5, listed the fitting result chart of two kinds of type signals.Fitting result chart when wherein Fig. 4 is the zero defect echo, the fractal property fitting result chart when Fig. 5 is defectiveness.

During from the angle analysis zero defect echoed signal of R/S method, also can estimate the fractal property of signal.And from the graftal of match as can be seen significantly fractal property exist, the linearity of its fitting a straight line is fine, as shown in Figure 4.The index Δ of weighing linearity quality hangs down nearly order of magnitude during than the defectiveness signal.Illustrating has fractal property to exist in the noise in the ultrasonic signal.

During from the angle analysis flaw echoes of R/S method, can estimate the Hurst index of signal.But, do not see that from the graftal of match tangible fractal property exists, linear very poor between match point, as shown in Figure 5.The index Δ of weighing its linear degree is bigger.So, can not find fractal property in the flaw indication with the R/S method.This mainly is that the flaw echoes general persistence is shorter, generally between a hundreds of data point (data point is spaced apart 0.01us) because the R/S method is to remove analytic signal from the angle of long-range correlativity.

The single FRACTAL DIMENSION numerical value of signal has randomness, has certain distribution range.The single kind that can't distinguish defective from the FRACTAL DIMENSION numerical value also is difficult to accurately distinguish flaw indication or noise signal.If combine with Hurst exponential quantity and fractal linearity index Δ, can distinguish flaw indication and noise signal effectively.

Listed the proper vector distributed area that ultrasonic sample signal is calculated among Fig. 3.Wherein, left area is the proper vector distributed area during signal when not having defective, and this interval corresponding signal mainly is a noise.As can be seen from the figure, this interval distribution is less, and characteristics of noise vector high conformity is described, wherein, the Hurst index mainly is distributed between 0.45～0.55, and the fractal index high conformity is described; And the linearity of weighing the fractal property quality is also less, mainly is distributed between 0.05～0.1, and fractal property preferably has been described.The distributed area of proper vector when zone, the right is defectiveness among the figure, this interval distribution range is bigger.This mainly is because the flaw indication fractal property is relatively poor, and Another reason causes owing to defect kind is different.Simultaneously, it can also be seen that from figure that two intervals do not have lap, this is very favourable to decision-making.

Table-1 checking result who has listed 100 samples.As can be seen from the table, the accuracy of differentiation is very high, can satisfy the requirement of actual test.

Table-1 pipe ultrasonic test sample book checking result

In order to verify the stability of this algorithm to noise, respectively 100 checking samples have been added the white Gaussian noise of different magnitudes, change the signal to noise ratio (S/N ratio) of signal.Under the situation of filtering signal is adjudicated not carrying out.Table-2 judgement accuracys rate of having listed under the different state of signal-to-noise.Therefrom as can be seen, noise is very little to the correctness influence of identification.Illustrate that the present invention has good stability to noise.

Table-2 noise levels are to judging nicety rate influence table

Claims (1)

1, the fractal elimination method of mass data during a kind of ultrasonic intelligent detects comprises following concrete steps:

The fractal characteristic extraction algorithm:

1) data to be analyzed is rejected the question resolves itself into binary identification problem, judge at first that promptly detection signal is defective echoed signal, does not still have the noise signal of flaw echo;

2) with data sementation amplitude to be analyzed normalization, utilize segmentation Hurst exponentiation algorithm based on (R/S), extract the Hurst index of ultrasonic signal; R is the discrete-time series cumulative departure, and S is the standard deviation of discrete-time series;

3) extract 2) in the fractal linearity index of respective signal, promptly based on cumulative departure based on the straight line y=kx+b of the segmentation Hurst exponentiation algorithm match of (R/S);

Set up data fractal characteristic staqtistical data base:

4) according to the practical study object, with above-mentioned 1)-3) method different types of known target body is carried out the fractal characteristic statistical study, try to achieve its statistics fractal characteristic, be stored in the property data base use when differentiating;

Measured data is judged classification:

5) to measured data by 1)-3) and method try to achieve its fractal property, utilize 4 again) in the fractal statistical feature data are classified; The zero defect signal, promptly noise signal is disallowable, and the defectiveness echoed signal is retained, and store.

Images